Articulo de referencia

Probit

Gráfico de la función probit En estadística , la función probit convierte una probabilidad (un número entre 0 y 1) en una puntuación. Esta puntuación indica cuántas desviaciones...

Gráfico de la función probit

En estadística , la función probit convierte una probabilidad (un número entre 0 y 1) en una puntuación. Esta puntuación indica cuántas desviaciones estándar se encuentra un valor con respecto a la media en una distribución normal estándar (o "curva de campana"). Por ejemplo, una probabilidad de 0,5 (50 %) representa el punto medio exacto de la distribución, por lo que su puntuación probit es 0. Una probabilidad menor, como 0,025 (2,5 %), se encuentra muy a la izquierda de la curva, lo que corresponde a una puntuación probit de aproximadamente -1,96.

La función se utiliza ampliamente en modelos probit , un tipo de análisis de regresión para resultados binarios (por ejemplo, éxito/fracaso o aprobado/suspenso). Se desarrolló inicialmente en toxicología para analizar relaciones dosis-respuesta, como la variación del porcentaje de plagas eliminadas por un plaguicida en función de su concentración. [ 1 ] La función probit también se utiliza para crear gráficos Q-Q , una herramienta gráfica para evaluar si un conjunto de datos sigue una distribución normal.

Matemáticamente, la función probit es la función cuantil (la inversa de la función de distribución acumulativa (FDA)) asociada con la distribución normal estándar. Si la FDA se denota por Φ(z){\displaystyle \Phi (z)} , entonces la función probit se define como: probit(pag)=Φ1(pag)parapag(0,1).{\displaystyle \operatorname {probit} (p)=\Phi ^{-1}(p)\quad {\text{para}}\quad p\in (0,1).} Esto significa que para cualquier probabilidadpag{\displaystyle p} , la función probit encuentra el valorz{\displaystyle z}de tal manera que el área bajo la curva normal estándar a la izquierda dez{\displaystyle z}es igual apag{\displaystyle p}.

Desarrollo conceptual

La idea de la función probit fue publicada por Chester Ittner Bliss en un artículo de 1934 en Science sobre cómo tratar datos como el porcentaje de una plaga muerta por un pesticida . [ 2 ] Bliss propuso transformar el porcentaje muerto en una "unidad de probabilidad" (o "probit") que estaba linealmente relacionada con la definición moderna (la definió arbitrariamente como igual a 0 para 0,0001 y 1 para 0,9999): [ 3 ] Incluyó una tabla para ayudar a otros investigadores a convertir sus porcentajes de muerte a su probit, que luego podrían graficar contra el logaritmo de la dosis y así, se esperaba, obtener una línea más o menos recta. Este llamado modelo probit sigue siendo importante en toxicología, así como en otros campos. El enfoque se justifica en particular si la variación de la respuesta puede racionalizarse como una distribución lognormal de tolerancias entre sujetos en prueba, donde la tolerancia de un sujeto particular es la dosis justo suficiente para la respuesta de interés.

El método introducido por Bliss fue continuado en Probit Analysis , un texto importante sobre aplicaciones toxicológicas de DJ Finney . [ 4 ] Los valores tabulados por Finney pueden derivarse de los probits definidos aquí sumando un valor de 5. Esta distinción es resumida por Collett . [ 5 ] La metodología probit, incluyendo la optimización numérica para el ajuste de funciones probit, se introdujo antes de la disponibilidad generalizada de la computación electrónica. Al usar tablas, era conveniente tener probits uniformemente positivos. Las áreas de aplicación comunes no requieren probits positivos.

Simetrías

Debido principalmente al teorema del límite central , la distribución normal estándar desempeña un papel fundamental en la teoría de la probabilidad y la estadística. Si consideramos el hecho conocido de que la distribución normal estándar sitúa el 95% de la probabilidad entre -1,96 y 1,96 y es simétrica respecto a cero, se deduce que Φ(1,96)=0,025=1Φ(1,96).{\displaystyle \Phi (-1,96)=0,025=1-\Phi (1,96).}

La función probit proporciona el cálculo "inverso", generando un valor de una variable aleatoria normal estándar, asociada con una probabilidad acumulada especificada. Continuando con el ejemplo, probit(0,025)=1,96=probit(0,975).{\displaystyle \operatorname {probit} (0.025)=-1.96=-\operatorname {probit} (0.975).}

En general, Φ(probit(pag))=pag,{\displaystyle \Phi (\operatorname {probit} (p))=p,} y probit(Φ(z))=z.{\displaystyle \operatorname {probit} (\Phi (z))=z.}

Diagnóstico de la desviación de una distribución respecto a la normalidad.

Además de servir de base para importantes tipos de regresión, la función probit es útil en el análisis estadístico para diagnosticar desviaciones de la normalidad, según el método de representación gráfica Q-Q. Si un conjunto de datos es en realidad una muestra de una distribución normal , la representación gráfica de los valores frente a sus puntuaciones probit será aproximadamente lineal. Se pueden diagnosticar desviaciones específicas de la normalidad, como asimetría , colas pesadas o bimodalidad , a partir de la detección de desviaciones específicas de la linealidad. Si bien la representación gráfica Q-Q puede utilizarse para comparar con cualquier familia de distribuciones (no solo la normal), la representación gráfica Q-Q normal es un procedimiento de análisis exploratorio de datos relativamente estándar, ya que la suposición de normalidad suele ser un punto de partida para el análisis.

Cálculo

La función de distribución acumulada (FDA) de la distribución normal y su inversa no están disponibles en forma cerrada , y su cálculo requiere el uso cuidadoso de procedimientos numéricos. Sin embargo, estas funciones están ampliamente disponibles en software para modelado estadístico y de probabilidad, y en hojas de cálculo. En entornos informáticos donde se dispone de implementaciones numéricas de la función de error inversa , la función probit se puede obtener como probit(pag)=2terreno1(2pag1).{\displaystyle \operatorname {probit} (p)={\sqrt {2}}\,\operatorname {erf} ^{-1}(2p-1).} Un ejemplo es MATLAB , donde erfinvestá disponible una función. El lenguaje Mathematica implementa InverseErf. Otros entornos implementan directamente la función probit, como se muestra en el siguiente código R.

> qnorm ( 0.025 ) [1] -1.959964 > pnorm ( -1.96 ) [1] 0.02499790

Wichura proporciona un algoritmo rápido para calcular la función probit con 16 decimales; este se utiliza en R para generar variables aleatorias para la distribución normal. [ 6 ]

Una ecuación diferencial ordinaria para la función probit

Otro método de cálculo se basa en la formación de una ecuación diferencial ordinaria (EDO) no lineal para probit, según el método de Steinbrecher y Shaw. [ 7 ] Abreviando la función probit comow(pag){\displaystyle w(p)} , la EDO es dwdpag=1F(w){\displaystyle {\frac {dw}{dp}}={\frac {1}{f(w)}}} dóndeF(w){\displaystyle f(w)}es la función de densidad de probabilidad de w .

En el caso de la gaussiana: dwdpag=2π miw22.{\displaystyle {\frac {dw}{dp}}={\sqrt {2\pi }}\ e^{\frac {w^{2}}{2}}.} Diferenciando de nuevo: d2wdpag2=w(dwdpag)2,{\displaystyle {\frac {d^{2}w}{dp^{2}}}=w\left({\frac {dw}{dp}}\right)^{2},} con las condiciones centrales (iniciales) w(1/2)=0,w(1/2)=2π.{\displaystyle {\begin{aligned}w\left(1/2\right)&=0,\\w'\left(1/2\right)&={\sqrt {2\pi }}.\end{aligned}}}

Esta ecuación puede resolverse mediante varios métodos, incluido el enfoque clásico de series de potencias. A partir de este, se pueden desarrollar soluciones de precisión arbitrariamente alta basadas en el enfoque de Steinbrecher para la serie de la función de error inversa. La solución de la serie de potencias viene dada por w(pag)=π2k=0dk(2k+1)(2pag1)(2k+1),{\displaystyle w(p)={\sqrt {\frac {\pi }{2}}}\sum _{k=0}^{\infty }{\frac {d_{k}}{(2k+1)}}(2p-1)^{(2k+1)},} donde los coeficientesdk{\displaystyle d_{k}}Satisfacer la recurrencia no lineal dk+1=π4j=0kdjdkj(j+1)(2j+1){\displaystyle d_{k+1}={\frac {\pi }{4}}\sum _{j=0}^{k}{\frac {d_{j}d_{kj}}{(j+1)(2j+1)}}} cond0=1{\displaystyle d_{0}=1} . En esta forma la razóndk+1/dk{\displaystyle d_{k+1}/d_{k}}enfoques1{\displaystyle 1}comok{\displaystyle k}se acerca al infinito.

Logit

Comparación de la función logit con un probit escalado (es decir, la función de distribución acumulada inversa de la distribución normal ), comparando logit( x ) vs. Φ −1 ( x ) / π/8 , lo que hace que las pendientes sean iguales en el origen.

Estrechamente relacionadas con la función probit (y el modelo probit ) están la función logit y el modelo logit . La inversa de la función logística viene dada por logit(pag)=ln(pag1pag).{\displaystyle \operatorname {logit} (p)=\ln \left({\frac {p}{1-p}}\right).}

De forma análoga al modelo probit, podemos suponer que dicha cantidad está relacionada linealmente con un conjunto de predictores, lo que da lugar al modelo logit , base en particular del modelo de regresión logística , la forma más común de análisis de regresión para datos de respuesta categórica. En la práctica estadística actual, los modelos de regresión probit y logit se suelen tratar como casos del modelo lineal generalizado .

Véase también

Referencias

  1. Bliss, CI (1934). "El método de los probits" . Science . 79 ( 2037): 38–39 . Bibcode : 1934Sci....79...38B . doi : 10.1126/science.79.2037.38 . JSTOR 1659792. PMID 17813446 .  
  2. Bliss (1934) , págs. 38–39.
  3. Bliss (1934) , p. 39 : "Estas unidades de probabilidad arbitrarias se han denominado 'probits'...".  
  4. Finney, DJ (1971) [1.ª ed. 1947]. Análisis probit (3.ª ed.). Cambridge University Press. ISBN  0-521-08041-XOCLC 174198382 
  5. Collett, D. (1991). Modelling Binary Data . Chapman and Hall / CRC Press. p. 55. La definición original de un probit [con 5 añadido] era principalmente para evitar tener que trabajar con probits negativos; ... Esta definición todavía se usa en algunos ámbitos, pero en los principales paquetes de software estadístico para lo que se conoce como análisis probit , los probits se definen sin la adición de 5.  
  6. Wichura, MJ (1988). "Algoritmo AS241: Los puntos porcentuales de la distribución normal" (PDF) . Applied Statistics . 37 (3). Blackwell Publishing: 477– 484. doi : 10.2307/2347330 . JSTOR 2347330 . 
  7. Steinbrecher, G.; Shaw, WT (2008). "Mecánica de cuantiles" (PDF) . European Journal of Applied Mathematics . 19 (2): 87– 112. doi : 10.1017/S0956792508007341 . S2CID 6899308 . 
  • ¿Qué función de enlace: Logit, Probit o Cloglog? 12.04.2023