<html><body><div style="color:#000; background-color:#fff; font-family:Courier New, courier, monaco, monospace, sans-serif;font-size:12pt"><div><div class="MsoNormal">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.</div>
<div class="MsoNormal"><o:p> </o:p></div>
<div class="MsoNormal">a<- 'a'; b<- 'b'</div>
<div class="MsoNormal">t<-c(a,a,a,a,a,a,a,a,b,b,b,b,b,b,b,b)</div>
<div class="MsoNormal">x<-c(0,0,10,10,30,30,45,45,0,0,10,10,30,30,45,45)</div>
<div class="MsoNormal">r<-c(2,3,4,5,7,9,11,13,5,4,7,8,8,9,10,11)</div>
<div class="MsoNormal">d<-data.frame(t,x,r);d</div>
<div class="MsoNormal">da<-subset(d, t=='a'); da</div>
<div class="MsoNormal"><span lang="EN-US">db<-subset(d,
t=='b'); db<o:p></o:p></span></div>
<div class="MsoNormal"><span lang="EN-US">m<-lm(r~x,
data=d)<o:p></o:p></span></div>
<div class="MsoNormal">ma<-lm(r~x, data=da)</div>
<div class="MsoNormal">mb<-lm(r~x, data=db)</div>
<div class="MsoNormal"><span style="font-size:16.0pt;line-height:115%;color:red">m;ma;mb<o:p></o:p></span></div>
<div class="MsoNormal"><span style="font-size:12.0pt;line-height:115%;font-family:
"Times New Roman","serif";color:red">deviance(m)/df.residual(m)<o:p></o:p></span></div>
<div class="MsoNormal"><span style="font-size:12.0pt;line-height:115%;font-family:
"Times New Roman","serif";color:red">deviance(ma)/df.residual(ma)<o:p></o:p></span></div>
<div class="MsoNormal"><span style="font-size:12.0pt;line-height:115%;font-family:
"Times New Roman","serif";color:red">deviance(mb)/df.residual(mb)<o:p></o:p></span></div></div></div></body></html>