这段代码运行良好,我得到的结果如下:
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))
resultscor<-apply(file_tot,c(1,2),function(x){cor(x[,1],x[,2],use = "na.or.complete")})
我想仅在P-value is lower than 0.05
. 这个功能将完成这项工作:
return_cor = function(x, y) {
z = cor.test(x,y)
if(z[[3]] < 0.05) {
return(z[[5]])
} else {
return(NA)
}
}
但是我收到了这个错误:
Error in cor.test.default(x, y) : not enough finite observations
这两个函数都工作得很好。我们如何将两个函数合并为一个函数,以便在 P 值为(某个值,阈值)时计算相关性,并且即使少于 3 对也进行计算。