R/colVars.R

In RBesT: R Bayesian Evidence Synthesis Tools

Defines functions colVars

Documented in colVars

#' Fast column-wise calculation of unbiased variances
#' @keywords internal
colVars <- function(a) {
    n <- dim(a)[[1]]
    c <- dim(a)[[2]]
    return(.colMeans(((a - matrix(.colMeans(a, n, c), nrow = n, ncol = c, byrow = TRUE)) ^ 2), n, c) * n / (n - 1))
}

Try the RBesT package in your browser