ラグ1と2で負の自己相関が発生する可能性はありますか?


8

私の統計ノートを通過するいくつかの心のゲームをやっているだけ...

ラグ1と2で負の値を持つACFを見てきました-ここで空白を持っているかもしれませんが、ラグ1で負のACが高くなると、(-1,1、-1,1、 ...)そのため、ACが正と負を交互に繰り返すことを期待しますか?

ここで私が完全に間違っている場合-ラグ1と2の両方に強い負のACがある簡単な構成例はありますか?

ありがとうございました!


5
(+1)これは、共分散行列の正定性に関する質問と密接に関連しています。stats.stackexchange.com/…を参照してください。その条件は、ACFがラグ1と2で同時に得ることができるマイナスに制約を課します
whuber

回答:


3

次のDGP、MA()プロセスは、ラグ1および2で負の自己相関があります。2

Yt=10.5あなたt1.25あなたt2+あなたt

DGPをシミュレートし、ACFを確認するためのRコードを次に示します。

library(stats)
library(forecast) # for the Acf() function

# number of "observations"
n<-500 
# initialization periods
j<-1000

# choose parameters
alpha<-10
theta<-c(-.5,-.25)
Q<-length(theta)

# generate iid disturbances
u<-rnorm(n+j,0,2)

# define the DGP and generate data series iteratively
y<-rep(alpha,n+j)
for(k in (Q+1):(n+j)){
  y[k]<-alpha + sum(theta*u[k-c(1:Q)]) + u[k] 
}

# get rid of the initialization periods
Y<-y[-c(1:j)]

# confirm the parameters
arima(Y,c(0,0,Q))

#   Call:
#   arima(x = Y, order = c(0, 0, Q))
#   
#   Coefficients:
#             ma1      ma2  intercept
#         -0.4763  -0.2546     9.9979
#   s.e.   0.0448   0.0485     0.0246
#   
#   sigma^2 estimated as 4.124:  log likelihood = -1064.03,  aic = 2134.05

# look at the ACF/PACF
par(mfrow=c(2,1))
Acf(Y)
pacf(Y)

ここに画像の説明を入力してください


1
(+1)良い例。いくつかのタイプミスはありますか?それがあるべきy<-rep(0,n+j)acfの代わりにAcf?またlibrary(stats)、誰かが知らない場合のために。
Matthew Gunn 2017年

1
マシューに感謝、私はライブラリを追加しました。acf()機能は、私はそれ省いグラフから0 = 1が遅れるので、使用したい予測パッケージからのものです。私は入れay<-rep(a,n+j)の初期値がそのようにすることを、一貫性を保つためにラインy(彼らはオーバーとにかく書き込まれますよその無条件の平均値に設定されたが、何のARの用語が存在しないので、値がループの前にそこに格納されているかは重要ではありません)。
Elon Plotkin 2017年

1
y<-rep(a,n+j)確かにタイプミスです!私はそれがそうであることを意味しましたy<-rep(alpha,n+j)
Elon Plotkin 2017年

1
y1=1y=0.522
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.