Я работаю в R. У меня есть матрица генной экспрессии со 100 пациентами и 10000 генами. Я делаю кластерный анализ пациентов и получаю разные группы пациентов. Получено 5 кластеров, и каждому пациенту присваивается число 1,2,3,4 или 5 в зависимости от того, к какому кластеру принадлежит пациент. У меня также есть кадр клинических данных в R всех 100 пациентов в матрице экспрессии генов. Клинические данные содержат время выживания в днях, если смерть наступила или нет (0 или 1). Затем я объединяю клинические данные и информацию о кластеризации, и новый фрейм данных выглядит следующим образом:
Я хочу построить графики Шоенфельда для каждого кластера. Это означает, что я хочу 5 участков Шенфельда.
Я использую следующий код, чтобы сделать это:
fitcox <- coxph(Surv(Survival,Death) ~ as.factor(cluster), data = dataX)
summary(fitcox)
ftest <- cox.zph(fitcox)
ftest
ggcoxzph(ftest)
Проблема в том, что я получаю только один график Шоенфельда вместо нескольких графиков (по одному для каждого кластера) в R графиках как на этом рисунке
Буду признателен, если кто-нибудь исправит этот код, чтобы получить графики для каждого кластера. Есть пять кластеров. Если в качестве эталона берется один кластер, то ожидается 4 участка.
Спасибо.