Tự động tương quan tiêu cực ở độ trễ 1 và 2 có thể xảy ra không?


8

Chỉ cần làm một vài trò chơi trí tuệ đi qua ghi chú thống kê của tôi ...

Tôi đã thấy các ACF xung quanh có giá trị âm ở độ trễ 1 và 2 - Tôi có thể có một tâm trí trống rỗng ở đây, nhưng sẽ không có AC âm cao ở độ trễ 1 ngụ ý một chuỗi như (-1,1, -1,1, ...), và như vậy, chúng ta sẽ mong đợi AC cũng xen kẽ giữa tích cực và tiêu cực?

Nếu tôi hoàn toàn sai ở đây - có một ví dụ dễ tạo ra khi chúng ta có AC âm mạnh cho cả độ trễ 1 và 2 không?

Cảm ơn bạn!


5
(+1) Điều này liên quan chặt chẽ đến các câu hỏi về tính tích cực của các ma trận hiệp phương sai: xem stats.stackexchange.com/ . Điều kiện đó đặt ra các ràng buộc về mức độ tiêu cực của ACF có thể nhận được đồng thời ở độ trễ 1 và 2.
whuber

Câu trả lời:


3

DGP sau đây, một quá trình MA ( ), có hiện tượng tự tương quan âm ở độ trễ 1 và 2:2

Yt= =10-.5bạnt-1-0,25bạnt-2+bạnt

Đây là một số mã R để mô phỏng DGP và tự mình xem ACF:

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)

nhập mô tả hình ảnh ở đây


1
(+1) Ví dụ hay. Một vài lỗi chính tả? Có nên y<-rep(0,n+j)acfthay vì Acf? Cũng library(stats)chỉ trong trường hợp ai đó không biết.
Matthew Gunn

1
Cảm ơn Matthew, tôi đã thêm các thư viện. Các acf()chức năng là từ các gói dự báo, mà tôi muốn sử dụng vì nó bỏ qua tụt 0 = 1 từ biểu đồ. Tôi đặt atrong y<-rep(a,n+j)dòng cho phù hợp, do đó các giá trị ban đầu của yđã được thiết lập để bình vô điều kiện của mình, nhưng vì không có điều kiện AR, nó không có vấn đề gì giá trị được lưu trữ ở đó trước khi vòng lặp (họ sẽ vượt qua bằng văn bản anyways ).
Elon Plotkin

1
y<-rep(a,n+j)thực sự là một lỗi đánh máy! Ý tôi là nó lày<-rep(alpha,n+j)
Elon Plotkin

1
y1= =1yn= =-(0,5)n-2n2
Khi sử dụng trang web của chúng tôi, bạn xác nhận rằng bạn đã đọc và hiểu Chính sách cookieChính sách bảo mật của chúng tôi.
Licensed under cc by-sa 3.0 with attribution required.