преобразовать журнал возвращает в фактическую цену прогноза временного ряда, используя R - PullRequest
1 голос
/ 28 апреля 2019

У меня есть эти смоделированные данные, и я установил модель ARMA-GARCH, используя пакет rugarch. Мой код, насколько это возможно,

    ar.sim<-arima.sim(model=list(ar=c(.9,-.2),ma=c(-.7,.1)),n=100)
    logr=diff(log(na.omit(ar.sim)))
require(rugarch)
    gar<-ugarchspec(variance.model = list(model = "sGARCH", garchOrder = c(2, 1)), 
                        mean.model = list(armaOrder = c(2, 1)), 
                        distribution.model = "norm");
    fitg=ugarchfit(spec = gar,data = ar.sim,solver = "hybrid");
    ugarchforecast(fitg,n.ahead =10) 

В этой модели я использовал возврат журналов. Таким образом, мой прогноз также основан на результатах журналов. Но мне нужна актуальная цена. Я гуглил, чтобы найти любую функцию R, которая преобразует этот журнал возврата в фактическую цену. Но я не смог найти ни одного.

Есть ли какая-либо функция в R для извлечения фактической цены из этого журнала возврата мне нужно сделать это вручную?

1 Ответ

2 голосов
/ 28 апреля 2019

Цена должна быть [начальная цена] * exp (кумулятивный возврат журналов). Например:

df <- data.frame(price = c(90, 108, 81, 105, 180))

df$log = log(df$price)
df$logr = c(NA, diff(df$log))
df$logr_na0 = ifelse(is.na(df$logr), 0, df$logr)
df$cuml_log= cumsum(df$logr_na0)
df$reconstructed_price_norm = exp(df$cuml_log)

initial_price <- 90
df$reconstructed_price = initial_price * df$reconstructed_price_norm

выход

> df
  price      log       logr   logr_na0   cuml_log reconstructed_price_norm reconstructed_price
1    90 4.499810         NA  0.0000000  0.0000000                 1.000000                  90
2   108 4.682131  0.1823216  0.1823216  0.1823216                 1.200000                 108
3    81 4.394449 -0.2876821 -0.2876821 -0.1053605                 0.900000                  81
4   105 4.653960  0.2595112  0.2595112  0.1541507                 1.166667                 105
5   180 5.192957  0.5389965  0.5389965  0.6931472                 2.000000                 180
...