При тестировании предположения параллельной регрессии в порядковой логистической регрессии я обнаружил, что существует несколько подходов. Я использовал как графический подход (как подробно описано в книге Харрелла), так и подход, подробно описанный с использованием порядкового пакета в R.
Однако я также хотел бы запустить тест Бранта (из Stata) как для отдельных переменных, так и для всей модели. Я посмотрел вокруг, но не могу найти его в R.
Есть ли реализация теста Бранта в R?
Ответы:
Я реализовал тест brant в R. Пакет и функция называются brant, и теперь они доступны в CRAN.
Тест Бранта был определен Роллином Брантом для проверки предположения параллельной регрессии (Brant R. (1990). Оценка пропорциональности в модели пропорциональных шансов для порядковой логистической регрессии. Biometrics , 46 , 1171–1178).
Вот пример кода:
В этом примере предположение о параллельной регрессии выполняется, поскольку все значения p превышают 0,05. Омнибус для всей модели, остальное для индивидуальных коэффициентов.
источник
Да, на самом деле, порядковый пакет, который вы связали, может это сделать (хотя они не называют это тестом Бранта). Взгляните на страницы 6 и 7 вашей ссылки, которые демонстрируют «тест отношения правдоподобия для равных угловых коэффициентов или предположения пропорциональных шансов», и это именно то, что вы ищете.
источник
Некоторые заметки по теме
Пакет R
VGAM
вCumulative
команде (Порядковая регрессия с кумулятивными вероятностями) позволяет изменять предположения о пропорциональных коэффициентах с помощью опцииparallel=FALSE
.Известно, что это общая проблема (из книги: Модели регрессии для категориально-зависимых переменных с использованием Stata, второе издание, Дж. Скотт Лонг, Джереми Фриз)
Эта статья углубляется в эту тему, будучи ясной и хорошо написанной, но она не учитывает пакет VGAM или «кумулятивную» команду: Порядковый логистический регресс в эпидемиологических исследованиях
источник
Этот урок о порядковой логистической регрессии в R охватывает проверку предположения о пропорциональных шансах.
источник