栅格corLocal有两个以上的图层

时间:2017-04-13 16:57:13

标签: r correlation raster r-raster cross-correlation

我想知道当我有超过r raster个图层时,我可以使用corLocal two函数来使用类似的方法来分析相关性。 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)

1 个答案:

答案 0 :(得分:0)

就我而言,我使用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)

对于基于像素的偏相关分析,我认为你可以采用与上述相同的方式。