Модели многочлена для порядковых ответов

Результатом переменной отклика может быть один из ограниченного набора возможных значений. Если существует только два возможных исхода, такие как штекер и розетка для пола, эти ответы называются бинарными ответами. Если существует несколько результатов, то они называются polytomous ответами. Некоторые примеры polytomous ответов включают уровни болезни (умеренный, средний, серьезный), предпочтенные районы, чтобы жить в городе, разновидностях для определенного цветочного типа, и так далее. Иногда среди категорий ответа может быть естественный порядок. Эти ответы называются порядковыми ответами.

Упорядоченное расположение может быть свойственно от выбора категории, таково как индивидуум, которому, не удовлетворяемый, удовлетворяют, или очень удовлетворено онлайновым обслуживанием клиентов. Упорядоченное расположение может также быть введено классификацией скрытой (непрерывной) переменной, такой как в случае индивидуума, находящегося в низком риске, среднем риске или группе высокого риска для разработки определенной болезни, на основе количественной медицинской меры, такой как кровяное давление.

Можно задать модель регрессии многочлена, которая использует естественное упорядоченное расположение среди категорий ответа. Эта порядковая модель описывает отношение между интегральными вероятностями категорий и переменных предикторов.

Различные функции ссылки могут описать это отношение с логитом и пробитом, являющимся наиболее используемым.

  • Логит: ссылка по умолчанию функционирует mnrfit использование для порядковых категорий является функцией ссылки логита. Это моделирует журнал совокупные разногласия. 'link','logit' пара "имя-значение" задает это в mnrfit. Совокупные разногласия журнала являются логарифмом отношения вероятности, что ответ принадлежит категории со значением, меньше чем или равным категории j, P (yc j), и вероятность, что ответ принадлежит категории со значением, больше, чем категория j, P (y> c j).

    Порядковые модели обычно основаны на предположении, что эффекты переменных предикторов являются тем же самым для всех категорий на логарифмическом масштабе. Таким образом, модель имеет различные прерывания, но общие наклоны (коэффициенты) среди категорий. Эта модель называется параллельной регрессией или пропорциональной моделью разногласий. Это - значение по умолчанию для порядковых ответов и 'interactions','off' пара "имя-значение" задает эту модель в mnrfit.

    Пропорциональная модель разногласий

    ln(P(yc1)P(y>c1))=ln(π1π2++πk)=α1+β1X1+β2X2++βpXp,ln(P(yc2)P(y>c2))=ln(π1+π2π3++πk)=α2+β1X1+β2X2++βpXp,ln(P(yck1)P(y>ck1))=ln(π1+π2++πk1πk)=αk1+β1X1+β2X2++βpXp,

    где πj, j = 1, 2..., k, является вероятностями категории.

    Например, для переменной отклика с тремя категориями, существуют 3 – 1 = 2 уравнения можно следующим образом:

    ln(π1π+2π3)=α1+β1X1+β2X2++βpXp,ln(π+1π2π3)=α2+β1X1+β2X2++βpXp.

    Под пропорциональным предположением разногласий частичным эффектом переменного предиктора 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*=βX+ε,

    где остаточный член ε имеет стандартное нормальное распределение. Предположим, что существует следующее отношение между скрытой переменной y* и наблюдаемой переменной y:

    y=c1 ifα0<y*α1,y=c2 ifα1<y*α2,y=ckifαk1<y*αk,

    где α 0 = – ∞ и α k = ∞. Затем интегральная вероятность y, находящегося в категории j или одна из более ранних категорий, P (yc j), равна

    P(ycj)=P(y*<αj)=P(βX+ε<αj)=P(ε<αjβX)=Φ(αjβX),

    где Φ является стандартной нормальной кумулятивной функцией распределения. Таким образом,

    Φ1(P(ycj))=αjβX,

    где αj соответствует точкам разделения скрытой переменной и прерывания в модели регрессии. Это только содержит под предположениями о нормальной скрытой переменной и параллельной регрессии. В более общем плане, для переменной отклика с категориями k и несколькими предикторами, упорядоченная нелинейная модель вероятности с нормальным распределением

    Φ1(P(yc1))=α1+β1X1++βpXp,Φ1(P(yc2))=α2+β1X1++βpXp,Φ1(P(yck1))=αk1+β1X1++βpXp,

    где P (yc j) = π 1 + π 2 +... + π j.

    Коэффициенты указывают на удар модульного изменения в переменном предикторе на вероятности состояния. Положительный коэффициент, β 1, например, указывает на увеличение базовой скрытой переменной с увеличением соответствующего переменного предиктора, X 1. Следовательно, это вызывает уменьшение в P (yc 1) и увеличение P (yc 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.

Смотрите также

| | | |

Похожие темы

Для просмотра документации необходимо авторизоваться на сайте