Я пытаюсь сделать так, чтобы ось X отражала дату дня, когда я строю результаты прогноза из модели autoarima. Я использовал autoplot, который, кажется, дает тот же стиль ggplot. Я нашел примеры, но они для ежемесячных или ежеквартальных. Мои данные в днях. Вот мой MWE:
x<-structure(list(A = c(0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L), B = c(77L, 62L, 75L,
25L, 28L, 23L, 11L, 20L, 10L, 17L, 13L, 4L, 2L, 1L, 12L, 17L,
11L, 6L, 12L, 7L), C = c(6L, 11L, 1L, 1L, 7L, 1L, 5L, 3L, 7L,
8L, 6L, 4L, 9L, 0L, 3L, 5L, 4L, 3L, 3L, 0L), D = c(0L, 1L, 1L,
1L, 1L, 0L, 1L, 1L, 0L, 4L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 0L,
0L), Date = c("2020-04-05", "2020-04-04", "2020-04-03", "2020-04-02",
"2020-04-01", "2020-03-31", "2020-03-30", "2020-03-29", "2020-03-28",
"2020-03-27", "2020-03-26", "2020-03-25", "2020-03-24", "2020-03-23",
"2020-03-22", "2020-03-21", "2020-03-20", "2020-03-19", "2020-03-18",
"2020-03-17")), row.names = c(NA, 20L), class = "data.frame")
код:
x[,"Date"] = as.Date(x[,"Date"], format="%Y-%m-%d")
str(dat)
#rearrange
x<-dplyr::arrange(x, Date)
# frequency here in days=365
tsdata<-ts(x$B,frequency = 365 )
autoplot(tsdata, geom = "point", ts.colour = ('dodgerblue3'))
autoarima1<-auto.arima(tsdata)
#h=25 DAYS
d.forecast <- forecast(autoarima1, level = c(95), h = 25)
autoplot(d.forecast,ts.colour = 'firebrick1', predict.colour = 'blue',
predict.linetype = 'dashed', conf.int = FALSE)
autoplot(d.forecast$residuals)
qqnorm(d.forecast$residuals)
autoplot(acf(d.forecast$residuals, plot = FALSE), conf.int.fill = '#0000FF', conf.int.value = 0.95, conf.int.type = 'ma')
autoplot(pacf(d.forecast$residuals, plot = FALSE), conf.int.fill = '#0000FF', conf.int.value = 0.95, conf.int.type = 'ma')
Я хочу получить что-то похожее на это здесь , Есть предложения?