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

Via nls que não altera as estimativas dos parâmetros você pode fazer: mg=nls(r ~ a[t] + b[t]*x,data=d,start=list(a=c(1,1),b=c(1,1))) ma=nls(r ~ a + b[t]*x,data=d,start=list(a=1,b=c(1,1))) anova(mg,ma) mb=nls(r ~ a[t] + b*x,data=d,start=list(a=c(1,1),b=1)) anova(mg,mb) (S,f,P) Allaman \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5076 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

Obrigado Ivan. ________________________________ De: Ivan Bezerra Allaman <ivanalaman@yahoo.com.br> Para: R Brasil <r-br@listas.c3sl.ufpr.br> Enviadas: Terça-feira, 18 de Outubro de 2011 15:51 Assunto: Re: [R-br] Identidade de modelos de regresão via lm Via nls que não altera as estimativas dos parâmetros você pode fazer: mg=nls(r ~ a[t] + b[t]*x,data=d,start=list(a=c(1,1),b=c(1,1))) ma=nls(r ~ a + b[t]*x,data=d,start=list(a=1,b=c(1,1))) anova(mg,ma) mb=nls(r ~ a[t] + b*x,data=d,start=list(a=c(1,1),b=1)) anova(mg,mb) (S,f,P) Allaman \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5076 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature} _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
participantes (2)
-
Ivan Bezerra Allaman
-
Marcelo Cardoso mello