時系列オブジェクトのプロット
時系列分析では、解析対象のデータを可視化することがとても有効です。ts クラスのデータであれば、plot() 関数には時間インデックス情報を自動で図に反映するメソッドがあります。
ここでは eu_stocks データセット（R では既定で EuStockMarkets として利用可能）を使います。このデータには 1991〜1998 年の主要な欧州株価指数の終値が日次で含まれます。具体的には、ドイツ（DAX）、スイス（SMI）、フランス（CAC）、イギリス（FTSE）です。市場が開いている日に観測されているため、週末や祝日の観測値はありません。本演習では、このデータが等間隔で観測され、4 次元の時系列であると近似します。
この章の締めくくりとして、これまでに学んだ関数をこの新しいデータセットに適用してみましょう。
この演習はコースの一部です
Rで学ぶ時系列分析
演習の手順
is.ts()を使って、eu_stocksが ts オブジェクトかどうかを確認します。start()、end()、frequency()をそれぞれ使って、eu_stocksの開始時点、終了時点、頻度を表示します。plot()コマンドでeu_stocksのシンプルなプロットを作成します。ts.plot()コマンドで、eu_stocksのより複雑な時系列プロットを作成します。あらかじめ用意されたコードにeu_stocksデータセットを入力し、他の引数はそのままにしてください。- 用意されたコードを使って、作成した時系列プロットに凡例を追加します。
実践的なインタラクティブ演習
このサンプルコードを完成させて、この演習に挑戦してみましょう。
# Check whether eu_stocks is a ts object
# View the start, end, and frequency of eu_stocks
# Generate a simple plot of eu_stocks
# Use ts.plot with eu_stocks
ts.plot(___, col = 1:4, xlab = "Year", ylab = "Index Value", main = "Major European Stock Indices, 1991-1998")
# Add a legend to your ts.plot
legend("topleft", colnames(eu_stocks), lty = 1, col = 1:4, bty = "n")