2017-02-15 40 views
5

survival paketini kullanarak R'de bir sağkalım analizi çalıştırmaya çalışıyorum. Yalnızca sansürlü verilerim var (etkinlik zamanı, sonra olay "1" veya sansür "0"). Surv işlevini kullanarak bir hayatta kalma nesnesi oluşturdum. Bununla birlikte, survdiff ve coxph hatalar vererek, bunların doğru sansürlü veriler olarak tanınmadığını düşündürmektedir. Herhangi bir tavsiye? survdiff içindeSağkalım sansürlenmemiş verileri tanımıyor

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

hatası (surv.dfs ~ grubu): coxph sağ sansürlenmiş veriler sadece

> coxph(surv.dfs ~ group) 

hatası (surv.dfs ~ grubu): Cox modeli doesn Destek "sağır" sağkalım verileri

+1

Welcome Stackoverflow için. [Tekrarlanabilir bir R örneği] sağlayabilir misiniz (http://stackoverflow.com/q/5963269/3250126)? dput (surv.dfs) 'buna yardımcı olmalıdır. – loki

cevap

7

sizin için bir faktör sınıfının mı? sayısal veya boole olmalıdır. senin Olay sınıfları 0 veya 1 eğer öyleyse, deneyin:

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

Teşekkürler, işe yaradı. Etkinliğim bir faktör olarak sınıflandırıldı. –