Welches Verhältnis der unabhängigen Verteilungen ergibt eine Normalverteilung?


12

Das Verhältnis zweier unabhängiger Normalverteilungen ergibt eine Cauchy-Verteilung. Die t-Verteilung ist eine Normalverteilung geteilt durch eine unabhängige Chi-Quadrat-Verteilung. Das Verhältnis von zwei unabhängigen Chi-Quadrat-Verteilungen ergibt eine F-Verteilung.

Ich suche ein Verhältnis unabhängiger stetiger Verteilungen, das eine normalverteilte Zufallsvariable mit mittlerem μ und Varianz ergibt ?σ2

Es gibt wahrscheinlich unendlich viele mögliche Antworten. Können Sie mir einige dieser möglichen Antworten geben? Ich würde es besonders begrüßen, wenn die beiden unabhängigen Verteilungen, deren Verhältnis berechnet wird, gleich sind oder zumindest eine ähnliche Varianz haben.


2
Der Wikipedia-Artikel zu Ratio Disributions enthält zwar keine Beispiele für den von Ihnen gesuchten Fall, ist aber eine interessante Lektüre.
Avraham

2
Ein eher spezieller Fall ist X eine Standardnormale und Y unabhängig ±1 jeweils mit Wahrscheinlichkeit 12 , dannX,YundXY hat den gleichen Mittelwert und die gleiche Varianz undXY ist normalverteilt.
Henry

1
" Das Verhältnis zweier unabhängiger Chi-Quadrat-Verteilungen ergibt eine F-Verteilung. " Es gibt eine Beta-Prime-Verteilung. Um ein F zu erhalten, müssen Sie jedes Chi-Quadrat mit seinem df skalieren.
Glen_b

2
Eine Reihe von Dingen überzeugt mich überhaupt nicht, dass es notwendig ist, alle Ihre Bedingungen zu erfüllen.
Glen_b

1
am Beispiel der Methode zur Erzeugung von Normalvariablen (z. B. Box-Muller) (die die Kreismethode verwendet) würde ich sagen, dass es keine Verhältnisse von Gleichverteilungen gibt , die eine Normalverteilung ergeben (vorausgesetzt, es werden Gleichverteilungen verlangt)
Nikos M.

Antworten:


5

Sei wobeiEeine Exponentialverteilung mit Mittelwert2σ2undZ=±1mit gleicher Wahrscheinlichkeit hat. SeiY2=1/Y1=ZEE2σ2Z=±1 wobeiBBeta(0,5,0,5) ist. UnterAnnahme(Z,E,B)voneinander unabhängig, dannY1unabhängig davon istY2undY1/Y2~Normale(0,σ2). Daher haben wirY2=1/BBBeta(0.5,0.5)(Z,E,B)Y1Y2Y1/Y2Normal(0,σ2)

  1. unabhängig von Y 2 ;Y1Y2
  2. Beide kontinuierlich; so dass
  3. .Y1/Y2Normal(0,σ2)

Ich habe nicht herausgefunden, wie man ein bekommt . Es ist schwieriger zu sehen, wie dies zu tun ist, da das Problem sich darauf reduziert, A und B zu finden, die so unabhängig sind, dass A - B μNormal(μ,σ2)AB das ist ziemlich viel härter alsmachenA/B~Normale(0,1)für unabhängigAundB.

ABμBNormal(0,1)
A/BNormal(0,1)AB

1
Wenn das wahr ist, ist das fantastisch.
Neil G

2
@NeilG es ist wahr; Das Produkt aus meinem Beta und Exponential ist ein Gamma mit der Form 1/2 (da Sie das Beta und ein unabhängiges Gamma mithilfe von Gammas erstellen können). Dann ist die Quadratwurzel davon halbnormal unter Verwendung der Tatsache, dass das Quadrat einer Normalen ein Chi-Quadrat ist.
Kerl

1
Wir hatten kürzlich eine Frage zu einem Produkt mit zwei normalverteilten Variablen (ich kann es nicht zurückfinden). Diese Frage hatte einen Kommentar oder eine Antwort in Bezug auf die Box-Muller-Transformation , die eine Normalverteilung (oder genauer eine bivariate Normalverteilung) aus dem Produkt zweier transformierter gleichverteilter Variablen berechnet. Diese Antwort hängt viel damit zusammen, nimmt aber die Umkehrung einer dieser Variablen in der Box-Muller-Transformation. cc: @kjetilbhalvorsen
Sextus Empiricus

1

Ich würde es besonders begrüßen, wenn die beiden unabhängigen Verteilungen, deren Verhältnis berechnet wird, gleich sind 

Es gibt keine Möglichkeit , dass eine normale Variable kann als ein Verhältnis von zwei unabhängigen Variablen mit der geschrieben wird gleiche Verteilung oder Verteilungsfamilie (wie die F-Verteilung , die das Verhältnis von zwei skaliert wird χ2 verteilte Variablen oder die Cauchy-Verteilung , die ist das Verhältnis zweier normalverteilter Variablen mit dem Mittelwert Null).

  • Angenommen, für jedes A,BF wobei F die gleiche Verteilung oder Verteilungsfamilie ist, gilt

    X=ABN(μ,σ2)

  • Wir müssen auch in der Lage sein, A und B umzukehren (wenn eine normale Variable als Verhältnis zweier unabhängiger Variablen mit derselben Verteilung oder Verteilungsfamilie geschrieben werden kann, kann die Reihenfolge umgekehrt werden)

    1X=BAN(μ,σ2)

  • Wenn jedoch XN(μ,σ2) ist, kann X1N(μ,σ2) nicht wahr sein (die Inverse einer normalverteilten Variablen ist keine andere normalverteilte Variable).

Umfassende Schlussfolgerung: Wenn die Variablen in einer Verteilungsfamilie FX als Verhältnis von Variablen in einer anderen Verteilungsfamilie FY können, muss die Familie FX unter Berücksichtigung des Kehrwerts geschlossen werden (dh für jede Variable, deren Verteilung in ist) FX die Verteilung des Kehrwerts wird auch in FX ) angegeben.

ZB ist das Inverse einer verteilten Cauchy-Variablen auch Cauchy-verteilt. Das Inverse einer F-verteilten Variablen ist ebenfalls F-verteilt.

  • Dieses "Wenn" ist kein "Wenn", das Gegenteil ist nicht wahr. Wenn sich X und 1/X in derselben Verteilungsfamilie befinden, kann es sein, dass es nicht immer möglich ist, eine Verhältnisverteilung mit einem Nenner und einem Nenner derselben Verteilungsfamilie zu schreiben.

    Gegenbeispiel: Wir können uns Verteilungsfamilien vorstellen, für die wir für jedes X in der Familie 1/X in derselben Familie haben, aber kein P(X=1)=0 . Dies widerspricht der Tatsache, dass für eine Verhältnisverteilung, bei der Nenner und Nenner die gleiche Verteilung haben,P(X=1)0 (und etwas Ähnliches kann für kontinuierliche Verteilungen wie das Integral entlang der Linie X / Y ausgedrückt werden = 1 in einem Streudiagramm von X hat Y eine Dichte ungleich Null, wenn X und Y die gleiche Verteilung haben und unabhängig sind).


Sieh es nicht. Scheint mir, dass nur weil und B / C normal sind, das A / D nicht machtA/DB/C normal. A/DB/C
Carl

Besser. Jetzt macht es Sinn.
Carl

1
Ich verstehe nicht, wie die zweite Aussage aus der ersten folgt. Wenn es etwas , deren Quotient normal ist, warum folgt daraus, dass deren Quotient in der anderen Reihenfolge ebenfalls normal sein sollte? Die Frage stellte keine Frage nach einer Verteilungsfamilie, so dass der Quotient aller Elementpaare normal ist. A,B
Neil G

1
Ich verstehe nicht, was du sagst. Im Idealfall ist Ihre Antwort ein schlüssiges Argument, ohne dass jemand die Änderungen lesen muss. Momentan scheint es so, als ob Ihre zweite Aussage ("wir müssen sie auch haben") nicht aus der ersten folgt.
Neil G

1
@kjetilbhalvorsen wie muss es überarbeitet werden? Ich habe den Teil der Frage beantwortet, in dem es heißt: "Ich wäre besonders dankbar, wenn die beiden unabhängigen Verteilungen, deren Verhältnis berechnet wird, gleich wären . " Ich verstehe nicht, wie sich die Antwort des Mannes darauf bezieht.
Sextus Empiricus

0

Nun, hier ist einer, aber ich werde ihn nicht beweisen, sondern nur in der Simulation zeigen.

Machen Sie zwei Beta-Verteilungen mit gleich großen Formparametern (hier n = 40 , 000 ), subtrahieren Sie 1/2 von xBeta(200,200)n=40,000x Werten von einem von ihnen und nennen Sie es "Zähler". Das gibt uns ein PDF mit einer maximalen Reichweite von , aber weil die Formparameter so groß sind, erreichen wir nie die Maximalwerte des Bereichs. Hier ist ein Histogramm einesn(12,12) „Zähler“ n=40,000enter image description here

Als nächstes nennen wir die zweite Beta-Verteilung "Nenner", ohne etwas zu subtrahieren, sodass sie den üblichen Beta-Verteilungsbereich von und einer davon so aussieht(0,1)

enter image description here

Da die Formen so groß sind, nähern wir uns auch hier nicht dem Maximalbereich mit den Werten. Als nächstes werden wir den Quotienten plotten als PDF mit der überlagerten Normalverteilung.numeratordenominator

enter image description here

In diesem Fall hat das Normalverteilungsergebnis und prüft auf Normalität, die so aussiehtμ0.0000204825,σ0.0501789

(StatisticP-ValueAnderson-Darling0.7997860.481181Baringhaus-Henze1.405850.0852017Cramér-von Mises0.1231450.482844Jarque-Bera ALM4.481030.106404Kolmogorov-Smirnov0.004523280.386335Kuiper0.007980630.109127Mardia Combined4.481030.106404Mardia Kurtosis1.538490.123929Mardia Skewness2.093990.147879Pearson χ2134.3530.571925Watson U20.1138310.211187)

Mit anderen Worten, wir können nicht beweisen, dass das Verhältnis nicht normal ist, auch wenn wir uns sehr bemühen, dies zu tun.

Jetzt, warum? Intuition meinerseits, die ich im Überfluss habe. Dem Leser überlassener Beweis, falls vorhanden (möglicherweise über die Begrenzung der Methode der Momente, aber auch dies ist nur Intuition).

Beta(20,20)Beta(20,20)12tμ0.000251208,σ0.157665,df33.0402

enter image description here

StatisticP-ValueAnderson-Darling0.2752620.955502Cramér-von Mises0.03511080.956524Kolmogorov-Smirnov0.003209360.804486Kuiper0.005565010.657146Pearson χ2145.0770.323168Watson U20.03510420.878202

Another hint N(0,1)N(10,1/1000) Student's t μ0.0000535722,σ0.0992765,df244.154

enter image description here

(StatisticP-ValueAnderson-Darling0.5016770.745102Cramér-von Mises0.06968240.753515Kolmogorov-Smirnov0.003556880.692225Kuiper0.006083820.501133Pearson χ2142.880.370552Watson U20.06032070.590369)

5
You are clearly very close to a normal distribution. However, that isn't at all the same thing as having a normal distribution, and I don't believe the ratio of a centered symmetric beta to an ordinary symmetric beta with the same parameters is ever to be actually normal. I'd be very interested in being wrong about this though.
Glen_b -Reinstate Monica

2
Your solution definitely is not Normal. You could generalize this approach: take any distribution that is approximately Normal and divide it by a distribution with its probability concentrated near a nonzero number. The result (obviously) will be close to Normal--but it still will not be Normal. Applying a bunch of tests is unconvincing because all it shows it that you didn't generate sufficiently large samples to demonstrate the non-Normality.
whuber

1
@whuber At 108 one could also show in normal machine precision that noise will cause anything to be not normal. I do not have a super computer to do that in extended precision. What you could show, and why I wanted you to look at this, was to prove or disprove these things mathematically, not just criticize around the edges with unachievable goals.
Carl

2
Let me get to the heart of the matter, then: (1) disproving normality is a simple exercise in integral approximation--no need to give the details here. You can, e.g., readily prove the 200th moment is infinite. (2) Your answer confuses distributions with samples. It is this fundamental confusion that I object to; it's the reason why I think this answer is more misleading than helpful. BTW, I did not write my last comment lightly: I performed that test. I did it not with a supercomputer, but with a decade-old PC workstation, and the whole process took just seconds.
whuber

1
@whuber Which approximation are you testing? The first, the second or the third? BTW, if they are only approximations, so be it. I suggest only that in the limiting case that they might be exact. All of statistics is an approximation so I do not share your apprehension.
Carl

-3

I imagine there are many possibilities. Here there is one that I can think of. It is known (Zolotarev) that, given X1G,X2G two standard normal distributed r.v., and XγC a Cauchy distributed r.v.

X1GX2G=XγC

Then, by Duality of the Stable distribution, we know that XγC1/X1/γC (where γ is the scale parameter of the Cauchy). So you get that the Normal distribution can be a result from a ratio between a Normal and a Cauchy:

X1G=X2G/X1/γC

for the desired μ I would just move both distributions to be centred there. (at μ). For the σ, in the mentioned wikipedia page about ratio distributions, there are the general formulas for the ratio of two normal distributions, you would just need to replace the scale factor of the Cauchy by its inverse value (γ1/γ).


4
Please test your hypothesis, either by explicit calculation of the ratio or via simulation. Either will show that your claim is incorrect. The error lies in assuming that distribution ratios can be "canceled" to "solve for" the numerator.
whuber

1
hm, it is true that the passage of the X2G to the right is shady. I'll check.
chuse
Durch die Nutzung unserer Website bestätigen Sie, dass Sie unsere Cookie-Richtlinie und Datenschutzrichtlinie gelesen und verstanden haben.
Licensed under cc by-sa 3.0 with attribution required.