2014-02-12 1 views
0

안녕하세요. 미리 감사드립니다. 나는 그것이으로 구별 된 후 시리즈 X의 부분 집합에 단위근 테스트를 실행하기 위해 노력하고있어R에서 차분을 계산 한 후 하위 집합의 NA 값을 사용해야하는 이유는 무엇입니까?

X_diff <- diff(X, differences = 1) 

urdfTest(subset(X_diff, dat[["Year"]] > 2001), lags = 4, type = c("c"), doplot = TRUE) 

나는 다음과 같은 오류가 X_diff : 'Error in na.fail.default(as.ts(x)) : missing values in object'

지금 내가 눈치가 나는 전체를 볼 수있는 경우 X_diff 변수 인 경우 NA 값을 얻지 못합니다. 그러나, 내가 X_diff 변수 (아래 코드)의 서브 세트를 보면 끝 부분에서 NA 값을 얻을 수 있으며, 이는 조건에 맞는 모든 연도에 해당됩니다.

TBG_diff[which(dat[["Year"]] > 2001)] 

NA가 나타나지 않고 어떻게이 NA 오류가없이으로 구별 시리즈의 부분 집합을 실행할 수 있습니까?

+0

'traceback() '의 결과는 무엇입니까? – AdamO

+0

'X_diff <- X_diff [! which (is.na (X_diff))]'를 시도한 후에도 여전히 문제가 발생합니까? – Jota

+1

@Frank 당신은 x_diff [! is.na (x_diff)]를 의미 했습니까? 여기에'~'이 필요 없습니다. '--which (..)'도 효과가 있었지만'! (...)'은 효과가 없었습니다. – Vincent

답변

1

차이점은 차이점이 있으므로 존재하지 않는 연도를 찾고 있다는 것입니다 (데이터가있는 첫 해 이전 연도 또는 마지막 연도 이후 연도). 년에 대한 데이터가 있습니다.)

솔루션 년

편집 좀 더 그것에 대해 생각

의 실제 세트 1 년 이상 짧은 집합에 대한 차이 것, 나는 당신이보다는 당신 전에 차이를 부분 집합해야한다고 생각 다른 방법으로. 비슷한 것

X_2001 <- subset(X, dat[["Year"]] > 2001) 
urdfTest(X_2001, MORE CODE HERE) 
+0

이 답변은 명확하게 시계열 차이에 대한 논리에 관한 것이며, R의 구현 세부 사항보다는 많은 구현에 적용될 것입니다. 이것이 CV에서 주제와 관련이없는 것으로 여겨지는 이유는 분명하지 않습니다. –

+0

전적으로 R 코드 인 것처럼 보입니다. –

+0

나는 그 질문이 그런 식으로 질문을 받았다는 것에 동의하지만 실제 문제는 통계적 분석 (차이 변환이 어떤 영향을 미칠지 생각하는 의미에서)의 하나 인 것처럼 보인다. 아마 나는 잘못 본 것이다. 그것은 처음이 아닐 것이다. –

관련 문제