Skip to contents

Estimate Signal Spread for a Given Population

Source: R/EstimateSpread.R
EstimateSpread.Rd

This function estimates the signal spread (variance components) for a specified population using the USERM result object. It computes weighted intercept and slope variances across detectors and fluorochromes based on the system matrix and signal intensities.

Usage

EstimateSpread(Userm, population_id)

Arguments

Userm

A list containing USERM analysis results. Must include:

  • Intensity_mtx: A matrix of signal intensities (fluorochrome x population).

  • detectors: A character vector of detector names.

  • fluors: A character vector of fluorochrome names.

  • A: A system matrix mapping detectors to fluorochromes.

  • Res: A named list of result objects for each fluorochrome, each containing interceptMtx and slopMtx.

population_id

A character string specifying the population ID to estimate spread for. Must match a column name in Userm$Intensity_mtx.

Value

A list containing:

  • population_id: The input population ID.

  • A: The system matrix used.

  • fluors: Fluorochrome names.

  • detectors: Detector names.

  • intensity_matrix: Signal intensity values for the specified population.

  • intercept_sigma2_col: Estimated intercept variance per channel.

  • slop_sigma2_col: Estimated slope variance per channel.

Examples

if (FALSE) { # \dontrun{
  EstimateSpread(Userm, "Population_1")
} # }