İki bağımsız rastgele değişkenin çarpımı


15

Yaklaşık 1000 değere sahip bir örnek var. Bu veriler adlı iki bağımsız rasgele değişkenin ürününden elde edilir . İlk rastgele değişkenin eşit dağılımı . İkinci rasgele değişkenin dağılımı bilinmemektedir. İkinci ( ) rastgele değişkenin dağılımını nasıl tahmin edebilirim ?ξ U ( 0 , 1 ) ψξ*ψξ~U(0,1)ψ


8
Bu, dekonvolüsyon problemi olarak adlandırılan bir versiyonudur: ürünün günlüğüne geçerseniz, terimlerin birinin dağılımını bildiğinizde toplamın tahmini dağılımını elde edersiniz. Vikipedi üzerinde kontrol edin .
Xi'an

3
Ayrıca , çapraz onaylanmış ilgili soruya bakın : günlük dönüşümünü uyguladıktan sonra sorun eşdeğerdir.
Xi'an,

3
@ Xi'an: Güzel bağlantılar. neredeyse kesinlikle ... umarım ayrıştırarak ve parçaları ayrı ayrı ele alarak bu durumun ölümcül bir şekilde ihlal edilmesinden . ψ = ψ + - ψ -ψ0ψ=ψ+ψ
kardinal

2
@cardinal Bazı veriler negatif olabileceği zaman tahmin sorununun nasıl ele alınacağını merak ediyorum. Ayrışma nasıl belirlenir? (Az veri atama sezgisel metod birden bileşen ve büyük veri için üstel ile evrişim eğiliminde olacaktır çünkü bana başka bakışlar Suboptimal için gelen değerleri açmak için nispeten büyük pozitif gözlemler içine bileşeni. ) Tahmin edicinin eşzamanlı olarak karışımın tanımını ve dekonvolüsyonu ele alması gerektiği gibi görünüyor ve bu da zor görünüyor. 1 ψ -11ψ
whuber

2
@ Açıklama için teşekkürler. Hayır, gürültü değil: logaritma açısından düşündüğüm için, negatif olmadığını unutmuştum . ξ
whuber

Yanıtlar:


3

pozitif gerçek çizgi üzerinde desteği olduğu varsayılarak , Burada ve , verilerin ampirik dağılımıdır. Elde ettiğimiz bu denklemin kaydını alarak,ξψX F n F n

ξψ=X
XFnFn

Log(ξ)+Log(ψ)=Log(X)

Böylece Levy'nin süreklilik teoremi ve ve karakteristik fonksiyonlarını alma bağımsızlığı ile : ψξψ

ΨLog(ξ)(t)ΨLog(ψ)(t)=ΨLog(X)

Şimdi, Böylece, , t h e r e f o r e - L O g ( ξ ) ~ e x s ( 1 ) Ψ L O g ( ξ ) ( - t ) = ( 1 + i t ) - 1ξUnif[0,1],thereforeLog(ξ)Exp(1)

ΨLog(ξ)(t)=(1+it)1

Verilen bu ile rasgele numune .X-1. . . X1000ln(X)Ψln(X)=1nk=11000exp(itXk),X1...X1000ln(X)

Artık 'nin karakteristik fonksiyonu aracılığıyla dağılımını tamamen belirleyebiliriz :Log(ψ)

(1+it)1ΨLog(ψ)(t)=1nk=11000exp(itXk)

nin moment üreten fonksiyonlarının var olduğunu ve olduğunu varsayarsak , yukarıdaki denklemi moment üreten fonksiyonlar açısından yazabiliriz:ln(ψ)t<1

MLog(ψ)(t)=1nk=11000exp(tXk)(1t)

Bu dağılımını elde etmek için moment kavramı ters sonra yeterli Bunun ve böyleceln(ϕ)ϕ


bunu R'deki bir örnekle açıklayabilir misiniz?
Andy

Elbette. Yarın göndermeye çalışacağım.
Drmanifold
Sitemizi kullandığınızda şunları okuyup anladığınızı kabul etmiş olursunuz: Çerez Politikası ve Gizlilik Politikası.
Licensed under cc by-sa 3.0 with attribution required.