使用 as.POSIXlt/ct 进行 R 预测

Question

美好的一天

我在此处的一篇文章中读到“函数 Forecast::plot.forecast 并非设计用于与 axis.Date 或 axis.POSIXct (在包 Forecast 中不使用)一起使用。”这可以在这里看到:custom axis labels plotting a forecast in R

尽管如此，他们还是设法使用预测包和一些代码来获取正确的轴标签。但是，此示例适用于季度数据。另外，这里使用“as.POSIXlt”的示例适用于每周数据:Forecasting time series data

我尝试过使用该代码，但无法使其适用于每月数据。所以我的轴标签仍然是错误的。我很困惑。请您帮助建议我如何使用预测包正确反射(reflect)轴标签？

示例

library(forecast)
headcount<-c(2475,2468,2452,2464,2500,2548,2536,2565,2590,2608,2625,2663)
date<-c("2013/01/31","2013/02/28","2013/03/31","2013/04/30","2013/05/31","2013/06/30",
"2013/07/31","2013/08/31","2013/09/30","2013/10/31","2013/11/30","2013/12/31")

x<-data.frame(headcount,date)
t<-ts(x$headcount,start=c(2013,1),end=c(2013,12),frequency=12)
fit<-forecast(t,h=12)
plot(forecast(fit))

通过这样做，轴标签将显示为 2013.0、2013.5、2014.5

我知道这只是一年的数据。我只是对如何修复每月数据的轴标签感兴趣，

亲切的问候d

Answer 1

这是使用提供的链接的可能解决方案

plot(forecast(fit), axes = FALSE)
a <- seq(as.Date(date[1]) + 1, by = "months", length = length(date) + 11)
axis(1, at = as.numeric(a)/365.3 + 1970, labels = format(a, format = "%m/%Y"), cex.axis = 0.9)
axis(2, cex.axis = 0.9)