Odwracanie transformacji Fouriera dla rozkładu Fishera

23

Charakterystyczną funkcją rozkładu Fishera jest: gdzie to zbieżna funkcja hipergeometryczna . Ja próbuje rozwiązać odwrotną transformatę Fouriera \ mathcal {f} _ {t, ^ {X} - 1} z n -convolution odzyskania gęstości zmiennej x , to jest: \ mathcal {C {T} _ , x} ^ {- 1} \ left (C (t) ^ n \ right) w celu uzyskania rozkładu sumy nC ( t ) = Γ ( α + 1F(1,α)U

C(t)=Γ(α+12)U(12,1α2,itα)Γ(α2)
U n x F - 1 t , x ( C ( t ) n ) nFt,x1nx
Ft,x1(C(t)n)
nZmienne losowe dystrybuowane przez Fishera. Zastanawiam się, czy ktoś ma jakiś pomysł, ponieważ wydaje się bardzo trudny do rozwiązania. Próbowałem wartości α=3 i n=2 bezskutecznie. Uwaga: dla n=2 przez splot otrzymuję pdf średniej (nie sumy):

3(12(x2+3)(5x23)x2+9(20x4+27x2+9)log(4x23+1)+23(x2+15)(4x2+3)x3tan1(2x3))π2x3(x2+3)3(4x2+3)
,

gdzie x jest średnią z 2 zmiennych. Wiem, że jest nieporęczny, ale chciałbym dowiedzieć się o przybliżeniu rozkładu basenu.

Nero
źródło
czy to pytanie jest żywe?
Brethlosze
1
Tak, wciąż jest otwarty.
Nero
1
Zakładam, że masz pakiet symboliczny, prawda?
Brethlosze

Odpowiedzi:

5

Nie ma gęstej postaci zamkniętej dla splotu statystyki F, więc próba analitycznego odwrócenia funkcji charakterystycznej prawdopodobnie nie doprowadzi do niczego użytecznego.

W statystyce matematycznej pochylone rozwinięcie Edgewortha (znane również jako przybliżenie punktu siodłowego) jest znaną i często stosowaną techniką aproksymacji funkcji gęstości, biorąc pod uwagę funkcję charakterystyczną. Przybliżenie punktu siodłowego, jeśli często niezwykle dokładne. Ole Barndorff-Nielsen i David Cox napisali podręcznik wyjaśniający tę technikę matematyczną.

Istnieją inne sposoby podejścia do problemu bez użycia funkcji charakterystycznej. Można się spodziewać, że rozkład splotu będzie kształtem przypominał rozkład F. Można wypróbować przybliżenie takie jak dla splotu, a następnie wybrać i aby poprawić pierwsze dwa momenty rozkładu. Jest to łatwe, biorąc pod uwagę znaną średnią i wariancję rozkładu F.aF(n,k)nak

αna=nk=α

Gordon Smyth
źródło