Question

我想计算p值并将其与相关栅格地图合并但是我不知道如何计算光栅p值并用相关栅格图来表示它。我的R编程代码如下所示

 r.stack <- stack
(y13,y14,y15,y16,y17,y18,y19,y20,y21,y22,y23,y24,y1,y2,y3,y4,y5,y6,y7,y8,y9,y10,y11,y12)
z1 <- r.stack[[1:12]]
z2 <- r.stack[[13:24]]
z3 <- stack(z1,z2)


r <- calc(z3, fun=function(x) cor(x[1:12], x[13:24], method='pearson'))

plot(r,col=colorRampPalette(c("red", "yellow", "blue"))(255))

# I have no idea how to get P- value from the Correlation #

而y1到y12是干旱（SEPI）光栅图像和 y 13到24是NDVI光栅图像

提前感谢您的帮助.....

Answer 1

我已经得到了答案，希望它也可以帮助其他人。我corLocal打算corr。代码如下所示，

r.stack <- stack (y13,y14,y15,y16,y17,y18,y19,y20,y21,y22,y23,y24 ,y3, y4,y5,y6, y7,y8,y9,y10,y11,y12,y1,y2)

z1 <- r.stack[[1:12]] z2 <- r.stack[[13:24]] z3 <- stack(z1,z2)

r <- corLocal(r.stack[[1:12]], r.stack[[13:24]], test=TRUE ) plot(r) # r and p-value map will be shown

# only correlation cells where the p-value < 0.05 will be ommited

r.mask <- mask(r[[1]], r[[2]] < 0.05, maskvalue=FALSE)

plot(r.mask)

如何从相关栅格中查找p值

1 个答案: