Articulo de referencia

Función generadora de momentos

En teoría de la probabilidad y estadística , la función generadora de momentos de una variable aleatoria de valor real es una función generadora que proporciona una especificaci...

En teoría de la probabilidad y estadística , la función generadora de momentos de una variable aleatoria de valor real es una función generadora que proporciona una especificación alternativa de la distribución de probabilidad de la variable aleatoria . De este modo, constituye la base de una vía alternativa para obtener resultados analíticos en comparación con el trabajo directo con funciones de densidad de probabilidad o funciones de distribución acumulativa . Existen resultados particularmente sencillos para las funciones generadoras de momentos de distribuciones definidas por sumas ponderadas de variables aleatorias. Sin embargo, no todas las variables aleatorias poseen funciones generadoras de momentos.

Como su nombre lo indica, la función generadora de momentos se puede utilizar para calcular los momentos de una distribución : el n -ésimo momento con respecto a 0 es la n -ésima derivada de la función generadora de momentos, evaluada en 0.

Además de las distribuciones univariadas de valores reales, las funciones generadoras de momentos también pueden definirse para variables aleatorias con valores vectoriales o matriciales, e incluso pueden extenderse a casos más generales.

La función generadora de momentos de una distribución de valores reales no siempre existe, a diferencia de la función característica . Existen relaciones entre el comportamiento de la función generadora de momentos de una distribución y las propiedades de la misma, como la existencia de momentos.

Definición

Dejarincógnita{\displaystyle X}sea ​​una variable aleatoria con función de distribución acumulada (CDF)Fincógnita{\displaystyle F_{X}}. La función generadora de momentos (mgf) deincógnita{\displaystyle X}(oFincógnita{\displaystyle F_{X}}), denotado porMETROincógnita(t){\displaystyle M_{X}(t)}, es

METROincógnita(t)=mi[mitincógnita]{\displaystyle M_{X}(t)=\operatorname {E} \left[e^{tX}\right]}

siempre que exista esta expectativat{\displaystyle t}en algún vecindario abierto de 0. Es decir, hay unh>0{\displaystyle h>0}de tal manera que para todost{\displaystyle t}satisfactorio h<t<h{\displaystyle -h<t<h}, mi[mitincógnita]{\displaystyle \operatorname {E} \left[e^{tX}\right]}existe. Si la esperanza no existe en un entorno abierto de 0, decimos que la función generadora de momentos no existe. [ 1 ]

En otras palabras, la función generadora de momentos de X es la esperanza de la variable aleatoria.mitincógnita{\displaystyle e^{tX}}. En términos más generales, cuandoincógnita=(incógnita1,,incógnitanorte)T{\displaystyle \mathbf {X} =(X_{1},\ldots ,X_{n})^{\mathrm {T} }}, unnorte{\displaystyle n}vector aleatorio de dimensión , yt{\displaystyle \mathbf {t} }es un vector fijo, uno usatincógnita=tTincógnita{\displaystyle \mathbf {t} \cdot \mathbf {X} =\mathbf {t} ^{\mathrm {T} }\mathbf {X} }en lugar de tincógnita{\displaystyle tX}:METROincógnita(t):=mi[mitTincógnita].{\displaystyle M_{\mathbf {X} }(\mathbf {t} ):=\operatorname {E} \left[e^{\mathbf {t} ^{\mathrm {T} }\mathbf {X} }\right].}

METROincógnita(0){\displaystyle M_{X}(0)}Siempre existe y es igual a  1. Sin embargo, un problema clave con las funciones generadoras de momentos es que los momentos y la función generadora de momentos pueden no existir, ya que las integrales no tienen por qué converger absolutamente. Por el contrario, la función característica o transformada de Fourier siempre existe (porque es la integral de una función acotada en un espacio de medida finita ) y, para algunos fines, puede utilizarse en su lugar.

La función generadora de momentos se denomina así porque se puede utilizar para encontrar los momentos de la distribución. [ 2 ] El desarrollo en serie demitincógnita{\displaystyle e^{tX}}es

mitincógnita=1+tincógnita+t2incógnita22¡+t3incógnita33¡++tnorteincógnitanortenorte¡+.{\displaystyle e^{tX}=1+tX+{\frac {t^{2}X^{2}}{2!}}+{\frac {t^{3}X^{3}}{3!}}+\cdots +{\frac {t^{n}X^{n}}{n!}}+\cdots .}

Por eso, METROincógnita(t)=mi[mitincógnita]=1+tmi[incógnita]+t2mi[incógnita2]2¡+t3mi[incógnita3]3¡++tnortemi[incógnitanorte]norte¡+=1+tmetro1+t2metro22¡+t3metro33¡++tnortemetronortenorte¡+,{\displaystyle {\begin{aligned}M_{X}(t)&=\operatorname {E} [e^{tX}]\\[1ex]&=1+t\operatorname {E} [X]+{\frac {t^{2}\operatorname {E} [X^{2}]}{2!}}+{\frac {t^{3}\operatorname {E} [X^{3}]}{3!}}+\cdots +{\frac {t^{n}\operatorname {E} [X^{n}]}{n!}}+\cdots \\[1ex]&=1+tm_{1}+{\frac {t^{2}m_{2}}{2!}}+{\frac {t^{3}m_{3}}{3!}}+\cdots +{\frac {t^{n}m_{n}}{n!}}+\cdots ,\end{aligned}}}

dóndemetronorte{\displaystyle m_{n}}es elnorte{\displaystyle n}-ésimo momento . DiferenciandoMETROincógnita(t){\displaystyle M_{X}(t)}i{\displaystyle i}tiempos con respecto at{\displaystyle t}y configuraciónt=0{\displaystyle t=0}, obtenemos eli{\displaystyle i}-momento sobre el origen,metroi{\displaystyle m_{i}}.

Siincógnita{\displaystyle X}es una variable aleatoria continua, la siguiente relación entre su función generadora de momentosMETROincógnita(t){\displaystyle M_{X}(t)}y la transformada de Laplace bilateral de su función de densidad de probabilidad.Fincógnita(incógnita){\displaystyle f_{X}(x)}contiene:

METROincógnita(t)=L{Fincógnita}(t),{\displaystyle M_{X}(t)={\mathcal {L}}\{f_{X}\}(-t),}

ya que la transformada de Laplace bilateral de la PDF se da como

L{Fincógnita}(s)=misincógnitaFincógnita(incógnita)dincógnita,{\displaystyle {\mathcal {L}}\{f_{X}\}(s)=\int _{-\infty }^{\infty }e^{-sx}f_{X}(x)\,dx,}

y la definición de la función generadora de momentos se expande (por la ley del estadístico inconsciente ) a METROincógnita(t)=mi[mitincógnita]=mitincógnitaFincógnita(incógnita)dincógnita.{\displaystyle M_{X}(t)=\operatorname {E} \left[e^{tX}\right]=\int _{-\infty }^{\infty }e^{tx}f_{X}(x)\,dx.}

Esto es consistente con la función característica deincógnita{\displaystyle X}ser una rotación de Wick deMETROincógnita(t){\displaystyle M_{X}(t)}cuando existe la función generadora de momentos, como la función característica de una variable aleatoria continuaincógnita{\displaystyle X}es la transformada de Fourier de su función de densidad de probabilidadFincógnita(incógnita){\displaystyle f_{X}(x)}y en general cuando una funciónF(incógnita){\displaystyle f(x)}es de orden exponencial , la transformada de Fourier deF{\displaystyle f}es una rotación de Wick de su transformada de Laplace bilateral en la región de convergencia. Consulte la relación entre las transformadas de Fourier y Laplace para obtener más información.

Ejemplos

Aquí se muestran algunos ejemplos de la función generadora de momentos y la función característica para su comparación. Se puede observar que la función característica es una rotación de Wick de la función generadora de momentos.METROincógnita(t){\displaystyle M_{X}(t)}cuando esto último existe.

Cálculo

La función generadora de momentos es la esperanza de una función de la variable aleatoria, y se puede escribir como:

Tenga en cuenta que para el caso dondeincógnita{\displaystyle X}tiene una función de densidad de probabilidad continuaF(incógnita){\displaystyle f(x)}, METROincógnita(t){\displaystyle M_{X}(-t)}es la transformada de Laplace de dos lados deF(incógnita){\displaystyle f(x)}.

METROincógnita(t)=mitincógnitaF(incógnita)dincógnita=(1+tincógnita+t2incógnita22¡++tnorteincógnitanortenorte¡+)F(incógnita)dincógnita=1+tmetro1+t2metro22¡++tnortemetronortenorte¡+,{\displaystyle {\begin{aligned}M_{X}(t)&=\int _{-\infty }^{\infty }e^{tx}f(x)\,dx\\[1ex]&=\int _{-\infty }^{\infty }\left(1+tx+{\frac {t^{2}x^{2}}{2!}}+\cdots +{\frac {t^{n}x^{n}}{n!}}+\cdots \right)f(x)\,dx\\[1ex]&=1+tm_{1}+{\frac {t^{2}m_{2}}{2!}}+\cdots +{\frac {t^{n}m_{n}}{n!}}+\cdots ,\end{aligned}}}

dóndemetronorte{\displaystyle m_{n}}es elnorte{\displaystyle n}momento .

Transformaciones lineales de variables aleatorias

Si la variable aleatoriaincógnita{\displaystyle X}tiene función generadora de momentosMETROincógnita(t){\displaystyle M_{X}(t)}, entoncesαincógnita+β{\displaystyle \alpha X+\beta }tiene función generadora de momentosMETROαincógnita+β(t)=miβtMETROincógnita(αt){\displaystyle M_{\alpha X+\beta }(t)=e^{\beta t}M_{X}(\alpha t)}

METROαincógnita+β(t)=mi[mi(αincógnita+β)t]=miβtmi[miαincógnitat]=miβtMETROincógnita(αt){\displaystyle M_{\alpha X+\beta }(t)=\operatorname {E} \left[e^{(\alpha X+\beta )t}\right]=e^{\beta t}\operatorname {E} \left[e^{\alpha Xt}\right]=e^{\beta t}M_{X}(\alpha t)}

Combinación lineal de variables aleatorias independientes

SiSnorte=i=1norteaiincógnitai{\textstyle S_{n}=\sum _{i=1}^{n}a_{i}X_{i}}donde las X i son variables aleatorias independientes y las a i son constantes, entonces la función de densidad de probabilidad para S n es la convolución de las funciones de densidad de probabilidad de cada una de las X i , y la función generadora de momentos para S n viene dada por

METROSnorte(t)=METROincógnita1(a1t)METROincógnita2(a2t)METROincógnitanorte(anortet).{\displaystyle M_{S_{n}}(t)=M_{X_{1}}(a_{1}t)M_{X_{2}}(a_{2}t)\cdots M_{X_{n}}(a_{n}t)\,.}

Variables aleatorias con valores vectoriales

Para variables aleatorias con valores vectorialesincógnita{\displaystyle \mathbf {X} }con componentes reales , la función generadora de momentos viene dada por

METROincógnita(t)=mi[mit,incógnita]{\displaystyle M_{X}(\mathbf {t} )=\operatorname {E} \left[e^{\langle \mathbf {t} ,\mathbf {X} \rangle }\right]}

dóndet{\displaystyle \mathbf {t} }es un vector y,{\displaystyle \langle \cdot ,\cdot \rangle }es el producto escalar .

Propiedades importantes

Las funciones generadoras de momentos son positivas y log-convexas , con M (0) = 1.

Una propiedad importante de la función generadora de momentos es que determina de forma única la distribución. En otras palabras, siincógnita{\displaystyle X}yY{\displaystyle Y}son dos variables aleatorias y para todos los valores de t , 

METROincógnita(t)=METROY(t),{\displaystyle M_{X}(t)=M_{Y}(t),} entonces Fincógnita(incógnita)=FY(incógnita){\displaystyle F_{X}(x)=F_{Y}(x)}

para todos los valores de x (o equivalentemente, X e Y tienen la misma distribución). Esta afirmación no es equivalente a la afirmación "si dos distribuciones tienen los mismos momentos, entonces son idénticas en todos los puntos". Esto se debe a que, en algunos casos, los momentos existen, pero la función generadora de momentos no, debido al límite.

límitenortei=0nortetimetroii¡{\displaystyle \lim _{n\to \infty }\sum _{i=0}^{n}{\frac {t^{i}m_{i}}{i!}}}

Puede que no exista. La distribución lognormal es un ejemplo de cuándo ocurre esto: sus momentos sonmi[incógnitanorte]=minorteμ+norte2σ2/2{\displaystyle \operatorname {E} [X^{n}]=e^{n\mu +n^{2}\sigma ^{2}/2}}y son todos finitos, pero su función generadora de momentosmi[mitincógnita]{\displaystyle \operatorname {E} \left[e^{tX}\right]}no está definida para ningún t positivo ya que la integral diverge y, por lo tanto, no está en un entorno de 0; existen otras distribuciones con los mismos momentos. [ 4 ]

Cálculos de momentos

La función generadora de momentos se denomina así porque si existe en un intervalo abierto alrededor de t = 0 , entonces es la función generadora exponencial de los momentos de la distribución de probabilidad :

metronorte=mi[incógnitanorte]=METROincógnita(norte)(0)=dnorteMETROincógnitadtnorte|t=0.{\displaystyle m_{n}=\operatorname {E} \left[X^{n}\right]=M_{X}^{(n)}(0)=\left.{\frac {d^{n}M_{X}}{dt^{n}}}\right|_{t=0}.}

Es decir, siendo n un entero no negativo, el n -ésimo momento alrededor de 0 es la n -ésima derivada de la función generadora de momentos, evaluada en t = 0 .

Otras propiedades

La desigualdad de Jensen proporciona una cota inferior simple para la función generadora de momentos: METROincógnita(t)miμt,{\displaystyle M_{X}(t)\geq e^{\mu t},} dóndeμ{\displaystyle \mu }es la media de X.

La función generadora de momentos se puede utilizar junto con la desigualdad de Markov para acotar la cola superior de una variable aleatoria real X. Esta afirmación también se denomina cota de Chernoff . Dado queincógnitamiincógnitat{\displaystyle x\mapsto e^{xt}}es monótonamente creciente parat>0{\displaystyle t>0}, tenemos Pr(incógnitaa)=Pr(mitincógnitamita)miatmi[mitincógnita]=miatMETROincógnita(t){\displaystyle \Pr(X\geq a)=\Pr(e^{tX}\geq e^{ta})\leq e^{-at}\operatorname {E} \left[e^{tX}\right]=e^{-at}M_{X}(t)} para cualquiert>0{\displaystyle t>0}y cualquier a , siempre queMETROincógnita(t){\displaystyle M_{X}(t)}existe. Por ejemplo, cuando X es una distribución normal estándar ya>0{\displaystyle a>0}, podemos elegirt=a{\displaystyle t=a}y recuerda queMETROincógnita(t)=mit2/2{\displaystyle M_{X}(t)=e^{t^{2}/2}}Esto daPr(incógnitaa)mia2/2{\displaystyle \Pr(X\geq a)\leq e^{-a^{2}/2}}, que está dentro de un factor de 1+ a del valor exacto.

Diversos lemas, como el lema de Hoeffding o la desigualdad de Bennett, proporcionan límites para la función generadora de momentos en el caso de una variable aleatoria acotada con media cero.

Cuandoincógnita{\displaystyle X}Si es no negativo, la función generadora de momentos proporciona una cota simple y útil para los momentos: mi[incógnitametro](metrotmi)metroMETROincógnita(t),{\displaystyle \operatorname {E} [X^{m}]\leq \left({\frac {m}{te}}\right)^{m}M_{X}(t),} Para cualquierincógnita,metro0{\displaystyle X,m\geq 0}yt>0{\displaystyle t>0}.

Esto se deduce de la desigualdad.1+incógnitamiincógnita{\displaystyle 1+x\leq e^{x}}en el que podemos sustituirincógnita=tincógnita/metro1{\displaystyle x'=tx/m-1}implicatincógnita/metromitincógnita/metro1{\displaystyle tx/m\leq e^{tx/m-1}}para cualquierincógnita,t,metroR{\displaystyle x,t,m\in \mathbb {R} }Ahora bien, sit>0{\displaystyle t>0}yincógnita,metro0{\displaystyle x,m\geq 0}, esto se puede reorganizar aincógnitametro(metro/(tmi))metromitincógnita{\displaystyle x^{m}\leq (m/(te))^{m}e^{tx}}. Tomando la esperanza en ambos lados se obtiene el límite enmi[incógnitametro]{\displaystyle \operatorname {E} [X^{m}]}en términos demi[mitincógnita]{\displaystyle \operatorname {E} [e^{tX}]}.

Como ejemplo, consideremosincógnitaChi cuadrado{\displaystyle X\sim {\text{Chi-Squared}}}conk{\displaystyle k}grados de libertad. Luego, a partir de los ejemplosMETROincógnita(t)=(12t)k/2{\displaystyle M_{X}(t)=(1-2t)^{-k/2}}. Recogiendot=metro/(2metro+k){\displaystyle t=m/(2m+k)}y sustituyendo en el límite: mi[incógnitametro](1+2metro/k)k/2mimetro(k+2metro)metro.{\displaystyle \operatorname {E} [X^{m}]\leq {\left(1+2m/k\right)}^{k/2}e^{-m}{\left(k+2m\right)}^{m}.} Sabemos que en este caso el límite correcto esmi[incógnitametro]2metroΓ(metro+k/2)/Γ(k/2){\displaystyle \operatorname {E} [X^{m}]\leq 2^{m}\Gamma (m+k/2)/\Gamma (k/2)}Para comparar los límites, podemos considerar el comportamiento asintótico para valores grandes de .k{\displaystyle k}. Aquí el límite de la función generadora de momentos eskmetro(1+metro2/k+O(1/k2)){\displaystyle k^{m}(1+m^{2}/k+O(1/k^{2}))}donde el límite real eskmetro(1+(metro2metro)/k+O(1/k2)){\displaystyle k^{m}(1+(m^{2}-m)/k+O(1/k^{2}))}Por lo tanto, la cota de la función generadora de momentos es muy fuerte en este caso.

Relación con otras funciones

Relacionadas con la función generadora de momentos, existen otras transformaciones comunes en la teoría de la probabilidad:

Función característica
La función característicaφincógnita(t){\displaystyle \varphi _{X}(t)}está relacionado con la función generadora de momentos a través deφincógnita(t)=METROiincógnita(t)=METROincógnita(it):{\displaystyle \varphi _{X}(t)=M_{iX}(t)=M_{X}(it):}La función característica es la función generadora de momentos de iX o la función generadora de momentos de X evaluada en el eje imaginario. Esta función también puede considerarse como la transformada de Fourier de la función de densidad de probabilidad , la cual, por lo tanto, puede deducirse mediante la transformada inversa de Fourier.
Función generadora de cumulantes
La función generadora de cumulantes se define como el logaritmo de la función generadora de momentos; algunos, en cambio, la definen como el logaritmo de la función característica , mientras que otros llaman a esta última la segunda función generadora de cumulantes.
Función generadora de probabilidad
La función generadora de probabilidad se define comoGRAMO(z)=mi[zincógnita].{\displaystyle G(z)=\operatorname {E} \left[z^{X}\right].}Esto implica inmediatamente queGRAMO(mit)=mi[mitincógnita]=METROincógnita(t).{\displaystyle G(e^{t})=\operatorname {E} \left[e^{tX}\right]=M_{X}(t).}

Véase también

Referencias

Citas

  1. Casella, George ; Berger, Roger L. (1990). Inferencia estadística . Wadsworth & Brooks/Cole. pág.  61. ISBN 0-534-11958-1.
  2. Bulmer, MG (1979). Principios de estadística . Dover. págs. 75–79 . ISBN  0-486-63760-3.
  3. Kotz et al.pág. 37 utilizando 1 como número de grados de libertad para recuperar la distribución de Cauchy
  4. Heyde, CC. (1963), "Sobre una propiedad de la distribución lognormal", Journal of the Royal Statistical Society, Serie B , vol. 25, pp. 392–393 , doi : 10.1007/978-1-4419-5823-5_6 , ISBN   978-1-4419-5822-8{{citation}}: Incompatibilidad de ISBN/Fecha ( ayuda )

Fuentes

  • Casella, George; Berger, Roger (2002). Inferencia estadística (2.ª  ed.). Thomson Learning. págs. 59–68 . ISBN  978-0-534-24312-8.