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

Результатом переменной отклика может быть один из ограниченного набора возможных значений. Если существует только два возможных исхода, такие как штекер и розетка для пола, эти ответы называются бинарными ответами. Если существует несколько результатов, то они называются 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.

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

| | | |

Похожие темы

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