Alguém já usou o procedimento Marascuilo para comparar múltiplas proporções?

O procedimento Marascuilo, conforme descrito aqui, parece ser um teste que aborda o problema de múltiplas comparações de proporções quando você deseja testar quais proporções específicas são diferentes umas das outras após rejeitar o nulo em um teste qui-quadrado geral.

No entanto, não estou muito familiarizado com este teste. Então, minhas perguntas:

1. Com que nuances (se houver) devo me preocupar ao usar este teste?

2. Conheço pelo menos duas outras abordagens (veja abaixo) para resolver o mesmo problema. Qual teste é a melhor abordagem ?:

   • Realizando o "quadrado qui particionado" consultado nesta resposta por @Brett Magill

   • Usando um método de Holm – Bonferroni para ajustar valores de p.

Apenas uma resposta parcial, porque nunca ouvi falar desse método. Pelo que li no link que você forneceu, parece ser um procedimento de etapa única (bem como Bonferroni, exceto que refazemos as estatísticas do teste em vez do valor-p) que provavelmente é muito conservador.

Em R, existe uma função pairwise.prop.test()que permite qualquer correção para comparações múltiplas (métodos FWER de etapa única ou reduzida ou baseada em FDR), mas é encerrada o que você já sugeriu (embora Bonferroni seja conservador demais, mas ainda muito usado na prática). Uma abordagem de reamostragem, usando permutação, também pode ser interessante. O coinpacote R fornece uma estrutura de teste bem estabelecida a esse respeito, consulte o §5 da Implementação de uma classe de testes de permutação: o pacote de moedas , mas nunca tive que lidar com testes de permutação de dados categóricos de maneira post-hoc.

Sobre a análise de tabelas de contingência subdivididas, geralmente considero associações específicas como um guia para desenvolver hipóteses adicionais (como para quaisquer comparações não planejadas), mas essa é outra questão. Geralmente, apenas uso ferramentas de visualização, como mosaico de Michael Friendly , os resíduos de Pearson, e se procuro explicar padrões específicos de associação, uso modelos log-lineares.

Eu gostaria de ver o procedimento Marascuilo usado com mais frequência. Frequentemente, vejo pessoas calculando o qui-quadrado em um subconjunto da tabela principal, ou seja, duas categorias ao mesmo tempo, mas sem realmente fazer o particionamento corretamente. A razão pela qual eles fazem dessa maneira, até onde eu entendi, é que eles não podem suportar agrupar as categorias, pois isso tornará a interpretação realmente difícil. No final do dia, isso depende do público, pois, se não o conhecem, podem recomendar a abordagem Bonferroni usual