
18 Out
2011
18 Out
'11
17:32
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)