Existe-t-il un moyen simple d'expliquer pourquoi la procédure de Benjamini et Hochberg (1995) contrôle réellement le taux de fausses découvertes (FDR)? Cette procédure est si élégante et compacte et pourtant la preuve de pourquoi elle fonctionne sous indépendance (apparaissant en annexe de leur article de 1995 ) n'est pas très accessible.
intuition
false-discovery-rate
teaching
Trisoloriansunscreen
la source
la source
Réponses:
Voici un
R
code pour générer une image. Il affichera 15 valeurs de p simulées tracées par rapport à leur ordre. Ils forment donc un motif de points ascendant. Les points sous les lignes rouges / violettes représentent des tests significatifs au niveau 0,1 ou 0,2. Le FDR est le nombre de points noirs sous la ligne divisé par le nombre total de points sous la ligne.J'espère que cela pourrait donner une idée de la forme de la distribution des valeurs de p ordonnées. Le fait que les lignes soient correctes et non par exemple une courbe en forme de parabole a à voir avec la forme des distributions d'ordre. Cela doit être calculé explicitement. En fait, la ligne n'est qu'une solution conservatrice.
la source
set.seed(<some number>)
et publier le chiffre obtenu pour les personnes qui ne lisent pas R?