仅绘制 ETS 模型的季节性成分 - R

Question

这可能是个简单的问题，但我的R技能还在学习阶段。我试图在 ETS 时间序列模型中获取季节性成分的图，我还希望 x 轴显示月份。下面的代码将生成下图，我想要的是将最后一部分（季节）制作成一个单一的图表，其中 x 轴不仅显示年份，还显示月份。

library(datasets)
fit <- ets(AirPassengers)
plot(fit)

我试图找到适合对象的正确组件，但未能成功隔离它。

Answer 1

试试这个：

library(forecast)
fit <- ets(AirPassengers)
plot.2=function (x, ...){
    if (!is.null(x$lambda)) 
 y <- BoxCox(x$x, x$lambda)
    else y <- x$x
    plot(cbind(season = x$states[,"s1"]), main = paste("Decomposition by", x$method, 
            "method"), ...)
}

plot.2(fit)

你可能想看看plot.ets

Answer 2

fit$states[,"s1"] 是季节的 ts.data。包 xts 帮助您处理日期格式。

library(datasets); library(xts); library(forecast)

fit <- ets(AirPassengers)
season <- fit$states[,"s1"]

plot(as.xts(season), major.format = "%Y-%m", auto.grid=F)

仅绘制 ETS 模型的季节性成分 - R

Plotting Just the Seasonal Component of ETS Model - R

r

time-series

forecasting