[R-br] Res: ajuste simultaneo

Gustavo Marcatti vgp.gustavo em yahoo.com.br
Quarta Abril 6 14:18:24 BRT 2011


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]
=============================================


      
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110406/3679872a/attachment.html>


Mais detalhes sobre a lista de discussão R-br