Probar la transformación integral de probabilidad sin suponer que el CDF está aumentando estrictamente


9

Sé que la prueba de la transformación integral de probabilidad se ha dado varias veces en este sitio. Sin embargo, las pruebas que encontré utilizan la hipótesis de que el CDF está aumentando estrictamente (junto, por supuesto, con la hipótesis de que es una variable aleatoria continua). Sé que en realidad la única hipótesis requerida es que es una variable aleatoria continua, y no se requiere una monotonicidad estricta. ¿Me puedes mostrar cómo?FX(X)XX

Como ya estoy aquí, también aprovecho la ocasión para pedir una aplicación simple de la transformación integral de probabilidad :) ¿puede mostrarme que, si tiene CDF e es el truncamiento de a , entonces se distribuye como donde ?XFX(X)YX[una,si]YFX-1(U)U[FX(una),FX(si)]


2
si fuera tan amable, en la prueba de su enlace, podría señalar dónde está el requisito de que FX(X)tiene que estar estrictamente en aumento. ¡Gracias!
Erosennin

1
@Erosennin, la prueba supone la existencia de la inversa de FX(X).
DeltaIV

¡Gracias! Pero, ¿hay alguna vez un CDF que no esté aumentando estrictamente? Probablemente ya hayas pensado en esto, sin embargo ...
Erosennin

2
Por supuesto que lo hay. La variable aleatoria cuyo pdf es igual a 1/2 en [0,0.5], 0 en [0.5,1] y 1/2 en [1,1.5], tiene un CDF que es continuo, pero no aumenta estrictamente.
DeltaIV

La parte difícil es tratar con la parte no absolutamente continua de F. La idea se aclara al considerar el caso extremo de discretoF. En stats.stackexchange.com/a/36246/919 doy un algoritmo que implementa la transformación integral de probabilidad en ese caso (además de proporcionar un código de trabajo). Emulando ese algoritmo para arbitrarioFcontestará tu pregunta
whuber

Respuestas:


9

En el enlace de wikipedia proporcionado por el OP, la transformación integral de probabilidad en el caso univariante se da de la siguiente manera

Supongamos que una variable aleatoria X tiene una distribución continua para la cual la función de distribución acumulativa (CDF) es FX. Entonces la variable aleatoriaY=FX(X)Tiene una distribución uniforme.
PRUEBA
Dada cualquier variable aleatoriaX, definir Y=FX(X). Entonces:

FY(y)=Problema(Yy)=Problema(FX(X)y)=Problema(XFX-1(y))=FX(FX-1(y))=y

FY es solo el CDF de un UnorteyoFormetro(0 0,1)variable aleatoria. Así,Y tiene una distribución uniforme en el intervalo [0 0,1].

El problema con lo anterior es que no está claro cuál es el símbolo FX-1representa. Si representara el inverso "habitual" (que existe solo para biyecciones), entonces la prueba anterior se mantendría solo para CDF continuos y estrictamente crecientes. Pero este no es el caso, ya que para cualquier CDF trabajamos con la función cuantil (que es esencialmente una inversa generalizada),

FZ-1(t)inf{z:FZ(z)t},t(0 0,1)

Según esta definición, la serie de igualdades de Wikipedia sigue siendo válida para los CDF continuos. La igualdad crítica es

Problema(XFX-1(y))=Problema(Xinf{X:FX(X)y})=Problema(FX(X)y)

que se cumple porque estamos examinando un CDF continuo. Esto en la práctica significa que su gráfica es continua (y sin partes verticales, ya que es una función y no una correspondencia). A su vez, esto implica que el infimum (el valor de inf {...}), lo denotaX(y), siempre será tal que FX(X(y))=y. El resto es inmediato.

Con respecto a los CDF de distribuciones discretas (o mixtas), no es (no puede ser) cierto queY=FX(X) sigue un uniforme U(0 0,1), pero aún es cierto que la variable aleatoria Z=FX-1(U) tiene función de distribución FX(por lo que el muestreo de transformación inversa todavía se puede utilizar). Una prueba se puede encontrar en Shorack, GR (2000). Probabilidad para los estadísticos . ch.7 .


3
+1 También se proporciona una prueba similar en la pág. 54 de Inferencia estadística de Casella y Berger, segunda edición.
StatsStudent

1
@ Analista1 Gracias, es bueno tener múltiples referencias.
Alecos Papadopoulos
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.