平均がゼロではない AR(1) 時系列をシミュレートする正しい方法が見つからないようです。53 個のデータ ポイントが必要です。rho = .8、平均 = 300 です。
ただし、arima.sim(list(order=c(1,0,0), ar=.8), n=53, mean=300, sd=21)
1500年代の値を教えてください。例えば：
1480.099 1480.518 1501.794 1509.464 1499.965 1489.545 1482.367 1505.103 (など)
私も試し arima.sim(n=52, model=list(ar=c(.8)), start.innov=300, n.start=1)
ましたが、次のようにカウントダウンします。
238.81775870 190.19203239 151.91292491 122.09682547 96.27074057 [6] 77.17105923 63.15148491 50.04211711 39.68465916 32.46837830 24.78357345 21.27437183 15.93486092 13.40199333 10.99762449 8.70208879 5.62264196 3.15086491 2.13809323 1.30009732
arima.sim(list(order=c(1,0,0), ar=.8), n=53,sd=21) + 300そして、どれが正しい答えを与えるように見えるかを試しました。例えば：
280.6420 247.3219 292.4309 289.8923 261.5347 279.6198 290.6622 295.0501 264.4233 273.8532 261.9590 278.0217 300.6825 291.4469 291.5964 293.5710 285.0330 274.5732 285.2396 298.0211 319.9195 324.0424 342.2192 353.8149 and so on..
しかし、私はこれが正しいことをしているのか疑問に思っていますか? それはまだ正しい数で自動相関していますか？