2017-02-15 6 views
5

Я пытаюсь провести анализ выживаемости в R, используя пакет survival. У меня есть только цензурные данные (время для события, затем событие «1» или цензор «0»). Я создал объект выживания, используя функцию Surv. Однако survdiff и coxph дают ошибки, предполагая, что они не признают его как цензурные данные. Любой совет?Выживание, не признающее право цензурированных данных

> surv.dfs <- Surv(DaysFromTx,Event) 
> surv.dfs 
[1] 99:1 334:1 1024+ 1341+ 210+ 1069+ 890+ 1242+ 255+ 228+ 349+ 300+ 717+  2+ 657+ 995+ 491+ 1544+ 265:1 440+ 362:1 845+ 
[23] 669+ 1176+ 718+ 768+ 1171+ 2276+ 1152+ 207+ 1138+ 1002+ 942+ 644+ 1110+ 179:1 1535+ 841+ 923+ 904+ 367:1 959+ 746+ 1256+ 
[45] 83:1 439:1 69+ 449+ 591+ 983+ 787+ 704+ 825+ 747+ 28+ 41+ 907+ 181:1 371+ 388+ 166:1 702+ 647+ 944+ 903+ 797+ 
[67] 1095+ 770:1 1118+ 63:1 1762+ 1662+ 127:1 634+ 312+ 483+ 
> survdiff(surv.dfs ~ group) 

Ошибка в survdiff (surv.dfs ~ группа): Право цензурированные данные только

> coxph(surv.dfs ~ group) 

Ошибка в coxph (surv.dfs ~ группа): Cox модель Безразлично 't поддерживать «честные» данные о выживании

+1

Добро пожаловать в Stackoverflow. Можете ли вы предоставить [воспроизводимый пример R] (http://stackoverflow.com/q/5963269/3250126)? 'dput (surv.dfs)' должен помочь с этим. – loki

ответ

7

- Ваше событие факторного класса? он должен быть числовым или логическим. поэтому, если ваши классы событий являются 0 или 1, попробуйте:

df$Event <- as.numeric(df$Event) 
+1

Спасибо, что сработало. Мое событие было классифицировано как фактор. –