Сравнение выживаемости в двух и более группах

Обзор и файл данных. Следующие четыре примера основаны на данных Crowley and Hu (1977) о выживаемости пациентов после операции по трансплантации сердца. Данные содержатся в файле Heart.sta.

Откройте этот файл данных через меню Файл - Открыть.

Скорее всего, он будет расположен в каталоге /Examples/Datasets системы STATISTICA.

Первые шесть переменных в этом множестве данных – даты, точнее, даты трансплантации сердца (в следующей последовательности: месяц-день-год), и даты, когда соответствующий пациент либо умер, либо был изъят из наблюдения (иными словами, цензурирован, например, с пациентом была утрачена связь).

Переменная Цензурировано - Censored является индикатором цензурирования с кодом, который показывает, является соответствующее наблюдение завершенным или цензурированным (0-завершенное; 1-цензурированное). Переменная Больница - Hospital представляет собой фиктивную группирующую переменную, которая показывает, к какой из трех больниц относится пациент.

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

Переменная Hospital в файле данных Heart.sta есть (фиктивная) группирующая переменная, которая идентифицирует больницу, в которую пациент был принят для проведения трансплантации сердца.

В этом исследовании представлены три больницы: Hillview (код 1), St_Andreas (код 2), and Biner (код 3). На основе сравнения функций выживания в этих трех больницах применением некоторых двувыборочных критериев из модуля Анализ выживаемости можно сделать выводы как о выживаемости пациентов, так и о возможностях этих тестов.

Спецификация анализа. Откройте файл данных Heart.sta.

Затем выберите Анализ выживаемости из меню Анализ - Углубленные методы анализа для отображения стартовой панели модуля Анализ выживаемости. Дважды щелкните по строке Сравнение нескольких выборок, чтобы открыть диалоговое окно Сравнение выживаемости в нескольких группах.

Нажмите кнопку Переменные, чтобы появилось стандартное окно Выберите переменные. Выберите первые 6 переменных в качестве Времен жизни (1, 2 или 6). STATISTICA будет интерпретировать первую и четвертую переменные в этом списке как месяцы, вторую и пятую - как дни, а третью и шестую - как год.

Далее необходимо определить переменную Censored как Цензурирующую переменную и переменную Hospital как Группирующую переменную.

Нажмите кнопку OK для возвращения в диалоговое окно Сравнение выживаемости в нескольких группах. Дважды щелкните по полю ввода Код для полных наблюдений, чтобы появилось диалоговое окно Переменная 7. Здесь выберите Complete и нажмите OK. Затем дважды щелкните по полю Код для цензурированных наблюдений и выберите Censored.

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

Диалоговое окно Сравнение выживаемости в нескольких группах теперь выглядит следующим образом:

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

Критерий хи-квадрат дает почти значимое отклонение в этом примере (p<.06). Отсюда мы можем сделать предварительное заключение, что имеются определенные различия между больницами.

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

Ясно видно, что функция выживания в больнице города Бинер медленнее уменьшается на начальном участке времени по сравнению с другими больницами.

Отсюда мы можем заключить, что по какой-то причине пациенты с трансплантированным сердцем в больнице Бинера имеют больше шансов выживания, в частности в течение критических первых 100 дней после трансплантации.

Вернитесь к диалоговому окну Результаты сравнения выживаемости в нескольких группах и нажмите на вкладку Дополнительно. Чтобы просмотреть таблицы времен жизни для этих трех больниц, нажмите кнопку Процент выживших по группам. Появятся отдельные электронные таблицы для каждой больницы. Данные таблицы подтверждают приведенное выше заключение.

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

Нажмите кнопку Отмена в диалоговом окне Результаты сравнения выживаемости в нескольких группах и в диалоговом окне Сравнение выживаемости в нескольких группах для возвращения в стартовую панель.

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

Чтобы сравнить группы, дважды щелкните по полю ввода Код для первой группы для отображения диалогового окна Переменная 11. Выберите здесь первую группу для сравнения (в нашем примере выбрана Hillview) и нажмите OK.

Потом дважды щелкните по полю ввода Код для второй группы и выберите Biner. Диалоговое окно Сравнение выживаемости в двух группах теперь выглядит следующим образом:

Теперь нажмите кнопку OK, чтобы перейти к диалоговому окну Двувыборочные критерии.

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

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

Связанные определения:
Гомперца распределение
Интенсивность отказов
Таблица времён жизни
Функция выживаемости
Цензурированные данные (неполные данные)

В начало

Содержание портала