Результатом переменной отклика может быть один из ограниченного набора возможных значений. Если существует только два возможных исхода, такие как штекер и розетка для пола, эти ответы называются бинарными ответами. Если существует несколько результатов, то они называются polytomous ответами. Некоторые примеры polytomous ответов включают уровни болезни (умеренный, средний, серьезный), предпочтенные районы, чтобы жить в городе, разновидностях для определенного цветочного типа, и так далее. Иногда среди категорий ответа может быть естественный порядок. Эти ответы называются порядковыми ответами.
Упорядоченное расположение может быть свойственно от выбора категории, таково как индивидуум, которому, не удовлетворяемый, удовлетворяют, или очень удовлетворено онлайновым обслуживанием клиентов. Упорядоченное расположение может также быть введено классификацией скрытой (непрерывной) переменной, такой как в случае индивидуума, находящегося в низком риске, среднем риске или группе высокого риска для разработки определенной болезни, на основе количественной медицинской меры, такой как артериальное давление.
Можно задать модель регрессии многочлена, которая использует естественное упорядоченное расположение среди категорий ответа. Эта порядковая модель описывает отношение между интегральными вероятностями категорий и переменных предикторов.
Различные функции ссылки могут описать это отношение с логитом и пробитом, являющимся наиболее используемым.
Логит: ссылка по умолчанию функционирует mnrfit использование для порядковых категорий является функцией ссылки логита. Это моделирует журнал совокупные разногласия. 'link','logit' пара "имя-значение" задает это в mnrfit. Совокупные разногласия журнала являются логарифмом отношения вероятности, что ответ принадлежит категории со значением, меньше чем или равным категории j, P (y ≤ c j), и вероятность, что ответ принадлежит категории со значением, больше, чем категория j, P (y> c j).
Порядковые модели обычно основаны на предположении, что эффекты переменных предикторов являются тем же самым для всех категорий на логарифмическом масштабе. Таким образом, модель имеет различные точки пересечения, но общие наклоны (коэффициенты) среди категорий. Эта модель называется параллельной регрессией или пропорциональной моделью разногласий. Это - значение по умолчанию для порядковых ответов и 'interactions','off' пара "имя-значение" задает эту модель в mnrfit.
Пропорциональная модель разногласий
где πj, j = 1, 2..., k, является вероятностями категории.
Например, для переменной отклика с тремя категориями, существуют 3 – 1 = 2 уравнения можно следующим образом:
Под пропорциональным предположением разногласий частичным эффектом переменного предиктора X является инвариантным к выбору категории переменной отклика, j. Например, если существует три категории, то коэффициенты описывают удар переменного предиктора на относительном риске или регистрируют разногласия значения отклика, находящегося в категории 1 по сравнению с категориями 2 или 3, или в категории 1 или 2 по сравнению с категорией 3.
Таким образом модульное изменение в переменной X 2 означало бы изменение в совокупных разногласиях значения отклика, находящегося в категории 1 по сравнению с категориями 2 или 3 или категорией 1 или 2 по сравнению с категорией 3 на коэффициент exp (β 2), учитывая все остальное равное.
Можно альтернативно подобрать модель с различной точкой пересечения и наклонами среди категорий при помощи 'interactions','on' аргумент пары "имя-значение". Однако с помощью этой опции для порядковых моделей, когда равная модель наклонов является истинными причинами снижение эффективности (вы теряете преимущество оценки меньшего количества параметров).
Пробит: 'link','probit' аргумент пары "имя-значение" использует функцию ссылки пробита, которая основана на нормально распределенном предположении скрытой переменной. Для порядковых переменных отклика это также называется упорядоченной нелинейной моделью вероятности с нормальным распределением. Рассмотрите модель регрессии, которая описывает отношение скрытого variable y* порядкового процесса и вектора из переменных предикторов, X,
где остаточный член ε имеет стандартное нормальное распределение. Предположим, что существует следующее отношение между скрытой переменной y* и наблюдаемой переменной y:
где α 0 = – ∞ и α k = ∞. Затем интегральная вероятность y, находящегося в категории j или одна из более ранних категорий, P (y ≤ c j), равна
где Φ является стандартной нормальной кумулятивной функцией распределения. Таким образом,
где αj соответствует точкам разделения скрытой переменной и точки пересечения в модели регрессии. Это только содержит под предположениями о нормальной скрытой переменной и параллельной регрессии. В более общем плане, для переменной отклика с категориями k и несколькими предикторами, упорядоченная нелинейная модель вероятности с нормальным распределением
где P (y ≤ c j) = π 1 + π 2 +... + π j.
Коэффициенты указывают на удар модульного изменения в переменном предикторе на вероятности состояния. Положительный коэффициент, β 1, например, указывает на увеличение базовой скрытой переменной с увеличением соответствующего переменного предиктора, X 1. Следовательно, это вызывает уменьшение в P (y ≤ c 1) и увеличение P (y ≤ c k).
После оценки коэффициентов модели с помощью mnrfit, можно оценить интегральные вероятности или совокупное число в каждой категории с помощью mnrval с 'type','cumulative' опция пары "имя-значение". mnrval принимает содействующие оценки и статистику модели mnrfit возвращается и оценивает категориальные вероятности или номер в каждой категории и их доверительных интервалах. Можно задать который категория или условные вероятности или числа, чтобы оценить путем изменения значения 'type' аргумент пары "имя-значение".
[1] Маккуллаг, P. и Дж. А. Нелдер. Обобщенные линейные модели. Нью-Йорк: Chapman & Hall, 1990.
[2] Долго, J. S. Модели регрессии для категориальных и ограниченных зависимых переменных. Мудрые публикации, 1997.
[3] Добсон, A. J. и А. Г. Барнетт. Введение в Обобщенные линейные Модели. Чепмен и Холл/CRC. Taylor & Francis Group, 2008.
fitglm | mnrfit | mnrval | glmfit | glmval
