[R-br] Packge Systemfit method SUR - ajuste simultâneo

[Thiago Wendling gonçalves de oliveira]

Prezados, boa tarde!!!
Meu nome é Thiago e sou doutorando em Eng. Florestal da UFPR
Estou fazendo um trabalho de regressão linear para estimativa de biomassa de folhas, galhos e fuste de uma espécie florestal.A ideia é utilizar o ajuste simultâneo das equações, pois a biomassa aérea é a soma das biomassas mencionadas anteriormente.As variáveis independentes seriam a Biomassa da folha, galhos, fuste e aerea.As variáveis independentes seriam DAP (diâmetro à altura do peito) e Ht (altura total).

O meu código do R foi o seguinte: 
install.packages("systemfit")library(systemfit)
#Equações:#Biomassa_folhas =b0 + b1*DAP+b2*(DAP^2)+b3*((DAP^2)*Ht)#Biomassa_galhos=b0+b1*DAP+b2*(DAP^2)#Biomassa_fuste=b0+b1*DAP+b2*((DAP^2)*Ht)#Biomassa_aerea = Biomassa_folhas+Biomassa_galhos+Biomassa_fuste# Ou -> Biomassa_aerea = (b0 + b1*DAP+b2*(DAP^2)+b3*((DAP^2)*Ht)) + (b0+b1*DAP+b2*(DAP^2)) + (b0+b1*DAP+b2*((DAP^2)*Ht))

#Biomassa folhas
eq1 <- tbl_araucaria$Bfolhas_kg~tbl_araucaria$DAP_cm+tbl_araucaria$dap2+tbl_araucaria$dap2Ht#Biomassa de galhoseq2 <- tbl_araucaria$Bgalhos_kg~tbl_araucaria$DAP_cm+tbl_araucaria$dap2#Biomassa do fusteeq3 <- tbl_araucaria$Bfuste_kg~tbl_araucaria$DAP_cm+tbl_araucaria$dap2Htsystem_eq <- list(eq1 = eq1, eq2 = eq2, eq3 = eq3)
#SUR
sur <- systemfit(system_eq, method = "SUR", data = tbl_araucaria)summary(sur)
Para eu relacionar o ajuste simultâneo com a biomassa aérea (que seria a soma das biomassas) acredito que meu código ainda esteja incompleto e não sei se falta algum procedimento no código systemfit.
Eu escrevi a equação para a biomassa aérea, considerando as minhas variáveis independentes dos três modelos (eq1,eq2,eq3):eq4 <- tbl_araucaria$Barea_kg~tbl_araucaria$DAP_cm+tbl_araucaria$dap2Ht+tbl_araucaria$dap2

Contudo a estimativa ainda não deu certo, pois acredito que as variáveis dependentes para a Biomassa_area ainda não estejam corretas.
Alguém que trabalha com ajuste simultâneo e fez algum procedimento parecido no R poderia me auxiliar?
Desde já agradeço.
Thiago Wendling Gonçalves de Oliveira(41) 98858-2031 Universidade Federal do Paraná
 
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20170927/60144092/attachment.html>