ggfortify не поддерживает переопределение нескольких коварируемых переменных?

Question

Следующий пример модифицируется из примера документа добавлением «age»

d.coxph <- (survfit(Surv(time, status) ~ sex+age, data = lung))
autoplot(d.coxph)

Я получу следующую ошибку:

Ошибка в levels<- (*tmp*, значение = if (nl == nL) as.character (метки), иначе paste0 (метки,: Уровень фактора [41] дублируется

Введите номер кадра или 0 для выхода

1: автоплот (d.coxph)> Ошибка в levels<- (*tmp*, значение = if (nl == nL) as.character (метки) else paste0 (метки,: Уровень фактора [41] дублируется

Введите номер кадра или 0 для выхода

1: автоплот (d.coxph) 2: autoplot.survfit (d.coxph) 3: укрепить 4: fortify.survfit (объект, Survivconnect = Survivconnect, весело = весело) 5: фактор (rep (groupID, модель $ страты), уровни = groupID)

2: autoplot.survfit (d.coxph) 3: укрепить 4: fortify.survfit (объект, Survivconnect = Survivconnect, весело = весело) 5: фактор (rep (groupID, модель $ страты), уровни = groupID)

Answer 1

Возможное решение - разделить непрерывную переменную age на категории и рассмотреть взаимодействие между полом и возрастом в формуле выживаемости:

library(survival)
data(lung)
lung$age2cat <- cut(lung$age,breaks=2)
lung$sex <- factor(lung$sex, labels=c("F","M"))
d.coxph <- survfit(Surv(time, status) ~ interaction(sex,age2cat), data = lung)

autoplot(d.coxph, conf.int=F, surv.size=1)