This function uses mean(values ^ 2) - mean(values) ^ 2
gene_vars(dat, means = NULL, rescale = FALSE)a matrix or dg[CT]Matrix of expression values with genes as rows and samples as columns.
(optional) a vector of means for each gene/row. If NULL (default), will be computed using gene_means().
a logical indicating whether or not to rescale using rescale_samples(). Default is FALSE.
a numeric vector of variance values for each gene
Results will differ from the var() function.