2017-02-15 14 views
5

Próbuję uruchomić analizę przeżycia w R za pomocą pakietu survival. Mam tylko dane z prawą cenzurą (czas do zdarzenia, następnie zdarzenie "1" lub cenzor "0"). Stworzyłem obiekt przeżycia przy użyciu funkcji Surv. Jednak survdiff i coxph zawierają błędy, co sugeruje, że nie rozpoznają go jako prawnie cenzurowanego. Jakakolwiek rada?Przetrwanie nie rozpoznaje poprawnie ocenzurowane dane

> 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) 

błędach w survdiff (surv.dfs ~ grupa): Prawo ucięta tylko

> coxph(surv.dfs ~ group) 

błędach w coxph (surv.dfs ~ grupa): modelu Cox robi "t poparcie" w porządku "dane o przeżyciu

+1

Witamy Stackoverflow. Czy możesz podać [odtwarzalny przykład R] (http://stackoverflow.com/q/5963269/3250126)? 'dput (surv.dfs)' powinno pomóc w tym. – loki

Odpowiedz

7

jest Twoim wydarzeniem w klasie czynnikowej? powinien to być numeryczny lub boolowski. więc jeśli klas zdarzeń są 0 lub 1, spróbuj:

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

Dziękuję, to zadziałało. Moje wydarzenie zostało zaklasyfikowane jako czynnik. –