具有两层以上的栅格 corLocal

Question

我想知道当我有多个two 层时，如何使用r rastercorLocal 函数来分析相关性。 ppcor 包有 pcor 可以处理的函数，但我很想知道如何使用 raster 包来实现它。这是一个示例，但显然它仅适用于 two layers，并且不允许您使用偏相关。

library (raster)
bio <- crop(raster::getData('worldclim', var='bio', res=10),extent(100, 120, 30, 40)) 
biocor <-  corLocal(bio[[1]],bio[[2]], method='kendall',test=F, exact=FALSE)

Answer 1

在我的例子中，我使用 corLocal 函数来查找基于像素的 r 和 p 值

首先，我堆叠了我的栅格图层

Stack1 <- stack  (x1,x2,x3,x4, y1,y2,y3,y4) 

那么，

Corr_coff <- corLocal(Stack1[[5:8]], Stack1[[1:4]], test=TRUE ) # test=T; P-value will be returned 
#TESTING!!!!!! Pearson's r  and p-value map MUST be showed.

plot(Corr_coff)

对于基于像素的偏相关分析，我想你可以按照上面的方法来做。