En estadística inferencial, específicamente en inferencia predictiva, un intervalo de predicción es una estimación de un intervalo de valores en el que se producirá una observación futura con determinada probabilidad, dado lo que ya se ha observado. Los intervalos de predicción se utilizan a menudo en análisis de la regresión.
Se utilizan tanto en estadística frecuentista como en estadística bayesiana: un intervalo de predicción guarda la misma relación con una observación futura que un intervalo de confianza frecuentista o intervalo de credibilidad bayesiano muestran con relación a un parámetro de población no observable: los intervalos de predicción estiman la distribución de futuros sucesos individuales según su valor, mientras que los intervalos de confianza y de credibilidad predicen la distribución de la media real de la población u otra cantidad de interés que no puede ser observada directamente (por ejemplo, cuando solo se poseen muestras de poblaciones muy grandes).
Introducción Por ejemplo, si se hace la suposición paramétrica de que la distribución subyacente de un determinado fenómeno observable es una distribución normal, y se dispone de un conjunto de muestras {X1, ..., Xn}, entonces los intervalos de confianza y credibilidad se pueden usar para estimar la media μ y la desviación típica σ de la población subyacente, mientras que los intervalos de predicción se pueden usar para estimar el valor de la siguiente variable de la muestra, Xn+1.
Alternativamente, en términos Bayesianos, un intervalo de predicción se puede describir como un intervalo de credibilidad para una variable por sí misma, en lugar de servir para estimar un parámetro de la distribución de la variable.
El concepto de intervalos de predicción no debe estar restringido a la inferencia sobre un determimado valor futuro de una muestra, sino que puede extenderse a casos más complicados. Por ejemplo, en el contexto de las inundaciones fluviales (donde los análisis a menudo se basan en los valores del mayor caudal de cada año), puede haber interés en hacer inferencias sobre la mayor inundación que probablemente se experimentará en los próximos 50 años.
Dado que los intervalos de predicción solo se refieren a observaciones pasadas y futuras, en lugar de a parámetros de población no observables, algunos estadísticos los recomiendan como un método mejor que los intervalos de confianza, como Seymour Geisser, siguiendo el enfoque sobre observables de Bruno de Finetti.
Distribución normal Dada una muestra con distribución normal, cuyos parámetros son desconocidos, es posible dar intervalos de predicción en el sentido frecuentista, es decir, determinar un intervalo [a, b] basado en estadísticas de la muestra mediante experimentos repetidos u observaciones, de forma que el evento Xn+1 se situará en el citado intervalo durante un determinado porcentaje del número de veces que se observe; recibiendo el nombre de "intervalos de confianza predictivos". Chapter 2: Non-Bayesian predictive approaches
Una técnica general de los intervalos de predicción frecuentista es encontrar y calcular una cantidad pivote de los observables X1, ..., Xn, Xn+1 (lo que significa que se establece una función para deducir el comportamiento de los sucesos observables y sus parámetros, y cuya distribución de probabilidad no depende de unos parámetros preestablecidos), que pueden invertirse para dar una probabilidad de que la observación futura Xn+1 caiga en algún intervalo calculado en función de los valores observados hasta el momento, . Cuando se trabaja con una cantidad pivote, dependiendo solo de los observables, su estudio se denomina estadística auxiliar. Confidence Regions El método usual para construir cantidades pivote es tomar la diferencia de dos variables que dependan de su ubicación, para que esa ubicación se cancele, y luego tomar la proporción de dos variables que dependan de la escala, haciendo que la relación de escala también se cancele.
La cantidad pivotal más familiar es la t de Student, que puede deducirse mediante este método y se utiliza en el apartado siguiente.
Media conocida, varianza conocida Un intervalo de predicción [ℓ, u] para una observación futura X en una distribución normal N (µ, σ2) con media y varianza conocidas se pueden calcular a partir de
donde , la unidad tipificada de X, se distribuye como estándar normal.
Por lo tanto
o
siendo z el cuantil en la distribución normal estándar para el cual:
o equivalentemente;
Intervalo dePredicción z 75% 1.15Tabla A2 en 90% 1.64 95% 1.96 99% 2.58
El intervalo de predicción se escribe convencionalmente como:
Por ejemplo, para calcular el intervalo de predicción con una probabilidad del 95% para una distribución normal con una media (µ) de 5 y una desviación estándar (σ) de 1, entonces z es aproximadamente 2 (véase la tabla adjunta). Por lo tanto, el límite inferior del intervalo de predicción es aproximadamente 5 - (2·1) = 3, y el límite superior es aproximadamente 5 + (2·1) = 7, dando así un intervalo de predicción de aproximadamente [3 , 7].
Estimación de los parámetros Para una distribución con parámetros desconocidos, un enfoque directo de la predicción es estimar los parámetros y luego usar la función de cuantil asociada; por ejemplo, se puede usar la media de la muestra como estimación para μ y la varianza 's' '2 como una estimación para' 'σ2. Tenga en cuenta que hay dos opciones naturales para s2 aquí: dividir por produce una estimación imparcial, mientras que dividir por n produce el máxima verosimilitud, y cualquiera de ellas puede ser utilizada. Luego se usa la función de cuantiles con estos parámetros estimados para dar un intervalo de predicción.
Este enfoque es utilizable, pero el intervalo resultante no tendrá la interpretación de muestreo repetida, no es un intervalo de confianza predictivo.
Para la secuela, usa la media de la muestra:
y la varianza de la muestra (imparcial):
Media desconocida, varianza conocida Dada Non-Bayesian Predictive Approaches es una distribución normal con una media desconocida μ pero varianza conocida 1, la media de la muestra de las observaciones tiene una distribución mientras que la futura observación tiene una distribución Tomando la diferencia de estos valores se cancelan la μ y se obtiene una distribución normal de varianza por lo tanto
La resolución de da la distribución de predicción a partir de la cual se pueden calcular intervalos como antes. Este es un intervalo de confianza predictivo en el sentido de que si se usa un rango de cuantiles de 100p%, luego, en las aplicaciones repetidas de este cálculo, la futura observación caerá en el intervalo predicho 100p% de las veces.
Téngase en cuenta que esta distribución de predicción es más conservadora que el uso de la media estimada y la varianza conocida 1, ya que utiliza la varianza , por lo que produce intervalos más amplios. Esto es necesario para que la propiedad del intervalo de confianza deseado se mantenga.
Media conocida, varianza desconocida Por el contrario, dada una distribución normal con media conocida 0, pero varianza desconocida , la varianza muestral de las observaciones tiene, prescindiendo de la escala, una distribución ; más precisamente:
mientras que la futura observación tiene distribución Al tomar la proporción de la observación futura y la desviación estándar de la muestra, se cancela el σ, lo que produce una distribución t de Student con n - 1 grados de libertad:
La resolución de da la distribución de predicción a partir de la que se pueden calcular intervalos como antes.
Téngase en cuenta que esta distribución de predicción es más conservadora que el uso de una distribución normal con la desviación estándar estimada y la media conocida 0, ya que utiliza la distribución t en lugar de la distribución normal, por lo que produce intervalos más amplios. Esto es necesario para que la propiedad del intervalo de confianza deseado se mantenga.
Media desconocida, varianza desconocida Combinando lo anterior para una distribución normal con μ y σ2 desconocidas se obtiene la siguiente estadística auxiliar: Confidence Regions
Esta combinación simple es posible porque la media muestral y la varianza muestral de la distribución normal son estadísticas independientes; esto solo es cierto para la distribución normal y, de hecho, caracteriza la distribución normal.
Resolviendo para se obtiene la distribución de predicción:
La probabilidad de que caiga en un intervalo dado es entonces:
donde Ta es el 100(1 - p/2)ésimo percentil de la distribución t de Student, con n & minus; 1 grados de libertad. Por lo tanto, los números:
son los puntos finales de un intervalo de predicción de 100(1 - p) para .
Métodos no paramétricos
Se pueden calcular intervalos de predicción sin ningún supuesto sobre la población. Formalmente, este es un método no paramétrico."Prediction Intervals ", Statistics @ SUNY Oswego
Supóngase que se extrae al azar una muestra de dos observaciones X1 y X2 de una población en la que se supone que los valores tienen una distribución de probabilidad continua.
¿Cuál es la probabilidad de que X2 > X1?
La respuesta es exactamente el 50%, independientemente de la población subyacente: la probabilidad de elegir 3 y luego 7 es la misma que la de seleccionar 7 y luego 3, independientemente de la probabilidad particular de elegir 3 o 7. Por lo tanto, si se elige un solo punto de muestra X1, entonces el 50% de las veces el siguiente punto de la muestra será mayor, lo que arroja (X1, +∞) como un intervalo de predicción del 50% para X2. De manera similar, el 50% de las veces será más pequeño, lo que produce otro intervalo de predicción del 50% para X2, a saber (−∞, X1). Téngase en cuenta que el supuesto de una distribución continua evita la posibilidad de que los valores puedan ser exactamente iguales; lo que complicaría las cosas.
De manera similar, si se tiene una muestra {X1, ..., Xn}, entonces la probabilidad de que la siguiente observación Xn+1 sea la más grande es de 1/(n + 1), ya que todas las observaciones tienen la misma probabilidad de ser el máximo. De la misma manera, la probabilidad de que Xn+1 sea el valor más pequeño es 1/(n + 1). El otro (n − 1)/(n + 1) de las veces, Xn+1 cae entre el máximo y el mínimo de la muestra {X1, ..., Xn}. Por lo tanto, al anotar el máximo y el mínimo de la muestra como M y m, esto produce un (n − 1)/(n + 1) intervalo de predicción de [m, M].
Por ejemplo, si n = 19, entonces [m, M] da un intervalo de predicción del 18/20= 90%, es decir, que el 90% de las veces, la vigésima observación debería caer entre la observación más pequeña y la más grande observada hasta el momento. Del mismo modo, n = 39 proporciona un intervalo de predicción del 95%, y n = 199 proporciona un intervalo de predicción del 99%.
Más generalmente, si X(j) y X(k) son observaciones ordenadas de la muestra con j < k y j + k = n + 1, entonces [X(j), X(k)] es un intervalo de predicción para Xn+1 con una cobertura de probabilidad (significación estadística) igual a .
Este hecho se puede visualizar dibujando los n puntos de la muestra sobre una recta, que divide la línea en n + 1 secciones (n − 1 segmentos entre muestras, y 2 intervalos que van hasta el infinito en ambos extremos), y observando que Xn+1 tiene la misma posibilidad de caer en cualquiera de estas n + 1 secciones. Por lo tanto, también se puede elegir cualquier k de estas secciones y dar un k / (n + 1) intervalo de predicción (o conjunto, si las secciones no son consecutivas). Por ejemplo, si n = 2, entonces la probabilidad de que X3 caiga entre las dos observaciones previas es de 1/3.