<html><head></head><body><div style="color:#000; background-color:#fff; font-family:Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:16px"><div id="yiv8175261886"><div id="yui_3_16_0_ym19_1_1506534880703_13828"><div style="background-color: rgb(255, 255, 255);" id="yui_3_16_0_ym19_1_1506534880703_13827"><div id="yiv8175261886yui_3_16_0_ym19_1_1506534880703_7280" style="color: rgb(0, 0, 0); font-family: "Helvetica Neue", Helvetica, Arial, "Lucida Grande", sans-serif; font-size: 16px;"><span></span></div><div style="color: rgb(0, 0, 0); font-family: "Helvetica Neue", Helvetica, Arial, "Lucida Grande", sans-serif; font-size: 16px;"></div><div id="yiv8175261886yui_3_16_0_ym19_1_1506534880703_7264" style="color: rgb(0, 0, 0); font-family: "Helvetica Neue", Helvetica, Arial, "Lucida Grande", sans-serif; font-size: 16px;"><span style="font-family: "Helvetica Neue", "Segoe UI", Helvetica, Arial, "Lucida Grande", sans-serif; font-size: 13px;" id="yui_3_16_0_ym19_1_1506534880703_13942">Prezados, boa tarde!!!</span></div><div id="yiv8175261886yui_3_16_0_ym19_1_1506534880703_7264" style="color: rgb(0, 0, 0); font-family: "Helvetica Neue", Helvetica, Arial, "Lucida Grande", sans-serif; font-size: 16px;"><span style="font-family: "Helvetica Neue", "Segoe UI", Helvetica, Arial, "Lucida Grande", sans-serif; font-size: 13px;"><br></span></div><div id="yiv8175261886yui_3_16_0_ym19_1_1506534880703_7264" style="color: rgb(0, 0, 0); font-family: "Helvetica Neue", Helvetica, Arial, "Lucida Grande", sans-serif; font-size: 16px;"><span style="font-family: "Helvetica Neue", "Segoe UI", Helvetica, Arial, "Lucida Grande", sans-serif; font-size: 13px;" id="yui_3_16_0_ym19_1_1506534880703_14057">Meu nome é Thiago e sou doutorando em Eng. Florestal da UFPR</span></div><div dir="ltr" style="color: rgb(0, 0, 0); font-family: "Helvetica Neue", "Segoe UI", Helvetica, Arial, "Lucida Grande", sans-serif; font-size: 13px;" id="yui_3_16_0_ym19_1_1506534880703_13943"><br id="yui_3_16_0_ym19_1_1506534880703_13944"></div><div dir="ltr" style="color: rgb(0, 0, 0); font-family: "Helvetica Neue", "Segoe UI", Helvetica, Arial, "Lucida Grande", sans-serif; font-size: 13px;" id="yui_3_16_0_ym19_1_1506534880703_13943">Estou fazendo um trabalho de regressão linear para estimativa de biomassa de folhas, galhos e fuste de uma espécie florestal.</div><div dir="ltr" style="color: rgb(0, 0, 0); font-family: "Helvetica Neue", "Segoe UI", Helvetica, Arial, "Lucida Grande", sans-serif; font-size: 13px;" id="yui_3_16_0_ym19_1_1506534880703_13943">A ideia é utilizar o ajuste simultâneo das equações, pois a biomassa aérea é a soma das biomassas mencionadas anteriormente.</div><div dir="ltr" style="color: rgb(0, 0, 0); font-family: "Helvetica Neue", "Segoe UI", Helvetica, Arial, "Lucida Grande", sans-serif; font-size: 13px;" id="yui_3_16_0_ym19_1_1506534880703_13943">As variáveis independentes seriam a Biomassa da folha, galhos, fuste e aerea.</div><div dir="ltr" style="color: rgb(0, 0, 0); font-family: "Helvetica Neue", "Segoe UI", Helvetica, Arial, "Lucida Grande", sans-serif; font-size: 13px;" id="yui_3_16_0_ym19_1_1506534880703_13943">As variáveis independentes seriam DAP (diâmetro à altura do peito) e Ht (altura total).</div><div dir="ltr" style="color: rgb(0, 0, 0); font-family: "Helvetica Neue", "Segoe UI", Helvetica, Arial, "Lucida Grande", sans-serif; font-size: 13px;" id="yui_3_16_0_ym19_1_1506534880703_13943"><br></div><div dir="ltr" style="color: rgb(0, 0, 0); font-family: "Helvetica Neue", "Segoe UI", Helvetica, Arial, "Lucida Grande", sans-serif; font-size: 13px;" id="yui_3_16_0_ym19_1_1506534880703_13943"><br></div><div dir="ltr" style="color: rgb(0, 0, 0); font-family: "Helvetica Neue", "Segoe UI", Helvetica, Arial, "Lucida Grande", sans-serif; font-size: 13px;" id="yui_3_16_0_ym19_1_1506534880703_13943">O meu código do R foi o seguinte: </div><div dir="ltr" style="color: rgb(0, 0, 0); font-family: "Helvetica Neue", "Segoe UI", Helvetica, Arial, "Lucida Grande", sans-serif; font-size: 13px;" id="yui_3_16_0_ym19_1_1506534880703_13943"><br></div><div dir="ltr" style="" id="yui_3_16_0_ym19_1_1506534880703_13943"><div dir="ltr" style="" id="yui_3_16_0_ym19_1_1506534880703_14287"><font face="Helvetica Neue, Segoe UI, Helvetica, Arial, Lucida Grande, sans-serif" id="yui_3_16_0_ym19_1_1506534880703_14288"><span style="font-size: 13px;" id="yui_3_16_0_ym19_1_1506534880703_14289">install.packages("systemfit")</span></font></div><div dir="ltr" style="" id="yui_3_16_0_ym19_1_1506534880703_14290"><font face="Helvetica Neue, Segoe UI, Helvetica, Arial, Lucida Grande, sans-serif" id="yui_3_16_0_ym19_1_1506534880703_14291"><span style="font-size: 13px;" id="yui_3_16_0_ym19_1_1506534880703_14292">library(systemfit)</span></font></div><div dir="ltr" style="" id="yui_3_16_0_ym19_1_1506534880703_14293"><font face="Helvetica Neue, Segoe UI, Helvetica, Arial, Lucida Grande, sans-serif" id="yui_3_16_0_ym19_1_1506534880703_14294"><span style="font-size: 13px;" id="yui_3_16_0_ym19_1_1506534880703_14295"><br id="yui_3_16_0_ym19_1_1506534880703_14296"></span></font></div><div dir="ltr" style="" id="yui_3_16_0_ym19_1_1506534880703_14293"><font face="Helvetica Neue, Segoe UI, Helvetica, Arial, Lucida Grande, sans-serif"><span style="font-size: 13px;">#Equações:</span></font></div><div dir="ltr" style="" id="yui_3_16_0_ym19_1_1506534880703_14293"><font face="Helvetica Neue, Segoe UI, Helvetica, Arial, Lucida Grande, sans-serif" id="yui_3_16_0_ym19_1_1506534880703_14711"><span style="font-size: 13px;" id="yui_3_16_0_ym19_1_1506534880703_14710">#Biomassa_folhas =b0 + b1*DAP+b2*(DAP^2)+b3*((DAP^2)*Ht)</span></font></div><div dir="ltr" style="" id="yui_3_16_0_ym19_1_1506534880703_14293"><font face="Helvetica Neue, Segoe UI, Helvetica, Arial, Lucida Grande, sans-serif" id="yui_3_16_0_ym19_1_1506534880703_14803"><span style="font-size: 13px;" id="yui_3_16_0_ym19_1_1506534880703_14802">#Biomassa_galhos=b0+b1*DAP+b2*(DAP^2)</span></font></div><div dir="ltr" style="" id="yui_3_16_0_ym19_1_1506534880703_14293"><font face="Helvetica Neue, Segoe UI, Helvetica, Arial, Lucida Grande, sans-serif" id="yui_3_16_0_ym19_1_1506534880703_14765"><span style="font-size: 13px;" id="yui_3_16_0_ym19_1_1506534880703_14766">#Biomassa_fuste=b0+b1*DAP+b2*((DAP^2)*Ht)</span></font></div><div dir="ltr" style="" id="yui_3_16_0_ym19_1_1506534880703_14293"><font face="Helvetica Neue, Segoe UI, Helvetica, Arial, Lucida Grande, sans-serif" id="yui_3_16_0_ym19_1_1506534880703_14801"><span style="font-size: 13px;" id="yui_3_16_0_ym19_1_1506534880703_14800">#Biomassa_aerea = Biomassa_folhas+Biomassa_galhos+Biomassa_fuste</span></font></div><div dir="ltr" style="" id="yui_3_16_0_ym19_1_1506534880703_14293"><font face="Helvetica Neue, Segoe UI, Helvetica, Arial, Lucida Grande, sans-serif"><span style="font-size: 13px;"># Ou -> Biomassa_aerea = (</span></font><span style="font-family: "Helvetica Neue", "Segoe UI", Helvetica, Arial, "Lucida Grande", sans-serif; font-size: 13px;" id="yui_3_16_0_ym19_1_1506534880703_14807">b0 + b1*DAP+b2*(DAP^2)+b3*((DAP^2)*Ht)) + (</span><span style="font-family: "Helvetica Neue", "Segoe UI", Helvetica, Arial, "Lucida Grande", sans-serif; font-size: 13px;" id="yui_3_16_0_ym19_1_1506534880703_14812">b0+b1*DAP+b2*(DAP^2)) + (</span><span style="font-family: "Helvetica Neue", "Segoe UI", Helvetica, Arial, "Lucida Grande", sans-serif; font-size: 13px;" id="yui_3_16_0_ym19_1_1506534880703_14828">b0+b1*DAP+b2*((DAP^2)*Ht))</span></div><div dir="ltr" style="" id="yui_3_16_0_ym19_1_1506534880703_14293"><span style="font-family: "Helvetica Neue", "Segoe UI", Helvetica, Arial, "Lucida Grande", sans-serif; font-size: 13px;"><br></span></div><div dir="ltr" style="" id="yui_3_16_0_ym19_1_1506534880703_14293"><span style="font-family: "Helvetica Neue", "Segoe UI", Helvetica, Arial, "Lucida Grande", sans-serif; font-size: 13px;"><br></span></div><div dir="ltr" style="" id="yui_3_16_0_ym19_1_1506534880703_14297"><span style="font-size: 13px; font-family: "Helvetica Neue", "Segoe UI", Helvetica, Arial, "Lucida Grande", sans-serif;" id="yui_3_16_0_ym19_1_1506534880703_14549">#Biomassa folhas</span><br></div><div dir="ltr" style="" id="yui_3_16_0_ym19_1_1506534880703_14310"><font face="Helvetica Neue, Segoe UI, Helvetica, Arial, Lucida Grande, sans-serif" id="yui_3_16_0_ym19_1_1506534880703_14311"><span style="font-size: 13px;" id="yui_3_16_0_ym19_1_1506534880703_14312">eq1 <- tbl_araucaria$Bfolhas_kg~tbl_araucaria$DAP_cm+tbl_araucaria$dap2+tbl_araucaria$dap2Ht</span></font></div><div dir="ltr" style="" id="yui_3_16_0_ym19_1_1506534880703_14313"><font face="Helvetica Neue, Segoe UI, Helvetica, Arial, Lucida Grande, sans-serif" id="yui_3_16_0_ym19_1_1506534880703_14314"><span style="font-size: 13px;" id="yui_3_16_0_ym19_1_1506534880703_14315">#Biomassa de galhos</span></font></div><div dir="ltr" style="" id="yui_3_16_0_ym19_1_1506534880703_14316"><font face="Helvetica Neue, Segoe UI, Helvetica, Arial, Lucida Grande, sans-serif" id="yui_3_16_0_ym19_1_1506534880703_14317"><span style="font-size: 13px;" id="yui_3_16_0_ym19_1_1506534880703_14318">eq2 <- tbl_araucaria$Bgalhos_kg~tbl_araucaria$DAP_cm+tbl_araucaria$dap2</span></font></div><div dir="ltr" style="" id="yui_3_16_0_ym19_1_1506534880703_14319"><font face="Helvetica Neue, Segoe UI, Helvetica, Arial, Lucida Grande, sans-serif" id="yui_3_16_0_ym19_1_1506534880703_14320"><span style="font-size: 13px;" id="yui_3_16_0_ym19_1_1506534880703_14321">#Biomassa do fuste</span></font></div><div dir="ltr" style="" id="yui_3_16_0_ym19_1_1506534880703_14322"><font face="Helvetica Neue, Segoe UI, Helvetica, Arial, Lucida Grande, sans-serif" id="yui_3_16_0_ym19_1_1506534880703_14323"><span style="font-size: 13px;" id="yui_3_16_0_ym19_1_1506534880703_14324">eq3 <- tbl_araucaria$Bfuste_kg~tbl_araucaria$DAP_cm+tbl_araucaria$dap2Ht</span></font></div><div dir="ltr" style="" id="yui_3_16_0_ym19_1_1506534880703_14325"><font face="Helvetica Neue, Segoe UI, Helvetica, Arial, Lucida Grande, sans-serif" id="yui_3_16_0_ym19_1_1506534880703_14326"><span style="font-size: 13px;" id="yui_3_16_0_ym19_1_1506534880703_14327">system_eq <- list(eq1 = eq1, eq2 = eq2, eq3 = eq3)</span></font></div><div dir="ltr" style="" id="yui_3_16_0_ym19_1_1506534880703_14328"><font face="Helvetica Neue, Segoe UI, Helvetica, Arial, Lucida Grande, sans-serif" id="yui_3_16_0_ym19_1_1506534880703_14329"><span style="font-size: 13px;" id="yui_3_16_0_ym19_1_1506534880703_14330"><br id="yui_3_16_0_ym19_1_1506534880703_14331"></span></font></div><div dir="ltr" style="" id="yui_3_16_0_ym19_1_1506534880703_14332"><span style="font-size: 13px; font-family: "Helvetica Neue", "Segoe UI", Helvetica, Arial, "Lucida Grande", sans-serif;">#SUR</span><br></div><div dir="ltr" style="" id="yui_3_16_0_ym19_1_1506534880703_14360"><font face="Helvetica Neue, Segoe UI, Helvetica, Arial, Lucida Grande, sans-serif" id="yui_3_16_0_ym19_1_1506534880703_14361"><span style="font-size: 13px;" id="yui_3_16_0_ym19_1_1506534880703_14362">sur <- systemfit(system_eq, method = "SUR", data = tbl_araucaria)</span></font></div><div dir="ltr" style="" id="yui_3_16_0_ym19_1_1506534880703_14363"><font face="Helvetica Neue, Segoe UI, Helvetica, Arial, Lucida Grande, sans-serif" id="yui_3_16_0_ym19_1_1506534880703_14364"><span style="font-size: 13px;" id="yui_3_16_0_ym19_1_1506534880703_14365">summary(sur)</span></font></div><div dir="ltr" style="" id="yui_3_16_0_ym19_1_1506534880703_14366"><font face="Helvetica Neue, Segoe UI, Helvetica, Arial, Lucida Grande, sans-serif" id="yui_3_16_0_ym19_1_1506534880703_14367"><span style="font-size: 13px;" id="yui_3_16_0_ym19_1_1506534880703_14368"><br id="yui_3_16_0_ym19_1_1506534880703_14369"></span></font></div><div dir="ltr" style="" id="yui_3_16_0_ym19_1_1506534880703_14366"><font face="Helvetica Neue, Segoe UI, Helvetica, Arial, Lucida Grande, sans-serif" id="yui_3_16_0_ym19_1_1506534880703_14424"><span style="font-size: 13px;" id="yui_3_16_0_ym19_1_1506534880703_14423">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.</span></font></div><div dir="ltr" style="" id="yui_3_16_0_ym19_1_1506534880703_14366"><font face="Helvetica Neue, Segoe UI, Helvetica, Arial, Lucida Grande, sans-serif"><span style="font-size: 13px;"><br></span></font></div><div dir="ltr" style="" id="yui_3_16_0_ym19_1_1506534880703_14366"><font face="Helvetica Neue, Segoe UI, Helvetica, Arial, Lucida Grande, sans-serif" id="yui_3_16_0_ym19_1_1506534880703_14479"><span style="font-size: 13px;" id="yui_3_16_0_ym19_1_1506534880703_14480">Eu escrevi a equação para a biomassa aérea, considerando as minhas variáveis independentes dos três modelos (eq1,eq2,eq3):</span></font></div><div dir="ltr" style="" id="yui_3_16_0_ym19_1_1506534880703_14366"><span style="font-family: "Helvetica Neue", "Segoe UI", Helvetica, Arial, "Lucida Grande", sans-serif; font-size: 13px;" id="yui_3_16_0_ym19_1_1506534880703_14469">eq4 <- tbl_araucaria$Barea_kg~tbl_araucaria$DAP_cm+tbl_araucaria$dap2Ht+tbl_araucaria$dap2</span><font face="Helvetica Neue, Segoe UI, Helvetica, Arial, Lucida Grande, sans-serif"><span style="font-size: 13px;"><br></span></font></div><div dir="ltr" style="" id="yui_3_16_0_ym19_1_1506534880703_14366"><font face="Helvetica Neue, Segoe UI, Helvetica, Arial, Lucida Grande, sans-serif"><span style="font-size: 13px;"><br></span></font></div><div dir="ltr" style="" id="yui_3_16_0_ym19_1_1506534880703_14366"><span style="font-size: 13px; font-family: "Helvetica Neue", "Segoe UI", Helvetica, Arial, "Lucida Grande", sans-serif;" id="yui_3_16_0_ym19_1_1506534880703_14564">Contudo a estimativa ainda não deu certo, pois acredito que as variáveis dependentes para a Biomassa_area ainda não estejam corretas.</span></div><div dir="ltr" style="" id="yui_3_16_0_ym19_1_1506534880703_14366"><br></div><div dir="ltr" style="" id="yui_3_16_0_ym19_1_1506534880703_14366"><font face="Helvetica Neue, Segoe UI, Helvetica, Arial, Lucida Grande, sans-serif" id="yui_3_16_0_ym19_1_1506534880703_14612"><span style="font-size: 13px;" id="yui_3_16_0_ym19_1_1506534880703_14611">Alguém que trabalha com ajuste simultâneo e fez algum procedimento parecido no R poderia me auxiliar?</span></font></div><div dir="ltr" style="" id="yui_3_16_0_ym19_1_1506534880703_14366"><font face="Helvetica Neue, Segoe UI, Helvetica, Arial, Lucida Grande, sans-serif"><span style="font-size: 13px;"><br></span></font></div><div dir="ltr" style="" id="yui_3_16_0_ym19_1_1506534880703_14366"><font face="Helvetica Neue, Segoe UI, Helvetica, Arial, Lucida Grande, sans-serif" id="yui_3_16_0_ym19_1_1506534880703_14879"><span style="font-size: 13px;" id="yui_3_16_0_ym19_1_1506534880703_14878">Desde já agradeço.</span></font></div></div><div dir="ltr" style="color: rgb(0, 0, 0); font-family: "Helvetica Neue", "Segoe UI", Helvetica, Arial, "Lucida Grande", sans-serif; font-size: 13px;" id="yui_3_16_0_ym19_1_1506534880703_13951"><br></div><div class="yiv8175261886signature" id="yiv8175261886yui_3_16_0_ym19_1_1506534880703_7279" style="color: rgb(0, 0, 0); font-family: "Helvetica Neue", Helvetica, Arial, "Lucida Grande", sans-serif; font-size: 16px;"><div id="yiv8175261886yui_3_16_0_ym19_1_1506534880703_7278"><b id="yiv8175261886yui_3_16_0_ym19_1_1506534880703_7277">Thiago Wendling Gonçalves de Oliveira</b></div><div style="font-family:arial, helvetica, clean, sans-serif;" id="yui_3_16_0_ym19_1_1506534880703_13833"><i id="yui_3_16_0_ym19_1_1506534880703_13832">(41) 98858-2031</i></div><div style="font-family:arial, helvetica, clean, sans-serif;font-style:italic;" id="yui_3_16_0_ym19_1_1506534880703_13834"> Universidade Federal do Paraná</div><div style="font-family:arial, helvetica, clean, sans-serif;font-style:italic;" id="yui_3_16_0_ym19_1_1506534880703_13835"><br></div><div style="font-family:arial, helvetica, clean, sans-serif;font-style:italic;" id="yui_3_16_0_ym19_1_1506534880703_13836"> </div></div></div></div></div></div></body></html>