Я пишу следующий код для вычисления оценки Каплана-Мейера в R.
library(survival)
data(veteran)
attach(veteran)
library(Rfast)
KM <- function(u){
t <- sort_unique(u)
Y <- sapply(t, function(x) sum(u>=x))
d <- sapply(t, function(x) sum(u==x & status == 1))
cumprod(1-d/Y)
}
KM(time) # Kaplan-Meier estimates value
Но я получаю неправильные результаты, когда запускаю следующий код для оценки интеграла 1
integrate(Vectorize(KM),10,999)
Было бы очень полезно, если бы кто-нибудь помог мне написать код R для оценки интеграла 1 в R.