我想知道为什么当我在一个只有0的数据帧上运行cor.test来获得相关相似性矩阵时它会返回NA . 在这种情况下，数据框中的0表示基于我们的测量值的小部件的实际丰度 . 所以我希望每个被关联的向量的相关性估计^ 2为1.我认为我没有很好地理解相关性是什么以及如何在cor.test中实现它 . 任何有关理解相关性和cor.test的建议或帮助都会很棒 . 代码如下 .
corpij <- function(i,j,data) {(cor.test(data[,i],data[,j])$estimate)^2}
corp <- Vectorize(corpij, vectorize.args=list("i","j"))
r2_scores <- outer(1:n,1:n,corp,ready_cor)
最好