Phylogenetic Pearson r using phyl.vcv function in “phytools”

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?