R/confint.mlePT.R

In tweeDEseq: RNA-seq data analysis using the Poisson-Tweedie family of distributions

confint.mlePT <- function(object, parm, level=0.95, ...)
 {
 
    cf <- object$par
    pnames <- names(cf)
    if (missing(parm)) 
        parm <- pnames
    else if (is.numeric(parm)) 
        parm <- pnames[parm]
    a <- (1 - level)/2
    a <- c(a, 1 - a)
    pct <- stats:::format.perc(a, 3)
    fac <- qnorm(a)
    ci <- array(NA, dim = c(length(parm), 2L), dimnames = list(parm, 
        pct))
    ses <- object$se
    ci[] <- cf[parm] + ses %o% fac
    ci
 }