CDF élevé au pouvoir?


15

Si est un CDF, il semble que ( ) soit également un CDF.FZFZ(z)αα>0

Q: Est-ce un résultat standard?

Q: Existe-t-il un bon moyen de trouver une fonction avec X \ equiv g (Z) st F_X (x) = F_Z (z) ^ \ alpha , où x \ equiv g (z)gXg(Z)FX(x)=FZ(z)αxg(z)

Fondamentalement, j'ai un autre CDF en main, FZ(z)α . Dans un sens réduit, j'aimerais caractériser la variable aléatoire qui produit ce CDF.

EDIT: Je serais heureux si je pouvais obtenir un résultat analytique pour le cas spécial . Ou du moins sachez qu'un tel résultat est insoluble.ZN(0,1)


2
Oui, c'est un résultat assez connu et facile à généraliser. (Comment?) Vous pouvez également trouver , au moins implicitement. C'est essentiellement une application de la technique de transformation inverse probablement couramment utilisée pour générer des variables aléatoires d'une distribution arbitraire. g
cardinal

2
@cardinal Veuillez répondre. L'équipe se plaint plus tard que nous ne nous battons pas avec un faible taux de réponse.

1
@mbq: Merci pour vos commentaires, que je comprends et respecte grandement. Veuillez comprendre que parfois des considérations de temps et / ou de lieu ne me permettent pas de poster une réponse, mais permettent un commentaire rapide qui peut amener le PO ou d'autres participants à démarrer. Soyez assuré qu'à l'avenir, si je peux poster une réponse, je le ferai. J'espère que ma participation continue par le biais de commentaires sera également satisfaisante.
cardinal

2
@cardinal Certains d'entre nous sont également coupables des mêmes, pour les mêmes raisons ...
whuber

2
@brianjd Oui, c'est un résultat bien connu qui a été utilisé pour produire industriellement des distributions "généralisées", voir . Il existe de nombreuses transformations comme celle-ci et les gens les utilisent à cet effet: ils trouvent une transformation paramétrique, l'appliquent à une distribution et voilá, vous avez un papier juste en calculant simplement ses propriétés. Et bien sûr, la normale est la première «victime».

Réponses:


11

J'aime les autres réponses, mais personne n'a encore mentionné ce qui suit. L'événement {Ut, Vt} se produit si et seulement si {max(U,V)t} , donc si U et V sont indépendants et W=max(U,V) , alors FW(t)=FU(t)FV(t) de manière àα un nombre entier positif (exemple,α=n ) prendreX=max(Z1,...Zn) où lesZ s » sont iid

Pour α=1/n nous pouvons changer pour obtenir FZ=FXn , donc X serait cette variable aléatoire telle que le maximum de n copies indépendantes a la même distribution que Z (et ce ne serait pas l'un de nos amis familiers , en général).

Le cas de α un nombre rationnel positif (disons, ) découle du précédent depuis α=m/n

(FZ)m/n=(FZ1/n)m.

Pour an irrationnel, choisissez une séquence de rationnels positifs convergeant vers ; alors la séquence (où nous pouvons utiliser nos astuces ci-dessus pour chaque ) convergera en distribution vers le souhaité.αakαXkkX

Ce n'est peut-être pas la caractérisation que vous recherchez, mais cela donne au moins une idée de la façon de penser pour convenablement. D'un autre côté, je ne sais pas vraiment à quel point cela peut vraiment être plus agréable: vous avez déjà le CDF, donc la règle de chaîne vous donne le PDF, et vous pouvez calculer les moments jusqu'au coucher du soleil ...? Il est vrai que la plupart des n'auront pas de familier pour , mais si je voulais jouer avec un exemple pour chercher quelque chose d'intéressant, je pourrais essayer distribué uniformément sur l'intervalle d'unité avec , .FZααZXα=2ZF(z)=z0<z<1


EDIT: J'ai écrit quelques commentaires dans la réponse @JMS, et il y avait une question sur mon arithmétique, donc je vais écrire ce que je voulais dire dans l'espoir que ce soit plus clair.

@cardinal correctement dans le commentaire à la réponse @JMS a écrit que le problème se simplifie en ou plus généralement lorsque n'est pas nécessairement , nous avons Mon point était que lorsque a une belle fonction inverse, nous pouvons simplement résoudre la fonction avec l'algèbre de base. J'ai écrit dans le commentaire que devrait être

g1(y)=Φ1(Φα(y)),
ZN(0,1)
x=g1(y)=F1(Fα(y)).
Fy=g(x)g
y=g(x)=F1(F1/α(x)).

Prenons un cas spécial, branchons les choses et voyons comment cela fonctionne. Soit une distribution Exp (1), avec CDF et CDF inverse Il est facile de tout brancher pour trouver ; une fois que nous avons terminé, nous obtenons Donc, en résumé , mon affirmation est que si et si nous définissons alors aura un CDF qui ressemble à Nous pouvons le prouver directement (regardezX

F(x)=(1ex), x>0,
F1(y)=ln(1y).
g
y=g(x)=ln(1(1ex)1/α)
XExp(1)
Y=ln(1(1eX)1/α),
Y
FY(y)=(1ey)α.
P(Yy) et utiliser l'algèbre pour obtenir l'expression, dans l'avant-dernière étape, nous avons besoin de la transformation intégrale de probabilité). Juste dans le cas (souvent répété) où je suis fou, j'ai fait quelques simulations pour vérifier que cela fonctionne, ... et c'est le cas. Voir ci-dessous. Pour faciliter le code, j'ai utilisé deux faits:
If XF then U=F(X)Unif(0,1).
If UUnif(0,1) then U1/αBeta(α,1).

Le tracé des résultats de la simulation suit.

ECDF and F to the alpha

Le code R utilisé pour générer le tracé (moins les étiquettes) est

n <- 10000; alpha <- 0.7
z <- rbeta(n, shape1 = alpha, shape2 = 1)
y <- -log(1 - z)
plot(ecdf(y))
f <- function(x) (pexp(x, rate = 1))^alpha
curve(f, add = TRUE, lty = 2, lwd = 2)

La coupe a l'air plutôt bien, je pense? Peut-être que je ne suis pas fou (cette fois)?


Voir mon commentaire dans la réponse @JMS. Pour la réponse est qui n'est pas fermée mais peut être calculé assez facilement. Et vous pouvez le rendre plus facile en reconnaissant que l'entrée du CDF inverse est une distribution bêta convenablement choisie. La réponse sera agréable dans les cas où le CDF inverse est agréable, et il y en a qui courent. ZN(0,1)g(z)=Φ1(Φ1/α(z))

Il serait bon de revérifier votre arithmétique.
cardinal

@cardinal errr ... OK, je l'ai fait, ... et c'est vrai? Pourriez-vous s'il vous plaît signaler l'erreur?

(+1) Toutes mes excuses. Je ne sais pas où était ma tête lorsque j'ai regardé cela pour la première fois. C'est évidemment (eh bien, aurait dû l'être!) Correct.
cardinal

@cardinal, pas de mal, pas de faute. J'avoue, cependant, que vous m'avez vraiment fait transpirer pendant une minute! :-)

14

Preuve sans mots

enter image description here

La courbe bleue inférieure est , la courbe rouge supérieure est F α (en typant le cas α < 1 ), et les flèches indiquent comment passer de z à x = g ( z ) .FFαα<1zx=g(z)


Une belle photo! Q: Qu'est-ce qui a été attiré? TikZ?
lowndrul

1
@brianjd: Si je me souviens bien, @whuber fait beaucoup de ses intrigues en utilisant Mathematica.
cardinal

3
@cardinal Vous avez raison. En fait, j'utilise tout ce qui est pratique et semble pouvoir faire du bon travail rapidement. FWIW, voici le code:Module[ {y, w, a = 0.1, z = 3.24, f = ChiDistribution[7.6], xmin=0, xmax=5}, y = CDF[f,z]; w = InverseCDF[f, y^(1/a)]; Show[ Plot[{CDF[f, x],CDF[f,x]^a} , {x, xmin, xmax}, Filling->{1->{2}}], Graphics[{ Dashed, Arrow[{{z,0}, {z,y}}], Arrow[{{z,y}, {w,y}}], Arrow[{{w,y}, {w,0}}] }] ] ]
whuber

6

Q1) Oui. Il est également utile pour générer des variables qui sont ordonnées stochastiquement; vous pouvez le voir sur la jolie photo de @ whuber :). permute l'ordre stochastique.α>1

Que ce soit un cdf valide est juste une question de vérification des conditions requises: doit être cadlag , non décroissant et limité à 1 à l'infini et 0 à l'infini négatif. F z a ces propriétés donc elles sont toutes faciles à montrer.Fz(z)α10Fz

Q2) Il semble que ce serait assez difficile analytiquement, à moins que soit spécialFZ


@JMS: Et ? ZN(0,1)
lowndrul

2
@brianjd: Je ne le crois pas. Soit une fonction continue strictement monotone (donc, ayant un inverse bien défini g - 1 ) qui satisfait vos conditions. Alors, il faut que Φ α ( u ) = P ( g ( Z ) u ) = P ( Z g - 1 ( u ) ) = Φ ( g - 1 ( u ) ) et donc g - 1gg1Φα(u)=P(g(Z)u)=P(Zg1(u))=Φ(g1(u)) . L'inverse est donc identifié de manière assez explicite, mais pas g lui-même. C'est ce que je voulais dire dans mon commentaire précédent à propos de ladécouverteimplicite de g . g1(u)=Φ1(Φα(u))gg
cardinal

@brianjd - Ce que @cardinal a dit :) Je ne pouvais même pas penser à un cas spécial pour où vous obtiendriez un formulaire fermé (pour ne pas dire qu'il n'y en a pas bien sûr). FZ
JMS

@JMS: serait un exemple positif. ZU[0,1]
cardinal

@cardinal Je n'aurais jamais pensé à une distribution aussi rare ... mais maintenant que vous le mentionnez, un devrait fonctionner en général, vous donnant en retour un B e t a ( a α , 1 ) . Beta(a,1)Beta(aα,1)
JMS
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.