次の2つのステップを実行したい:
- 与えられた時系列に基づいて、AR(1) プロセスを調整したい、つまりパラメータを推定したい。
- 推定されたパラメータに基づいて、AR(1) プロセスをシミュレートしたいと考えています。
これが私のアプローチでした:
set.seed(123)
#Just generate random AR(1) time series; based on this, I want to estimate the parameters
ts_AR <- arima.sim(n=10000, list(ar=c(0.5)))
#1. Estimate parameters with arima()
model_AR <- arima(ts_AR, order=c(1,0,0))
#Looks actually good
model_AR
Series: ts_AR
ARIMA(1,0,0) with non-zero mean
Coefficients:
ar1 intercept
0.4891 -0.0044
s.e. 0.0087 0.0195
sigma^2 estimated as 0.9974: log likelihood=-14176.35
AIC=28358.69 AICc=28358.69 BIC=28380.32
#2. Simulate based on model
arima.sim(model=model_AR, n = 100)
Error in arima.sim(model = model_AR, n = 100) :
'ar' part of model is not stationary
私は時系列の最大の専門家ではありませんが、持続性パラメーターが 1 未満の AR(1) プロセスは定常モデルになるはずだと確信しています。ただし、エラーメッセージは何か違うことを教えてくれます。だから私はここで愚かなことをしますか?もしそうなら、推定されたパラメータに基づいて AR(1) プロセスをシミュレートするために、なぜ、何をすべきか。arima
または、の出力をモデルの入力として に渡すことはできませんarima.sim
か? しかし、そのようなエラーメッセージがどのように表示されるのかわかりません... 「モデル入力を読み取れません。次のようになるはずです...」のようなものを期待します。