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).
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 galhos
eq2 <- tbl_araucaria$Bgalhos_kg~tbl_araucaria$DAP_cm+tbl_araucaria$dap2
#Biomassa do fuste
eq3 <- tbl_araucaria$Bfuste_kg~tbl_araucaria$DAP_cm+tbl_araucaria$dap2Ht
system_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.