Como obter previsões estritamente positivas?

15

Estou trabalhando em uma série temporal cujos valores são estritamente positivos . Trabalhando com vários modelos, incluindo AR, MA, ARMA, etc, não consegui encontrar uma maneira fácil de obter previsões estritamente positivas.

Estou usando R para fazer minhas previsões, e tudo o que pude encontrar foi forecast.hts {hts} que possui um parâmetro positivo descrito aqui:

Prever uma série temporal hierárquica ou agrupada, pacote 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

Alguma sugestão para séries temporais não hierárquicas? E a generalização sobre o uso de outras restrições, como mínimo, máximo, etc?

Mesmo se não implementado em R, sugestões de artigos, modelos ou transformações gerais úteis de variáveis ​​seriam apreciadas.

Ho1
fonte
3
Uma das coisas mais fáceis, mas nem sempre corretas, nesse caso, é simplesmente prever o log da variável.
Mvctas
4
Para ecoar parcialmente o @mpiktas, uma abordagem é trabalhar na escala de log. Na prática, isso geralmente melhora vários aspectos do modelo de uma só vez. Embora os intervalos de previsão se transformem perfeitamente, é necessário ter cuidado com as previsões médias (se a normalidade for razoável nos logs, você poderá obter uma estimativa para a média do normal do log que normalmente é razoável se o tamanho da amostra for grande). Uma alternativa que às vezes pode funcionar para alguns modelos simples de séries temporais é usar um modelo Gamma.
Glen_b -Reinstala Monica

Respostas:

13

Com o forecastpacote para R, basta definir lambda=0ao montar um modelo. Por exemplo:

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

lambdalambdaλ=0 0lambda=0

Veja http://www.otexts.org/fpp/2/4 para mais discussões.

Rob Hyndman
fonte
Obrigado Prof Hyndman por sua ajuda. Acho que deveria reler esse capítulo seriamente! Você acha que mencionar isso no capítulo 2-4 pode ajudar? Acho que sim! :-) Algumas perguntas permanecem para mim: Algum tipo de transformação pode ser usado para valores mínimos (ou máximos) possíveis? Estou tentando fazer isso com uma função baseada em log, mas afinal o intervalo de confiança resultante é matematicamente correto?
• Ho1
11
Faça a pergunta min / max separadamente. Sim, os intervalos de previsão estão corretos quando retrotraduzidos.
amigos estão
11
@ Ho1 Análise de séries temporais aplicada à previsão gerencial da NELSON; O Holden-Day 1973, pp162-165, discute isso em detalhes ... com uma opinião diversa #
IrishStat
Infelizmente não funcionou como esperado, uma vez que mudou o método e, em vez de um agradável esperado variação y em previsões, ele só fez uma linha reta em torno da média
Diego Duarte