Por que controlar o FDR é menos rigoroso do que controlar o FWER?


16

Eu li que controlar o FDR é menos rigoroso do que controlar o FWER, como na Wikipedia :

Os procedimentos de controle de FDR exercem um controle menos rigoroso sobre a descoberta falsa em comparação com os procedimentos da taxa de erro familiar (FWER) (como a correção de Bonferroni). Isso aumenta o poder ao custo de aumentar a taxa de erros do tipo I, ou seja, rejeitar a hipótese nula de nenhum efeito quando deve ser aceita.

Mas eu queria saber como isso é mostrado matematicamente?

Existe alguma relação entre FDR e FWER?


4
Você leu o artigo original? É tudo o que se poderia esperar de um artigo estatístico: uma única idéia fundamental, uma história clara e concisa para contar, um exemplo útil e (breves!) Provas precisas.
cardeal

Respostas:


10

De fato, @cardinal está certo de que o artigo é o mais claro possível. Então, pelo que vale a pena, caso você não tenha acesso ao artigo, aqui está uma versão um pouco elaborada de como Benjamini – Hochberg argumenta:

O FDR é o valor esperado da proporção de falsas rejeições v para todas as rejeições r . Agora, r é, obviamente, a soma de rejeições falsas e corretas; chame o último s .Qevrrs

Em resumo, (usando letras maiúsculas para variáveis ​​aleatórias e letras minúsculas para valores realizados),

Qe=E(VR)=E(VV+S)=:E(Q).

Toma-se se R = 0 .Q=0R=0

Agora, existem duas possibilidades: ou todos os nulos são verdadeiros ou apenas m 0 < m deles são verdadeiros. No primeiro caso, não pode haver rejeições corretas, portanto r = v . Assim, se houver alguma rejeição ( r 1 ), q = 1 , caso contrário, q = 0 . Conseqüentemente,mm0<mr=vr1q=1q=0

FDR=E(Q)=1P(Q=1)+0P(Q=0)=P(Q=1)=P(V1)=FWER

Portanto, neste caso, de modo que qualquer procedimento que controla o F D R trivialmente também controla o F W E R e vice-versa.FDR=FWERFDRFWER

No segundo caso em que , se v > 0 (portanto, se houver pelo menos uma rejeição falsa), obviamente temos (sendo esta uma fração com também v no denominador) que v / r 1 . Isto implica que a função de indicador que assume o valor 1 se houver pelo menos uma falsa rejeição, 1 V 1 nunca será inferior a Q , 1 V 1Q . Agora, tenha expectativa de ambos os lados da desigualdade, que pela monotonicidade de Em0<mv>0vv/r11V1Q1V1QE deixa a desigualdade intacta,

E(1V1)E(Q)=FDR

O valor esperado de uma função indicadora é a probabilidade do evento no indicador, temos , que mais uma vez é o F W E R .E(1V1)=P(V1)FWER

Assim, quando temos um procedimento que controla no sentido de F W E Rα , devemos ter esse F D Rα .FWERFWERαFDRα

Por outro lado, tendo controlo em algum α pode vir com um substancialmente maior F W E R . Intuitivamente, aceitar uma fração esperada diferente de zero de rejeições falsas ( F D R ) de um total potencialmente grande de hipóteses testadas pode implicar uma probabilidade muito alta de pelo menos uma rejeição falsa ( F W E R ).FDRαFWERFDRFWER

Portanto, um procedimento deve ser menos rigoroso quando apenas o controle é desejado, o que também é bom para energia. É a mesma idéia que em qualquer teste de hipótese básica: quando você testa no nível de 5%, rejeita com mais frequência (nulos corretos e falsos) do que quando testa no nível de 1% simplesmente porque possui um valor crítico menor.FDR


2
(+1) Boa exposição. Obviamente, no primeiro caso, também podemos dizer que o controle FWER implica o controle FDR (que é o problema em questão). Além disso, pode ser interessante ressaltar que essa propriedade não possui premissas distributivas (por exemplo, independência) nas estatísticas de teste, diferentemente do procedimento fornecido no artigo original para controle do FDR.
cardeal
Ao utilizar nosso site, você reconhece que leu e compreendeu nossa Política de Cookies e nossa Política de Privacidade.
Licensed under cc by-sa 3.0 with attribution required.