이상 감지를위한 결 측값이있는 시계열의 STL


12

관측치가 누락 된 일련의 기후 데이터에서 이상 값을 감지하려고합니다. 웹을 검색하면서 사용 가능한 많은 접근 방법을 찾았습니다. 이 중 stl 분해는 추세와 계절 성분을 제거하고 나머지를 연구한다는 의미에서 매력적입니다. 읽기 STL : 황토에 계절 - 트렌드 분해 절차 기반은 , stl누락 된 값에도 불구하고 적용 할 이상치의 영향을받지 가능한 가변성을 할당하는 설정을 결정 유연하게 나타납니다. 그러나 그것을 적용하려고 R에 따라 모든 매개 변수를 관찰 4 년으로하고, 정의 http://stat.ethz.ch/R-manual/R-patched/library/stats/html/stl.html , I 만남 오류:

"time series contains internal NAs"(때 na.action=na.omit),
"series is not periodic or has less than two periods"(때 na.action=na.exclude)

주파수가 올바르게 정의되어 있는지 다시 확인했습니다. 블로그에서 관련 질문을 보았지만이를 해결할 수있는 제안을 찾지 못했습니다. 결 stl측값이있는 시리즈 에 적용 할 수 없습니까? 나는 유물을 도입 (및 결과적으로 감지 ...)하고 싶지 않기 때문에 그것들을 보간하는 것을 매우 꺼려합니다. 같은 이유로 ARIMA 접근법을 대신 사용하는 것이 얼마나 좋은지 알 수 없습니다 (그리고 결 측값이 여전히 문제가 될 경우).

stl측값이있는 시리즈 에 적용하는 방법을 알고 있거나 내 선택이 방법 론적으로 적절하지 않다고 생각되거나 더 나은 제안이있는 경우 공유하십시오. 나는 그 분야에서 아주 새롭고 관련 정보의 더미에 압도당했습니다.


같은 문제가 있습니다. STL 메소드는 LOESS를 통해 누락 값을 채울 수 있어야하지만이 R 구현은 그렇게하지 않는 것 같습니다.
jf328

답변:


5

ARIMA 모델은 더미 변수를 쉽게 통합하여 결 측값을 처리합니다. 이것을 펄스 표시기라고합니다. 방법론은 간단하고 http://www.unc.edu/~jbhill/tsay.pdf에 문서화되어 있습니다. 일반적으로이 방법은 펄스, 레벨 시프트, 계절 펄스 및 로컬 시간 추세에 관한 현재 잔차 시리즈 정보를 추출합니다.


제안과 매우 흥미로운 기사에 감사드립니다!
effie

3

이에 @Julius에 따르면 사용하는 것이 가능하다 stl함께 na.approx에서 zoo사용하는 패키지 stl(x, na.action = na.approx, ...). 이것은 일종의 보간수행 합니다.


0

Hafen의 stlplus패키지 를 확인하고 싶을 수도 있습니다 . 자세한 내용은 그의 논문에 있습니다.

CRAN에서 다음과 같이 설치할 수 있습니다.

install.packages("stlplus")

또는 github에서 직접

devtools::install_github("hafen/stlplus").
당사 사이트를 사용함과 동시에 당사의 쿠키 정책개인정보 보호정책을 읽고 이해하였음을 인정하는 것으로 간주합니다.
Licensed under cc by-sa 3.0 with attribution required.