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

###