<div dir="ltr"><div><div><div><div><div><div><div>Boa tarde,<br><br></div>Gostaria de comparar a variação no acúmulo de N foliar entre duas comunidades vegetais utilizando regressões lineares filogenéticas (PGLS ou OLS). Tal análise foi realizada por Penũelas et al 2011 (tabela 2) ao comparar o acúmulo de nutrientes floiares entre duas comunidades vegetais, uma em Borneo e a outra no Hawaii (<a href="http://www.esajournals.org/doi/abs/10.1890/ES10-00185.1">http://www.esajournals.org/doi/abs/10.1890/ES10-00185.1</a>).<br>
</div>O pacote CAPER possui a função 'pgls' capaz de rodar esse tipo de análise, entretanto, eu só encontrei exemplos sobre comparações de duas variáveis dentro da mesma comunidade:<br><br>require(caper)<br>data(shorebird)<br>
shorebird <- comparative.data(shorebird.tree, shorebird.data, Species, vcv=TRUE)<br>mod <- pgls(log(Egg.Mass) ~ log(M.Mass), shorebird)<br></div>summary(mod)<br><br><br></div>Alguém poderia me dar uma luz sobre como fazer isso, mas comparando duas populações?<br>
</div>Agradeço imensamente por qualquer ajuda.<br></div>Abraço<br><br></div>Marcelo<br></div>