boot.clos: Bootstrap LOS

In n8thangreen/HESmanip: Hospital length of stay data preprocessing and analysis

Description

Bootstrap standard error estimate for length of stay.

Usage

boot.clos(data, state.names, tra, cens.name, s = 0, nboot)

Arguments

data
state.names
tra
cens.name
s
nboot

See Also

boot.clos.dtable

Examples

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (data, state.names, tra, cens.name, s = 0, nboot) 
{
    res <- double(nboot)
    for (i in seq_len(nboot)) {
        index.inf <- sample(x = unique(data$id[data$to == 1]), 
            replace = TRUE)
        index.mix <- sample(x = unique(data$id[data$from == 0 & 
            data$to == 2]), replace = TRUE)
        index <- c(index.inf, index.mix)
        linds <- sapply(index, function(x) which(x == data$id))
        indrep <- sapply(linds, function(y) length(y))
        inds <- unlist(linds)
        new.id <- rep(seq_along(index), indrep)
        dboot <- cbind(data[inds, ], new.id)
        dboot$id <- dboot$new.id
        dboot$new.id <- NULL
        tr.prob <- etm::etm(dboot, state.names, tra, cens.name, 
            s, cova = FALSE)
        res[i] <- etm::clos(tr.prob)$e.phi
    }
    res
  }

n8thangreen/HESmanip documentation built on March 21, 2020, 12:20 a.m.