Articulo de referencia

Corrección de Bonferroni

La corrección de Bonferroni es un método para contrarrestar el problema de las comparaciones múltiples en estadística . Recibe su nombre del matemático Carlo Emilio Bonferroni ....

La corrección de Bonferroni es un método para contrarrestar el problema de las comparaciones múltiples en estadística . Recibe su nombre del matemático Carlo Emilio Bonferroni .

Fondo

La prueba de hipótesis estadística se basa en rechazar la hipótesis nula cuando la probabilidad de los datos observados sería baja si la hipótesis nula fuera verdadera. Cuando se prueban múltiples hipótesis, la probabilidad de observar un evento raro aumenta, por lo que la probabilidad de rechazar incorrectamente una hipótesis nula (es decir, cometer un error de tipo I ) aumenta si varias hipótesis nulas son verdaderas. [ 1 ]

La corrección de Bonferroni compensa ese aumento al probar cada hipótesis individual a un nivel de significancia deα/metro{\displaystyle \alpha /m}, dóndeα{\displaystyle \alpha }es el nivel alfa general deseado ymetro{\displaystyle m}es el número de hipótesis. [ 2 ] Por ejemplo, si se está probando un ensayometro=20{\displaystyle m=20}hipótesis con un objetivo general deseadoα=0,05{\displaystyle \alpha =0.05}, el nivel alfa corregido por Bonferroni seríaα=0,05/20=0,0025{\displaystyle \alpha =0,05/20=0,0025}El método recibe su nombre por el uso de las desigualdades de Bonferroni . [ 3 ]

La corrección de Bonferroni también puede aplicarse como ajuste del valor p: con este método, en lugar de ajustar el nivel alfa, cada valor p se multiplica por el número de pruebas (los valores p ajustados que superan 1 se reducen a 1), y el nivel alfa permanece sin cambios. Las decisiones de significación con este método serán las mismas que con el ajuste del nivel alfa.

Definición

DejarH1,,Hmetro{\displaystyle H_{1},\ldots ,H_{m}}sea ​​una familia de hipótesis nulas y dejemos quepag1,,pagmetro{\displaystyle p_{1},\ldots ,p_{m}}sean sus correspondientes valores p .metro{\displaystyle m}Sea el número total de hipótesis nulas, y seametro0{\displaystyle m_{0}}sea ​​el número de hipótesis nulas verdaderas (que presumiblemente es desconocido para el investigador). La tasa de error familiar (FWER) es la probabilidad de rechazar al menos una hipótesis verdadera.Hi{\displaystyle H_{i}}, es decir, de cometer al menos un error de tipo I. La corrección de Bonferroni rechaza la hipótesis nula para cadapagiαmetro{\displaystyle p_{i}\leq {\frac {\alpha }{m}}}, controlando así el FWER enα{\displaystyle \leq \alpha }La prueba de este control se deduce de la desigualdad de Boole , como sigue:

FWER=PAG{i=1metro0(pagiαmetro)}i=1metro0{PAG(pagiαmetro)}metro0αmetroα.{\displaystyle {\text{FWER}}=P\left\{\bigcup _{i=1}^{m_{0}}\left(p_{i}\leq {\frac {\alpha }{m}}\right)\right\}\leq \sum _{i=1}^{m_{0}}\left\{P\left(p_{i}\leq {\frac {\alpha }{m}}\right)\right\}\leq m_{0}{\frac {\alpha }{m}}\leq \alpha .}

Este control no requiere ninguna suposición sobre la dependencia entre los valores p ni sobre cuántas de las hipótesis nulas son verdaderas. [ 4 ]

Extensiones

Generalización

En lugar de probar cada hipótesis en elα/metro{\displaystyle \alpha /m}nivel, las hipótesis pueden ser probadas en cualquier otra combinación de niveles que sumenα{\displaystyle \alpha }, siempre que el nivel de cada prueba se decida antes de examinar los datos. [ 5 ] Por ejemplo, para dos pruebas de hipótesis, un generalα{\displaystyle \alpha }Se podría mantener un valor de 0,05 realizando una prueba a 0,04 y la otra a 0,01.

Intervalos de confianza

El procedimiento propuesto por Dunn [ 6 ] puede utilizarse para ajustar los intervalos de confianza . Si se establecemetro{\displaystyle m}intervalos de confianza, y desea tener un nivel de confianza general de1α{\displaystyle 1-\alpha }, cada intervalo de confianza individual puede ajustarse al nivel de1αmetro{\displaystyle 1-{\frac {\alpha }{m}}}. [ 6 ]

Problemas continuos

Al buscar una señal en un espacio de parámetros continuo, también puede presentarse el problema de las comparaciones múltiples o el efecto de búsqueda en otro lugar. Por ejemplo, un físico podría estar buscando descubrir una partícula de masa desconocida considerando un amplio rango de masas; este fue el caso durante la detección del bosón de Higgs , que le valió el Premio Nobel . En tales casos, se puede aplicar una generalización continua de la corrección de Bonferroni empleando la lógica bayesiana para relacionar el número efectivo de ensayos.metro{\displaystyle m}, a la relación de volumen anterior a posterior. [ 7 ]

Alternativas

Existen métodos alternativos para controlar la tasa de error familiar . Por ejemplo, el método de Holm-Bonferroni y la corrección de Šidák son procedimientos universalmente más potentes que la corrección de Bonferroni, lo que significa que siempre son al menos igual de potentes. Sin embargo, a diferencia del procedimiento de Bonferroni, estos métodos no controlan el número esperado de errores de tipo I por familia (la tasa de error de tipo I por familia). [ 8 ]

Crítica

Con respecto al control de la tasa de error familiar (FWER) , la corrección de Bonferroni puede ser conservadora si hay un gran número de pruebas y/o los estadísticos de prueba están correlacionados positivamente, debido a la flexibilidad del límite de unión para conjuntos superpuestos. [ 9 ]

Las correcciones para pruebas múltiples, incluido el procedimiento de Bonferroni, aumentan la probabilidad de errores de tipo II cuando las hipótesis nulas son falsas, es decir, reducen la potencia estadística . [ 10 ] [ 9 ]

Referencias

  1. Mittelhammer, Ron C.; Judge , George G .; Miller, Douglas J. (2000). Fundamentos econométricos . Cambridge University Press. págs. 73–74 . ISBN  978-0-521-62394-0.
  2. Miller, Rupert G. (1966). Inferencia estadística simultánea . Springer. ISBN 9781461381228.
  3. ^ Bonferroni, CE, Teoria statistica delle classi e calcolo delle probabilità, Pubblicazioni del R Istituto Superiore di Scienze Economiche e Commerciali di Firenze 1936
  4. Goeman, Jelle J.; Solari, Aldo (2014). "Pruebas de hipótesis múltiples en genómica". Statistics in Medicine . 33 (11): 1946– 1978. doi : 10.1002/sim.6082 . PMID 24399688 . S2CID 22086583 .  
  5. Neuwald, AF; Green, P (1994). "Detección de patrones en secuencias de proteínas". J. Mol. Biol . 239 (5): 698– 712. doi : 10.1006/jmbi.1994.1407 . PMID 8014990 . 
  6. 1 2 Dunn, Olive Jean (1961). "Comparaciones múltiples entre medias" (PDF) . Journal of the American Statistical Association . 56 (293): 52– 64. Bibcode : 1961JASA...56...52D . CiteSeerX 10.1.1.309.1277 . doi : 10.1080/01621459.1961.10482090 . 
  7. Bayer, Adrian E.; Seljak, Uroš (2020). "El efecto de mirar hacia otro lado desde una perspectiva bayesiana y frecuentista unificada" . Journal of Cosmology and Astroparticle Physics . 2020 (10): 009. arXiv : 2007.13821 . Bibcode : 2020JCAP...10..009B . doi : 10.1088/1475-7516/2020/10/009 . S2CID 220830693 . 
  8. Frane, Andrew (2015). "¿Son relevantes las tasas de error de tipo I por familia en las ciencias sociales y del comportamiento?" . Journal of Modern Applied Statistical Methods . 14 (1) jmasm.eP1855: 12– 23. doi : 10.22237/jmasm/1430453040 .
  9. 1 2 Moran, Matthew (2003). "Argumentos para rechazar el método Bonferroni secuencial en estudios ecológicos". Oikos . 100 (2): 403– 405. Bibcode : 2003Oikos.100..403M . doi : 10.1034/j.1600-0706.2003.12010.x .
  10. Nakagawa, Shinichi (2004). "Adiós a Bonferroni: los problemas de la baja potencia estadística y el sesgo de publicación" . Behavioral Ecology . 15 (6): 1044– 1045. doi : 10.1093/beheco/arh107 .
  • Bonferroni, calculadora en línea Sidak