<html><body><div style="color:#000; background-color:#fff; font-family:HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:12pt"><div>Olá pessoal,</div><div>Gostaria da opinião de vocês de como proceder para gerar uma matriz de coeficientes para dados agrupados utilizando apenas a função lm, considerando os grupos como fatores.</div><div>Segue um RCMR de exemplificação com os respectivos comentários. Em seguida apresento uma solução via "nlme".</div><div>De ante mão já agradeço a ajuda de todos.</div><div>Forte abraço.</div><div>Samuel</div><div><br></div><div>#RMCR</div><div><br></div><div>data(Orange)</div><div>head(Orange)</div><div>class(Orange$Tree)</div><div>Orange$Tree <- as.numeric(Orange$Tree)</div><div>ml <- lm(circumference~age*factor(Tree), data=Orange)</div><div>summary(ml)</div><div>coef(ml)</div><div>matrix(coef(ml),length(unique(Orange$Tree)),2) # Nesta matriz que preciso garantir a
integridade dos parametros ao seu respectivo grupo</div><div><br></div><div><br></div><div># Solucao via nlme</div><div>require(nlme)</div><div>ml_list <- lmList(circumference~age|factor(Tree), data=Orange)</div><div>coef_list <- data.frame(Tree=rownames(coef_list), b0=coef(ml_list)[,1], b1=coef(ml_list)[,2])</div><div>coef_list # a ideia é gerar uma saida semelhante a esta porem apenas com a função lm</div><div><br></div><div><span></span></div><div>###</div><div><br></div><div></div><div> </div><div style="font-family: courier, monaco, monospace, sans-serif;"><font size="3"><b>====================================</b></font></div><div style="font-family: courier, monaco, monospace, sans-serif;"><font size="3"><b>Samuel P. C. Carvalho<br></b><font color="#000000" size="2">Engº Florestal [UFLA]</font></font></div><div style="font-family: courier, monaco, monospace, sans-serif; color: rgb(0, 0, 0); font-size: 16px; background-color:
transparent; font-style: normal;"><font size="3"><font color="#000000" size="2">Mestre em Ciências Florestais [UFLA]</font></font></div><div style="font-family: courier, monaco, monospace, sans-serif;"><font color="#000000" size="2">Doutor em Recursos Florestais [ESALQ/USP]</font></div><div style="font-family: courier, monaco, monospace, sans-serif;"><font size="2">=============================================</font></div></div></body></html>