Snažím se získat Bonferroniho korekci správně.
Předpokládejme, že mám datovou sadu se 100 hodnotami p, 20 významné minimálně na 10% úrovni (hodnota p < 0,1) a 80 nevýznamné (hodnota p> = 0,1).
Při použití Bonferroniho korekce s R pomocí p.adjust() jaká je správná n hodnota?
Mám uvést n = 100 (všechny hodnoty p) nebo n = 20 (pouze významné hodnoty p)?
Děkujeme za jakoukoli pomoc
Odpovědět
Neznamená to Nezáleží na tom, zda jsou vaše hodnoty p významné nebo ne. Vaše n je počet testů (= hodnot p), které chcete opravit. Pokud vás zajímají výsledky všech 100 testů, pak n = 100.
Komentáře
- takže použití n = 100 by mělo platit také v případě, že mám zájem opravit pouze hodnoty p klobouk je významný?
- Ano, neměli byste opravovat na základě hodnot p. Před spuštěním testů byste měli kontrolovat míru chyb typu 1 u více testů. Nejprve určete, kolik testů chcete skutečně spustit. Za provádění zbytečných testů budete potrestáni.
- Děkuji. Ale s mými skutečnými daty z 1 000 testů dostanu po korekci Bonferroni pouze 2 významné hodnoty p. To je podle mého názoru příliš přísné.
- Existují další opravy, které by vám mohly poskytnout větší sílu k detekci významných účinků. Důležité body jsou: 1) proč provádíte 100 (1 000?) Testů? 2) proč si myslíte, že vaším cílem je najít význam? To by ' tak nemělo být, měli byste se snažit provádět správné závěry a opravy – Bonferroni nebo jiné – vám pomohou kontrolovat míru chyb typu 1, což pomůže v tomto směru.
- Neznám ' tvé otázky týkající se oboru, designu nebo výzkumu, takže je mimo mě, zda všechny takové testy opravdu potřebuješ. Pokud se jedná o průzkumný výzkum, možná zvažte jiné alternativy než hodnoty NHST ap.