Usualmente utilizo o EVIEWS para ajustar o modelo de Clutter aos meus dados, porém a fim de curiosidade estou tentando ajusta-lo usando o R.

Vi que existe um pacote chamado "systemfit", o qual estou utilizando, porém, não tenho tido sucesso.

library("systemfit")

eqAB <- LNB2 ~ offset(X1) + X2 + X3 -1

eqVol <- LNVol2 ~ INVI2 + S1 + LNB2

system <- list( AB = eqAB, Vol = eqVol )

inst <- ~ X1 + X2 + X3 + INVI2 + S1

instlist <- inst

fit2sls2 <- systemfit( system, "2SLS", inst = instlist, data = dados )

print( fit2sls2 )

__________________________________________________

Gustavo Martins Soares

Doutorando em Silvicultura e Manejo Florestal - ESALQ/USP

Mestre em Silvicultura e Manejo Florestal - ESALQ/USP

Engenheiro Florestal - Universidade Federal de Viçosa (UFV)

CREA: 16633

(19)98179-4424

Laboratório de Ecofisiologia Florestal e Silvicultura - (19)3447-6689

___________________________________________________