2016-07-08 11 views
2

これはおそらく簡単な質問ですが、私のRスキルはまだ学習段階にあります。私はETS時系列モデルで季節成分のプロットを取得しようとしており、x軸に月を表示したいと考えています。下のコードは、次のグラフを生成します。私が望むのは、最後のセクション(シーズン)をとり、x軸に年だけでなく月を示す単一のグラフにすることです。ETSモデルのちょうど季節成分をプロットする - R

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

ETS Plot

私はフィット、オブジェクトの右側のコンポーネントを見つけることを試みたが、私はそれを単離するのに失敗しました。

答えて

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) 

enter image description here

+0

ありがとう、これはちょうど私が必要としたものでした。私はそれがフィット$状態の "s1"部分だと思っていましたが、私はそれに関する文書を見つけることができず、私の一般的な知識はちょうど危険なほどです。 – user1723699

2

が、これは試してみてください:

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) 

あなたはplot.ets

+0

おかげで、これはだけでなく、私は答えとしてマークされたものを作品、その他が最初だったように見えます(正直に言うことができない)が、私はあなたに与えましたあなたの助けをupvote、私はあなたに両方のチェックを与えることができれば願っています。再度、感謝します。 – user1723699

関連する問題