R/epoisCensored.jackknife.R

In EnvStats: Package for Environmental Statistics, Including US EPA Guidance

epoisCensored.jackknife <-
function (x, censored, censoring.side, est.fcn, ci.type, conf.level) 
{
    N <- length(x)
    jack.vec <- numeric(N)
    new.N <- N - 1
    for (i in 1:N) {
        new.x <- x[-i]
        new.censored <- censored[-i]
        new.n.cen <- sum(new.censored)
        if (new.n.cen == 0) 
            jack.vec[i] <- mean(new.x)
        else jack.vec[i] <- do.call(est.fcn, list(x = new.x, 
            censored = new.censored, censoring.side = censoring.side, 
            ci = FALSE))$parameters[1]
    }
    jack.vec
}