R/medIntensities.R

In MarioniLab/cydar: Using Mass Cytometry for Differential Abundance Analyses

#' @export
#' @importFrom flowCore markernames
#' @importFrom SummarizedExperiment assayNames assay<- assayNames<-
#' @importFrom S4Vectors metadata
medIntensities <- function(x, markers)
# Computes the median intensity for each sample, in each hypersphere,
# for each marker. The idea is to fit this to a linear model.
#
# written by Aaron Lun
# created 2 December 2016
{
    .check_cell_data(x)
    all.markers <- markernames(x, mode="all")
    selected <- all.markers[.chosen_markers(markers, all.markers)]
    
    all.leftovers <- markernames(x, mode="unused")
    if (!all(selected %in% all.leftovers)) {
        stop("markers used for counting cannot be used for computing median intensities")
    }
    keep <- all.leftovers %in% selected

    unused <- .raw_unusedIntensities(x)[keep,,drop=FALSE]
    out <- median_int_by_sample(t(unused), .raw_cellAssignments(x), .raw_sample_id(x), ncol(x))
    
    chosen.leftovers <- all.leftovers[keep]
    old.names <- assayNames(x)
    for (j in seq_along(chosen.leftovers)) { 
        assay(x, length(old.names)+j, withDimnames=FALSE) <- out[[j]]
    }
    assayNames(x) <- c(old.names, sprintf("med.%s", chosen.leftovers))

    x
}