2017-04-13 17 views
0

나는 비슷한 방법을 사용하여 two 개 이상의 레이어가있을 때 r rastercorLocal 함수를 사용하여 상관 관계를 분석 할 수 있습니다. ppcor 패키지는 pcor 기능을 처리 할 수 ​​있습니다. 그러나 raster 패키지를 사용하여 구현할 수있는 방법을 알고 싶습니다. 여기에 예제가 있지만 분명히 two layers과 함께 작동하며 부분 상관 관계를 사용할 수 없습니다. 제 경우2 개 이상의 레이어가있는 래스터 corLocal

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) 

답변

0

, 난 내 래스터 층들

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

누적 화소 기반 R 및 p- 값

먼저 찾아 corLocal 함수를 사용

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) 

픽셀 기반 부분 상관 분석의 경우 위와 같은 방법으로 처리 할 수 ​​있다고 생각합니다.