R – Bonferronin korjaus p.adjust () -toiminnolla, mikä ' on oikea n-arvo?

Yritän saada Bonferroni-korjauksen oikein.

Oletetaan, että minulla on tietojoukko, jossa on 100 p-arvoa, 20 merkitsevä vähintään 10%: n tasolla (p-arvo < 0,1) ja 80 ei merkitsevä (p-arvo> = 0,1).

Bonferroni-korjausta sovellettaessa R: n kanssa käyttämällä p.adjust() mitä oikea n -arvo on?

Laitetaanko n = 100 (kaikki p-arvot) vai n = 20 (vain merkittävät p-arvot)?

Kiitos avusta

Vastaa

Ei ”t riippumatta siitä, ovatko p-arvosi merkitseviä vai eivät. N on testien lukumäärä (= p-arvot), jotka haluat korjata. Jos olet kiinnostunut kaikkien 100 testin tuloksista, n = 100.

Kommentit

  • joten n = 100: n käyttöä tulisi soveltaa myös, jos olen kiinnostunut korjaamaan vain p-arvoja t hatut ovat merkittäviä?
  • Kyllä, sinun ei pitäisi korjata p-arvojen perusteella. Sinun tulisi hallita tyypin 1 virhesuhdetta useissa testeissä ennen testien suorittamista. Määritä ensin, kuinka monta testiä haluat todella suorittaa. Saat rangaistuksen tarpeettomien testien suorittamisesta.
  • kiitos. Mutta todellisilla tiedoillani 1000 testistä saan vasta 2 merkittävää p-arvoa Bonferroni-korjauksen jälkeen. Se on mielestäni liian tiukka.
  • On olemassa muita korjauksia, jotka voivat antaa sinulle enemmän voimaa havaita merkittäviä vaikutuksia. Tärkeitä seikkoja ovat: 1) miksi suoritat 100 (1000?) Testiä? 2) miksi kuulostaa siltä, että tavoitteesi on löytää merkitys? Tämän ei pitäisi olla ', sinun pitäisi pyrkiä tekemään oikeat johtopäätökset ja korjaukset – Bonferroni tai muu – auttaa hallitsemaan tyypin 1 virhesuhdetta, mikä auttaa tähän suuntaan.
  • En tiedä ' tietämystäsi alastasi, suunnittelustasi tai tutkimuskysymyksistäsi, joten tarvitsen todella kaikkia tällaisia testejä. Jos tämä on tutkimustutkimusta, harkitse ehkä muita vaihtoehtoja kuin NHST- ja p-arvot.

Vastaa

Sähköpostiosoitettasi ei julkaista. Pakolliset kentät on merkitty *