The phyl.vcv function in 'phytools' R package computes a phylogenetic trait variance-covariance matrix between two variables.
Can I use this matrix to compute a phylogenetic Pearson r value? If so, can this r-value be used in a t test (with n-2 df) to test the significance of the correlation?