R/rowvars.r

In perishky/dmrff: Identifying differentially methylated regions efficiently with power and control

row.vars <- function(x, na.rm=F) {
    x <- x - rowMeans(x, na.rm=na.rm)
    if (na.rm)
        n <- rowSums(!is.na(x))
    else
        n <- ncol(x)
    rowSums(x^2,na.rm=na.rm)/(n-1)
    }
row.sds <- function(...) {
    sqrt(row.vars(...))
}