[R-br] matriz de coeficientes do modelo linear

Samuel Carvalho samukajm em yahoo.com.br
Terça Novembro 26 13:05:34 BRST 2013


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]
=============================================
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20131126/a8343929/attachment.html>


Mais detalhes sobre a lista de discussão R-br