Qualcuno ha usato la procedura Marascuilo per confrontare più proporzioni?

La procedura Marascuilo come qui descritta sembra essere un test che affronta il problema di confronti multipli per le proporzioni quando si desidera verificare quali proporzioni specifiche sono diverse l'una dall'altra dopo aver rifiutato il valore nullo in un test chi-quadro complessivo.

Tuttavia, non ho molta familiarità con questo test. Quindi, le mie domande:

1. Di quali sfumature dovrei preoccuparmi quando uso questo test?

2. Conosco almeno altri due approcci (vedi sotto) per affrontare lo stesso problema. Quale test è l'approccio "migliore" ?:

   • L'esecuzione di "chi square partizionato" è stata confermata in questa risposta da @Brett Magill

   • Utilizzo di un metodo Holm-Bonferroni per regolare i valori p.

Solo una risposta parziale perché non ho mai sentito parlare di questo metodo. Da quanto ho letto nel link che hai fornito, sembra essere una procedura a singolo passaggio (molto simile a Bonferroni, tranne per il fatto che rielaboriamo le statistiche dei test anziché il valore p) che è probabilmente troppo conservativo.

In R, c'è una funzione pairwise.prop.test()che consente qualsiasi correzione per confronti multipli (metodi FWER a singolo passaggio o step-down o basati su FDR), ma è uscito da ciò che hai già suggerito (anche se Bonferroni è di gran lunga troppo conservativo, ma comunque molto usato in pratica). Anche un approccio di ricampionamento, usando la permutazione, potrebbe essere interessante. Il coinpacchetto R fornisce un framework di test ben consolidato a questo proposito, vedere §5 di Implementazione di una classe di test di permutazione: il pacchetto di monete , ma non ho mai avuto a che fare con test di permutazione su dati categorici in un modo post-hoc.

Per quanto riguarda l'analisi delle tabelle di contingenza suddivise, in genere considero specifiche associazioni come una guida per sviluppare ulteriori ipotesi (come per qualsiasi confronto non pianificato), ma questa è un'altra domanda. In genere utilizzo solo strumenti di visualizzazione, come mosaicplot di Michael Friendly , i residui di Pearson, e se cerco di spiegare specifici schemi di associazione uso modelli log-linear.

Mi piacerebbe vedere la procedura Marascuilo utilizzata più spesso. Molto spesso vedo persone che calcolano il chi-quadrato su un sottoinsieme della tabella principale, cioè due categorie alla volta, ma senza effettivamente eseguire correttamente il partizionamento. Il motivo per cui lo fanno in questo modo per quanto ho capito è che non riescono a sopportare il raggruppamento delle categorie in quanto ciò renderà l'interpretazione davvero difficile. Alla fine della giornata dipende anche dal pubblico perché se non lo conoscono potrebbero semplicemente raccomandare il solito approccio Bonferroni