R – Bonferroni korrekció a p.adjust () segítségével, mi ' s a helyes n érték?

A Bonferroni korrekciót próbálom helyrehozni.

Tegyük fel, hogy van egy 100 p-értékű, 20-as adatkészletem legalább 10% -os (p-érték < 0,1) és 80 nem szignifikáns (p-érték> = 0,1).

A Bonferroni-korrekció alkalmazásakor R-vel p.adjust() használatával mi a helyes n érték?

Tegyem n = 100 (minden p-érték) vagy n = 20 (csak jelentős p-értékek)?

Köszönöm a segítséget

Válasz

Nem “t” mindegy, hogy a p értékei szignifikánsak-e vagy sem. Az n az a tesztek száma (= p érték), amelyet ki akar javítani. Ha érdekli mind a 100 teszt eredménye, akkor n = 100.

Megjegyzések

  • ezért az n = 100 használatát akkor is alkalmazni kell, ha csak a p a kalap jelentős?
  • Igen, a p értékek alapján nem szabad korrigálni. A tesztek futtatása előtt több tesztnél ellenőriznie kell az 1. típusú hibaarányt. Először határozza meg, hogy hány tesztet szeretne futtatni. Szükségtelen tesztek futtatásáért büntetést kapsz.
  • köszönöm. De az 1000 tesztből származó valós adataimmal csak 2 jelentős p-értéket kapok a Bonferroni-korrekció után. Ez szerintem túl szigorú.
  • Vannak más korrekciók, amelyek nagyobb hatalmat nyújthatnak a jelentős hatások észleléséhez. Fontos szempontok: 1) miért futtat 100 (1000?) Tesztet? 2) Miért hangzik úgy, hogy célja a jelentőség megtalálása? Nem szabad, hogy ' legyen a helyzet, arra kell törekednie, hogy helyes következtetéseket és javításokat tegyen – a Bonferroni vagy más – segít az 1. típusú hibaarány szabályozásában, ami ebben az irányban segít.
  • Nem tudom, hogy ' nem ismeri az Ön szakterületét, tervezési vagy kutatási kérdéseit, így rajtam kívül áll, hogy valóban szüksége van-e ilyen tesztekre. Ha ez feltáró kutatás, akkor fontoljon meg más alternatívákat, mint az NHST és a p értékeket.

Vélemény, hozzászólás?

Az email címet nem tesszük közzé. A kötelező mezőket * karakterrel jelöltük