当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])})