Справочная информация : в течение полугода наблюдения в течение 4 лет пациенты могут переключаться на другую группу препаратов. Чтобы учесть это, я преобразовал данные о выживании в форму подсчета. Я хочу сравнить кривые выживаемости для групп препаратов A, B и C. Я использую расширенную модель Кокса, но хочу сделать парные сравнения каждой функции опасности или провести стратифицированные логарифмические тесты. pairwise_survdiff
выдает ошибку из-за формы моих данных, я думаю.
Пример данных :
x<-data.frame(tstart=rep(seq(0,18,6),3),tstop=rep(seq(6,24,6),3), rx = rep(c("A","B","C"),4), death=c(rep(0,11),1))
x
Задача
При использовании survdiff
в пакете survival
,
survdiff(Surv(tstart,tstop,death) ~ rx, data = x)
Я получаю ошибку:
Error in survdiff(Surv(tstart, tstop, death) ~ rx, data = x) :
Right censored data only
Я думаю, что это связано с формой процесса подсчета, так как я не могу найти в Интернете пример, который сравнивает кривые выживаемости для изменяющихся во времени ковариат.
Вопрос : есть ли быстрое решение этой проблемы? Или есть альтернативный пакет / функция с той же универсальностью для сравнения кривых выживания, а именно с использованием разных методов? Как я могу реализовать стратифицированные логарифмические тесты, используя survidff
при подсчете данных формы процесса?
ПРИМЕЧАНИЕ : это было помечено как известная проблема в пакете Survminer, смотрите проблему с GitHub здесь, но обновление SurminMiner не решило мою проблему, и, используя один интервал времени, tstop-tstart не будет правильно, так как это оставит, например, несколько записей на 6 месяцев, а не на фактический интервал риска.