r - 栅格的相关性和 p 值

Question

我有 4 个栅格，我用这段代码计算了相关矩阵；但我需要知道 p 值，但我做不到。你能帮助我吗？

a<-stack(FS2,FP2,W2,P)

jnk=layerStats(a, 'pearson', na.rm=T)
corr_matrix=jnk$'pearson correlation coefficient'

标签： rr-rasterp-valuepearson-correlation

Answer 1

您可以使用 raster 包中的 corLocal 函数：

library(raster)
b <- stack(system.file("external/rlogo.grd", package="raster"))
b <- aggregate(b, 2, mean)

#flip one layer so values differ:
set.seed(0)
b[[2]] <- flip(b[[2]], 'y') + runif(ncell(b))
b[[1]] <- b[[1]] + runif(ncell(b))
plot(b)

#run correlation
x <- corLocal(b[[1]], b[[2]], test=TRUE )

#view pearson correlation and p-value
plot(x)