Supongamos que


18

¿Cuál es la forma más fácil de ver que la siguiente afirmación es verdadera?

Supongamos que Y1,,YniidExp(1) . Mostrar i=1n(YiY(1))Gamma(n1,1) .

Y(1)=min1inYi

Por , esto significa que f_ {X} (x) = \ dfrac {1} {\ beta} e ^ {- x / \ beta} \ cdot \ mathbf {1} _ {\ {x> 0 \}} .XExp(β)fX(x)=1βex/β1{x>0}

Es fácil ver que Y(1)Exponential(1/n) . Además, también tenemos que i=1nYiGamma(α=n,β=1) bajo la parametrización

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

Solución dada la respuesta de Xi'an : usando la notación en la pregunta original:

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)]nY(1)=i=1n(ni+1)[Y(i)Y(i1)]nY(1)=i=2n(ni+1)[Y(i)Y(i1)]+nY(1)nY(1)=i=2n(ni+1)[Y(i)Y(i1)].
De esto, obtenemos que i=2n(ni+1)[Y(i)Y(i1)]Gamma(n1,1)1) .

1
@MichaelChernick 1) No estoy seguro de si mi prueba de independencia es correcta, y 2) No estoy seguro de si el resultado que supuse anteriormente sobre la diferencia de las distribuciones Gamma es incluso correcto. Esto parece contradecir lo que se da aquí , pero ¿tal vez esta situación es diferente ya que esta diferencia incluye una de las estadísticas del pedido? No estoy seguro.
Clarinetista

1
@ Clarinetista, no estoy seguro. Tal vez intente trabajar con , que claramente es igual a la suma con la que está trabajando. La respuesta aquí puede ser útil: math.stackexchange.com/questions/80475/…i=2n(Y(i)Y(1))
jugador

3
¿Has intentado demostrar que cada - excepto uno , para el cual y, luego, usando el hecho de que la suma de iid Variables exponenciales se distribuirá Gamma? i Y i - Y ( 1 ) = 0 ( n - 1 )(YiY(1))Expon(1)iYiY(1)=0(n1)
Marcelo Ventura el

1
@jbowman Tenemos y condicionado a , lo dividimos por , dando , por lo tanto tenemos . Ahora, esto es lo que me molestó acerca de esta prueba: consideraba como una constante. Pero no es una constante. ¿Por qué funcionaría esto? Y ia e - a e - z i ( Z iY iA ) Exp ( 1 ) A Y ( 1 )
FZyo(zyo)=FYyo(zyo+un)=mi-(zyo+un)
Yyounmi-unmi-zyo(ZyoYyoUN)Exp(1)UNY(1)
Clarinetista

1
El punto es que no importa lo que es! ¡La distribución es siempre ! Notable, ¿no es así? Y de esto puede concluir que la distribución de es siempre para , independientemente del valor real de . Exp ( 1 ) Y i - Y [ 1 ] Exp ( 1 ) i > 1 Y [ 1 ]unExp(1)Yyo-Y[1]Exp(1)yo>1Y[1]
jbowman

Respuestas:


15

La prueba se da en la Madre de Todos los Libros de Generación Aleatoria, Generación de Variación Aleatoria No Uniforme de Devroye , en la p.211 (¡y es muy elegante!):

Teorema 2.3 (Sukhatme, 1937) Si definimos entonces los espacios exponenciales normalizados derivan de las estadísticas de orden de una muestra exponencial iid de tamaño son en sí mismas variables exponenciales iid( n - i + 1 ) ( E ( i ) - E ( i - 1 ) ) E ( 1 )E ( n ) nE(0)=0

(ni+1)(E(i)E(i1))
E(1)E(n)n

Prueba. Desde la la densidad conjunta de la estadística de orden escribe como Configurando , el cambio de variables de a tiene un constante jacobiano [¡casualmente igual apero esto no necesita ser calculado] y, por lo tanto, la densidad 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!(mi(1),...,mi(norte)) Y i = ( E ( i ) - E ( i - 1 ) ) ( E ( , , Y n ) 1 /
F(mi)=norte!Exp{-yo=1nortemi(yo)}=norte!Exp{-yo=1norte(norte-yo+1)(mi(yo)-mi(yo-1))}
Yyo=(mi(yo)-mi(yo-1))( Y 1(mi(1),...,mi(norte))(Y1,...,Ynorte)( S 1 , ... , S n )1/ /norte!(Y1,...,Ynorte) es proporcional a que establece el resultado. QED
Exp{-yo=1norteyyo}

Una alternativa que me sugirió Gérard Letac es verificar que tenga la misma distribución que (en virtud de la propiedad sin memoria), que hace la derivación de sencillo.( E 1

(mi(1),...,mi(norte))
n k=1(Ek-E(1)) n - 1 k=1Ek
(mi1norte,mi1norte+mi2norte-1,...,mi1norte+mi2norte-1+...+minorte1)
k=1norte(mik-mi(1))k=1norte-1mik

1
Gracias por esta respuesta! Me gustaría completar algunos detalles para cualquiera que lea esto en el futuro: son valores observados de , y la forma más fácil de ver que es escribir término por término. Debido a que la densidad de es proporcional a , separe el para ver que la densidad es proporcional a , por tanto .miyomiyoyo=1nortemi(yo)=yo=1norte(norte-yo+1)(mi(yo)-mi(yo-1))=yo=1nortemi(yo)yo=1norte(norte-yo+1)(mi(yo)-mi(yo-1))(Y1,...,Ynorte)Exp(-yo=1norteyyo)yyoyo=1nortemi-yyoY1,...,YnorteiidExp(1)
Clarinetista

5

Presento aquí lo que se ha sugerido en los comentarios de @jbowman.

Deje una constante . Deje que siga un y considere . Luegoun0 0YyoExp (1)Zyo=Yyo-un

Pr(ZyozyoYyoun)=Pr(Yyo-unzyoYyoun)

Pr(Yyozyo+unYyoun)=Pr(Yyozyo+un,Yyoun)1-Pr(Yyoun)

Pr(unYyozyo+un)1-Pr(Yyoun)=1-mi-zyo-un-1+mi-unmi-un=1-mi-zyo

cual es la función de distribución de .Exp (1)

Describamos esto: la probabilidad de que un rv caiga en un intervalo específico (el numerador en la última línea), dado que excederá el límite inferior del intervalo (el denominador), depende solo del longitud del intervalo y no en donde este intervalo se coloca en la línea real. Exp (1)Esta es una encarnación de la propiedad " falta de memoria " de la distribución exponencial, aquí en un entorno más general, libre de interpretaciones de tiempo (y se aplica a la distribución exponencial en general)

Ahora, al condicionar , a a ser no negativo y, de manera crucial, el resultado obtenido se mantiene . Entonces podemos decir lo siguiente: {Yyoun}ZyounR+

Si , entonces . YyoExp (1)Q0 0:Zyo=Yyo-Q0 0 ZyoExp (1)

¿Podemos encontrar un que sea libre de tomar todos los valores reales no negativos y para el cual la desigualdad requerida siempre se mantenga (casi seguramente)? Si podemos, entonces podemos prescindir del argumento condicionante. Q0 0

Y de hecho podemos. Es la estadística de orden mínimo , , . Entonces hemos obtenidoQ=Y(1)Pr(YyoY(1))=1

YyoExp (1)Yyo-Y(1)Exp (1)

Esto significa que

Pr(Yyo-Y(1)yyo-y(1))=Pr(Yyoyyo)

Entonces, si la estructura probabilística de permanece sin cambios si restamos la estadística de orden mínimo, se deduce que las variables aleatorias y donde independientes, también son independientes ya que el posible vínculo entre ellos, no tiene un efecto en la estructura probabilística.YyoZyo=Yyo-Y(1)Zj=Yj-Y(1)Yyo,YjY(1)

Entonces la suma contiene iid variables aleatorias (y un cero), y asíyo=1norte(Yyo-Y(1))norte-1 Exp (1)

yo=1norte(Yyo-Y(1))Gama(norte-1,1)
Al usar nuestro sitio, usted reconoce que ha leído y comprende nuestra Política de Cookies y Política de Privacidad.
Licensed under cc by-sa 3.0 with attribution required.