1. Из представленной выборки создать выборку данных, содержащую 40% пациентов, принадлежащих больнице
BINER.
HILLVIEW
ST_AND
2. Создать две случайные выборки, содержащие по 20% (30, 40 %) пациентов общей выборки.
Лабораторная работа № 3. Сравнение выживаемости в двух и более группах
Пример выполнения
Данные о выживаемости пациентов после операции по трансплантации сердца содержатся в файле Heart.sta (файл с данными из предыдущей лабораторной работы).
Откройте этот файл данных через меню Файл – Открыть.
Рисунок 31
Первые шесть переменных в этом множестве данных – даты, точнее, даты трансплантации сердца (в следующей последовательности: месяц-день-год), и даты, когда соответствующий пациент либо умер, либо был изъят из наблюдения (иными словами, цензурирован, например, с пациентом была утрачена связь).
Переменная Цензурировано – Censored является индикатором цензурирования с кодом, который показывает, является соответствующее наблюдение завершенным или цензурированным (0-завершенное; 1-цензурированное). Переменная Больница – Hospital представляет собой фиктивную группирующую переменную, которая показывает, к какой из трех больниц относится пациент.
Вместо записи данных в шести переменных (месяц, день и год начала наблюдения и месяц, день и год окончания наблюдения) мы можем ввести продолжительность жизни с помощью одной переменной или ввести даты начала и окончания как значения двух различных переменных.
В этом исследовании представлены три больницы: Hillview (код 1), St_AND (код 2), and Biner (код 3). На основе сравнения функций выживания в этих трех больницах применением некоторых двувыборочных критериев из модуля Анализ выживаемости можно сделать выводы как о выживаемости пациентов, так и о возможностях этих тестов.
Спецификация анализа.
Откройте файл данных Heart.sta.
Затем выберите Анализ выживаемости из меню Анализ – Углубленные методы анализа для отображения стартовой панели модуля Анализ выживаемости, дважды щелкните по строке Сравнение нескольких выборок, чтобы открыть диалоговое окно Сравнение выживаемости в нескольких группах.
Нажмите кнопку Переменные, чтобы появилось стандартное окно Выберите переменные. Выберите первые 6 переменных в качестве Времен жизни (1, 2 или 6). STATISTICA будет интерпретировать первую и четвертую переменные в этом списке как месяцы, вторую и пятую - как дни, а третью и шестую - как год.
Далее необходимо определить переменную Censored как Цензурирующую переменную и переменную Hospital как Группирующую переменную.
Нажмите кнопку OK для возвращения в диалоговое окно Сравнение выживаемости в нескольких группах. Дважды щелкните по полю ввода Код для полных наблюдений, чтобы появилось диалоговое окно Переменная 7. Здесь выберите Complete и нажмите OK. Затем дважды щелкните по полю Код для цензурированных наблюдений и выберите Censored.
Чтобы ввести коды для группирующей переменной, нажмите кнопку Коды (для групп) и потом в диалоговом окне сначала нажмите кнопку Все (чтобы отобрать все коды), а потом щелкните OK (чтобы ввести эти коды и завершить этот диалог).
Диалоговое окно Сравнение выживаемости в нескольких группах теперь выглядит следующим образом:
Рисунок 32
Теперь, чтобы продолжить, достаточно щелкнуть OK, и откроется диалоговое окно Результаты сравнения выживаемости в нескольких группах.
Рисунок 33
Критерий хи-квадрат дает почти значимое отклонение в этом примере (p <0.06). Отсюда мы можем сделать предварительное заключение, что имеются определенные различия между больницами.
Чтобы увидеть эти различия, мы можем построить графики функций выживания для различных групп, для чего достаточно нажать кнопку График функции выживаемости по группам и получим график – Кумулятивная доля выживших:
Рисунок 34
Видно, что функция выживания в больнице города Бинер медленнее уменьшается на начальном участке времени по сравнению с другими больницами.
Отсюда мы можем заключить, что по какой-то причине пациенты с трансплантированным сердцем в больнице Бинера имеют больше шансов выживания, в частности в течение критических первых 100 дней после трансплантации.
Вернитесь к диалоговому окну Результаты сравнения выживаемости в нескольких группах и нажмите на вкладку Дополнительно. Чтобы просмотреть таблицы времен жизни для этих трех больниц, нажмите кнопку Процент выживших по группам. Появятся отдельные электронные таблицы для каждой больницы.
Рисунок 35
Данные таблицы подтверждают приведенное выше заключение.