R/pred.int.norm.Modified.CA.on.r.K.R

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

pred.int.norm.Modified.CA.on.r.K <-
function (n, df = n - 1, n.mean = 1, r = 1, delta.over.sigma = 0,
    conf.level = 0.95, K.tol = .Machine$double.eps^(1/2),
    integrate.args.list = NULL)
{
    if (!is.vector(n, mode = "numeric") || length(n) != 1 ||
        !is.vector(df, mode = "numeric") || length(df) != 1 ||
        !is.vector(n.mean, mode = "numeric") || length(n.mean) !=
        1 || !is.vector(r, mode = "numeric") || length(r) !=
        1 || !is.vector(delta.over.sigma, mode = "numeric") ||
        length(delta.over.sigma) != 1 || !is.vector(conf.level,
        mode = "numeric") || length(conf.level) != 1)
        stop(paste("'n', 'df', 'r', 'delta.over.sigma',", "and 'conf.level' must be numeric scalars"))
    if (n < 2 || n != trunc(n))
        stop("'n' must be an integer greater than or equal to 2")
    if (df < 1 || df != trunc(df))
        stop("'df' must be an integer greater than or equal to 1")
    if (n.mean < 1 || n.mean != trunc(n.mean))
        stop("'n.mean' must be an integer greater than or equal to 1")
    if (r < 1)
        stop("'r' must be greater than or equal to 1")
    if (!is.finite(delta.over.sigma))
        stop("'delta.over.sigma' must be finite")
    if (conf.level <= 0 || conf.level >= 1)
        stop("'conf.level' must be between 0 and 1")
    fcn.to.min <- function(K, n.weird, df.weird, n.mean, r.weird,
        delta.over.sigma, conf.level, integrate.args.list) {
        (conf.level - pred.int.norm.Modified.CA.on.r.prob(n = n.weird,
            df = df.weird, n.mean = n.mean, K = K, delta.over.sigma = delta.over.sigma,
            r = r.weird, integrate.args.list = integrate.args.list))^2
    }
    K <- nlminb(start = 1, objective = fcn.to.min, lower = 0,
        control = list(x.tol = K.tol), n.weird = n, df.weird = df,
        n.mean = n.mean, r.weird = r, delta.over.sigma = delta.over.sigma,
        conf.level = conf.level, integrate.args.list = integrate.args.list)$par
    K
}