boot

Running examples in ‘boot-Ex.R’ failed
The error most likely occurred in:

> base::assign(".ptime", proc.time(), pos = "CheckExEnv")
> ### Name: censboot
> ### Title: Bootstrap for Censored Data
> ### Aliases: censboot cens.return
> ### Keywords: survival
> 
> ### ** Examples
> 
> library(survival)

Attaching package: ‘survival’

The following object is masked from ‘package:boot’:

    aml

> # Example 3.9 of Davison and Hinkley (1997) does a bootstrap on some
> # remission times for patients with a type of leukaemia.  The patients
> # were divided into those who received maintenance chemotherapy and 
> # those who did not.  Here we are interested in the median remission 
> # time for the two groups.
> data(aml, package = "boot") # not the version in survival.
> aml.fun <- function(data) {
+      surv <- survfit(Surv(time, cens) ~ group, data = data)
+      out <- NULL
+      st <- 1
+      for (s in 1:length(surv$strata)) {
+           inds <- st:(st + surv$strata[s]-1)
+           md
...[truncated]...
arguments to min; returning Inf
> 
> # Now for an example where a cox regression model has been fitted
> # the data we will look at the melanoma data of Example 7.6 from 
> # Davison and Hinkley (1997).  The fitted model assumes that there
> # is a different survival distribution for the ulcerated and 
> # non-ulcerated groups but that the thickness of the tumour has a
> # common effect.  We will also assume that the censoring distribution
> # is different in different age groups.  The statistic of interest
> # is the linear predictor.  This is returned as the values at a
> # number of equally spaced points in the range of interest.
> data(melanoma, package = "boot")
> library(splines)# for ns
> mel.cox <- coxph(Surv(time, status == 1) ~ ns(thickness, df=4) + strata(ulcer),
+                  data = melanoma)
> mel.surv <- survfit(mel.cox)
Error in a$specials[[i]] : subscript out of bounds
Calls: survfit ... survfit.coxph -> model.matrix.coxph -> delete.response -> off
Execution halted