R Under development (unstable) (2024-08-27 r87062 ucrt) -- "Unsuffered Consequences" Copyright (C) 2024 The R Foundation for Statistical Computing Platform: x86_64-w64-mingw32/x64 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > ## Reported by Duncan Murdoch (from StackOverflow) 2021-02-04 > > library(survival) > library(boot) Attaching package: 'boot' The following object is masked from 'package:survival': aml > statMeanSurv <- function(data, var) { + surv <- survfit(Surv(time, cens) ~ 1, data = data) + mean(surv$surv) + var + } > ## both have datasets aml > res <- censboot(boot::aml, statMeanSurv, R = 5L, var = 1, + parallel = "multicore", ncpus = 2L) > res$t [,1] [1,] 1.436232 [2,] 1.517110 [3,] 1.485934 [4,] 1.520575 [5,] 1.548173 > > proc.time() user system elapsed 1.23 0.20 1.40