2013-05-22 2 views
1

시계열 데이터 y[t],을 가지고 있는데, H(B) = (1 - \phi*B)/(1 + \theta*B) 연산자를 적용하려고합니다. 여기서 phitheta은 상수이고 B은 지연 연산자입니다. R에서 이것을 어떻게 할 수 있습니까? filter 명령을 사용하면 분모가 아닌 분자를 적용 할 수 있습니다.R에 전송 함수 필터를 적용하려면 어떻게해야합니까?

+0

재현 가능한 예를 제공하여 지금까지 시도한 것을 보여줄 수 있습니까? – Nishanth

+0

무슨 뜻인지 확실하지 않습니다. 필터 (y, c (1, -phi))를 사용하여 (1- \ phi B) y_t = y_t - \ phi y_ {t-1}을 적용 할 수 있지만 지연 연산자의 역함수는 적용 할 수 없습니다. 어떻게해야합니까 (1/(1- \ phi B)) y_t? SARIMA 추정의 특정 계수를 시계열 데이터 y_t에 적용해야하기 때문에 필요합니다. 이것은 SARIMA이므로 예측 라이브러리 함수를 사용할 수 없습니다. http://stats.stackexchange.com/questions/43370/filtering-using-arma-model-in-r – anqif

답변

2

filter 함수에 대해 method 인수를 확인하십시오.

method = "recursive"을 각각 사용하여 filter 함수를 두 번 적용하십시오.

0

분모의 테일러 확장은 시리즈 1/(1+x) = 1-x+x^2-x^3+...을 제공합니다. 따라서 \theta<1의 경우 \theta에 원하는 순서로 근사하기위한 적절한 방법을 얻을 수 있습니다. 희망이 도움이됩니다.

면책 조항 : 더 나은 방법이있을 수 있습니다. 나는 시계열 전문가가 아니다.

관련 문제