Supposons que

18

Quelle est la façon la plus simple de vérifier que l'énoncé suivant est vrai?

Supposons que . Afficher .Y1,,YniidExp(1)i=1n(YiY(1))Gamma(n1,1)

Notez que .Y(1)=min1inYi

Par XExp(β) , cela signifie que fX(x)=1βex/β1{x>0} .

Il est facile de voir que Oui(1)Exponentiel(1/n) . De plus, nous avons également i=1nYiGamma(α=n,β=1) sous le paramétrage

fY(y)=1Γ(α)βαxα1ex/β1{x>0}α,β>0.

Solution donnée par la réponse de Xi'an : en utilisant la notation dans la question d'origine:

i=1n[YiY(1)]=i=1n[Y(i)Y(1)]=i=1nY(i)nY(1)=i=1n{Y(i)Y(i1)+Y(i1)Y(1)+Y(1)}nY(1)=i=1nj=1i{Y(j)Y(j1)}nY(1) where Y(0)=0=j=1ni=jn{Y(j)Y(j1)}nY(1)=j=1n(nj+1)[Y(j)Y(j1)]-nOui(1)=je=1n(n-je+1)[Oui(je)-Oui(je-1)]-nOui(1)=je=2n(n-je+1)[Oui(je)-Oui(je-1)]+nOui(1)-nOui(1)=je=2n(n-je+1)[Oui(je)-Oui(je-1)].
De cela, nous obtenons que i=2n(ni+1)[Y(i)Y(i1)]Gamma(n1,1)1) .
Clarinettiste
la source
1
@MichaelChernick 1) Je ne sais pas si ma preuve d'indépendance est correcte, et 2) je ne sais pas si le résultat que j'ai deviné ci-dessus concernant la différence des distributions Gamma est même correct. Cela semble contredire ce qui est donné ici , mais peut-être que cette situation est différente puisque cette différence inclut l'une des statistiques de commande? Je ne suis pas sûr.
Clarinettiste
1
@ Clarinettiste, je n'en suis pas certain. Essayez peut-être de travailler avec , ce qui équivaut clairement à la somme avec laquelle vous travaillez. La réponse pourrait être utile ici: math.stackexchange.com/questions/80475/...je=2n(Oui(je)-Oui(1))
gammer
3
Avez-vous essayé de prouver que chacun - sauf pour un , pour lequel , puis, en utilisant le fait que la somme de iid Variantes exponentielles sera distribuée Gamma? i Y i - Y ( 1 ) = 0 ( n - 1 )(YiY(1))Expon(1)iYiY(1)=0(n1)
Marcelo Ventura
1
@jbowman Nous avons et conditionnés à , nous divisons cela par , donnant , donc nous avons . Maintenant, voici ce qui m'a dérangé sur cette preuve: je considérais comme une constante. Mais n'est pas une constante. Pourquoi cela fonctionnerait-il? Y ia e - a e - z i ( Z iY iA ) Exp ( 1 ) A Y ( 1 )
fZi(zi)=fYi(zi+a)=e(zi+a)
Yiaeaezi(ZiYiA)Exp(1)AY(1)
Clarinettiste
1
Le fait est que peu importe ce c'est! La distribution est toujours ! Remarquable, non? Et à partir de cela, vous pouvez conclure que la distribution de est toujours pour , quelle que soit la valeur réelle de . Exp ( 1 ) Y i - Y [ 1 ] Exp ( 1 ) i > 1 Y [ 1 ]uneExp(1)Ouije-Oui[1]Exp(1)je>1Oui[1]
jbowman

Réponses:

15

La preuve est donnée dans Mother of All Random Generation Books, Devroye's Non-uniform Random Variate Generation , p.211 (et c'est très élégant!):

Théorème 2.3 (Sukhatme, 1937) Si nous définissons alors les espacements exponentiels normalisés dérivés de statistiques de commande d'un échantillon exponentiel iid de taille sont eux-mêmes des variables exponentielles iid( n - i + 1 ) ( E ( i ) - E ( i - 1 ) ) E ( 1 )E ( n ) nE(0)=0

(n-je+1)(E(je)-E(je-1))
E(1)E(n)n

Preuve. Depuis le la densité conjointe de la statistique d'ordre s'écrit comme Réglage , le changement de variables de à a une constante jacobienne [accessoirement égale àmais cela n'a pas besoin d'être calculé] et donc la densité de (E(1),,E(n

i=1nei=i=1ne(i)=i=1nj=1i(e(j)e(j1))=j=1ni=jn(e(j)e(j1))=j=1n(nj+1)(e(j)e(j1))
f( e )=n!(E(1),,E(n)) Y i = ( E ( i ) - E ( i - 1 ) ) ( E ( 1 ) , , E ( n ) ) ( Y 1
f(e)=n!exp{i=1ne(i)}=n!exp{i=1n(ni+1)(e(i)e(i1))}
Yi=(E(i)E(i1))(E(1),,E(n))1 / n ! ( Y 1 , , Y n )(Y1,,Yn)1/n!(Y1,,Yn) est proportionnel à qui établit le résultat. QED
exp{i=1nyi}

Une alternative proposée par Gérard Letac est de vérifier que a la même distribution que (en vertu de la propriété sans mémoire), ce qui rend la dérivation de simple.( E 1

(E(1),,E(n))
n k=1(Ek-E(1)) n - 1 k=1Ek
(E1n,E1n+E2n1,,E1n+E2n1++En1)
k=1n(EkE(1))k=1n1Ek
Xi'an
la source
1
Merci pour cette réponse! Je voudrais compléter certains détails pour quiconque lira ceci à l'avenir: sont des valeurs observées de , et la façon la plus simple de voir que consiste à écrire terme par terme. Étant donné que la densité de est proportionnelle à , séparez pour voir que la densité est proportionnelle à , d'où .E i n i = 1 e ( i ) = n i = 1 ( n - i + 1 ) ( e ( i ) - e ( i - 1 ) ) = n i = 1 e ( i )n i = 1 ( n - i + 1 )eiEii=1ne(i)=i=1n(ni+1)(e(i)e(i1))=i=1ne(i)( Y 1 , , Y n ) exp ( - n i = 1 y i ) y i n i = 1 e - y i Y 1 , , Y n iid Exp ( 1 )i=1n(ni+1)(e(i)e(i1))(Y1,,Yn)exp(i=1nyi)yii=1neyiY1,,YniidExp(1)
Clarinettiste
5

Je présente ici ce qui a été suggéré dans les commentaires de @jbowman.

Soit une constante . Soit suivre un et considérons . alorsune0OuijeExp (1)Zje=Ouije-une

Pr(ZjezjeOuijeune)=Pr(Ouije-unezjeOuijeune)

Pr(Ouijezje+uneOuijeune)=Pr(Ouijezje+une,Ouijeune)1-Pr(Ouijeune)

Pr(uneOuijezje+une)1-Pr(Ouijeune)=1-e-zje-une-1+e-unee-une=1-e-zje

qui est la fonction de distribution de .Exp (1)

Décrivons ceci: la probabilité qu'un rv tombe dans un intervalle spécifique (le numérateur sur la dernière ligne), étant donné qu'il dépassera la borne inférieure de l'intervalle (le dénominateur), dépend uniquement de la la longueur de l'intervalle et non sur l'endroit où cet intervalle est placé sur la ligne réelle. Exp (1)Ceci est une incarnation de la propriété " sans mémoire " de la distribution exponentielle, ici dans un cadre plus général, libre d'interprétations temporelles (et cela vaut pour la distribution exponentielle en général)

Maintenant, en conditionnant sur on force à être négative, et surtout, le résultat obtenu est titulaire . Nous pouvons donc affirmer ce qui suit: {Ouijeune}ZjeuneR+

Si , alors . OuijeExp (1)Q0:Zje=Ouije-Q0 ZjeExp (1)

Pouvons-nous trouver un libre de prendre toutes les valeurs réelles non négatives et pour lequel l'inégalité requise tient toujours (presque sûrement)? Si nous le pouvons, alors nous pouvons nous passer de l'argument du conditionnement. Q0

Et en effet, nous le pouvons. C'est la statistique d'ordre minimum , , . Nous avons donc obtenuQ=Oui(1)Pr(OuijeOui(1))=1

OuijeExp (1)Ouije-Oui(1)Exp (1)

Cela signifie que

Pr(Ouije-Oui(1)yje-y(1))=Pr(Ouijeyje)

Donc, si la structure probabiliste de reste inchangée si l'on soustrait la statistique d'ordre minimum, il s'ensuit que les variables aléatoires et où indépendant, sont également indépendants puisque le lien possible entre eux, n'a pas d'effet sur la structure probabiliste.OuijeZje=Ouije-Oui(1)Zj=Ouij-Oui(1)Ouije,OuijOui(1)

Alors la somme contient iid variables aléatoires (et un zéro), et ainsije=1n(Ouije-Oui(1))n-1 Exp (1)

je=1n(Ouije-Oui(1))Gamma(n-1,1)
Alecos Papadopoulos
la source