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á