我可以在 "phytools" 中使用 phyl.vcv 函数来计算系统发育 Pearson r

Can I use phyl.vcv function in "phytools" to compute a phylogenetic Pearson r

'phytools' R 包中的 phyl.vcv 函数计算两个变量之间的系统发育性状方差-协方差矩阵。

我可以使用这个矩阵来计算系统发育 Pearson r 值吗?如果是这样,这个 r 值是否可以用于 t 检验(n-2 df)以检验相关性的显着性?

Liam Revell 在他的博客上使用 R 代码发布了一个解决方案:http://blog.phytools.org/2017/08/pearson-correlation-with-phylogenetic.html