Dlaczego w funkcji gęstości rozkładu beta występuje -1?

18

Dystrybucja beta pojawia się w dwóch parametryzacjach (lub tutaj )

(1)f(x)xα(1x)β

lub ten, który wydaje się być używany częściej

(2)f(x)xα1(1x)β1

Ale dlaczego dokładnie jest „ ” w drugiej formule?1

Pierwsze sformułowanie wydaje się intuicyjnie bardziej bezpośrednio odpowiadać rozkładowi dwumianowemu

(3)g(k)pk(1p)nk

ale „widziana” z perspektywy „sp . Jest to szczególnie wyraźne w modelu beta-dwumianowym, w którym można rozumieć jako wcześniejszą liczbę sukcesów, a to wcześniejszą liczbę niepowodzeń.αβ

Dlaczego więc druga forma zyskała popularność i jakie jest jej uzasadnienie ? Jakie są konsekwencje zastosowania dowolnej z parametryzacji (np. W połączeniu z rozkładem dwumianowym)?

Byłoby wspaniale, gdyby ktoś mógł dodatkowo wskazać pochodzenie takiego wyboru i wstępne argumenty przemawiające za nim, ale nie jest to dla mnie koniecznością.

Tim
źródło
3
W tej odpowiedzi podano głęboki powód : równa się stosunku do miary . To ogranicza twoje pytanie do „ dlaczego ta konkretna miara? Uznanie, że miarą jest sugeruje „właściwy” sposób zrozumienia tych rozkładów jest zastosowanie transformacji logistycznej: wówczas warunki „ ” znikną.x α ( 1 - x ) β d μ = d x / ( ( x ( 1 - x ) ) d μ = d ( log ( xfxα(1x)βdμ=dx/((x(1x))-1
dμ=d(log(x1x))
1
whuber
1
Myślę, że faktyczny powód, dla którego się to wydarzyło, jest historyczny - ponieważ pojawia się w ten sposób w funkcji beta, od której pochodzi nazwa dystrybucji. Jeśli chodzi o to , dlaczego ma to w potędze, spodziewam się, że ostatecznie będzie to związane z powodem, o którym wspomina whuber (choć historycznie nie ma to nic wspólnego z miarą, a nawet prawdopodobieństwem). 1
Glen_b
2
@Glen_b To coś więcej niż historyczne: istnieją głębokie powody. Wynika to z intymnego połączenia funkcji Beta i Gamma, redukując pytanie, dlaczego wykładnik w jest a nie . Jest tak, ponieważ jest sumą Gaussa . Odpowiednio, „słuszne” jest postrzeganie jako całki multiplikatywnego homomorfizmu razy addytywnego znaku stosunku do miary Haar w grupie . s - 1 s Γ Γ t t s t e - t d t / t R ×Γ(s)=0ts1etdts1sΓΓttstetdt/tR×
Whuber
1
@wh To dobry powód, dla którego funkcja gamma powinna być wybrana właśnie w ten sposób (i już zasugerowałem, że taki powód istniał powyżej i akceptuję jakąś formę rozumowania podobną do tego - ale koniecznie z innym formalizmem - do wyboru Eulera); wraz z gęstością występują odpowiednio ważne powody; ale to nie dowodzi, że tak naprawdę był to powód wyboru (dlaczego forma została wybrana w obecnej formie), tylko że jest to dobry powód, aby to zrobić. Forma funkcji gamma ... ctd
Glen_b -Reinstate Monica
1
Sam ctd ... może z łatwością być wystarczającym powodem, aby wybrać tę formę dla gęstości i aby inni poszli w jej ślady. [Często dokonuje się wyborów z prostszych powodów, które możemy później zidentyfikować, a następnie często potrzeba istotnych powodów, aby zrobić cokolwiek innego. Czy wiemy, że to właśnie dlatego początkowo wybrane] - można wyjaśnić jasno, że istnieje powód, dlaczego należy wybrać gęstość być w ten sposób, a nie dlatego, że jest w ten sposób. Wiąże się to z sekwencją ludzi dokonujących wyborów (korzystania z nich w ten sposób i podążania za nimi) oraz ich powodów w wybranym przez siebie momencie.
Glen_b

Odpowiedzi:

9

To opowieść o stopniach swobody i parametrach statystycznych oraz o tym, dlaczego miło jest, że oba mają bezpośrednie proste połączenie.

Historycznie terminy „ ” pojawiły się w badaniach Eulera dotyczących funkcji Beta. Używał tej parametryzacji do 1763 roku, podobnie jak Adrien-Marie Legendre: ich użycie ustanowiło kolejną konwencję matematyczną. Ta praca poprzedza wszystkie znane zastosowania statystyczne.1

Współczesna teoria matematyczna dostarcza wielu wskazówek, poprzez bogactwo zastosowań w analizie, teorii liczb i geometrii, że terminy „ ” faktycznie mają pewne znaczenie. Naszkicowałem niektóre z tych powodów w komentarzach do pytania.1

Bardziej interesujące jest to, jaka powinna być „właściwa” parametryzacja statystyczna. To nie jest tak jasne i nie musi być takie samo jak konwencja matematyczna. Istnieje ogromna sieć powszechnie używanych, znanych, powiązanych ze sobą rodzin rozkładów prawdopodobieństwa. Zatem konwencje stosowane do nazwania (to znaczy parametryzacji) jednej rodziny zwykle implikują powiązane konwencje do rodzin związanych z nazwami. Zmień jedną parametryzację, a będziesz chciał je wszystkie zmienić. Możemy zatem spojrzeć na te relacje w poszukiwaniu wskazówek.

Mało kto nie zgadza się, że najważniejsze rodziny dystrybucyjne wywodzą się z rodziny Normalnej. Przypomnijmy, że zmienna losowa jest „rozkładem normalnym”, gdy ma gęstość prawdopodobieństwa proporcjonalną do . Gdy i , mówi się, że ma standardowy rozkład normalny.( X - μ ) / σ f ( x ) exp ( - x 2 / 2 ) σ = 1 μ = 0 XX(Xμ)/σf(x)exp(x2/2)σ=1μ=0X

Wiele zestawów danych jest badanych przy użyciu stosunkowo prostych statystyk obejmujących racjonalne kombinacje danych i niskie moce (zwykle kwadraty). Gdy dane te są modelowane jako losowe próbki z rozkładu normalnego - tak, że każdy jest postrzegany jako realizacja zmiennej normalnej , wszystkie mają wspólny rozkład i są niezależne - rozkłady tych statystyk są określone przez ten rozkład normalny. Najczęściej pojawiają się w praktycex i X i X ix1,x2,,xnxiXiXi

  1. t ν = n - 1 t = ˉ Xtν , rozkład Studentat z „stopniami swobody”. Jest to rozkład statystyki nazwa gdzie modeluje średnią danych i jest standardowym błędem średniej. Dzielenie przez pokazuje, że musi wynosić lub więcej, skąd jest liczbą całkowitąν=n1ˉ X =(X1+X2++Xn)/nse(X)=(1/

    t=X¯se(X)
    X¯=(X1+X2++Xn)/n n-1n2ν1se(X)=(1/n)(X12+X22++Xn2)/(n1)X¯2n1n2ν1albo lepszy. Formuła, choć pozornie trochę skomplikowana, jest pierwiastkiem kwadratowym z racjonalnej funkcji danych stopnia drugiego: jest względnie prosta.
  2. χ 2 ν ν χ 2 1 / ν χ 2χν2 , (chi-kwadrat) rozkładu z "stopni swobody" (DF). Jest to rozkład sumy kwadratów niezależnych standardowych zmiennych normalnych . Rozkład średniej kwadratów tych zmiennych będzie zatem skalowanym przez : będę to określać jako „znormalizowany” .χ2ννχ21/νχ2

  3. F ( ν 1 , ν 2 ) χ 2 ν 1 ν 2Fν1,ν2 , Współczynnik podziału parametrów to stosunek dwóch niezależnych znormalizowane rozkładów i stopni swobody.F(ν1,ν2)χ2ν1ν2

Obliczenia matematyczne pokazują, że wszystkie trzy z tych rozkładów mają gęstości. Co ważne, gęstość jest proporcjonalna do całki w całkowej definicji Eulera funkcji Gamma ( ). Porównajmy je: Γχν2Γ

fχν2(2x)xν/21ex;fΓ(ν)(x)xν1ex.

To pokazuje, że dwukrotnie zmienna ma rozkład gamma z parametrem . Współczynnik połowy jest dość uciążliwy, ale odjęcie znacznie pogorszy związek. To już dostarcza przekonującej odpowiedzi na pytanie: jeśli chcemy, aby parametr zliczał liczbę kwadratowych zmiennych normalnych, które go wytwarzają (do ), to wykładnik gęstości funkcja musi być o połowę mniejsza niż ta liczba. v / 2 1 χ 2 1 / 2χν2ν/21χ21/2

Dlaczego współczynnik mniej kłopotliwy niż różnica ? Powodem jest to, że czynnik pozostanie spójny, gdy dodamy rzeczy. Jeśli suma kwadratów niezależnych normalnych norm jest proporcjonalna do rozkładu gamma z parametrem (razy jakiś czynnik), wówczas suma kwadratów niezależnych normalnych norm jest proporcjonalna do rozkładu gamma z parametrem (razy ten sam współczynnik) , skąd suma kwadratów wszystkich zmiennych jest proporcjonalna do rozkładu gamma z parametrem (wciąż razy ten sam współczynnik). 1 n n m m n + m m + n1/21nnmmn+mm+nBardzo pomocne jest to, że tak ścisłe dodanie parametrów naśladuje dodanie zliczeń.

Gdybyśmy jednak usunęli te matematycznie wyglądające „ ” ze wzorów matematycznych, te miłe relacje stałyby się bardziej skomplikowane. Na przykład, jeśli zmieniliśmy parametryzację rozkładów gamma tak, aby odnosiła się do rzeczywistej mocy we wzorze, tak że byłby powiązany z rozkładem „Gamma ” (ponieważ moc w jego PDF to ), wówczas suma trzech rozkładów musiałaby być nazwana rozkładem „Gamma ”. Krótko mówiąc, bliski związek addytywny między stopniami swobody a parametrem w rozkładach gamma zostałby utracony przez usunięciex χ 2 1 ( 0 ) x 1 - 1 = 0 χ 2 1 ( 2 ) - 11xχ12(0)x11=0χ12(2)1 ze wzoru i wchłaniając go w parametrze.

Podobnie funkcja prawdopodobieństwa rozkładu współczynnika jest ściśle związana z rozkładami Beta. Rzeczywiście, gdy ma rozkład współczynnika , rozkład ma Beta . Jego funkcja gęstości jest proporcjonalna doY C Z = ν 1 T / ( ν 1 Y + ν 2 ) ( ν 1 / 2 , ν 2 / 2 )FYFZ=ν1Y/(ν1Y+ν2)(ν1/2,ν2/2)

fZ(z)zν1/21(1z)ν2/21.

Ponadto - przyjmując te idee do pełnego koła - kwadrat rozkładu Studenta z df ma rozkład współczynnika z parametrami . Jeszcze raz widać, że utrzymanie konwencjonalnej parametryzacji utrzymuje wyraźny związek z podstawowymi liczbami, które przyczyniają się do stopni swobody.ν F ( 1 , ν )tνF(1,ν)

Z statystycznego punktu widzenia byłoby więc najbardziej naturalne i najprostsze zastosowanie odmiany konwencjonalnych parametryzacji matematycznych rozkładów i Beta: powinniśmy preferować nazywanie rozkładu a „ dystrybucja ”i dystrybucja Beta należy nazwać„ dystrybucją Beta ”. W rzeczywistości już to zrobiliśmy: właśnie dlatego nadal używamy nazw „rozkład chi-kwadrat” i „ współczynnik ” zamiast „gamma” i „beta”. Niezależnie od tego w żadnym wypadku nie chcielibyśmy usunąć „ΓΓ ( 2 α )Γ(α)Γ(2α)( 2 α , 2 β ) F - 1(α,β)(2α,2β)F1„terminy, które pojawiają się we wzorach matematycznych dla ich gęstości. Gdybyśmy to zrobili, stracilibyśmy bezpośredni związek między parametrami w gęstościach i liczeniem danych, z którymi są one powiązane: zawsze bylibyśmy wyłączeni o jeden.

Whuber
źródło
1
Dzięki za odpowiedź (już I + 1d). Mam tylko małe pytanie uzupełniające: może coś mi brakuje, ale czy nie poświęcamy bezpośredniej relacji z dwumianem za pomocą parametryzacji -1?
Tim
Nie jestem pewien, do którego „bezpośredniego związku z dwumianem” masz na myśli, Tim. Na przykład, gdy Beta dystrybucja służy jako sprzężoną przed dla próbki dwumianowego, wyraźnie parametry są dokładnie te właściwe do użycia: dodać (nie ) do liczby sukcesów i (nie ) do liczby awarii. a a - 1 b b - 1(a,b)aa1bb1
whuber
1

Notacja wprowadza Cię w błąd. Istnieje „ukryta ” w formule , ponieważ w , i musi być większa niż (drugi człon, który podałeś w swoim pytaniu mówi to wyraźnie). „S ” te same parametry, y w obu wzorach nie są; mają różne zakresy: w , oraz w , . Te zakresy dla i( 1 ) ( 1 ) α β - 1 α β ( 1 ) α , β > - 1 ( 2 ) α , β > 0 α β (1(1)(1)αβ1αβ(1)α,β>1(2)α,β>0αβsą niezbędne, aby zagwarantować, że całka gęstości nie będzie się różnić. Aby to zobaczyć, weź pod uwagę w przypadek (lub mniej) i , a następnie spróbuj zintegrować gęstość (jądro) między a . Odpowiednio, wypróbuj to samo w dla (lub mniej) i .α = - 1 β = 0 0 1 ( 2 ) α = 0 β = 1(1)α=1β=001(2)α=0β=1

Zen
źródło
2
Kwestia zakresu definicji dla i wydaje się znikać, gdy interpretowana jest całka, jak zrobił to Pochhammer w 1890 r., Jako specyficzna całka konturu. W takim przypadku można go przyrównać do wyrażenia, które określa funkcję analityczną dla wszystkich wartości i wszystkich złożonych. To rzuca światło na wątpliwości w pytaniu: dlaczego dokładnie przyjęto tę specyficzną parametryzację, skoro istnieje wiele innych możliwych parametryzacji, które wydają się być równie dobre? β α βαβαβ
whuber
1
Dla mnie wątpliwości PO wydają się znacznie bardziej podstawowe. Jest trochę zdezorientowany co do „-1” w (2), ale nie w (1) (oczywiście nieprawda). Wygląda na to, że twój komentarz odpowiada na inne pytanie (nawiasem mówiąc, o wiele bardziej interesujące).
Zen
2
Dziękuję za twój wysiłek i odpowiedź, ale wciąż nie odpowiada na moją główną obawę: dlaczego wybrano -1? Zgodnie z logiką, w zasadzie można wybrać dowolną wartość, zmieniając dowolną dolną granicę na coś innego. Nie rozumiem, dlaczego -1 lub 0 może być lepszą lub gorszą dolną granicą dla wartości parametrów, poza tym, że 0 jest „estetycznie” ładniejsze. Z drugiej strony Beta (0, 0) byłby miły „domyślny” dla równomiernego rozkładu przy użyciu pierwszego formularza. Tak, są to bardzo subiektywne komentarze, ale to moja główna uwaga: czy są jakieś niearbitralne powody takiego wyboru?
Tim
1
Zen, zgadzam się, że było pytanie, jak interpretować oryginalny post. Tim, dziękuję za wyjaśnienia.
whuber
1
Cześć Tim! Nie widzę żadnego definitywnego powodu, chociaż czyni to bardziej bezpośrednim związek z faktem, że dla , jeśli i są niezależne, następnie to , a gęstość jest proporcjonalna do . Ale wtedy możesz zakwestionować parametryzację rozkładu gamma ...U G a m m a ( α , 1 ) V G a m m aα,β>0UGamma(α,1)X = U / ( U + V ) B e t a ( α , β ) X x α - 1 ( 1 - x ) β - 1VGamma(β,1)X=U/(U+V)Beta(α,β)Xxα1(1x)β1
Zen
0

Dla mnie istnienie -1 w wykładniku wiąże się z rozwojem funkcji gamma. Motywacją funkcji Gamma jest znalezienie gładkiej krzywej do połączenia punktów silni. Ponieważ nie można obliczyćbezpośrednio, jeśli nie jest liczbą całkowitą, pomysłem było znalezienie funkcji dla dowolnego która spełnia relację powtarzalności zdefiniowaną przez silnię, a mianowiciex!x!xx0

f(1)=1f(x+1)=xf(x).

Rozwiązanie polegało na zbieżności całki. Dla funkcji zdefiniowanej jako

f(x+1)=0txexdt,

integracja przez części zapewnia:

f(x+1)=0txexdt=[txex]0+0xtx1exdt=limx(txex)0e0+x0tx1exdt=00+x0tx1exdt=xf(x).

Tak więc powyższa funkcja spełnia tę właściwość, a -1 w wykładniku wywodzi się z procedury całkowania przez części. Zobacz artykuł w Wikipedii https://en.wikipedia.org/wiki/Gamma_function .

Edycja: przepraszam, jeśli mój post nie jest w pełni jasny; Próbuję tylko wskazać, że moim zdaniem istnienie -1 w rozkładzie beta wynika z uogólnienia silni za pomocą funkcji gamma. Istnieją dwa warunki: oraz . Mamy, dlatego spełnia. Ponadto mamy . Jeśli chodzi o rozkład beta z parametrami , uogólnienie współczynnika dwumianowego tof(1)=1f(x+1)=xf(x)Γ(x)=(x1)!Γ(x+1)=xΓ(x)=x(x1)!=x!Γ(1)=(11)!=0!=1α,βΓ(α+β)Γ(α)Γ(β)=(α+β1)!(α1)!(β1)!. Tam mamy -1 w mianowniku dla obu parametrów.

teatr
źródło
Nie ma to sensu, ponieważ funkcja rekurencji spełniana przez silnię nie jest tym, co podajesz:(x+1)!xx!.
whuber
Funkcją spełniającą relację powtarzalności jest Gamma: . Tak to jest zdefiniowane. Γ ( x + 1 ) = x Γ ( x )f(x)Γ(x+1)=xΓ(x)
aatr
Tak: ale podana motywacja opiera się na funkcji czynnikowej , a nie na gamma.
whuber
Ważne jest, aby przypomnieć związek między gamma i silnią:. Γ(x)=(x1)!
aatr
Niestety, jest to logika kołowa: zaczynasz od silni, charakteryzujesz Gamma jako interpolując ją, a następnie dochodzisz do wniosku, że istnieje -1. W rzeczywistości twój post pokazuje -1, jakby wypadł przypadkowo, myląc Gamma z silnią. Niewielu uzna to za pouczające lub przekonujące.
whuber