[R-br] Res: ajuste simultaneo

[Benilton Carvalho]

O unico problema com essa solucao e' a inducao de vicio..... b

2011/4/6 Gustavo Marcatti <vgp.gustavo em 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 em yahoo.com.br>
> Para: r-br <r-br em 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 em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>
>