Hiệp phương sai biến ngẫu nhiên


12

Tôi có hai biến ngẫu nhiên X>0Y>0 .

Cho rằng tôi có thể ước tính

Cov(X,Y),
làm thế nào tôi có thể ước tính
Cov(log(X),log(Y))?

3
Câu hỏi trước đây hỏi về tương quan thay vì hiệp phương sai, nhưng nó có liên quan: stats.stackexchange.com/questions/35941/ Kẻ
Douglas Zare

Câu trả lời:


16

Người ta có thể thực hiện phương pháp mở rộng Taylor:

http://en.wikipedia.org/wiki/Taylor_Exansions_for_the_moments_of_fiances_of_random_variables

Biên tập:

Lấy , V = log ( Y ) .U=log(X)V=log(Y)

Sử dụng khai triển Taylor đa biến để tính xấp xỉ (theo cách tương tự với ví dụ ở cuối "Khoảnh khắc đầu tiên" trong liên kết thực hiện trường hợp đơn giản hơn của E ( X .1 / Y ) ) và sử dụng mở rộng đơn biến để tính các xấp xỉ cho E ( U )E ( V ) (như được đưa ra trong phần đầu tiên của cùng một phần) với độ chính xác tương tự. Từ những điều đó, tính hiệp phương sai (gần đúng).E(UV)E(X.1/Y))E(U)E(V)

Mở rộng ra mức độ gần đúng tương tự như ví dụ trong liên kết, tôi nghĩ rằng bạn kết thúc với các thuật ngữ về giá trị trung bình và phương sai của từng biến (chưa được dịch) và hiệp phương sai của chúng.

Chỉnh sửa 2:

Nhưng đây là một mẹo nhỏ có thể tiết kiệm một số nỗ lực:

Lưu ý rằng X = exp ( U )Y = exp ( V ) .E(XY)=Cov(X,Y)+E(X)E(Y)X=exp(U)Y=exp(V)

Với ta có E(exp(U))exp(μU)+exp(μU)

E[f(X)]f(μX)+f(μX)2σX2
E(exp(U))exp(μU)+exp(μU)2σU2exp(μU+12σU2)

Chỉnh sửa: Đó là bước cuối cùng sau từ Taylor xấp xỉ , đó là tốt cho nhỏ b (có tính b = 1exp(b)1+bb ).b=12σU2

(xấp xỉ đó là chính xác cho , V bình thường: E ( exp ( U ) ) = exp ( μ U + 1UV)E(exp(U))=exp(μU+12σU2)

Đặt W=U+V

E(XY)=E(exp(U).exp(V))=E(exp(W))

exp(μW)+exp(μW)2σW2exp(μW+12σW2)

và đưa ra , sau đóVar(W)=Var(U)+Var(V)+2Cov(U,V)

(Biên tập:)

exp(μW+ 1

1+Cov(X,Y)E(X)E(Y)=E(XY)E(X)E(Y)
exp(μU+μV+ 1
exp(μW+12σW2)exp(μU+12σU2).exp(μV+12σV2)
exp[Cov(U,V)]
exp(μU+μV+12(σU2+σV2+2Cov(U,V)))exp(μU+12σU2).exp(μV+12σV2)
exp[Cov(U,V)]

Do đó . Điều này phải chính xác chogaussianU,Vbivariate.Cov(U,V)log(1+Cov(X,Y)E(X)E(Y))U,V

Nếu bạn đã sử dụng xấp xỉ đầu tiên chứ không phải thứ hai, bạn sẽ nhận được một xấp xỉ khác ở đây.


Bạn có thể cho biết thêm một chút chi tiết xin vui lòng? Dù sao, thx cho đề xuất
user7064

Chỉnh sửa để biết chi tiết.
Glen_b -Reinstate Monica

Cảm ơn @Glend_b. Tôi sẽ chấp nhận khi chi tiết sẽ được thêm vào. Trong khi đó, +1 :-)
user7064

Đừng lo lắng; Lúc đó tôi đang bận, rồi hoàn toàn quên mất. Bây giờ đã sửa
Glen_b -Reinstate Monica

UVXY là nhỏ).
Glen_b -Reinstate Monica

8

XYCov(X,Y)XYCov(log(X),log(Y))log(X)log(Y)

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.