Caros(as)
Estou ajustando um mesmo modelo linear para "x" parcelas diferentes e para isso gostaria de gerar um objeto para cada uma destas parcelas e ajustar este modelo.
O que consigo, é fazer este ajuste inserindo o efeito da parcela na composição do modelo linear. Abaixo segue um CMR
dados <- data.frame(parcela=rep(1:3,each=5), x=rnorm(15), y=rnorm(15))
mod1 <- lm(y ~ x*factor(parcela), data=dados) #inserindo o efeito da parcela no modelo
ou
mod1 <- lm(y ~ x, subset(dados, parcela==1)) # este passo que gostaria de fazer de uma maneira mais automatizada
mod2 <- lm(y ~ x, subset(dados, parcela==2))
mod3<- lm(y ~ x, subset(dados, parcela==3))
Pela atenção
Obrigado
====================================
Samuel P. C. Carvalho
Mestre em Ciências Florestais [UFLA]
Doutorando em Recursos Florestais [ESALQ/USP]
=============================================