Articulo de referencia

False coverage rate

In statistics , a false coverage rate (FCR) is the average rate of false coverage , i.e. not covering the true parameters, among the selected intervals. The FCR gives a simultan...

In statistics, a false coverage rate (FCR) is the average rate of false coverage, i.e. not covering the true parameters, among the selected intervals.

The FCR gives a simultaneous coverage at a (1  α)×100% level for all of the parameters considered in the problem. The FCR has a strong connection to the false discovery rate (FDR). Both methods address the problem of multiple comparisons, FCR from confidence intervals (CIs) and FDR from P-value's point of view.

FCR was needed because of dangers caused by selective inference. Researchers and scientists tend to report or highlight only the portion of data that is considered significant without clearly indicating the various hypothesis that were considered. It is therefore necessary to understand how the data is falsely covered. There are many FCR procedures which can be used depending on the length of the CI – Bonferroni-selected–Bonferroni-adjusted, Adjusted BH-Selected CIs (Benjamini and Yekutieli 2005[1]). The incentive of choosing one procedure over another is to ensure that the CI is as narrow as possible and to keep the FCR. For microarray experiments and other modern applications, there are a huge number of parameters, often tens of thousands or more and it is very important to choose the most powerful procedure.

The FCR was first introduced by Daniel Yekutieli in his PhD thesis in 2001.[2]

Definitions

Not keeping the FCR means FCR>q{\displaystyle {\text{FCR}}>q} when q=VR=αm0R{\displaystyle q={\frac {V}{R}}={\frac {\alpha m_{0}}{R}}}, where m0{\displaystyle m_{0}} is the number of true null hypotheses, R{\displaystyle R} is the number of rejected hypothesis, V{\displaystyle V} is the number of false positives, and α{\displaystyle \alpha } is the significance level. Intervals with simultaneous coverage probability 1q{\displaystyle 1-q} can control the FCR to be bounded by q{\displaystyle q}.

Classification of multiple hypothesis tests

The following table defines the possible outcomes when testing multiple null hypotheses. Suppose we have a number m of null hypotheses, denoted by: H1, H2, ..., Hm. Using a statistical test, we reject the null hypothesis if the test is declared significant. We do not reject the null hypothesis if the test is non-significant. Summing each type of outcome over all Hi  yields the following random variables:

En m pruebas de hipótesis de las cualesmetro0{\displaystyle m_{0}}son hipótesis nulas verdaderas, R es una variable aleatoria observable y S , T , U y V son variables aleatorias no observables .

Los problemas abordados por FCR

Selección

La selección provoca una cobertura promedio reducida. La selección puede presentarse como un condicionamiento a un evento definido por los datos y puede afectar la probabilidad de cobertura de un IC para un solo parámetro . De manera equivalente, el problema de la selección cambia el sentido básico de los valores p . Los procedimientos FCR consideran que el objetivo de la cobertura condicional siguiendo cualquier regla de selección para cualquier conjunto de valores (desconocidos) para los parámetros es imposible de lograr. Una propiedad más débil en lo que respecta a los IC selectivos es posible y evitará declaraciones de cobertura falsas. FCR es una medida de la cobertura del intervalo después de la selección. Por lo tanto, aunque un IC 1  α no ofrece cobertura selectiva ( condicional ), la probabilidad de construir un IC sin cobertura es como máximo α , donde 

Pr[θdoI, CI construido]Pr[θdoI]α{\displaystyle \Pr[\theta \not \in \mathrm {CI} ,\ {\text{CI constructed}}]\leq \Pr[\theta \not \in \mathrm {CI} ]\leq \alpha }

Selección y multiplicidad

Al enfrentarse a la multiplicidad (inferencia sobre múltiples parámetros) y a la selección , la proporción esperada de cobertura sobre los parámetros seleccionados en 1−α no solo no es equivalente a la proporción esperada de no cobertura en α, sino que esta última ya no puede garantizarse mediante la construcción de intervalos de confianza marginales para cada parámetro seleccionado. Los procedimientos FCR resuelven esto tomando la proporción esperada de parámetros no cubiertos por sus intervalos de confianza entre los parámetros seleccionados, donde la proporción es 0 si no se selecciona ningún parámetro. Esta tasa de falsa cobertura (FCR) es una propiedad de cualquier procedimiento que se define por la forma en que se seleccionan los parámetros y la forma en que se construyen los múltiples intervalos.

Procedimientos de control

Procedimiento de Bonferroni (seleccionado por Bonferroni y ajustado por Bonferroni) para IC simultáneo

Intervalos de confianza simultáneos con el procedimiento de Bonferroni cuando tenemos m parámetros, cada intervalo de confianza marginal construido al nivel 1 − α/m. Sin selección, estos intervalos de confianza ofrecen cobertura simultánea, en el sentido de que la probabilidad de que todos los intervalos de confianza cubran sus respectivos parámetros es al menos 1 − α. Desafortunadamente, incluso una propiedad tan sólida no garantiza la propiedad de confianza condicional después de la selección.

FCR para IC simultáneo seleccionado mediante Bonferroni y ajustado mediante Bonferroni

El procedimiento Bonferroni-Bonferroni no puede ofrecer cobertura condicional, sin embargo, controla el FCR en <α De hecho, lo hace demasiado bien, en el sentido de que el FCR está demasiado cerca de 0 para valores grandes de θ. La selección de intervalos se basa en la prueba de Bonferroni, y luego se construyen IC de Bonferroni. El FCR se estima como, la proporción de intervalos que no cubren sus parámetros respectivos entre los IC construidos se calcula (estableciendo la proporción en 0 cuando no se selecciona ninguno). Donde la selección se basa en pruebas individuales no ajustadas y se construyen IC no ajustados.

IC seleccionados por BH ajustados por FCR

En el procedimiento BH para FDR después de ordenar los valores p P (1) ≤  • • •  P ( m ) y calcular R = max{ j  : P ( j ) ≤ jq / m }, las R hipótesis nulas para las cuales P ( i ) ≤ Rq / m son rechazadas. Si la prueba se realiza utilizando el procedimiento de Bonferroni, entonces el límite inferior del FCR puede caer muy por debajo del nivel deseado q , lo que implica que los intervalos son demasiado largos. Por el contrario, aplicar el siguiente procedimiento, que combina el procedimiento general con la prueba de control de FDR en el procedimiento BH, también produce un límite inferior para el FCR, q /2 ≤ FCR. Este procedimiento es preciso en el sentido de que para algunas configuraciones, el FCR se aproxima a q .

1. Ordene los valores p utilizados para probar las m hipótesis con respecto a los parámetros, P (1) ≤ • • • ≤ P ( m ).

2. Calcular R = max{ i  : P ( i ) ≤ iq / m }.

3. Seleccione los parámetros R para los cuales P ( i ) ≤ Rq / m , correspondientes a las hipótesis rechazadas.

4. Construya un IC de 1 Rq / m para cada parámetro seleccionado.  

Véase también

Referencias

Notas a pie de página

  1. Benjamini, Yoav; Yekutieli, Daniel (marzo de 2005). "Intervalos de confianza múltiples ajustados por la tasa de falsos descubrimientos para parámetros seleccionados" (pdf) . Journal of the American Statistical Association . 100 (469): 71– 93. doi : 10.1198/016214504000001907 .
  2. Resultados teóricos necesarios para aplicar la tasa de falsos descubrimientos en problemas estadísticos . Abril de 2001 (Sección 3.2, página 51)

Otras fuentes

  • Zhao, Zhigen; Hwang, JT Gene (2012). "Intervalos de confianza de control de la tasa de cobertura falsa de Bayes empírico" (pdf) . Journal of the Royal Statistical Society, Serie B. doi : 10.1111 /j.1467-9868.2012.01033.x .