E, se vc estiver se referindo a modelos com erros nas medidas, faca uma revisao da literatura basica (Cook '94; Carroll '96) e veja o pacote simex. b 2011/4/6 Walmes Zeviani <walmeszeviani@gmail.com>:
Samuel,
Adoraria ajudar mais não entendi muito bem. Você quer algo automático, de forma que todas as variáveis que hora estão do lado esquerdo (independentes) sejam usadas do lado direito com variável resposta/dependente? Meu tiro no escuro é esse
da <- as.data.frame(matrix(rnorm(6*50), ncol=6)) formulas <- sapply(1:ncol(da), function(i){ paste(names(da)[i], paste(names(da)[-i], collapse="+"), sep="~") }) formulas
ajustes <- lapply(formulas, function(f){ m0 <- lm(f, data=da) m0 }) lapply(ajustes, summary)
À disposição. Walmes.
========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
Em 6 de abril de 2011 13:11, Samuel Carvalho <samukajm@yahoo.com.br> escreveu:
Pessoal Alguem saberia me dizer como se faz um ajuste simultaneo no R? Algo tipo ### mod1 <- Y ~ X1 + X2 + X3 mod2 <- W ~ Z1 + Z2 + Y ### Valeu Samuel ==================================== 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
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br