ARMA モデルのシミュレーション
動画で見たように、任意の定常時系列はホワイトノイズの線形結合として表せます。さらに、どの ARMA モデルもこの形を持つため、定常時系列のモデリングに適した選択肢です。
R には、ARMA モデルからデータを生成するためのシンプルな関数 arima.sim() があります。たとえば、パラメータ 0.9 の MA(1) から 100 観測を生成する構文は arima.sim(model = list(order = c(0, 0, 1), ma = .9 ), n = 100) です。order = c(0, 0, 0) を使えばホワイトノイズも生成できます。
この演習では、さまざまな ARMA モデルからデータを生成します。各コマンドで観測数は200とし、結果をプロットしてください。
この演習はコースの一部です
Rで学ぶARIMAモデル
演習の手順
arima.sim()とplot()を使って、ホワイトノイズを生成してプロットします。arima.sim()とplot()を使って、パラメータ 0.9 の MA(1) を生成してプロットします。arima.sim()とplot()を使って、パラメータ 1.5 と -0.75 の AR(2) を生成してプロットします。
実践的なインタラクティブ演習
このサンプルコードを完成させて、この演習に挑戦してみましょう。
# Generate and plot white noise
WN <-
# Generate and plot an MA(1) with parameter .9
MA <-
# Generate and plot an AR(2) with parameters 1.5 and -.75
AR <-