Memperkirakan probabilitas kelangsungan hidup di R

14

Berdasarkan sampel dari kali bertahan hidup, saya ingin memperkirakan probabilitas bertahan waktu , untuk beberapa tertentu , menggunakan penduga Kaplan-Meier. Apakah mungkin untuk melakukan ini ? Harap dicatat, belum tentu waktu acara.nttRt

pengguna7064
sumber
1
Tentu saja: lihat fungsi survfit () dari paket survival [ketik bantuan (package = "survival")]
Stéphane Laurent
3
@Stephane Laurent: Fungsi surfit () menampilkan perkiraan kelangsungan hidup pada waktu kejadian. Tapi saya ingin memiliki prosedur otomatis untuk menghitung kelangsungan hidup itu kapan saja t. Terima kasih ...
user7064
Kemudian gunakan kira-kira ()
Stéphane Laurent
dapatkah saya memiliki contoh?
user7064

Jawaban:

23

Anda dapat menggunakan output survfitfungsi dari survivalpaket dan memberikannya kepada stepfun.

km <- survfit(Surv(time, status)~1, data=veteran)
survest <- stepfun(km$time, c(1, km$surv))

Sekarang survestadalah fungsi yang dapat dievaluasi kapan saja.

> survest(0:100)
  [1] 1.0000000 0.9854015 0.9781022 0.9708029 0.9635036 0.9635036 0.9635036
  [8] 0.9416058 0.9124088 0.9124088 0.8978102 0.8905109 0.8759124 0.8613139
 [15] 0.8613139 0.8467153 0.8394161 0.8394161 0.8175182 0.8029197 0.7883212
 [22] 0.7737226 0.7664234 0.7664234 0.7518248 0.7299270 0.7299270 0.7225540
 [29] 0.7225540 0.7151810 0.7004350 0.6856890 0.6856890 0.6783160 0.6783160
 [36] 0.6709430 0.6635700 0.6635700 0.6635700 0.6635700 0.6635700 0.6635700
 [43] 0.6561970 0.6488240 0.6414510 0.6340780 0.6340780 0.6340780 0.6267050
 [50] 0.6193320 0.6193320 0.5972130 0.5750940 0.5677210 0.5529750 0.5529750
 [57] 0.5456020 0.5456020 0.5456020 0.5382290 0.5382290 0.5308560 0.5308560
 [64] 0.5234830 0.5234830 0.5234830 0.5234830 0.5234830 0.5234830 0.5234830
 [71] 0.5234830 0.5234830 0.5161100 0.5087370 0.5087370 0.5087370 0.5087370
 [78] 0.5087370 0.5087370 0.5087370 0.4939910 0.4939910 0.4866180 0.4866180
 [85] 0.4791316 0.4791316 0.4791316 0.4716451 0.4716451 0.4716451 0.4640380
 [92] 0.4640380 0.4564308 0.4564308 0.4564308 0.4412164 0.4412164 0.4412164
 [99] 0.4412164 0.4257351 0.4179945
Brian Diggs
sumber
Surreal mengalami masalah, mencari jawaban di StackExchange, menemukan pertanyaan yang sama, dan menyadari bahwa Andalah yang menjawabnya hampir 7 tahun yang lalu ...
Brian Diggs
7

Parameter waktu dapat diteruskan ke fungsi ringkasan objek survfit:

summary(km, times=100)

Vektor juga dapat dikirimkan:

summary(km, times=0:100)
bensentropi
sumber