¿Por qué controlar FDR es menos estricto que controlar FWER?

He leído que controlar FDR es menos estricto que controlar FWER, como en Wikipedia :

Los procedimientos de control de FDR ejercen un control menos estricto sobre el descubrimiento falso en comparación con los procedimientos de tasa de error familiar (FWER) (como la corrección de Bonferroni). Esto aumenta la potencia a costa de aumentar la tasa de errores de tipo I, es decir, rechazar la hipótesis nula de ningún efecto cuando debería aceptarse.

Pero me preguntaba cómo se demuestra que es matemáticamente cierto.

¿Hay alguna relación entre FDR y FWER?

hypothesis-testing multiple-comparisons false-discovery-rate StackExchange para todos
fuente

¿Leíste el periódico original? Es casi todo lo que uno podría esperar en un documento de estadísticas: una sola idea fundamental, una historia clara y concisa para contar, un ejemplo útil y (¡breve!) Pruebas precisas.

cardenal

De hecho, @cardinal tiene toda la razón en que el documento es tan claro como parece. Entonces, para lo que vale, en caso de que no tenga acceso al documento, aquí hay una versión ligeramente elaborada de cómo argumenta Benjamini-Hochberg:

El FDR es el valor esperado de la proporción de falsos rechazos a todos los rechazos . Ahora, es, obviamente, la suma de los rechazos falsos y correctos; llamar a este último . $Q_e$ $v$ $r$ $r$ $s$

En resumen, (usando letras mayúsculas para variables aleatorias y letras minúsculas para valores realizados),

Q_{e} = E (\frac{V}{R}) = E (\frac{V}{V + S}) =: E (Q) .

$Q_e=E\left(\frac{V}{R}\right)=E\left(\frac{V}{V+S}\right)=:E\left(Q\right).$

Uno toma si . $Q=0$ $R=0$

Ahora, hay dos posibilidades: o bien todos los valores nulos son verdaderas o simplemente de ellos son verdaderas. En el primer caso, no puede haber rechazos correctos, entonces . Por lo tanto, si hay rechazos ( ), , de lo contrario . Por lo tanto, $m$ $m_0<m$ $r=v$ $r\geq 1$ $q=1$ $q=0$

F D R = E (Q) = 1 \cdot P (Q = 1) + 0 \cdot P (Q = 0) = P (Q = 1) = P (V \geq 1) = F W E R

$\newcommand{\FDR}{\mathrm{FDR}}\newcommand{\FWER}{\mathrm{FWER}}\FDR=E(Q)=1\cdot P(Q=1)+0\cdot P(Q=0)=P(Q=1)=P(V \geq 1)=\FWER$

Entonces, en este caso, de modo que cualquier procedimiento que controle el trivialmente también controla el y viceversa. $\FDR=\FWER$ $\FDR$ $\FWER$

En el segundo caso en el que , si (por lo tanto, si hay al menos un falso rechazo), obviamente tenemos (esto es una fracción con también en el denominador) que . Esto implica que la función indicadora que toma el valor 1 si hay al menos un falso rechazo, nunca será menor que , . Ahora, espere a ambos lados de la desigualdad, que por la monotonicidad de deja la desigualdad intacta, $m_0<m$ $v>0$ $v$ $v/r\leq 1$ $\mathbf 1_{V\geq 1}$ $Q$ $\mathbf 1_{V\geq 1}\geq Q$ $E$

E (1_{V \geq 1}) \geq E (Q) = F D R

$E(\mathbf 1_{V\geq 1})\geq E(Q)=\FDR$

El valor esperado de una función de indicador es la probabilidad del evento en el indicador, tenemos , que nuevamente es . $E(\mathbf 1_{V\geq 1})=P(V\geq 1)$ $\FWER$

Por lo tanto, cuando tenemos un procedimiento que controla el en el sentido de que , debemos tener ese . $\FWER$ $\FWER\leq \alpha$ $\FDR\leq\alpha$

Por el contrario, tener un control en algún puede venir con un sustancialmente mayor . Intuitivamente, aceptar una fracción esperada distinta de cero de falsos rechazos ( ) de un total potencialmente grande de hipótesis probadas puede implicar una probabilidad muy alta de al menos un falso rechazo ( ). $\FDR$ $\alpha$ $\FWER$ $\FDR$ $\FWER$

Por lo tanto, un procedimiento tiene que ser menos estricto cuando solo se desea el control , que también es bueno para el poder. Esta es la misma idea que en cualquier prueba de hipótesis básica: cuando prueba al nivel del 5%, rechaza con mayor frecuencia (tanto nulos correctos como falsos) que cuando prueba al nivel del 1% simplemente porque tiene un valor crítico menor. $\FDR$

Christoph Hanck
fuente

(+1) Buena exposición. Obviamente, en el primer caso también podemos decir que el control FWER implica control FDR (que es el asunto en cuestión). Además, puede valer la pena señalar que esta propiedad viene sin supuestos de distribución (por ejemplo, independencia) en las estadísticas de prueba, a diferencia del procedimiento dado en el documento original para el control del FDR.

cardenal

¿Por qué controlar FDR es menos estricto que controlar FWER?

Respuestas: