Jaki stosunek niezależnych dystrybucji daje rozkład normalny?

13

Stosunek dwóch niezależnych rozkładów normalnych daje rozkład Cauchy'ego. Rozkład t jest rozkładem normalnym podzielonym przez niezależny rozkład chi-kwadrat. Stosunek dwóch niezależnych rozkładów chi-kwadrat daje rozkład F.

Szukam współczynnika niezależnych ciągłych rozkładów, który daje normalnie rozłożoną zmienną losową o średniej i wariancji σ 2 ?μσ2)

Prawdopodobnie istnieje nieskończony zestaw możliwych odpowiedzi. Czy możesz podać mi niektóre z tych możliwych odpowiedzi? Byłbym szczególnie wdzięczny, gdyby dwa niezależne rozkłady, których stosunek jest obliczany, są takie same lub przynajmniej mają podobną wariancję.

Remi.b
źródło
2
Chociaż artykuł w Wikipedii o rozkładach proporcji nie zawiera przykładów przypadku, dla którego szukasz, jest to ciekawa lektura.
Avraham
2
Dość szczególnym przypadkiem jest normalna norma, a Y niezależnie ± 1 z prawdopodobieństwem 1XY±1 , a następnieX,YiX12)XY mają tę samą średnią i wariancję iXXY jest zwykle dystrybuowane. XY
Henry
1
Stosunek dwóch niezależnych rozkładów chi-kwadrat daje rozkład F ” --- cóż, niezupełnie. Daje rozkład beta-prime. Aby uzyskać F, musisz przeskalować każdy chi-kwadrat o jego df.
Glen_b
2
Wiele rzeczy sprawia, że ​​wcale nie jestem przekonany, że spełnienie wszystkich twoich warunków jest koniecznie.
Glen_b
1
biorąc przykład generowania metody zmiennych normalnych (np. Box-Muller) jako przykład (która wykorzystuje metodę koła) powiedziałbym, że nie ma stosunków rozkładów jednorodnych, które dają rozkład normalny (przy założeniu, że proszone są rozkłady jednolite)
Nikos M.

Odpowiedzi:

7

Niech gdzieEma rozkład wykładniczy ze średnią2σ2iZ=±1z jednakowym prawdopodobieństwem. NiechY2=1/Y1=Zmimi2)σ2)Z=±1 gdzieBBeta(0,5,0,5). Zakładając, że(Z,E,B)są wzajemnie niezależne, wówczasY1jest niezależne odY2iY1/Y2Normalne(0,σ2). Stąd mamyY2)=1/bbBeta(0,5,0,5)(Z,mi,b)Y1Y2)Y1/Y2)Normalna(0,σ2))

  1. niezależne od Y 2 ;Y1Y2)
  2. Oba ciągłe; takie, że
  3. .Y1/Y2)Normalna(0,σ2))

Nie wymyśliłem, jak uzyskać . Trudniej jest zrozumieć, jak to zrobić, ponieważ problem sprowadza się do znalezienia A i B, które są niezależne od A - B μNormalna(μ,σ2))ZAb , który jest trochę trudniejsze niż wykonywanieA/B~Normalny(0,1)dla niezależnegoAiB.

ZA-bμbNormalna(0,1)
ZA/bNormalna(0,1)ZAb
chłopak
źródło
1
Jeśli to prawda, to jest niesamowite.
Neil G
2
@NeilG to prawda; produktem mojej beta i wykładniczej jest gamma o kształcie 1/2 (z powodu tego, jak można zbudować beta i niezależną gamma za pomocą gamma). Zatem pierwiastek kwadratowy tego jest w połowie normalny, biorąc pod uwagę fakt, że kwadrat normy jest chi-kwadrat.
facet
1
Niedawno zadaliśmy pytanie o iloczyn dwóch zmiennych o rozkładzie normalnym (nie mogę go znaleźć z powrotem). To pytanie zawierało komentarz lub odpowiedź dotyczącą transformacji Boxa-Mullera, która oblicza rozkład normalny (a ściślej dwuwymiarowy rozkład normalny) z iloczynu dwóch transformowanych jednorodnych zmiennych rozkładowych. Ta odpowiedź bardzo się z tym wiąże, ale bierze odwrotność jednej z tych zmiennych w transformacji Boxa-Mullera. cc: @kjetilbhalvorsen
Sextus Empiricus
2

Byłbym szczególnie wdzięczny, gdyby dwa niezależne rozkłady, których stosunek jest obliczany, są takie same 

Jest ma możliwość, że normalna zmienna może być zapisane jako stosunek dwóch zmiennych niezależnych z tego samego rozkładu lub rodziny rozprowadzania (takich jak F-rozprowadzający, który jest stosunkiem dwóch skalowane χ2) zmiennych o rozkładzie lub Cauchy- rozprowadzający, który jest stosunek dwóch normalnych zmiennych rozkładowych ze średnią zerową).

  • Załóżmy, że: dla dowolnej ZA,bfa gdzie fa jest tą samą dystrybucją lub rodziną dystrybucji, mamy

    X=ZAbN.(μ,σ2))

  • Musimy także być w stanie odwrócić ZA i b (jeśli zmienna normalna może być zapisana jako stosunek dwóch niezależnych zmiennych o tym samym rozkładzie lub rodzinie rozkładów, wówczas kolejność można odwrócić)

    1X=bZAN.(μ,σ2))

  • Ale jeśli XN.(μ,σ2)) to X-1N.(μ,σ2)) nie może być prawdziwe (odwrotność normalnej zmiennej rozproszonej nie jest inną normalną zmienną rozkładową).

Szerszy wniosek: jeśli zmienne w dowolnej rodzinie rozkładów faX można zapisać jako stosunek zmiennych w innej rodzinie rozkładów faY to musi być tak, że rodzina faX jest zamknięta przyjmując odwrotność (tj. Dla każdej zmiennej, której rozkład jest w faX rozkład jego wzajemności będzie również wyrażony w faX ).

Np. Odwrotność zmiennej rozproszonej Cauchy'ego jest również rozproszona Cauchy'ego. Odwrotność zmiennej rozproszonej F. jest również rozproszona F.

  • To „jeśli” nie jest „iff”, odwrotność nie jest prawdą. Gdy X i 1/X należą do tej samej rodziny rozkładów, nie zawsze może być możliwe zapisanie jako rozkład stosunku z mianownikiem i mianownikiem z tej samej rodziny rozkładów.

    Przeciwprzykład: możemy sobie wyobrazić rodziny dystrybucji, dla których dla dowolnego X w rodzinie mamy 1/X w tej samej rodzinie, ale nie mamy P.(X=1)=0 . Jest to sprzeczne z faktem, że dla rozkładu stosunku, w którym mianownik i mianownik mają ten sam rozkład, musimy mieć P.(X=1)0 (i coś podobnego można wyrazić dla ciągłych rozkładów, takich jak całka wzdłuż linii X / Y = 1 na wykresie rozrzutu X, Y ma pewną niezerową gęstość, gdy X i Y mają ten sam rozkład i są niezależne).

Sextus Empiricus
źródło
Nie widzę tego Wydaje mi się, że tylko dlatego, że i B / C są normalne, to nie czyni A / DZA/reb/do normalny. ZA/reb/do
Carl
1
Nie rozumiem, w jaki sposób drugie zdanie wynika z pierwszego. Jeśli istnieją pewne takie, że ich iloraz jest normalny, dlaczego wynika z tego, że ich iloraz w innej kolejności również powinien być normalny? Pytanie nie wymagało takiej rodziny dystrybucji, że iloraz wszystkich par elementów jest normalny. ZA,b
Neil G
1
Nie rozumiem co mówisz. Idealnie byłoby, gdyby twoja odpowiedź była spójnym argumentem bez konieczności odczytywania zmian. W tej chwili wydaje się, że twoje drugie stwierdzenie („musimy też mieć”) nie wynika z pierwszego.
Neil G,
1
@kjetilbhalvorsen, jak należy to zmienić? Odpowiedziałem na tę część pytania, która mówi: „Byłbym szczególnie wdzięczny, gdyby dwa niezależne rozkłady, których stosunek jest obliczany, są takie same” . Nie rozumiem, jak odnosi się do tego odpowiedź faceta.
Sextus Empiricus
2
+1: proste, jasne, eleganckie i wnikliwe.
whuber
1

Cóż, tutaj jest jeden, ale go nie udowodnię, pokaż tylko w symulacji.

Dokonać dwóch rozkładów beta z równym dużej parametry kształtu (w tym przypadku n = 40 , 000 ), odejmowanie od 1/2 x -values jednego z nich i nazywają to „licznik”. To daje nam plik PDF, który ma zakres od maksymalnej ( - 1Beta(200,200)n=40,000x, ale ponieważ parametry kształtu są tak duże, nigdy nie osiągamy maksymalnych wartości zakresu. O to histogramn=40,000„licznik” (-12),12))n=40,000wprowadź opis zdjęcia tutaj

Następnie nazywamy drugi mianownik dystrybucji beta „od mianownikiem” bez odejmowania czegokolwiek, więc ma on zwykły zakres dystrybucji beta i jeden z nich wygląda następująco(0,1)

wprowadź opis zdjęcia tutaj

Ponownie, ponieważ kształty są tak duże, nie zbliżamy się do maksymalnego zakresu wartości. Następnie narysujemy licznik ilorazowy w formacie PDF z nałożonym rozkładem normalnym.licznik ułamkamianownik

wprowadź opis zdjęcia tutaj

Teraz w tym przypadku wynik rozkładu normalnego ma i testy normalności, które wyglądają takμ-0,0000204825,σ0,0501789

(StatystycznyWartość pAnderson-Darling0,7997860,481181Baringhaus-Henze1,405850,0852017Cramér-von Mises0,1231450,482844Jarque-Bera ALM4,481030,106404Kołmogorow-Smirnov0,004523280,38635Kuiper0,007980630,109127Połączone Mardia4,481030,106404Mardia Kurtosis1,538490,123929Skośność Mardii2.093990,147879osoba χ2)134,3530,571925Watson U2)0,1313310,211187)

Innymi słowy, nie możemy udowodnić, że stosunek nie jest normalny, nawet bardzo się starając.

Teraz dlaczego? Intuicja z mojej strony, którą mam w nadmiarze. Dowód pozostawiony czytelnikowi, jeśli taki istnieje (może przez limit metody chwil, ale znowu to tylko intuicja).

Beta(20,20)Beta(20,20)-12)tμ-0,000251208,σ0,157665,df33,0402

wprowadź opis zdjęcia tutaj

StatystycznyWartość pAnderson-Darling0,2752620,955502Cramér-von Mises0,03511080,956524Kołmogorow-Smirnov0,003209360,804486Kuiper0,005565010,657146osoba χ2)145,0770,323168Watson U2)0,03510420,878202

N.(0,1)N.(10,1/1000)t μ-0,0000535722,σ0,0992765,df244,154

wprowadź opis zdjęcia tutaj

(StatystycznyWartość pAnderson-Darling0,5016770,745102Cramér-von Mises0,06968240,753515Kołmogorow-Smirnov0,003556880,692225Kuiper0,006083820,501133osoba χ2)142,880,370552Watson U2)0,06032070,590369)
Carl
źródło
6
Jesteś wyraźnie bardzo blisko normalnej dystrybucji. Jednak to wcale nie jest to samo, co normalny rozkład i nie sądzę, aby stosunek wyśrodkowanej symetrycznej beta do zwykłej symetrycznej beta o tych samych parametrach kiedykolwiek był normalny. Byłbym bardzo zainteresowany tym, że się mylę.
Glen_b
3
Twoje rozwiązanie zdecydowanie nie jest normalne. Możesz uogólnić to podejście: weź dowolny rozkład, który jest w przybliżeniu Normalny i podziel go przez rozkład z prawdopodobieństwem skoncentrowanym w pobliżu liczby niezerowej. Wynik (oczywiście) będzie zbliżony do Normalnego - ale nadal nie będzie Normalny. Zastosowanie szeregu testów nie jest przekonujące, ponieważ wszystko to pokazuje, że nie wygenerowano wystarczająco dużych próbek, aby wykazać nienormalność.
whuber
1
108
2
Przejdźmy teraz do sedna sprawy: (1) obalenie normalności jest prostym ćwiczeniem w integralnym przybliżeniu - nie trzeba tutaj podawać szczegółów. Możesz np. Łatwo udowodnić, że 200. moment jest nieskończony. (2) Twoja odpowiedź myli rozkłady z próbkami. Sprzeciwiam się temu fundamentalnemu zamieszaniu; to jest powód, dla którego uważam, że ta odpowiedź jest bardziej myląca niż pomocna. BTW, nie napisałem lekko mojego ostatniego komentarza: wykonałem ten test. Nie zrobiłem tego z superkomputerem, ale z dziesięcioletnią stacją roboczą na PC, a cały proces zajął zaledwie kilka sekund.
whuber
2
@whuber Pozwolę sobie na ostatni punkt. Odpowiedziałem tutaj, aby wzbudzić zainteresowanie, co zrobiłem. Bez mojego zainteresowania ta rozmowa, która była nieaktywna od 2014 r., Nadal pozostanie bez odpowiedzi. Widzę moją rolę jako pobudzanie myśli, a nie bycie doskonałym czy chodzenie po wodzie.
Carl
-3

X1sol,X2)solXγdo

X1solX2)sol=Xγdo

Xγdo1/X1/γdoγ

X1sol=X2)sol/X1/γdo

μμσγ1/γ

chuse
źródło
4
Sprawdź swoją hipotezę, albo przez wyraźne obliczenie współczynnika, albo przez symulację. Albo pokaże, że twoje roszczenie jest nieprawidłowe. Błąd polega na założeniu, że współczynniki dystrybucji można „anulować”, aby „rozwiązać” licznik.
whuber
1
X2)sol