Var (X) đã biết, làm thế nào để tính Var (1 / X)?


13

Nếu tôi chỉ có Var(X) , làm sao tôi có thể tính Var(1X)?

Tôi không có bất kỳ thông tin về sự phân bố của X , vì vậy tôi không thể sử dụng chuyển đổi, hoặc bất kỳ phương pháp khác mà sử dụng phân bố xác suất của X .


Tôi nghĩ rằng điều này có thể giúp bạn.
Christoph_J

Câu trả lời:


18

Nó là không thể.

Hãy xem xét một chuỗi Xn của các biến ngẫu nhiên, trong đó

P(Xn=n1)=P(Xn=n+1)=0.5

Sau đó:

Var(Xn)=1for all n

Nhưng tiến đến 0 khi đi đến vô cùng:nVar(1Xn)n

Vmộtr(1Xn)= =(0,5(1n+1-1n-1))2

Ví dụ này sử dụng thực tế là là bất biến dưới các bản dịch của , nhưng thì không.X V a r ( 1Vmộtr(X)XVmộtr(1X)

Nhưng ngay cả khi chúng tôi giả sử , chúng tôi không thể tính toán : HãyV a r ( 1E(X)= =0Vmộtr(1X)

P(Xn= =-1)= =P(Xn= =1)= =0,5(1-1n)

P(Xn= =0)= =1ncho n>0

Sau đó tiếp cận 1 khi chuyển sang vô cùng, nhưng cho tất cả .n V a r ( 1Vmộtr(Xn)nnVmộtr(1Xn)= =n


20

Bạn có thể sử dụng chuỗi Taylor để lấy xấp xỉ các khoảnh khắc thứ tự thấp của một biến ngẫu nhiên được chuyển đổi. Nếu phân phối khá 'chặt chẽ' xung quanh giá trị trung bình (theo một nghĩa cụ thể), thì phép tính gần đúng có thể khá tốt.

Ví dụ

g(X)= =g(μ)+(X-μ)g'(μ)+(X-μ)22g"(μ)+Giáo dục

vì thế

Var[g(X)]=Var[g(μ)+(Xμ)g(μ)+(Xμ)22g(μ)+]=Var[(Xμ)g(μ)+(Xμ)22g(μ)+]=g(μ)2Var[(Xμ)]+2g(μ)Cov[(Xμ),(Xμ)22g(μ)+]+Var[(Xμ)22g(μ)+]

thường chỉ có thuật ngữ đầu tiên được thực hiện

Var[g(X)]g(μ)2Var(X)

Trong trường hợp này (giả sử tôi không mắc lỗi), với , . Var[1g(X)=1XVar[1X]1μ4Var(X)

Wikipedia: Taylor mở rộng cho các khoảnh khắc chức năng của các biến ngẫu nhiên

---

Một số ví dụ để minh họa điều này. Tôi sẽ tạo hai mẫu (phân phối gamma) trong R, một mẫu có phân phối 'không quá chặt chẽ' về ​​giá trị trung bình và một mẫu chặt chẽ hơn một chút.

 a <- rgamma(1000,10,1)  # mean and variance 10; the mean is not many sds from 0
 var(a)
[1] 10.20819  # reasonably close to the population variance

Giá trị gần đúng cho thấy phương sai của phải gần với( 1 / 10 ) 4 × 10 = 0,0011/một(1/10)4×10= =0,001

 var(1/a)
[1] 0.00147171

Tính toán đại số cho thấy phương sai dân số thực tế là1/6480,0054

Bây giờ cho người chặt chẽ hơn:

 a <- rgamma(1000,100,10) # should have mean 10 and variance 1
 var(a)
[1] 1.069147

Giá trị gần đúng cho thấy phương sai của phải gần với( 1 / 10 ) 4 × 1 = 0,00011/một(1/10)4×1= =0,0001

 var(1/a)
[1] 0.0001122586

Tính toán đại số cho thấy phương sai dân số của đối ứng là .102992×980,000104


1
Lưu ý rằng trong trường hợp này, một giả thuyết khá yếu dẫn đến kết luận rằng sẽ không tồn tại trung bình (phương sai whence) cho , nghĩa là gần đúng trong câu trả lời sẽ khá sai lệch. :-) Một giả thuyết ví dụ là có mật độ liên tục trong một khoảng quanh 0 và sao cho . Kết quả sau đó là do mật độ sẽ bị giới hạn từ 0 trong một khoảng thời gian . Tất nhiên, giả thuyết được đưa ra không phải là yếu nhất có thể. X f f ( 0 ) 0 [ - ϵ , ϵ ]1/XXff(0)0[-ε,ε]
Đức hồng y

Lý do chuỗi đối số Taylor sau đó không thành công là vì ẩn phần còn lại (lỗi), trong trường hợp này là và điều này hoạt động không tốt xung quanh . R ( x , μ ) = ( x + μ ) ( x - μ ) 2x = 0
R(x,μ)= =(x+μ)(x-μ)2xμ,
x= =0
Đức hồng y

Người ta thực sự phải cẩn thận về hành vi của mật độ gần 0. Lưu ý rằng trong các ví dụ gamma ở trên, phân bố nghịch đảo là gamma nghịch đảo, trong đó có một giá trị trung bình hữu hạn đòi hỏi ( là tham số hình dạng của gamma chúng ta đang đảo ngược). Hai ví dụ có và . Mặc dù vậy (với các bản phân phối "đẹp" để đảo ngược), việc bỏ qua các thuật ngữ cao hơn có thể đưa ra sự thiên vị đáng chú ý. α>1αα= =10α= =100
Glen_b -Reinstate Monica

điều này có vẻ đúng hướng, của một phân phối bình thường thay đổi qua lại thay vì phân phối bình thường tiêu chuẩn đối ứng: en.wikipedia.org/wiki/
Kẻ
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.