Plotten nur die saisonale Komponente des ETS-Modells - R

Question

Dies ist wahrscheinlich eine einfache Frage, aber meine R-Fähigkeiten sind immer noch in den Lernphasen. Ich versuche, nur eine Darstellung der saisonalen Komponente in einem ETS-Zeitreihenmodell zu erhalten und möchte auch, dass die x-Achse die Monate anzeigt. Der folgende Code erzeugt das folgende Diagramm. Ich möchte den letzten Abschnitt (Saison) nehmen und es zu einem einzigen Diagramm machen, wobei auch auf der x-Achse nicht nur die Jahre, sondern auch die Monate angezeigt werden.

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

Ich habe versucht, die richtige Komponente des fit-Objekt zu finden, aber ich war zu isolieren sie erfolglos.

Answer 1

fit$states[,"s1"] haben wollen, ist ts.data der Saison. Paket xts hilft Ihnen mit dem Format des Datums.

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

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

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

Answer 2

Versuchen Sie folgendes:

library(forecast) 
fit <- ets(AirPassengers) 
plot.2=function (x, ...){ 
    if (!is.null(x$lambda)) 
\t 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) 

Sie können einen Blick auf plot.ets