Correzione di Bonferroni con correlazione di Pearson e regressione lineare


9

Sto eseguendo statistiche su 5 IV (5 tratti di personalità, estroversione, gradevolezza, coscienza, nevroticismo, apertura) contro 3 DV Attitude to PCT, Attitude to CBT, Attitude to PCT vs CBT. Ho anche aggiunto età e sesso per vedere quali altri effetti ci sono.

Sto testando per vedere se i tratti della personalità possono prevedere gli atteggiamenti dei DV.

Inizialmente ho usato la correlazione di Pearson per tutte le variabili (45 test).

La scoperta principale è stata che l'estroversione era correlata all'atteggiamento della PCT a p = 0,05. Ma mentre stavo eseguendo 45 test ho fatto una correzione Bonferroni di alfa = 0,05 / 45 = 0,001, rendendo quindi questa scoperta insignificante.

Ho quindi eseguito una semplice regressione lineare su tutte le variabili, di nuovo l'estroversione è stata significativa con l'atteggiamento nei confronti del PCT. Se eseguo la correzione Bonferroni, questa risulta di nuovo insignificante.

Domande:

  1. Devo correggere Bonferroni alla correlazione di Pearson?
  2. Se lo faccio, e quindi rendere insignificante l'estroversione con un atteggiamento nei confronti della PCT, ha ancora senso fare una regressione lineare?
  3. Se eseguo una regressione lineare, devo eseguire anche la correzione Bonferroni per questo?
  4. Devo segnalare solo valori corretti o valori non corretti e corretti?

4
La domanda 1 è correlata a Look e dovresti trovare (una correlazione) e la domanda 3 a. La regolazione dei valori p in una regressione multipla per confronti multipli è una buona idea? . Più in generale i risultati di questa query potrebbero essere di interesse.
chl

Risposte:


4

Penso che Chl ti abbia indicato molti buoni materiali e riferimenti senza rispondere direttamente alla domanda. La risposta che do potrebbe essere un po 'controversa perché so che alcuni statistici non credono nell'adeguamento della molteplicità e molti bayesiani non credono nel valore p. In effetti, una volta ho sentito Don Berry dire che l'uso dell'approccio bayesiano in particolare nei progetti adattivi che controllano l'errore di tipo I non è un problema. Lo ha ripreso più tardi dopo aver visto quanto sia praticamente importante per la FDA assicurarsi che le droghe nocive non arrivino sul mercato.

La mia risposta è sì e no. Se esegui 45 test, devi sicuramente adattarti alla molteplicità, ma no a Bonferroni perché potrebbe essere troppo conservativo. L'inflazione dell'errore di tipo I quando si estraggono i dati per la correlazione è chiaramente un problema che ha attirato l'attenzione sul post citato "guarda e troverai la correlazione". Tutti e tre i collegamenti forniscono ottime informazioni. Quello che penso manchi è l'approccio di ricampionamento all'adeguamento del valore di p così sviluppato da Westfall e Young. Puoi trovare esempi nel mio libro bootstrap o dettagli completi nel loro libro di ricampionamento. La mia raccomandazione sarebbe quella di prendere in considerazione i metodi bootstrap o di permutazione per la regolazione del valore p e forse considerare il tasso di rilevamento falso rispetto al rigoroso tasso di errore familiare.

Link a Westfall and Young: http://www.amazon.com/Resampling-Based-Multiple-Testing-Adjustment-Probability/dp/0471557617/ref=sr_1_1?s=books&ie=UTF8&qid=1343398751&sr=1-1&keywords=peter+ Westfall

Libro recente di Bretz et al su molteplici confronti: http://www.amazon.com/Multiple-Comparisons-Using-Frank-Bretz/dp/1584885742/ref=sr_1_2?s=books&ie=UTF8&qid=1343398796&sr=1-2&keywords= peter + Westfall

Il mio libro con materiale nella sezione 8.5 e tonnellate di riferimenti bootstrap: http://www.amazon.com/Bootstrap-Methods-Practitioners-Researchers-Probability/dp/0471756210/ref=sr_1_2?s=books&ie=UTF8&qid=1343398953&sr=1 -2 & parole chiave = michael + Chernick


+1 La riproduzione della Griglia statistica di Munchausen di Graham Martin alla fine di Westfall & Young dice tutto in un modo molto avvincente. Puoi leggere questo nella funzione "look inside" di Amazon. (È quasi altrettanto divertente vedere Amazon offrire un prezzo di permuta di $ 7 per questo libro da $ 150.)
whuber

@whuber Penso di aver visto una volta un cartone animato che mostrava il Barone tirarsi fuori da un lago dai suoi stivali. Efron potrebbe essere stato saggio chiamarlo bootstrap poiché molti sono scettici sul fatto che si possa fare in statistica proprio come molti sono scettici sulla leggenda del Barone!
Michael R. Chernick,

4

Mi sembra che questa sia ricerca esplorativa / analisi dei dati , non di conferma. Cioè, non sembra che tu abbia iniziato con una teoria secondo cui solo l' estroversione dovrebbe essere correlata al PCT per qualche motivo. Quindi non mi preoccuperei troppo delle regolazioni alfa, poiché la ritengo più correlata al CDA, né ritengo che la tua scoperta sia necessariamente vera. Invece, ci penserei come qualcosa che potrebbesii vero e gioca con queste idee / possibilità alla luce di ciò che so sugli argomenti in questione. Avendo visto questo risultato, suona vero o sei scettico? Cosa significherebbe per le attuali teorie se fosse vero? Sarebbe interessante? Sarebbe importante? Vale la pena eseguire un nuovo studio (di conferma) per determinare se è vero, tenendo conto del tempo, degli sforzi e delle spese potenziali che ciò comporta? Ricorda che il motivo delle correzioni di Bonferroni è che ci aspettiamo che appaia qualcosa quando ci sono così tante variabili. Quindi penso che un euristico possa essere 'se questo studio fosse sufficientemente informativo, anche se la verità risulta essere no'? Se decidi che non ne vale la pena, questa relazione rimane nella categoria "potere" e vai avanti, ma se vale la pena farlo, provalo.


Se capisce veramente cos'è l'analisi dei dati esplorativi e non prende troppo sul serio le grandi correlazioni, sarei d'accordo con te. Ma le persone ammetteranno che stanno solo facendo analisi esplorative per filtrare le correlazioni bizzarre ma tuttavia si eccitano eccessivamente quando vedono qualcosa di promettente. Questo fa parte della natura umana. penso che fare l'aggiustamento usando FDR come criterio sia un modo ragionevole di mettere sotto controllo l'eccitazione.
Michael R. Chernick,

1
@MichaelChernick, non sono necessariamente in disaccordo con te. Volevo solo esprimere un'altra opinione e spesso mi piace dare una visione d'insieme, semi-filosofica, di cosa si tratta. Molti praticanti possono impantanarsi in dettagli che sembrano arcani per loro e rimangono senza comprensione radicata.
gung - Ripristina Monica

1
Non c'è disaccordo qui e capisco il tuo punto. Voglio solo aggiungere che se potessimo essere dispiaciuti e accettare il principio statistico e non attaccarci personalmente alla nostra ricerca con un interesse acquisito nel risultato, potremmo fare esattamente quello che dici. Ma è così difficile da fare. Immagina di lavorare per un'azienda farmaceutica che ha speso milioni per la ricerca clinica per un particolare farmaco e che fallisce. Il direttore medico ti chiederà di cercare 20 diversi sottogruppi e di trovarne uno che funzioni.
Michael R. Chernick,

1
L'analisi dei sottogruppi è uno degli aspetti più controversi della ricerca clinica. Senza aggiustamento della molteplicità non c'è modo di legittimarlo e farlo post hoc rende difficile vendere alla FDA. Questo è solo un esempio della mia esperienza degli ultimi anni che mi rende sensibile ai suggerimenti di ignorare la molteplicità.
Michael R. Chernick,

-1

Secondo quanto segue: http://birnlab.psychiatry.wisc.edu/resources/fMRI_TestRetest_Documentation.pdf

Quando si determina il significato, correggere il valore p per confronti multipli. Ad esempio, un valore p corretto da Bonferroni è il valore p diviso per il numero totale di confronti, che in questo caso sono le connessioni univoche m (m - 1) / 2.

Ad esempio, il valore p di cutoff per una correlazione è 0,05 e supponiamo che la tabella di correlazione sia 100 * 100. Quindi il tuo valore p dovrebbe essere regolato su 0,05 / (100 * 99/2).

La regressione lineare applica la correzione di Bonferroni in modo simile a quanto sopra.

So che la risposta sembra non essere correlata a ciò che chiedi. In tal caso, per favore fatemi sapere e farò del mio meglio per chiarire. Spero che sia d'aiuto.

Utilizzando il nostro sito, riconosci di aver letto e compreso le nostre Informativa sui cookie e Informativa sulla privacy.
Licensed under cc by-sa 3.0 with attribution required.