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