JackKnifeAcc: Calculation of JeckKnife acceleration parameter for BCA...

In MaciejDanko/BootstrapLT: CI for pash measures

Description

Calculation of JeckKnife acceleration parameter for BCA method.

Internal function

Usage

JackKnifeAcc(JackKnifeMat, dx)

Arguments

JackKnifeMat	JackKnife matrix with one row per each considered pash measure. See JackKnife_dx.
dx	A vector with integer counts (deaths).

Author(s)

Maciej J. Danko <[e-mail]> <[e-mail]>

References

Efron, B., & Tibshirani, R. J. (1993). An introduction to the bootstrap. New York: Chapman & Hall. Page 186.

Examples

# ************************** REMOVE THIS ******************************
# Some functions will not work as they are not in the namespace
# Before integrating with pash please add "pash:::" in front of them. 
# To test them just run all functions in the bootpash.r 
# ***************************** END ***********************************
## Not run: 

#Get some data
population.size <- 10000
obj <- Inputlx(x = australia_10y$x, lx = australia_10y$lx,nax = australia_10y$nax, nx = australia_10y$nx, last_open = TRUE)
dx <- lx2dx(round(object$lt$lx*population.size))
x <- obj$lt$x

#Acceleration parameters calculated by slow method:
times <- dx2age(dx,x)
gp <- function(times, x, pash.parent) getpash(dx = age2dx(times = times, x=x), x=x, pash.parent = pash.parent)
J1 <- as.matrix(sapply(seq_along(times), function(k) gp(times = times[-k], x = x, pash.parent = obj)))
a.slow <- JackKnifeAcc(J1)

#Acceleration parameters calculated by fast method:
J2 <- JackKnife_dx(dx = org.dx, x = org.x, pash.parent = obj)
a.fast <- JackKnifeAcc(J2,dx)

#Compare results
cbind(a.slow = a.slow, a.fast = a.fast, diff = round(a.slow - a.fast, 12))


## End(Not run)

MaciejDanko/BootstrapLT documentation built on May 14, 2019, 2:40 p.m.