R/calc_specificity.R

In bschilder/scKirby: Automated ingestion and conversion of various single-cell data formats.

calc_specificity <- function(X) {
  #### From EWCE::generate_celltype_data() #####
  normalised_meanExp <- Matrix::t(Matrix::t(X) * (1 / colSums(X)))
  spec <- normalised_meanExp / (apply(normalised_meanExp, 1, sum) + 1e-12)
  return(spec)
}