0
각각 60 개의 래스터가있는 두 래스터 스택에서 픽셀 단위의 Granger 인과 테스트를 수행하려고합니다.Raster Stack에 대한 Granger 인과 테스트 R
이library(raster)
library(lmtest)
r <- raster(ncol=10, nrow=10)
r[]=1:ncell(r)
S <- stack(r,r,r,r,r,r,r,r,r,r,r,r)
R <- stack(r,r,r,r,r,r,r,r,r,r,r,r)
FNO2<-stack(S,R)
"lmtest"패키지 번째 사용하여 원래의 기능은 다음과 같습니다 : 아래의 예는 20 래스터를 가지고
D<- grangertest(degg ~ dchick, order=4)
여기에 내가 래스터 스택에 원래 grangertest 기능을 실행했던 수정입니까?
funG <- function(x) { if (is.na(x[1])){ NA } else { grangertest(x[13:24] ~ x[1:12],order=1)}}
granger<-calc(FNO2,funG)
여기서 FNO2는 두 래스터 스택의 스택입니다. 아래 오류가 발생합니다.
Error in `colnames<-`(`*tmp*`, value = c("x", "y", "x_1", "y_1")) :
length of 'dimnames' [2] not equal to array extent
래스터를 위해이 기능을 어떻게 수정합니까?
감사 @RobertH. 이것은 효과가 있었다. –