matriz de coeficientes do modelo linear

Olá pessoal, 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. Segue um RCMR de exemplificação com os respectivos comentários. Em seguida apresento uma solução via "nlme". De ante mão já agradeço a ajuda de todos. Forte abraço. Samuel #RMCR data(Orange) head(Orange) class(Orange$Tree) Orange$Tree <- as.numeric(Orange$Tree) ml <- lm(circumference~age*factor(Tree), data=Orange) summary(ml) coef(ml) matrix(coef(ml),length(unique(Orange$Tree)),2) # Nesta matriz que preciso garantir a integridade dos parametros ao seu respectivo grupo # Solucao via nlme require(nlme) ml_list <- lmList(circumference~age|factor(Tree), data=Orange) coef_list <- data.frame(Tree=rownames(coef_list), b0=coef(ml_list)[,1], b1=coef(ml_list)[,2]) coef_list # a ideia é gerar uma saida semelhante a esta porem apenas com a função lm ### ==================================== Samuel P. C. Carvalho Engº Florestal [UFLA] Mestre em Ciências Florestais [UFLA] Doutor em Recursos Florestais [ESALQ/USP] =============================================

Mude a declaração do modelo para ir para a parametrização desejada ml <- lm(circumference~-1+factor(Tree)/age, data=Orange) summary(ml) À disposição. Walmes.
participantes (2)
-
Samuel Carvalho
-
walmes .