Preprocess an expression matrix for CLAMP

Source: R/utils.R

Filters genes by mean expression and variance, returning the filtered matrix and per-gene statistics.

preprocessCLAMP(Y, mean_cutoff = 0, var_cutoff = 0)

Arguments

Y

Numeric matrix of gene expression (rows = genes, cols = samples)

mean_cutoff

Numeric. Minimum row-mean required to keep a gene (default 0).

var_cutoff

Numeric. Minimum row-variance required to keep a gene (default 0).

Value

A list with components:

  • Y_filtered: filtered matrix (genes x samples)

  • rowStats: data.frame with columns mean and variance for each kept gene

  • kept_rows: integer vector of the original row indices that were kept

Examples

# construct a small example matrix
mat <- matrix(
    c(
        1, 5, 10,
        2, 6, 11,
        3, 7, 12,
        4, 8, 13
    ),
    nrow = 4, byrow = FALSE,
    dimnames = list(paste0("gene", seq_len(4)), paste0("sample", seq_len(3)))
)

# keep genes with mean >= 6 and variance >= 2
res <- preprocessCLAMP(mat, mean_cutoff = 6, var_cutoff = 2)