Làm thế nào để kiểm tra nếu một ma trận hiệp phương sai không khác không?


11

Nền tảng của nghiên cứu của tôi :

Trong một mẫu Gibbs trong đó chúng ta lấy mẫu (biến sở thích) và từ và , trong đó và là các vectơ ngẫu nhiên -chiều. Chúng tôi biết rằng quá trình này thường được chia thành hai giai đoạn:XYP(X|Y)P(Y|X)XYk

  1. Thời kỳ Burn-in, nơi chúng tôi loại bỏ tất cả các mẫu. Suy ra các mẫu là và .X1XtY1Yt
  2. Khoảng thời gian "Sau khi đốt cháy", trong đó chúng tôi lấy trung bình các mẫu là kết quả mong muốn cuối cùng của chúng tôi.X¯=1ki=1kXt+i

Tuy nhiên, các mẫu trong chuỗi "sau khi ghi" không được phân phối độc lập. Do đó, nếu tôi muốn kiểm tra phương sai của kết quả cuối cùng, nó sẽ trở thànhXt+1Xt+k

Var[X¯]=Var[i=1kXt+i]=1k2(i=1kVar[Xt+i]+i=1k1j=i+1kCov[Xt+i,Xt+j])

Ở đây thuật ngữ là ma trận hiệp phương sai áp dụng cho bất kỳ với .Cov[Xt+i,Xt+j]k×k(i,j)i<j

Ví dụ, tôi có

Xt+1=(1,2,1)Xt+2=(1,0,2)Xt+3=(1,0,0)Xt+4=(5,0,1)

sau đó tôi có thể ước tính ma trận hiệp phương sai vớiCov[Xt+i,Xt+i+1]

13i=13(Xt+iμt+i)(Xt+i+1μt+i+1)

Bây giờ tôi quan tâm đến việc ước tính kết quả là khác không đáng kể vì vậy tôi cần đưa nó vào ước tính phương sai của tôi là .Var[X¯]

Vì vậy, đây là câu hỏi của tôi :

  1. Chúng tôi lấy mẫu từ . Vì đang thay đổi, tôi nghĩ và không cùng phân phối, vì vậy không giống với . Là tuyên bố này đúng?Xt+iP(Xt+i|Yt+i)Yt+iXt+iXt+i+1Cov[Xt+i,Xt+j]Cov[Xt+i,Xt+i]
  2. Giả sử tôi có đủ dữ liệu để ước tính (các mẫu lân cận trong chuỗi), có cách nào để kiểm tra xem ma trận hiệp phương sai có đáng kể không Ma trận khác không? Nói rộng ra, tôi quan tâm đến một chỉ số hướng dẫn tôi một số ma trận hiệp phương sai có ý nghĩa cần được đưa vào ước lượng phương sai cuối cùng của tôi.Cov[Xt+i,Xt+i+1]

4
Trên thực tế, bây giờ điều này có vẻ như là một câu hỏi hay; Tôi nghĩ rằng một số người khác sẽ có vị trí tốt hơn để đưa ra câu trả lời tốt hơn tôi, vì vậy tôi muốn quảng bá điều này (đặt tiền thưởng cho nó) khi nó đủ điều kiện trong thời gian ngắn. [Câu trả lời ngắn: 1. Hai hiệp phương sai này khác nhau. 2. Bạn không cần kiểm tra xem các biến thiên liên tiếp có tương quan hay không (trong tất cả các trường hợp tầm thường nhất; thuật toán hoạt động bằng cách tạo ra các biến thiên phụ thuộc) - thú vị hơn để đo lường mối tương quan so với kiểm tra nó;] ... nếu những câu trả lời hay không xuất hiện Tôi sẽ mở rộng những bình luận ngắn đó thành một câu trả lời đầy đủ
Glen_b -Reinstate Monica

4
Có vẻ như câu hỏi của bạn rộng hơn nhiều so với câu hỏi tiêu đề của bạn. Cụ thể giải quyết câu hỏi tiêu đề của bạn, có bài kiểm tra tính toàn cầu của Bartlett cho phép kiểm tra xem ma trận hiệp phương sai mẫu có chéo không. Bạn có thể cần phải điều chỉnh nó theo kịch bản hiệp phương sai ("ma trận hiệp phương sai" của bạn thực sự không phải là ma trận hiệp phương sai, đó là một ma trận hiệp phương sai, đó là một khối chéo của ma trận hiệp phương sai đầy đủ của cả X_t và X_ { t + 1} với nhau). CC đến @Glen_b.
amip nói rằng Phục hồi lại

2
Tôi sẽ thêm rằng hiệp phương sai có xu hướng phân rã ít nhiều về mặt hình học (ngày càng tăng khi bạn di chuyển xa nhau hơn); các giá trị cách xa nhau về thời gian có xu hướng tương quan rất thấp ( không phải bằng 0 nhưng phần lớn không thể biết được ) trong khi những giá trị gần nhau đôi khi có thể khá phụ thuộc.
Glen_b -Reinstate Monica

1
@Tom 1. Tuy nhiên, với loạt văn phòng phẩm, ở độ trễ rất xa (4 không xa!), Điều gì xảy ra với ACF? 2. Bạn biết điều gì đó về cách các giá trị được tạo từ MCMC hoạt động mà bạn không thể nói về chuỗi thời gian tùy ý ... chúng là Markovian . Bạn sẽ lưu ý rằng những nhận xét trước đây của tôi không cho rằng độ trễ gần nhất phải thể hiện sự phân rã hình học (ví dụ: tôi không nói rằng không thể thấy mối tương quan cao hơn ở độ trễ 4 hơn 3). Bạn vẫn sẽ nhận được (nếu một số điều kiện nhất định) có xu hướng phân rã hình học trong ACF khi bạn di chuyển xa nhau.
Glen_b -Reinstate Monica

2
Nếu thời gian lấy mẫu của bạn quá ngắn, bạn không có ước tính chính xác cao về hiệp phương sai, bạn có thể phải đối phó với thực tế là ước tính của bạn về các điều khoản hiệp phương sai có lỗi tiêu chuẩn rất lớn. Với sự hiểu biết hiện tại của tôi, tôi thậm chí còn mạnh mẽ hơn để khẳng định lại sự phản đối của mình để kiểm tra các mối tương quan. Kiểm tra giả thuyết cho tương quan 0 và không khác không giải quyết vấn đề của bạn ở đây.
Glen_b -Reinstate Monica

Câu trả lời:


1
  1. Chúng tôi lấy mẫu từ . Vì đang thay đổi, tôi nghĩ và không cùng phân phối [...]Xt+iP(Xt+i|Yt+i)Yt+iXt+iXt+i+1

Bạn đang nhầm lẫn phân phối có điều kiện và vô điều kiện ở đây, xem thêm nhận xét tiếp theo của tôi. Có điều kiện trên và , . Nhưng toàn bộ điểm xây dựng sampler Gibbs của bạn là mẫu từ phân phối văn phòng phẩm của và . Nói một cách đại khái, nếu bạn đã chạy chuỗi của mình đủ lâu và để tuân theo phân phối cố định, thì bạn có thể nói có nghĩa là phân phối vô điều kiện của cũng là bất biến. Nói cách khác, nhưYt+i=y1Yt+i+1=y2P(Xt+i|Yt+i=y1)P(Xt+i+1|Yt+i+1=y2)XY{Yt}

P(Xt)=YP(Xt|Yt)dP(Yt),
Xtt và chúng tôi hội tụ đến các bản phân phối cố định, , vì và sẽ được rút ra một cách bất thường từ (cùng!) phân phối văn phòng phẩm . Mặt khác và như trước đây, một khi chúng ta tạo điều kiện cho và , điều này sẽ không giữ được nữa, bất kể lớn đến mức nào .P(Xt+i|Yt+i)=P(Xt+i+1|Yt+i+1)Yt+iYt+i+1P(Yt)Yt+i=y1Yt+i+1=y2t

[...] vì vậy không giống với . Là tuyên bố này đúng?Cov[Xt+i,Xt+j]Cov[Xt+i,Xt+i]

Đúng, điều này là chính xác - mặc dù , tức là và có cùng phân phối văn phòng phẩm. Tôi biết điều này có thể gây nhầm lẫn, nhưng chịu đựng tôi. Xác định với . Bằng cách thay thế lặp đi lặp lại, người ta có thể chỉ ra rằng và vì (các số vô hạn) của các quy tắc vẫn bình thường, nên nó giữ và do đó . Rõ ràng, vàXt+1XtXtXt+1Yt=0.8Yt1+εtεtiidN(0,1)Yt=i=0t0.8iεtiVar(Yt)=i=0t0.82i=110.82YtiidN(0,110.82)YtYt+1vẫn sẽ tương quan, nhưng chúng cũng sẽ đến từ cùng một phân phối ( ). Một tình huống tương tự giữ cho của bạn .Yt+1YtXt

  1. Giả sử tôi có đủ dữ liệu để ước tính (các mẫu lân cận trong chuỗi), có cách nào để kiểm tra xem ma trận hiệp phương sai có đáng kể không Ma trận khác không? Nói rộng ra, tôi quan tâm đến một chỉ số hướng dẫn tôi một số ma trận hiệp phương sai có ý nghĩa cần được đưa vào ước lượng phương sai cuối cùng của tôi.Cov[Xt+i,Xt+i+1]

Chà, nếu bạn có vô số quan sát, tất cả chúng sẽ có ý nghĩa cuối cùng. Rõ ràng, bạn không thể làm điều này trong thực tế, nhưng có nhiều cách 'băm nhỏ' việc mở rộng sau một số điều khoản, xem câu trả lời xuất sắc được chấp nhận ở đây. Về cơ bản, bạn xác định hạt nhân phân rã thành và gán trọng số cho ma trận hiệp phương sai đầu tiên mà bạn có thể tính toán. Nếu bạn muốn chọn theo cách nguyên tắc, bạn sẽ phải đào sâu một chút vào tài liệu, nhưng bài đăng tôi liên kết cung cấp cho bạn một số tài liệu tham khảo tốt để làm chính xác điều đó.k()0lTlT

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.