Как с помощью R можно разработать показатель индекса для прогнозирования общей выживаемости пациентов (ОВ)?
У меня есть короткий список из 4 предикторов-кандидатов, которые показали связь с ОС. Они получены в результате многомерной регрессии Кокса (прогон с coxph()
). Предикторами являются уровни белка, следовательно, все они являются непрерывными переменными.
Таблица данных выглядит примерно так (здесь показано только n=10):
days Status Prot1 Prot13 Prot7 Prot21
Subj_1 115.69 0 2.284498 6.319168 6.070115 8.457412
Subj_2 72.30 1 2.473034 6.066573 6.140178 8.225987
Subj_3 1.08 1 2.662481 6.212845 6.971018 8.128949
Subj_4 69.63 1 2.761391 5.902610 6.433883 7.876319
Subj_5 78.41 1 3.038122 6.355257 6.852981 7.500973
Subj_6 42.90 1 2.058549 6.020681 7.231307 8.164025
Subj_7 31.00 1 2.305096 5.415107 8.126941 8.566320
Subj_8 51.12 1 2.931978 5.574601 7.503275 7.529957
Subj_9 11.01 1 2.218814 6.270222 6.710297 8.193895
Subj_10 27.68 1 2.821947 6.132379 6.911071 8.428218
Возникает вопрос: Как я могу создать формулу, которая способна разделить этих пациентов на 2 группы: группу, в которой расчетная выживаемость составляет ‹60% в течение 1 года, и другую, в которую войдут пациенты с расчетной выживаемостью. > 60 % за тот же период времени?
Будет ли какой-нибудь function()
в R, который занимается этим?
Заранее большое спасибо.