Calcul exact de Benjamini-Hochberg (FDR)?Salut la communauté JMP,
Aussi loin que je me souvienne, j'ai utilisé une approximation du calcul du FDR Benjamini-Hochberg pour tenir compte de plusieurs tests:
> Sur une liste triée de valeurs nominales P (de plus petite à plus grande): BH p Val = Nom P Val * N Row () / Row ()
Cependant, je suis conscient de certaines limites avec cette approche:
- Pour les grandes valeurs nominales p, la formule décrite ci-dessus renvoie des valeurs> 1 (pas de beaucoup mais cela me dit que quelque chose ne va pas)
- Pour les valeurs p nominales qui sont égales (ou très proches les unes des autres), la formule renvoie des valeurs légèrement différentes qui, à mon avis, ne devraient pas être le cas
Par conséquent, y a-t-il une meilleure façon de calculer le FDR?
Merci de votre aide.
Meilleur,
TS
AnalytiqueAutreRe: Calcul exact de Benjamini-Hochberg (FDR)?Salut la communauté JMP,
Avec un peu de fouille, j'ai pu trouver un complément développé par John Sall de JMP qui calcule correctement la valeur p FDR False Discovery Rate PValue
Meilleur,
TS
Calcul exact de Benjamini-Hochberg (FDR) ?Bonjour la communauté JMP,
Aussi loin que je me souvienne, j'ai utilisé une approximation du calcul FDR de Benjamini-Hochberg pour tenir compte des tests multiples :
> Sur une liste triée des valeurs p nominales (de la plus petite à la plus grande) : BH p Val = Nom P Val * N Row()/Row ()
Cependant, je suis conscient de certaines limites de cette approche :
- Pour des valeurs p nominales élevées, la formule décrite ci-dessus renvoie des valeurs > 1 (de peu, mais cela m'indique que quelque chose ne va pas).
- Pour des valeurs p nominales égales (ou très proches), la formule renvoie des valeurs légèrement différentes, ce qui me semble anormal.
Existe-t-il donc une meilleure façon de calculer le FDR ?
Merci pour votre aide.
Meilleur,
TS
Études de marché et de consommationIngénierie de la qualité et des procédés
