엄격하게 긍정적 인 예측을 달성하는 방법?


15

나는 값이 엄격하게 양수인 시계열을 연구하고 있습니다. AR, MA, ARMA 등의 다양한 모델을 사용하여 긍정적 인 예측을 얻는 쉬운 방법을 찾을 수 없었습니다.

나는 예측을 위해 R 을 사용 하고 있으며 내가 찾을 수있는 모든 것은 여기에 설명 된 긍정적 인 매개 변수 가있는 predict.hts {hts}입니다 .

계층 적 또는 그룹화 된 시계열 예측, 패키지 HTS

## S3 method for class 'gts':
forecast((object, h,
  method = c("comb", "bu", "mo", "tdgsf", "tdgsa", "tdfp", "all"),
  fmethod = c("ets", "rw", "arima"), level, positive = FALSE,
    xreg = NULL, newxreg = NULL, ...))

positive
    If TRUE, forecasts are forced to be strictly positive

http://www.inside-r.org/packages/cran/hts/docs/forecast.gts

비 계층 적 시계열에 대한 제안 사항이 있습니까? 최소, 최대 등과 같은 다른 제약 조건 사용에 대한 일반화는 어떻습니까?

R로 구현되지 않더라도 기사, 모델 또는 유용한 일반 변수 변환에 대한 제안을 부탁드립니다.


3
이 경우 가장 쉬운 방법은 아니지만 항상 올바른 방법 중 하나는 단순히 변수의 로그를 예측하는 것입니다.
mpiktas

4
@mpiktas를 부분적으로 반향하려면 한 가지 방법은 로그 스케일을 사용하는 것입니다. 실제로 이것은 종종 모델의 여러 측면을 한 번에 개선합니다. 예측 구간은 양호하게 변환되지만 평균 예측은주의해야합니다 (로그에서 정규성이 합리적인 경우 표본 크기가 큰 경우 일반적으로 합리적인 로그 정규 평균의 추정치를 얻을 수 있음). 때때로 간단한 시계열 모델에서 작동 할 수있는 대안은 감마 모델을 사용하는 것입니다.
Glen_b-복지 주 모니카

답변:


13

forecastR 용 패키지를 사용하여 lambda=0모델을 피팅 할 때 간단히 설정하십시오 . 예를 들면 다음과 같습니다.

fit <- auto.arima(x, lambda=0)
forecast(fit)

lambdalambdaλ=0lambda=0

자세한 내용은 http://www.otexts.org/fpp/2/4 를 참조 하십시오 .


친절한 도움을 주신 Hyndman 교수님 께 감사드립니다. 나는 그 장을 진지하게 다시 읽어야한다고 생각합니다! 2-4 장에서 이것을 언급하면 ​​도움이 될 것이라고 생각하십니까? 나도 그렇게 생각해! :-) 몇 가지 질문이 남아 있습니다. 가능한 최소값 (또는 최대 값)으로 어떤 종류의 변환을 사용할 수 있습니까? 로그 기반 함수를 사용 하여이 작업을 수행하려고하지만 결과 신뢰 구간이 수학적으로 정확합니까?
Ho1

1
최소 / 최대 질문은 별도로 문의하십시오. 예, 역변환 될 때 예측 구간이 정확합니다.
Rob Hyndman

1
NELSON의 경영 예측을위한 @ Ho1 응용 시계열 분석; Holden-Day 1973 pp162-165에서 이에 대해 자세히 논의합니다 ... 다양한 의견으로
IrishStat

안타깝게도 방법이 변경되어 예상대로 작동하지 않았으며 예측에 대한 예상되는 변화가 아니라 평균 주위에 평평한 선을 만들었습니다.
Diego Duarte
당사 사이트를 사용함과 동시에 당사의 쿠키 정책개인정보 보호정책을 읽고 이해하였음을 인정하는 것으로 간주합니다.
Licensed under cc by-sa 3.0 with attribution required.