automatização de ajuste de modelos de regressao
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] =============================================
Tente assim: mod <- lapply(split(dados, dados$parcela), lm, formula = x ~ y) 2011/8/10 Samuel Carvalho <samukajm@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@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
participantes (2)
-
Henrique Dallazuanna -
Samuel Carvalho