仅绘制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)