Rozkład splotu kwadratowych zmiennych normalnych i chi-kwadratowych?


14

następujący problem pojawił się ostatnio podczas analizy danych. Jeśli zmienna losowa X podąża za rozkładem normalnym, a Y za (z n dof), jak rozkłada się ? Do tej pory wymyśliłem pdf : χn2Z=X2+Y2Y2

ψn2(x)=F(x)x=(0xtn/21et/22n/2Γ(n/2)dt)x=12n/2Γ(n/2)(x)n/21ex/2(x)x=12n/21Γ(n/2)xn/41ex/2

a także pewne uproszczenia dla całki splotowej ( ma pdf z m dof):χ 2 mX2χm2

Kmn(t):=(χm2ψn2)(t)=0tχm2(x)ψn2(tx)dx=(2(n+m)2+1Γ(m2)Γ(n2))10t(tx)n41xm21exp((tx+x)/2)dx

Czy ktoś widzi dobry sposób obliczenia tej całki dla dowolnego rzeczywistego t, czy też trzeba ją obliczyć numerycznie? A może brakuje mi dużo prostszego rozwiązania?


1
Gdyby nie było podniesione do kwadratu, miałbym jakąś konkretną radę. Nie sądzę, aby ten był wykonalny (niekoniecznie szczególnie pouczający, nawet gdyby okazał się wykonalny). Kusiłoby mnie, aby przyjrzeć się podejściom obliczeniowym, takim jak splot numeryczny lub symulacja, w zależności od tego, co dokładnie chcesz zrobić z wynikiem. Y
Glen_b

Moim zdaniem jest bardzo mało prawdopodobne, aby można było wykonać całkę.
Dave31415

1
@ Dave31415 Dla i nawet całka może być jawnie obliczane dla dodatnich wartości zintegrowanego z a . Będzie równa liniowej kombinacji funkcji wykładniczej i funkcji błędu ze współczynnikami, które są wielomianami w . Oceny można dokonać przez podstawienie . Na przykład przy otrzymujemy . m n m nmnm x=t-u2n=2,m=4 1tx=tu2n=2,m=414e18(2t+1)2(et2(2π(4t+3)(erfi(2t122)+erfi(122))+4e18)4et2+18(2t+1))
Whuber

Ładny. W przypadku liczb nieparzystych można to prawdopodobnie przybliżyć w przybliżeniu do średniej wyniku dla liczb parzystych ograniczających? Albo może nie.
Dave31415

1
Dziękuję za odpowiedzi! W niektórych parzystych przypadkach uzyskałem podobny wynik dotyczący funkcji Dawsona, ale wygląda na to, że będę musiał popracować nad ogólnym rozwiązaniem ...
Leo Szilard

Odpowiedzi:


2

Jeśli to pomaga, zmienna jest uogólnioną zmienną losową gamma (patrz np. Stacy 1962). Twoje pytanie dotyczy rozkładu sumy losowej zmiennej chi-kwadrat i uogólnionej zmiennej losowej gamma. Według mojej wiedzy gęstość wynikowej zmiennej nie ma wyrażenia w postaci zamkniętej. Dlatego uzyskany splot jest całką bez rozwiązania w postaci zamkniętej. Myślę, że utkniesz w numerycznym rozwiązaniu tego problemu.Y2


Stacy, EW (1962). Uogólnienie rozkładu gamma. Annals of Mathematical Statistics 33 (3) , s. 1187–1192.


0

To tylko wskazówka. Pearson typu III może mieć kwadrat chi. Czasami splot można znaleźć, splatając coś ze sobą. Udało mi się to zrobić dla zwojów ND i GD , dla których sam z sobą splotłem Pearson III. Jak to działa z ND i Chi-Squared, nie jestem pewien. Ale poprosiłeś o podpowiedzi, a to jest ogólna wskazówka. Mam nadzieję, że to powinno wystarczyć.2


2
Czy możesz wyjaśnić, jak to odpowiada na pytanie? Nie wydaje się to bezpośrednio powiązane.
whuber

Można dokonać splotu Pearsona typu III z samym sobą. Z jakiegoś powodu splot jednej rzeczy ze sobą jest łatwiejszy do rozwiązania niż splot jednej rzeczy z drugą. Na przykład rozwiązałem splot Pearsona typu III i uzyskałem splot ND z GD, pokrewnym problemem.
Carl

Wygląda na to, że nie pomogło, wkrótce zostanie usunięte.
Carl
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.