Konwergencja w dystrybucji \ CLT


9

Biorąc pod uwagę, że , warunkowe rozbieżność. o jest . ma marginalne rozproszenie. Poissona ( ) jest dodatnią stałą.N=nYχ2(2n)Nθθ

Pokaż, że jako , w rozkładzie.θ  (YE(Y))/Var(Y)N(0,1)

Czy ktoś mógłby zasugerować strategie rozwiązania tego problemu. Wygląda na to, że musimy użyć CLT (Central Limit Theorem), ale ciężko jest uzyskać wszelkie informacje na temat na własną rękę. Czy istnieje rv, który można wprowadzić, aby pobrać próbkę i wygenerować ?YY

To zadanie domowe, więc wskazówki są mile widziane.


Dla mnie też wygląda to na clt. Może jest to już dla ciebie oczywiste, ale jako theta-> Infinity, co stanie się z N?
PeterR

Czy powinienem patrzeć na rozkład N? Jeśli się z tym bawię, wygląda na to, że pdf zawsze będzie wynosił 0. Co mogę z tego wywnioskować?
user42102,

co oznacza zmienna losowa poissona (theta)?
PeterR

Zmieszałem N w tym pytaniu z rozmiarem próbki n w definicji CLT. Więc . Widzimy więc, że oczekiwana wartość N zbliża się do nieskończoności. Nie jestem jednak pewien, dokąd się udać. E(N)=θ
user42102,

1
Powinieneś przyjrzeć się niecentralnemu rozkładowi chi-kwadrat. Udowodnienie, że limit jest normalny, będzie bardziej skomplikowane niż prosta aplikacja CLT, której się obawiam.
caburke

Odpowiedzi:


3

rozwiązanie oparte na właściwościach charakterystycznych funkcji, które są zdefiniowane następująco Wiemy, że rozkład jest jednoznacznie zdefiniowany przez funkcję charakterystyczną, więc udowodnię, że a następnie wynika z pożądanej zbieżności.

ψX(t)=Eexp(itX).
ψ(YEY)/Var(Y)ψN(0,1)(t), when θ,

W tym celu będę musiał obliczyć średnią i wariancję , dla której stosuję prawo całkowitych oczekiwań / wariancji - http://en.wikipedia.org/wiki/Law_of_total_expectation . Użyłem, że średnia i wariancja rozkładu Poissona to i średnia i wariancja to i . Teraz pojawia się rachunek z charakterystycznymi funkcjami. Najpierw przepisuję definicję jakoY

EY=E{E(Y|N)}=E{2N}=2θ
Var(Y)=E{Var(Y|N)}+Var{E(Y|N)}=E{4N}+Var(2N)=4θ+4Var(N)=8θ
EN=Var(N)=θχ2n2E(Y|N=n)=2nVar(Y|N=n)=4nY
Y=n=1Z2nI[N=n], where Z2nχ2n2
Teraz używam twierdzenia, które stwierdza, że Charakterystyczna funkcja to , który pochodzi stąd: http://en.wikipedia.org/wiki/Characteristic_function_(probability_theory)
ψY(t)=n=1ψZ2n(t)P(N=n)
χ2n2ψZ2n(t)=(12it)n

Teraz obliczamy funkcję charakterystyczną dla za pomocą rozszerzenia Taylora dla Na koniec korzystamy z właściwości charakterystycznych funkcji Przeskoczyłem nad rachunkiem, ponieważ jest on teraz zbyt długi ...Yexp(x)

ψY(t)=n=1ψZ2n(t)P(N=n)=n=1(12it)nθnn!exp(θ)=n=1(θ(12it))n1n!exp(θ)=exp(θ12it)exp(θ)=exp(2itθ12it)
ψ(YEY)/Var(Y)(t)=exp(iEYVarY)ψY(t/VarY)=exp(t22)exp(1+2it8θ)exp(t22)=ψN(0,1)(t), when θ

1

Można to wykazać poprzez związek z niecentralnym rozkładem chisquared. Jest dobry artykuł w Wikipedii na ten temat, do którego będę się swobodnie odwoływał! https://en.wikipedia.org/wiki/Noncentral_chi-squared_distribution

Podałeś, że jest rozkładem chisquare z stopniami swobody, dla . Tutaj ma rozkład Poissona z oczekiwaniem .Y|N=n2nn=0,1,,Nθ

Mamy zatem, że funkcję gęstości (bezwarunkowo) można zapisać, stosując zasadę całkowitego prawdopodobieństwa, jako Która jest prawie gęstością niecentralnej zmiennej chisquared, z wyjątkiem parametru stopnia swobody , co jest naprawdę niezdefiniowane. (podano to w sekcji definicji artykułu w Wikipedii).Y

fY(y;0,θ)=i=0eθθii!fχ22i(y)
k=0

Aby więc uzyskać coś dobrze zdefiniowanego, zastępujemy powyższą formułę która jest gęstością niecentralnej zmiennej chisquared z stopni swobody i parametrem niecentralności . Tak więc w naszej analizie musimy pamiętać o przekroczeniu limitu, gdy po wzięciu limitu . Jest to bezproblemowe, ponieważ na granicy prawdopodobieństwo

fY(y;k,θ)=i=0eθθii!fχ22i+k(y)
k2θk0θθN=0 idzie do zera, więc masa punktowa przy zeru znika (zmienna krzyżowa o zerowych stopniach swobody musi być interpretowana jako masa punktowa przy zeru, więc nie mają funkcji gęstości).

Teraz, dla każdego ustalonego , użyj wyniku w wiki, rozkładach związanych z sekcjami, normalnych przybliżeniach, co daje poszukiwany standardowy normalny limit dla każdego . Następnie weź limit, gdy osiągnie zero, co daje wynik. kkk

Korzystając z naszej strony potwierdzasz, że przeczytałeś(-aś) i rozumiesz nasze zasady używania plików cookie i zasady ochrony prywatności.
Licensed under cc by-sa 3.0 with attribution required.