Ao testar a suposição de regressão paralela na regressão logística ordinal, percebo que existem várias abordagens. Eu usei a abordagem gráfica (conforme detalhado no livro de Harrell) e a abordagem detalhada usando o pacote ordinal em R.
No entanto, eu também gostaria de executar o teste Brant (da Stata) para as variáveis individuais e também para o modelo total. Eu olhei em volta, mas não consigo encontrá-lo implementado em R.
Existe uma implementação do teste de Brant em R?
Respostas:
Eu implementei o teste brant em R. O pacote e a função são chamados brant e agora estão disponíveis no CRAN.
O teste de Brant foi definido por Rollin Brant para testar a suposição de regressão paralela (Brant, R. (1990) Avaliando a proporcionalidade no modelo de chances proporcionais para a regressão logística ordinal. Biometrics , 46 , 1171-1178).
Aqui está um exemplo de código:
No exemplo, a suposição de regressão paralela é válida, porque todos os valores de p estão acima de 0,05. O Omnibus é para todo o modelo, o restante para os coeficientes individuais.
fonte
Sim - na verdade, o pacote ordinal que você vinculou pode fazê-lo (embora não o chame de teste de Brant). Dê uma olhada nas páginas 6 e 7 do seu link, que demonstram "um teste da razão de verossimilhança das pistas iguais ou da suposição de chances proporcionais", exatamente o que você está procurando.
fonte
Algumas notas sobre o tópico
O pacote R
VGAM
noCumulative
comando (Ordinal regressão com probabilidades cumulativas) permite alterar as odds suposições proporcionais, com a opçãoparallel=FALSE
.Sabe-se que é um problema comum (do livro: Modelos de regressão para variáveis dependentes categóricas usando Stata, segunda edição, por J. Scott Long, Jeremy Freese)
Este artigo é aprofundado neste tópico, sendo claro e bem escrito, mas não considera o pacote VGAM ou o comando "cumulativo": Regressão logística ordinal em estudos epidemiológicos
fonte
Este tutorial sobre regressão logística ordinal em R cobre o teste da premissa de chances proporcionais.
fonte