Latvian | Russian
   
Самоучитель по SPSS
Литература по SPSS
Скачать SPSS

PSPP

LimeSurvey survey application

Необходимый файл:
hkarz.sav
 

Бинарная логистическая регрессия

Back Menu Forward

16.4. Бинарная логистическая регрессия

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

Как правило, в случае с дихотомическими переменными речь идёт о некотором событии, которое может произойти или не произойти; бинарная логистическая регрессия в таком случае рассчитывает вероятность наступления события в зависимости от значений независимых переменных.

Вероятность наступления события для некоторого случая рассчитывается по формуле

где z = b1 • X1 + b2 • Х2 + ... + bn • Xn + a,

X1 — значения независимых переменных, b1 — коэффициенты, расчёт которых является задачей бинарной логистической регрессии, а — некоторая константа.

Если для р получится значение меньшее 0,5, то можно предположить, что событие не наступит; в противном случае предполагается наступление события.

В качестве примера рассмотрим два диагностических теста из области медицины на предмет обнаружения карциномы (злокачественной опухоли) мочевого пузыря: подсчет количества (типизация) Т-клеток и тест LAI. Результатами первого теста являются значения, принадлежащие к интервальной шкале, а тест LAI дает дихотомический результат: "положительно" или "отрицательно".

Оба теста были проведены со здоровыми людьми и заведомо больными пациентами. Результаты представлены в следующей таблице:

Коллектив Типизация t-клеток LAI Коллектив Типизация t-клеток LAI
болен 48.5 положительно болен 73.5 положительно
болен 55.5 положительно здоров 61.1 положительно
болен 57.5 положительно здоров 62.5 отрицательно
болен 58.5 положительно здоров 63.5 отрицательно
болен 61.0 положительно здоров 64.5 положительно
болен 61.5 положительно здоров 69.5 положительно
болен 61.5 положительно здоров 70.0 отрицательно
болен 62.0 положительно здоров 70.0 отрицательно
болен 62.0 положительно здоров 71.0 положительно
болен 62,0 положительно здоров 71,5 положительно
болен 62.5 положительно здоров 71.5 отрицательно
болен 63.0 положительно здоров 72.0 отрицательно
болен 63.5 положительно здоров 73.0 отрицательно
болен 65.0 положительно здоров 76.0 отрицательно
болен 65.0 отрицательно здоров 72.5 отрицательно
болен 66.5 отрицательно здоров 73.0 отрицательно
болен 66.5 отрицательно здоров 73.5 отрицательно
болен 66.5 положительно здоров 74.0 отрицательно
болен 68.5 положительно здоров 75.0 отрицательно
болен 69.0 отрицательно здоров 77.0 отрицательно
болен 71.0 положительно здоров 77.0 отрицательно
болен 71.0 положительно здоров 78.5 отрицательно
болен 71.0 положительно

Если сначала посмотреть на результаты типизации Т-клеток, то можно заметить, что здесь для здоровых людей значения в среднем выше, чем для больных. Следовательно, исходя из значений, получившихся при типизации Т-клеток, можно попытаться, вывести вероятность наличия карциномы мочевого пузыря.

Приведенные в таблице данные находятся в файле hkarz.sav. Больным присвоена кодировка 1, а здоровым 2; для теста LA1 кодировка 0 соответствует положительному результату, а 1 отрицательному.

  • Откройте файл hkarz.sav.

  • Выберите в меню Analyze... (Анализ) ► Regression... (Регрессия) ► Binary logistic... (Бинарная логистическая). Открывается диалоговое окно Logistic Regression (Логистическая регрессия).

  • Поместите переменную gruppe (группа), содержащую информацию о принадлежности к одному или второму коллективу (больным или здоровым), в поле для зависимых переменных, а переменную tzell — в поле ковариат. Результаты теста LAI сначала мы не будем использовать в расчёте.

Диалоговое окно Logistic Regression (Логистическая регрессия)

Рис. 16.16: Диалоговое окно Logistic Regression (Логистическая регрессия).

В качестве метода использования переменных в вычислениях предварительно установлен метод Enter (Вложение), при котором в расчёт одновременно вовлекаются все переменные объявленные ковариатами. Альтернативой здесь являются прогрессивная и обратная селекции. В случае наличия лишь одной ковариаты, как в указаном примере, для расчёта подходит только предварительно установленный метод.

Кнопка Select» (Выбрать) предоставляет возможность отбора определённых случаев для дальнейшего анализа.

Используя кнопку Categorical... (Категориальные) Вы можете подготовить для расчета категориальные переменные (то есть переменные, принадлежащие к номинальной шкале и имеющих более 2 значений). На этом мы остановимся более подробно, рассматривая второй пример.

При помощи кнопки Save... (Сохранить) Вы можете добавить в файл дополнительные переменные; активируйте к примеру в разделе Predicted Values (Спрогнозированные значения) предварительные установки Probabilities (Вероятности) и Принадлежность к группе.

Нажав на кнопку Options... (Опции), Вы сможете организовать вывод дополнительных статистических характеристик, различных диаграмм и произвести некоторые дополнительные установки. В данном расчёте мы этого делать не будем.

  • Начните расчёт нажатием ОК.

Наиболее важные результаты приведены в нижеследующей таблице, причём в 10 версии SPSS они уже выводятся в новой табличной форме.

Omnibus Tests of Model Coefficients (Универсальный критерий коэффициентов модели)

Chi-square (Хи-квадрат) Df Sig. (Значимость)
Step 1 (Шаг 1) Step (Шаг) 18,789 1 ,000
Block (Блок) 18,789 1 ,000
Model (Модель) 18,789 1 ,000

Model Summary (Сводная таблица модели)

Step (Шаг) -2 Log Likelihood (-2 логарифмическое правдоподобие) Сох & Snell R Square (R-квадрат Кокса и Шнела) R Square Nadelkerkes (R-квадрат Наделькеркеса)
1 43,394 ,341 ,456

Качество приближения регрессионной модели оценивается при помощи функции подобия. Мерой правдоподобия служит отрицательное удвоенное значение логарифма этой функции (-2LL). В качестве начального значения для -2LL применяется значение, которое получается для регрессионной модели, содержащей только константы. После добавления переменной влияния tzell значение -2LL равно 43,394; это значение на 18,789 меньше, чем начальное. Подобное снижение величины означает улучшение; разность обозначается как величина хи-квадрат и является очень значимой.

Это означает, что начальная модель после добавления переменной tzell претерпела значительное улучшение. Если при наличии некоторого количества независимых переменных анализ производится не при помощи метода вложения, а пошаговым образом, то получающиеся изменения отображаются в разделах "Блок" и "Шаг". При этом, если Вы производили ввод переменных в блочной форме, то показатель в разделе "Блок" приобретает особое значение.

Два других выведенных показателя, названные именами Кокса & Шела и Наделькеркеса, являются мерами определённости. Они также как и при линейной регрессии указывают на ту часть дисперсии, которую можно объяснить с помощью логистической регрессии. Мера определённости по Коксу и Шелу имеет тот недостаток, что значение равное 1 является теоретически не достижимым; этот недостаток устранен благодаря модификации данной меры по методу Наделькеркеса. Часть дисперсии, объяснимой с помощью логистической регрессии, в данном примере составляет 45,6 %.

Далее приводится классификационная таблица, в которой наблюдаемые показатели принадлежности к группе (1 = болен, 2 = здоров) противопоставляются предсказанным на основе рассчитанной модели.

Classification Table (Классификационная таблица) а

Observed (Наблюдаемый показатель) Predicted (Спрогнозировано)
GRUPPE (Группа) Percentage Correct (Процентный показатель верных показателей)
Krank (болен) Gesund (здоров)
Шаг 1 GRUPPE (Группа) Krank (болен) 18 6 75,0
Gesund (здоров) 4 17 81,0
Overall Percentage (Суммарный процентный показатель) 77,8

a. The cut value is ,500 (Разделительное значение равно ,500)

Из таблицы можно сделать вывод о том, что из общего числа больных, равного 24, тестом были признаны таковыми только 18 (в медицинской диагностике в таких случаях говорят о "строго положительных" результатах). Остальных 6 называют "ложно отрицательными"; они были признаны тестом здоровыми, хотя и являются больными. Из общего числа здоровых, равного 21, тестом были признаны таковыми только 17 ("строго отрицательные"), 4 признаны больными, хотя они и являются здоровыми ("ложно положительные"). В общем, правильно были распознаны 35 случаев из 45, это составляет 77,8 %.

В заключении выводятся результаты о рассчитанных коэффициентах и проверке их значимости:

Variables in the Equation (Переменные в уравнении)

В (Коэффициент регрессии В) S.E. (Стандартная ошибка) Wald (Вальд) df Sig. (Значимость) Ехр (В)
Step 1 (Шаг 1)а TZELL ,278 ,082 11,599 1 ,001 1,321
Constant (Константа) -19,005 5,587 11,571 1 ,001 ,000

a. Variable(s) entered on step 1: TZELL (Переменные, введенные на шаге 1: TZELL)

Проверка значимости отличия коэффициентов от нуля, проводится при помощи статистики Вальда, использующей распределение хи-квадрат, которая представляет собой квадрат отношения соответствующего коэффициента к его стандартной ошибке.

В приведенном примере получились сверх значимые коэффициенты а = -19,005 b1 = 0,278. При помощи этих двух значений коэффициентов мы можем для каждого значения Т-типизации рассчитать вероятность р. К примеру, для некоего обследуемого со значением Т-типизации 72 получим:

z = -19,005 + 0,278 • 72 = 1,018

и таким образом

Вероятность выздоровления

Рассчитанная вероятность р всегда указывает на исполнение предсказании, которое соответствует большей из двух кодировок зависимых переменных, в данном случае — на исполнение предсказания "здоров". Следовательно, рассматриваемый человек является здоровым с вероятностью 0,735.

Рассчитанная вероятность для всех случаев и связанная с ней принадлежность к группе кодировка 1 для болен и 2 для здоров) добавлены к файлу под именами рrе_1 и pgr_l.

Теперь подключим к нашему анализу тест LAI. Дополнительно к переменной tzell теперь в поле ковариат поместите и переменную lai.

Расчёт выдаст сначала заметно снизившееся значение -2LL (хи-квадрат = 25,668) и следующую классификационную таблицу. Доля правильно спрогнозированных диагнозов незначительно выросла (с 77,8 % до 80,0 %).

Classification Table (Классификационная таблица)а

Observed (Наблюдаемый показатель) Predicted (Спрогнозировано)
GRUPPE (Группа) Percentage Correct (Процентный показатель верных показателей)
Krank (болен) Gesund (здоров)
Шаг 1 GRUPPE (Группа) Krank (болен) 20 4 83,3
Gesund (здоров) 5 16 76,2
Overall Percentage (Суммарный процентный показатель) 80,0

а. The cut value is ,500 (Разделительное значение равно ,500)

Количество ложно отрицательных диагнозов снизилось на 2, а количество ложно положительных повысилось на 1. Для коэффициентов получим:

Variables in the Equation (Переменные в уравнении)

В (Коэффициент регрессии В) S.E. (Стандартная ошибка) Wald (Вальд) df Sig. (Значимость) Ехр (В)
Step 1 (Шаг 1) TZELL ,201 ,094 4,574 1 0,32 1,222
LAI 2,205 ,877 6,324 1 ,012 9,074
Constant (Константа) -14,645 6,328 5,356 1 ,021 ,000

a. Variable(s) entered on step 1: TZELL, LAI. (Переменные, вводимые на шаге 1: TZELL, LAI)

Для обследуемого с типизированным числом Т-клеток равным 72 получилась вероятность оказаться здоровым р = 0,735. Если в дополнении к этому и тест LAI отрицателен (кодировка 1), то эта же вероятность рассчитывается следующим образом:

Вероятность оказаться здоровым

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

Ещё один пример из области медицины, теперь уже с большим количеством независимых переменных, должен помочь нам разобраться в пошаговом методе анализа. Кроме того, в состав независимых переменных будет включена категориальная переменная.

Для данного примера в некоторой клинике со специальными автоматизированными методиками лечения были накоплены данные о пациентах с тяжёлыми (или даже смертельными) повреждениями лёгких. Из большого количества переменных были выбраны следующие:

Имя переменной Расшифровка
out Исход (0 = скончался, 1 = выздоровел)
alter Возраст
bzeit Время проведения искусственного дыхания в часах
kob Концентрация кислорода в воздушной массе для искусственного дыхания
agg Интенсивность искусственного дыхания
geschl Пол (1 = мужской, 2 = женский)
gr Рост
ursache (причина) Причина повреждения лёгких (1 = несчастный случай, 2 = воспаление лёгких, 3 = прочее)

Наряду с переменной out (исход), имеются переменные, при первом же взгляде на которые можно понять, что они с ней связаны. Причина повреждения лёгких является категориальной переменной, которая перед проведением анализа должна быть преобразована в несколько дихотомических переменных (к примеру, несчастный случай: да — нет).

Вопрос, на который нам предстоит найти ответ, звучит так: какое влияние на вероятность выздоровления оказывают отобранные переменные.

  • Откройте файл lunge.sav.

  • После выбора соответствующего меню в диалоговом окне Logistic Regression (Логистическая регрессия) переменной out присвойте статус независимой переменной, а всем остальным (кроме nr) присвойте статус ковариат. Здесь, как и при множественной линейной регрессии, ввод ковариат Вы можете производить по блокам.

Из-за вовлечения в анализ большого количества переменных компьютер должен решить, какие из них в конечном случае будут отобраны для использовании в уравнении вероятности. Поэтому здесь должен быть выбран не метод вложения, который включает в расчёт все переменные, а один из пошаговых методов.

Метод прямой селекции начинается с использования одних лишь констант на стартовом этапе, а затем последовательно подключаются переменные, которые демонстрируют сильную корреляцию с зависимыми переменными. Далее опять следует проверка того, какие переменные должны быть исключены, причём в качестве критерия проверки выбирается либо статистика Вальдовского (Wald), либо функция правдоподобия, либо один из вариантов, называемых "условной статистикой" (которые, однако, не рекомендуются). Метод обратной селекции сначала берёт в расчёт все переменные, а затем в обратном порядке происходит исключение малозначимых переменных.

  • Выберите в качестве метода Forward: LR (Прямой:LR) и щёлкните на кнопке Categorical... (Категориальные), чтобы поместить переменную ursache в поле, предусмотренное для категориальных ковариат.

Количество образовываемых "фиктивных" дихотомических переменных должно быть всегда на 1 меньше, чем число количество заданных категорий. Категория, оказавшаяся лишней, называется эталонной категорией и, в соответствии с предварительными установками, является последней категорией. При помощи поля контрастов (Contrast) Вы можете управлять особенностями вовлечения в анализ образованных Фиктивных переменных; при контрасте равном Deviation (Отклонение) все категории кроме эталонной будут проверяются относительно суммарного эффекта.

  • Установите контраст Deviation (Отклонение) и при помощи щелчка на Continue (Далее) вернитесь в исходное диалоговое окно.

  • Начните расчёт нажатием ОК.

Вы можете проследить, какие переменные вовлекаются в анализ и как улучшается вероятность прогноза после вовлечения каждой новой переменной. На завершающей стати анализа присутствуют четыре переменные, а именно: возраст, время проведения искусственного дыхания, рост и концентрация кислорода в воздушной массе для искусственного дыхания.

Точность исполнения прогноза, которая достигается при использовании этих четыpex переменных, составляет 71,0%; её можно увидеть в нижеследующей классификанионной таблице.

Classification Table (Классификационная таблица)а

Observed (Наблюдаемый показатель) Predicted (Спрогнозировано)
Outcome (Исход) Percentage Correct (Процентный показатель верных прогнозов)
gestorben (скончался) ueberlebt (выздоровел)
Step 1 (Шаг 1) Outcome (Исход) gestorben (скончался) 29 34 46,0
ueberlebt (выздоровел) 14 54 79,4
Overall Percentage (Суммарный процентный показатель) 63,4
Step 2 (Шаг 2) Outcome (Исход) gestorben (скончался) 32 31 50,8
ueberlebt (выздоровел) 16 52 76,5
Overall Percentage (Суммарный процентный показатель) 64,1
StepS (Шаг 3) Outcome (Исход) gestorben (скончался) 33 30 52,4
ueberlebt (выздоровел) 19 49 72,1
Overall Percentage (Суммарный процентный показатель) 62,6
Step 4 (Шаг 4) Outcome (Исход) gestorben (скончался) 37 26 58,7
ueberlebt (выздоровел) 12 56 82,4
Overall Percentage (Суммарный процентный показатель) 71,0

a. The cut value is ,500 (Разделительное значение равно ,500)

Прогноз оправдался для 58,7% умерших пациентов и для 82,4% выздоровевших. Значения коэффициента b. и константы а для расчёта вероятности (выздоровления) находятся в следующей таблице:

Variables in the Equation (Переменные в уравнении)

В (Коэффициент регрессии В) S.E. (Стандартная ошибка) Wald
(Вальд)
df Sig. (Значимость) Ехр(В)
Шаг 1а BZEIT -,081 ,028 8,482 1 ,004 ,922
Константа 1,104 ,385 8,205 1 ,004 3,017
Шаг 2b GR ,038 ,017 5,109 1 ,024 1,039
BZEIT -,073 ,028 6,688 1 ,010 ,930
Константа -5,460 2,924 3,487 1 ,062 ,004
Шаг 3c KОВ -2,678 1,264 4,489 1 ,034 ,069
GR ,037 ,017 4,622 1 ,032 1,038
BZEIT -,077 ,029 6,866 1 ,009 ,926
Константа -2,995 3,192 ,880 1 ,348 ,050
Шаг 4d ALTER -,037 ,017 4,653 1 ,031 ,963
КОВ -3,028 1,302 5,410 1 ,020 ,048
GR ,044 ,017 6,650 1 ,010 1,045
BZEIT -,062 ,029 4,639 1 ,031 ,940
Константа -2,884 3,079 ,877 1 ,349 ,056

a. Variable(s) entered on step 1: BZEIT. (Переменные, вводимые на шаге 1: BZEIT.)
b. Variable(s) entered on step 2: GR. (Переменные, вводимые на шаге 2: GR.)
с. Variable(s) entered on step 3: КОВ. (Переменные, вводимые на шаге 3: КОВ.)
d. Variable(s) entered on step 4: ALTER. (Переменные, вводимые на шаге 4: ALTER.)

Если мы рассмотрим случай с 30-тилешим пациентом, с ростом 180 см, которому делали искусственное дыхание в течении 10 часов при концентрации кислорода в смеси равной 0,7, то исходя из соотношения

z = -2,884 - 0,037 • 30 - 0,062 • 10 + 0,044 • 180 -3,028 • 0,7 = 1,126 получим вероятность выздоровления

следовательно, вероятность выздоровления пациента равна 0,755 или 75,5%.


Back Menu Forward


Top.LV Latvijas Reitingi e-TIRGUS.LV on-line.lv Izglītība Education Яндекс.Метрика
Клиенты

Roche

TNS Latvija

Latvijas Pilsoniskā alianse

LU FSI

Sendigo

Prime Mail

Valodu mācību centrs

Latvijas infektoloģijas centrs
Карта сайта
www.citariga.lv || www.limesurvey.ru || www.exsobalt.lv