Comment expliqueriez-vous la fonction de génération de moment (MGF) en termes simples?


15

Qu'est-ce qu'une fonction de génération de moment (MGF)?

Pouvez-vous l'expliquer en termes simples et avec un exemple simple et facile?

Veuillez limiter autant que possible les notations mathématiques formelles.


1
Vous voulez un exemple simple, facile ... mais sans notation mathématique? Je ne suis pas sûr qu'une telle chose serait très facile à faire - du moins pas sans risquer de donner une impression trompeuse de ce que vous traitez. Je suppose que l'on pourrait donner le mgf d'une variable aléatoire dégénérée qui est toujours sans avoir besoin de beaucoup de notation mathématique, mais ça ne sera pas éclairant si vous voulez vraiment comprendre mgfs. 0
Glen_b -Reinstate Monica

1
Je ne sais pas s'il existe une manière intuitive de la comprendre, vous pourriez simplement la considérer comme un moyen de "coder" une distribution (au moins lorsqu'elle existe, cette idée fonctionne un peu mieux avec des fonctions caractéristiques).
dsaxton

1
Une fonction de génération de moment - lorsqu'elle existe - est un moyen de coder tous les moments non négatifs entiers d'une variable aléatoire dans une fonction, et dont ils peuvent être extraits à nouveau; mgfs peut être utilisé pour effectuer des calculs particuliers qui ne sont parfois pas si faciles à faire autrement. Je ne m'attends pas à ce que cela aide beaucoup.
Glen_b -Reinstate Monica

1
Je suis sûr que vous avez vu Joe Blitztein répondre à la même question sur Quora
Antoni Parellada

Réponses:


18

Supposons qu'une intuition sans équation ne soit pas possible, et insistons toujours pour résumer les mathématiques à l'essentiel pour avoir une idée de ce qui se passe: nous essayons d'obtenir les moments statistiques , qui, après la référence obligatoire à la physique , nous définissons comme la valeur attendue d'une puissance d'une variable aléatoire. Pour une variable aléatoire continue, le ème moment brut est par LOTUS :k

(1)E[Xk]=Xkpdfdx

La fonction de génération de moment , est un moyen de contourner cette intégrale (Eq.1) en effectuant à la place:

MX(t):=E[etX],

(2)E[etX]=etXpdfdx

Pourquoi? Parce que c'est plus facile et il y a une propriété fantastique du MGF qui peut être vue en étendant la série Maclaurin de etX

etX=1+X1!t+X22!t2+X33!t3+

Prenant l'attente des deux côtés de cette série de puissance:

MX(t)=E[etX](3)=1+E[X]1!t+E[X2]2!t2+E[X3]3!t3+

les moments apparaissent "perchés" sur cette "corde à linge" polynomiale, prêts à être éliminés en différenciant simplement k fois et en évaluant à zéro une fois que nous avons traversé l'intégration plus facile (en éq. (2)) une seule fois pour tous les moments! Le fait qu'il s'agisse d'une intégration plus facile est plus apparent lorsque le pdf est exponentiel.

Pour récupérer le k ème moment:

MX(k)(0)=dkdtkMX(t)|t=0

Le fait que finalement il y ait un besoin de différenciation en fait un déjeuner pas gratuit - en fin de compte, c'est une transformation Laplace bilatérale du pdf avec un signe changé dans l'exposant:

L{pdf(x)}(s)=esxpdf(x)dx

tel que

(4)MX(t)=L{pdf(x)}(s).

Ceci, en effet, nous donne une voie physique vers l'intuition. La transformée de Laplace agit sur le pdf et le décompose en moments. La similitude avec une transformée de Fourier est inévitable : un FT mappe une fonction à une nouvelle fonction sur la ligne réelle, et Laplace mappe une fonction à une nouvelle fonction sur le plan complexe. La transformée de Fourier exprime une fonction ou un signal sous la forme d'une série de fréquences, tandis que la transformée de Laplace résout une fonction en ses moments . En fait, une autre manière d'obtenir des moments consiste à utiliser une transformée de Fourier ( fonction caractéristique ). Le terme exponentiel dans la transformée de Laplace est en général de la forme est avecs=σ+iω , correspondant auxexponentielles réelles et sinusoïdales imaginaires, et donnant des tracés tels quecelui-ci:


[ Extrait du Guide du scientifique et de l'ingénieur sur le traitement du signal par Steven W. Smith ]


Par conséquent, la fonction MX(t) décompose le pdf quelque sorte en ses "fréquences constituantes" lorsque σ=0. D'après l'éq. (4):

MX(t)=E[esX]=esxpdf(x)dx=e(σ+iω)xpdf(x)dx=eσxeiωxpdf(x)dx

ce qui nous laisse avec l'intégrale impropre de la partie de l'expression en rouge, correspondant à la transformée de Fourier du pdf.

En général, l' intuition des pôles de transformation de Laplace d'une fonction serait qu'ils fournissent des informations sur les composantes exponentielle (décroissance) et fréquentielle de la fonction (dans ce cas, le pdf).


Xketx(1)(2)k(2)k(3)0


2
E[etX]E[Xk]

2
Je souhaite que les profanes qui comprennent cette réponse soient mes élèves :)
Aksakal

3

M(t)=etμ+1/2σ2t2
ddtM(t)|t=0=μ+σ2t|t=0=μ

De plus, puisque ce MGF code tout sur la distribution, si vous savez manipuler la fonction, vous pouvez appliquer des opérations sur toutes les caractéristiques de la distribution à la fois! Pourquoi n'utilisons-nous pas toujours MGF? Tout d'abord, ce n'est pas dans toutes les situations que le MGF est l'outil le plus simple. Deuxièmement, MGF n'existe pas toujours.

Au-dessus du profane

f(x)=12πex2/2

Vous pouvez calculer son moment tel que la moyenne et l'écart-type, et l'utiliser sur des variables transformées et des fonctions sur des normales aléatoires, etc.

Vous pouvez considérer le MGF de la distribution normale comme une alternative au PDF. Il contient la même quantité d'informations. J'ai déjà montré comment obtenir la moyenne.

σ2=x212πex2/2dx=?
M(t)=et2/2
σ2=d2dt2M(t)|t=0=ddtt|t=0=1


1
Pouvez-vous développer le "tout" qu'il code sur la distribution?
ColorStatistics

2
Pour apprécier l'argument avancé par @ColorStatistics, veuillez consulter stats.stackexchange.com/questions/25010 .
whuber

@whuber: Merci, whuber. J'étudierai cette référence. C'est un sujet que je cherche à mieux comprendre.
ColorStatistics

Comment prouver que MGF & PDF contiennent la même quantité d'informations?
Aerin
En utilisant notre site, vous reconnaissez avoir lu et compris notre politique liée aux cookies et notre politique de confidentialité.
Licensed under cc by-sa 3.0 with attribution required.