У меня есть ежедневные временные ряды, например:
45, 89, 56, 87, 56, 88, 45, 90, 54, 87, 65, 77, 97 , 57, 85, 86, 58, 93, 77, 44, 82, 78
Дата начала - 23 марта 2020 г.
Как я могу получить график, полученный путем выполнения следующее (источник: http://jkunst.com/highcharter/):
library("forecast")
airforecast <- forecast(auto.arima(AirPassengers), level = 95)
hchart(airforecast)
Для моих данных я проделал сложный путь:
library("forecast")
y = c(45, 89, 56, 87, 56, 88, 45, 90, 54, 87, 65, 77, 97, 57, 85, 86, 58, 93, 77, 44, 82, 78)
TS = ts(y)
date = as.Date("2020-03-23", "%Y-%m-%d") + 0:21
H = 3
fc <- forecast(auto.arima(TS, stationary = FALSE, seasonal = FALSE), level = 95, h = H)
fdate = last(date) + 1:H
n = length(TS)
alldate = c(date, fdate, fdate)
L = c(rep(NA, n), fc$lower, fc$upper) # 95% Confidence Limits
M = c(rep(NA, n), fc$mean, rep(NA, H))
y = c(y, rep(NA, 2*H))
df = data.frame(date = alldate, y, M, L)
hc <- highchart() %>%
hc_xAxis(categories = df$date,
title = list(text="Date"),
showFirstLabel = TRUE,
showLastLabel = TRUE
) %>%
hc_yAxis(title = list(text="Count"),
showFirstLabel = TRUE,
showLastLabel = TRUE) %>%
hc_add_series(name = "Observed", type = "line", data = df$y) %>%
hc_add_series(name = "95% confidence bound", type = "scatter", data = df$L, dashStyle = "Dot", color = "red") %>%
hc_add_series(name = fc$method, type = "scatter", data = df$M, dashStyle = "Dot", color = "black") %>%
hc_title(text="Predicted")
hc
Затем два доверительных диапазона не выровнены. Почему?
Спасибо!