光柵與兩層以上的corlocal

Question

我想知道如何使用r rastercorLocal函數爲類似的方法來分析我有多於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) 

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) 

對於基於像素的偏相關分析，我認爲你可以用與上面相同的方法。