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