[R-br] automatização de ajuste de modelos de regressao

Henrique Dallazuanna wwwhsd em gmail.com
Quarta Agosto 10 14:11:19 BRT 2011


Tente assim:

mod <- lapply(split(dados, dados$parcela), lm, formula = x ~ y)

2011/8/10 Samuel Carvalho <samukajm em yahoo.com.br>:
> 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]
> =============================================
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> 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.
>



-- 
Henrique Dallazuanna
Curitiba-Paraná-Brasil
25° 25' 40" S 49° 16' 22" O


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