r - 当 P 小于 0.05 时，我们如何强制我们的代码计算相关性？

Question

我正在计算两个数据集之间的相关图。我得到的结果如下：

dir1 <- list.files("D:thly", "*.bin", full.names = TRUE)
dir2 <- list.files("D:002", "*.envi", full.names = TRUE)
file_tot <- array(dim = c(1440, 720, 11, 2))

}

我只想在 P 值低于 0.05 时计算相关性。

Answer 1

您的问题不是非常（即根本不是）可重复的，但这应该让您开始。诀窍是计算所有相关性，但仅绘制/返回高于阈值的相关性。该函数cor.test执行相关性测试。

return_cor = function(x, y)  {
    z = cor.test(x,y)
    if(z[[3]] < 0.05) {
        return(z[[4]])
    } else { 
        return(NA) 
    }
}

和一个测试：

x = rnorm(10); y = rnorm(10)
R> return_cor(x, x)
cor 
  1 
R> return_cor(x, y)
[1] NA

所以

resultsr = apply(file_tot, c(1,2), function(x) {return_cor(x[,1],x[,2])})