R/sumCI.R

In ybkamaleri/nbReg: Resultattjeneste for Barnediabetesregisteret (BDR)

##' Bregne Konfidens Interval
##'
##' Funksjonen brukes til å bregne Konfidens Interval (CI)
##'
##' @param data Datasettet
##' @param maalvar Utvalgte variabel
##' @param gpvars Gruppe
##' @param conf.Int Konfidence nivå f.eks .95
##'
##' @import dplyr
##' @import tidyr
##' 
##' @export
##' 

sumCI <- function(data, maalvar, gpvars, conf.Int=.95) {

    summ = data %>%
        filter_(.dots = lazyeval::interp(~!is.na(var), var = as.name(maalvar))) %>%
        group_by_(.dots = gpvars) %>%
        summarise_(.dots = list(
                       N = lazyeval::interp(~n()),
                       Mean =  lazyeval::interp(~mean(var, na.rm = TRUE), var = as.name(maalvar)),
                       SD =  lazyeval::interp(~sd(var, na.rm = T), var = as.name(maalvar))))
    ## For sample
    ## ciMult = qt(conf.Int/2 + .5, summ$N-1) #konvertere f.eks CI=95% blir .975

    ## For total populasjon
    ciMult = qt(conf.Int/2 + .5, summ$N) #konvertere f.eks CI=95% blir .975

    summ$SE = summ$SD/sqrt(summ$N)
    summ$CI = ciMult * summ$SE

    summ.na.zero <- data %>%
       expand_(gpvars) %>%
       left_join(summ)

    summ.na.zero[is.na(summ.na.zero)] <- 0

    return(invisible(summ.na.zero))
    
}