Форум врачей-аспирантов

Здравствуйте, гость ( Вход | Регистрация )

> Анализ выживаемости в R
kont
сообщение 5.09.2016 - 16:24
Сообщение #1





Группа: Пользователи
Сообщений: 149
Регистрация: 11.02.2014
Пользователь №: 26005



Подскажите, пожалуйста, как в R мне разделить график Каплана -Мейерса по группам, чтобы каждая группа имела свой цвет.
Например
# Kaplan-Meier non-parametric analysis by group
kmsurvival1 <- survfit(Surv(kaplan17$time, kaplan17$event) ~ kaplan17$group)
plot(kmsurvival1, xlab="Time", ylab="Survival Probability")

График выглядит не информативно.
http://funkyimg.com/view/2gk7a
Как сделать, чтобы первая группа имела зеленый цвет, а вторая красный.
Вернуться в начало страницы
 
+Ответить с цитированием данного сообщения
 
 
Открыть тему
Ответов
DrgLena
сообщение 8.09.2016 - 12:17
Сообщение #2





Группа: Пользователи
Сообщений: 1325
Регистрация: 27.11.2007
Пользователь №: 4573



К-М строится по неполным данным, если у вас есть для каждого временного интервала полные данные, жив - умер, то вам не нужны К-М. Для неполных данных вводятся некоторые понятия для каждого интервала есть понятия введенных в интервал и подверженных риску. число подверженых риску=число введенных в интервал минус половина не умерших. С сензурированными данными знакомы? Постройте таблицы выживаемости и укажите 8 интервалов и получите эти цифры. В R я это не делала. Кокс регрессия с предиктором группа странно выглядит, там есть инфекция, вроде, в данных.
Вернуться в начало страницы
 
+Ответить с цитированием данного сообщения
 

Сообщений в этой теме


Добавить ответ в эту темуОткрыть тему