De la distribution uniforme à la distribution exponentielle et vice-versa


20

Ceci est probablement une question triviale, mais ma recherche a été infructueuse jusqu'à présent, y compris cet article wikipedia , et le « Recueil des distributions » le document .

Si a une distribution uniforme, cela signifie-t-il que suit une distribution exponentielle?e XXeX

De même, si suit une distribution exponentielle, cela signifie-t-il que suit une distribution uniforme?l n ( Y )Yln(Y)


3
Pourquoi vous attendriez-vous à ce qu'il en soit ainsi? À cause du nom? Consultez en.wikipedia.org/wiki/… pour voir comment les autres distributions sont liées à l'exponentielle. Aussi ...exp(X)[0,)
Tim

Non, je pense que je suis en train de suivre les analogies avec les transformations de fonctions standard, en oubliant qu'avec les distributions, les choses sont différentes.
luchonacho

Réponses:


25

Ce n'est pas le cas que l'exponentialisation d'une variable aléatoire uniforme donne une exponentielle, ni la prise du log d'une variable aléatoire exponentielle ne donne un uniforme.

Soit uniforme sur et soit .( 0 , 1 ) X = exp ( U )U(0,1)X=exp(U)

FX(x)=P(Xx)=P(exp(U)x)=P(Ulnx)=lnx,1<x<e

Donc .fx(x)=ddxlnx=1x,1<x<e

Ce n'est pas une variable exponentielle. Un calcul similaire montre que le log d'une exponentielle n'est pas uniforme.

Soit une exponentielle standard, donc .F Y ( y ) = P ( Y y ) = 1 - e - yYFOui(y)=P(Ouiy)=1-e-y,y>0

Soit . Alors .F V ( v ) = P ( V v ) = P ( lnV=lnOuiFV(v)=P(Vv)=P(lnOuiv)=P(Ouiev)=1-e-ev,v<0

Ce n'est pas un uniforme. (En effet, est une variable aléatoire distribuée par Gumbel , vous pouvez donc appeler la distribution de un Gumbel inversé .)V-VV

Cependant, dans chaque cas, nous pouvons le voir plus rapidement en considérant simplement les limites des variables aléatoires. Si est uniforme (0,1), il se situe entre 0 et 1, donc se situe entre et ... donc ce n'est pas exponentiel. De même, pour exponentielle, est sur , donc cela ne peut pas être uniforme (0,1), ni même aucun autre uniforme.X = exp ( U ) 1 e Y ln Y ( - , )UX=exp(U)1eOuilnOui(-,)

Nous pourrions également simuler, et le voir à nouveau immédiatement:

Tout d'abord, exponentier un uniforme -

histogramme d'uniforme exponentialisé avec la densité théorique superposée

[la courbe bleue est la densité (1 / x sur l'intervalle indiqué) que nous avons élaborée ci-dessus ...]

Deuxièmement, le journal d'un exponentiel:

histogramme du log d'une variable exponentielle

Ce que nous pouvons voir est loin d'être uniforme! (Si nous différencions le cdf que nous avons élaboré auparavant, ce qui donnerait la densité, il correspond à la forme que nous voyons ici.)

En effet, la méthode cdf inverse indique que prendre le négatif du log d'un varié uniforme (0,1) donne un varié exponentiel standard, et inversement, exponentier le négatif d'un exponentiel standard donne un uniforme. [Voir aussi transformation intégrale de probabilité ]

Cette méthode nous dit que si , . Si nous appliquons l'inverse du cdf comme transformation sur , un uniforme standard, la variable aléatoire résultante a la fonction de distribution .Y =U=FOui(Oui)Oui=F-1(U)UFOui

Si nous laissons uniforme (0,1), alors . Soit . (Notez que est également uniforme sur (0,1), vous pouvez donc laisser , mais nous suivons la méthode cdf inverse en entier ici)P ( U u ) = u Y = - ln ( 1 - UUP(Uu)=uOui=-ln(1-U)1-UOui=-lnU

Alors , qui est le cdf d'une exponentielle standard.P(Ouiy)=P(-ln(1-U)y)=P(1-Ue-y)=P(U1-e-y)=1-e-y

[Cette propriété de la transformation cdf inverse est la raison pour laquelle la transformation est réellement requise pour obtenir une distribution exponentielle, et la transformation intégrale de probabilité est la raison pour laquelle l'exponentiation du négatif d'une exponentielle négative revient à un uniforme.]Journal


Très bonne réponse! Merci. Je le vois maintenant. J'ai calculé le CDF dans les deux cas, et j'obtiens le négatif du log dans le premier cas et la valeur absolue d'un inverse, dans le second. Je pense que ma confusion est sur la pensée en termes de transformations de fonctions standard, qui ne suit pas quand il s'agit de distributions. +1 pour les graphiques!
luchonacho

6

Vous l'avez presque en arrière. Tu as demandé:

  • XeX

  • Ouiln(Oui)

En réalité

  • X[0,1]-Journale(X)1
  • Oui1e-Oui[0,1]

Plus généralement, vous pourriez dire:

  • X[une,b]-1kJournale(X-uneb-une)k
  • Ouike-kOui[0,1]une+(b-une)e-kOui[une,b]
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.