[R-br] Identidade de modelos de regresão via lm

[Marcelo Cardoso mello]

Como posso comparar os dois modelo (ma e mb) na CMR abaixo?
Tipo um teste de identidade de modelos em lm? Como testar se os tratamentos 'a' e 'b' tem (estatisticamente) a mesma curva de regressão? Obrigado por qualquer
dica.
 
a<- 'a'; b<- 'b'
t<-c(a,a,a,a,a,a,a,a,b,b,b,b,b,b,b,b)
x<-c(0,0,10,10,30,30,45,45,0,0,10,10,30,30,45,45)
r<-c(2,3,4,5,7,9,11,13,5,4,7,8,8,9,10,11)
d<-data.frame(t,x,r);d
da<-subset(d, t=='a'); da
db<-subset(d,
t=='b'); db
m<-lm(r~x,
data=d)
ma<-lm(r~x, data=da)
mb<-lm(r~x, data=db)
m;ma;mb
deviance(m)/df.residual(m)
deviance(ma)/df.residual(ma)
deviance(mb)/df.residual(mb)
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111018/0d8af983/attachment.html>