Ciao comunità JMP,
Da che ho memoria, ho sempre utilizzato un'approssimazione del calcolo FDR di Benjamini-Hochberg per tenere conto dei test multipli:
> Su un elenco ordinato di valori P nominali (dal più piccolo al più grande): BH p Val = Nom P Val * N Row()/Row ()
Tuttavia, sono consapevole di alcuni limiti di questo approccio:
- Per valori p nominali elevati, la formula descritta sopra restituisce valori > 1 (non di molto, ma mi dice che qualcosa non va).
- Per valori p nominali uguali (o molto simili tra loro), la formula restituisce valori leggermente diversi, cosa che a mio avviso non dovrebbe accadere.
Esiste quindi un modo migliore per calcolare l'FDR?
Grazie per l'aiuto.
Migliore,
TS
Salve JMP Community,
Con un po 'di ricerca sono riuscito a trovare un add-in sviluppato da John Sall di JMP che calcola correttamente il valore p di FDR
https://community.jmp.com/t5/JMP-Add-Ins/False-Discovery-Rate-PValue/ta-p/21353
Migliore,
TS
