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]
=============================================
_______________________________________________
R-br mailing list
[email protected]
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.
