Ước tính xác suất sống sót trong R


14

Dựa trên một mẫu của lần sống sót, tôi muốn ước tính xác suất sống sót của thời gian , đối với một số cụ thể , sử dụng công cụ ước tính Kaplan-Meier. Có thể làm điều này trong ? Xin vui lòng, lưu ý rằng không nhất thiết phải là một thời gian sự kiện.nttRt


1
Tất nhiên: xem chức năng Survfit () của gói sinh tồn [loại trợ giúp (gói = "sống sót")]
Stéphane Laurent

3
@Stephane Laurent: Hàm Surfit () cho ra tỷ lệ sống ước tính tại các thời điểm sự kiện. Nhưng tôi muốn có một quy trình tự động để tính toán sự sống sót đó bất cứ lúc nào. Cảm ơn ...
user7064

Sau đó sử dụng khoảngfun ()
Stéphane Laurent

Tôi có thể có một ví dụ?
dùng7064

Câu trả lời:


23

Bạn có thể sử dụng đầu ra của survfithàm từ survivalgói và cung cấp cho stepfun.

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

Bây giờ survestlà một chức năng có thể được đánh giá bất cứ lúc nào.

> 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

Surreal đang gặp vấn đề, tìm kiếm câu trả lời trên StackExchange, tìm ra câu hỏi tương tự và nhận ra bạn là người đã trả lời nó gần 7 năm trước ...
Brian Diggs

7

Một tham số thời gian có thể được chuyển đến chức năng tóm tắt của đối tượng Survfit:

summary(km, times=100)

Một vectơ cũng có thể được thông qua:

summary(km, times=0:100)
Khi sử dụng trang web của chúng tôi, bạn xác nhận rằng bạn đã đọc và hiểu Chính sách cookieChính sách bảo mật của chúng tôi.
Licensed under cc by-sa 3.0 with attribution required.