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

O unico problema com essa solucao e' a inducao de vicio..... b 2011/4/6 Gustavo Marcatti <vgp.gustavo@yahoo.com.br>:

Samuel, acho que você está querendo o ajuste por minimos quadrados em dois estagios. Pela sua área e pelo modelo de exemplo acho que você está interessado em ajustar o modelo de Clutter. Existe o pacote sem e systemfit para MQ2S. Inclusive um dos autores de systemfit é da nossa área de Ciência Florestal. Confesso que nunca utilizei esses pacotes, eu ajusto via lm mesmo. Observe a função: clutter <- function(data){         est1 <- lm(lnB2 ~ offset(x1) + x2 + x3 -1, data = data)         data$lnB2fit <- fitted(est1)         est2 <- lm(lnVtCC2 ~ invI2 + valorS + lnB2fit, data = data)         coefs <- list( a0=coef(est1)[[1]],a1=coef(est1)[[2]],b1=coef(est2)[[1]],              b2=coef(est2)[[2]], b3=coef(est2)[[3]], b4=coef(est2)[[4]] )         unlist(coefs) } Ajusta o modelo de Clutter em uma das parametrizações mais utilizadas.

Gustavo Marcatti

Eng. Florestal UFV ________________________________ De: Samuel Carvalho <samukajm@yahoo.com.br> Para: r-br <r-br@listas.c3sl.ufpr.br> Enviadas: Quarta-feira, 6 de Abril de 2011 13:11:43 Assunto: [R-br] ajuste simultaneo

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