Often used as a QC parameter for genes in scRNA-seq data, usually once the LQ cells have been removed. The matrix is converted to TPM space from logTPM (more specifically, from log2(TPM/10 + 1)), gene averages are calculated and finally the values are converted back to log2(TPM + 1), without dividing by 10 since the data is now simulated bulk.
aggr_gene_expr(m, na.rm = T, isBulk = FALSE)matrix of expression values to test