¿Qué proporción de distribuciones independientes da una distribución normal?

12

La relación de dos distribuciones normales independientes da una distribución de Cauchy. La distribución t es una distribución normal dividida por una distribución chi-cuadrado independiente. La relación de dos distribuciones chi-cuadrado independientes da una distribución F.

¿Estoy buscando una relación de distribuciones continuas independientes que proporcione una variable aleatoria normalmente distribuida con media y varianza ?σ 2μσ2

Probablemente hay un conjunto infinito de posibles respuestas. ¿Me puede dar algunas de estas posibles respuestas? Apreciaría particularmente si las dos distribuciones independientes cuya relación se calcula son iguales o al menos tienen una varianza similar.

Remi.b
fuente
2
Si bien el artículo de Wikipedia sobre distribuciones de razón no proporciona ejemplos del caso que busca, es una lectura interesante.
Abraham
2
Un caso bastante especial es X normal normal e Y independientemente ±1 cada uno con probabilidad 12 , luegoX,YyXY tiene la misma media y varianza yXY se distribuye normalmente.
Henry
1
" La relación de dos distribuciones chi-cuadrado independientes da una distribución F " --- bueno, no del todo. Da una distribución beta-prime. Para obtener una F necesitas escalar cada chi-cuadrado por su df.
Glen_b -Reinstalar Monica
2
Varias cosas me hacen no estar del todo convencido de que es necesariamente posible cumplir con todas sus condiciones.
Glen_b -Reinstale a Monica el
1
tomando el método de generación de variables normales (por ejemplo, Box-Muller) como ejemplo (que usa el método de círculo), diría que no hay proporciones de distribuciones uniformes que den una distribución normal (suponiendo que se soliciten distribuciones uniformes)
Nikos M.

Respuestas:

5

Deje dondeEtiene una distribución exponencial con media2σ2yZ=±1con igual probabilidad. DejeY2=1/Y1=Zmimi2σ2Z=±1 dondeBBeta(0.5,0.5). Suponiendo que(Z,E,B)son mutuamente independientes, entoncesY1es independiente deY2eY1/Y2Normal(0,σ2). Por lo tanto tenemosY2=1/ /sisiBeta(0.5 0.5,0.5 0.5)(Z,mi,si)Y1Y2Y1/ /Y2Normal(0 0,σ2)

  1. independiente de Y 2 ;Y1Y2
  2. Ambos continuos; tal que
  3. .Y1/ /Y2Normal(0 0,σ2)

No he descubierto cómo obtener un . Es más difícil ver cómo hacerlo, ya que el problema se reduce a encontrar A y B, que son independientes de tal manera que A - B μNormal(μ,σ2)UNsi que es bastante más difícil que hacer queA/BNormal(0,1)paraAyBindependientes.

UN-siμsiNormal(0 0,1)
UN/ /siNormal(0 0,1)UNsi
chico
fuente
1
Si esto es cierto, esto es asombroso.
Neil G
2
@NeilG es cierto; El producto de mi beta y exponencial es una gamma con forma 1/2 (debido a cómo puedes construir la beta y una gamma independiente usando gammas). Entonces la raíz cuadrada de eso es medio normal usando el hecho de que el cuadrado de una normal es chi-cuadrado.
chico
1
Recientemente tuvimos una pregunta pidiendo un producto de dos variables que se distribuye normalmente (no puedo encontrarlo de nuevo). Esa pregunta tenía un comentario o respuesta relacionada con la transformación de Box-Muller que calcula una distribución normal (o más precisamente una distribución normal bivariada) del producto de dos variables distribuidas uniformes transformadas. Esta respuesta se relaciona mucho con eso, pero toma el inverso de una de esas variables en la transformación de Box-Muller. cc: @kjetilbhalvorsen
Sextus Empiricus
1

Apreciaría especialmente si las dos distribuciones independientes cuya proporción se calcula son las mismas 

No hay posibilidad de que una variable normal pueda escribirse como una razón de dos variables independientes con la misma distribución o familia de distribución (como la distribución F, que es la razón de dos variables distribuidas χ2 escaladas o la distribución Cauchy, que es la razón de dos variables distribuidas normales con media cero).

  • Supongamos que: para cualquier UN,siF donde F es la misma distribución o familia de distribución tenemos

    X=UNsinorte(μ,σ2)

  • También debemos poder revertir UN y si (si una variable normal puede escribirse como una razón de dos variables independientes con la misma distribución o familia de distribución, entonces el orden puede revertirse)

    1X=siUNnorte(μ,σ2)

  • Pero si Xnorte(μ,σ2) entonces X-1norte(μ,σ2) no puede ser cierto (lo inverso de una variable distribuida normal no es otra variable distribuida normal).

Conclusión más amplia: si las variables en cualquier familia de distribución FX pueden escribirse como una relación de variables en otra familia de distribución FY entonces debe ser que la familia FX está cerrada tomando el recíproco (es decir, para cualquier variable cuya distribución está en FX la distribución de su recíproco también estará en FX ).

Por ejemplo, el inverso de una variable distribuida de Cauchy también está distribuido por Cauchy. La inversa de una variable distribuida en F también está distribuida en F.

  • Este 'si' no es un 'iff', lo contrario no es cierto. Cuando X y 1/ /X están en la misma familia de distribución, entonces no siempre se puede escribir como una distribución de razón con nominador y denominador de la misma familia de distribución.

    Contraejemplo: podemos imaginar familias de distribución para las cuales para cualquier X en la familia tenemos 1/ /X en la misma familia pero no tenemos PAG(X=1)=0 0 . Esto es contradictorio con el hecho de que para una distribución de razón donde el denominador y el nominador tienen la misma distribución, debemos tener PAG(X=1)0 0 (y algo similar se puede expresar para distribuciones continuas como la integral a lo largo de la línea X / Y = 1 en un diagrama de dispersión de X, Y tiene una densidad distinta de cero cuando X e Y tienen la misma distribución y son independientes).

Sexto empírico
fuente
No lo veo Me parece que solo porque y B / C son normales, eso no hace que A / DUN/ /resi/ /C normal. UN/ /resi/ /C
Carl
Mejor. Ahora tiene sentido.
Carl
1
No entiendo cómo la segunda declaración se sigue de la primera. Si existe algún tal que su cociente sea normal, ¿por qué se deduce que su cociente en el otro orden también debería ser normal? La pregunta no pedía una familia de distribución tal que el cociente de todos los pares de elementos sea normal. UN,si
Neil G
1
No entiendo lo que dices. Idealmente, su respuesta sería un argumento coherente sin requerir que alguien lea las ediciones. En este momento, parece que su segunda declaración ("también debemos tener") no se sigue de la primera.
Neil G
1
@kjetilbhalvorsen, ¿cómo debe revisarse? He respondido la parte de la pregunta que especifica "Agradecería especialmente si las dos distribuciones independientes cuya relación se calcula son las mismas" . No veo cómo la respuesta de chico se relaciona con eso.
Sextus Empiricus
0

Bueno, aquí hay uno, pero no lo probaré, solo lo mostraré en simulación.

Haga dos distribuciones beta con los mismos parámetros de forma grande (aquí, n = 40 , 000 ), reste 1/2 de los valores x de uno de ellos y llámelo "numerador". Eso nos da un PDF que tiene un rango máximo de ( - 1Beta(200,200)norte=40,000X, pero debido a que los parámetros de forma son tan grandes, nunca llegamos a los valores máximos del rango. Aquí hay un histograma de unn=40,000"numerador" (-12,12)norte=40,000ingrese la descripción de la imagen aquí

A continuación, llamamos al segundo denominador de distribución beta "sin restar nada", por lo que tiene el rango de distribución beta habitual de y uno de esos se ve así(0 0,1)

ingrese la descripción de la imagen aquí

Nuevamente, debido a que las formas son tan grandes, no nos acercamos al rango máximo con los valores. A continuación, graficamos el numerador del cociente como PDF con la distribución normal superpuesta.numeradordenominador

ingrese la descripción de la imagen aquí

Ahora, en este caso, el resultado de distribución normal tiene y pruebas de normalidad que se ven asíμ-0.0000204825,σ0.0501789

(EstadísticaValor pAnderson-Darling0.7997860.481181Baringhaus-Henze1.405850.0852017Cramér-von Mises0.1231450.482844Jarque-Bera ALM4.481030.106404Kolmogorov-Smirnov0.004523280.386335Kuiper0.007980630.109127Mardia combinada4.481030.106404Kurtosis de Mardia1.538490.123929Torcedura de Mardia2.093990.147879Pearson χ2134,3530.571925Watson U20.1138310.211187)

En otras palabras, no podemos probar que la relación no sea normal, incluso si tratamos de hacerlo con mucha fuerza.

¿Ahora por qué? Intuición de mi parte, que tengo en exceso. Prueba dejada al lector, si existe alguna (tal vez a través del límite del método de los momentos, pero de nuevo eso es solo intuición).

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

ingrese la descripción de la imagen aquí

EstadísticaValor pAnderson-Darling0.2752620.955502Cramér-von Mises0,03511080.956524Kolmogorov-Smirnov0.003209360,804486Kuiper0.005565010.657146Pearson χ2145,0770.323168Watson U20.03510420.878202

norte(0 0,1)norte(10,1/ /1000)t μ-0.0000535722,σ0.0992765,df244,154

ingrese la descripción de la imagen aquí

(EstadísticaValor pAnderson-Darling0.5016770.745102Cramér-von Mises0.06968240,753515Kolmogorov-Smirnov0.003556880.692225Kuiper0.006083820.501133Pearson χ2142,880.370552Watson U20,06032070.590369)
Carl
fuente
55
Está claramente muy cerca de una distribución normal. Sin embargo, eso no es lo mismo que tener una distribución normal, y no creo que la relación de una beta simétrica centrada a una beta simétrica ordinaria con los mismos parámetros sea realmente normal. Sin embargo, me interesaría mucho estar equivocado sobre esto.
Glen_b: reinstala a Monica el
2
Su solución definitivamente no es normal. Podría generalizar este enfoque: tome cualquier distribución que sea aproximadamente Normal y divídala por una distribución con su probabilidad concentrada cerca de un número distinto de cero. El resultado (obviamente) estará cerca de Normal, pero aún no será Normal. Aplicar un montón de pruebas no es convincente porque todo lo que demuestra es que no generó muestras suficientemente grandes para demostrar la no normalidad.
whuber
1
108
2
Permítanme llegar al meollo del asunto, entonces: (1) refutar la normalidad es un simple ejercicio de aproximación integral: no es necesario dar los detalles aquí. Puede, por ejemplo , demostrar fácilmente que el momento 200 es infinito. (2) Su respuesta confunde las distribuciones con muestras. Es esta confusión fundamental a la que me opongo; Es la razón por la que creo que esta respuesta es más engañosa que útil. Por cierto, no escribí mi último comentario a la ligera: realicé esa prueba. No lo hice con una supercomputadora, sino con una estación de trabajo para PC de hace una década, y todo el proceso tomó solo unos segundos.
whuber
1
@whuber ¿Qué aproximación estás probando? ¿El primero, el segundo o el tercero? Por cierto, si son solo aproximaciones, que así sea. Solo sugiero que, en el caso límite, podrían ser exactos. Todas las estadísticas son aproximadas, así que no comparto su aprensión.
Carl
-3

X1sol,X2solXγC

X1solX2sol=XγC

XγC1/ /X1/ /γCγ

X1sol=X2sol/ /X1/ /γC

μμσγ1/ /γ

chuse
fuente
44
Pruebe su hipótesis, ya sea mediante el cálculo explícito de la relación o mediante la simulación. Cualquiera de los dos mostrará que su reclamo es incorrecto. El error radica en suponer que las relaciones de distribución se pueden "cancelar" para "resolver" el numerador.
whuber
1
X2sol