##' 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))
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.