Function qc_metrics

src/memory/statistics/qc.rs:448–472

Convenient function to calculate standard QC metrics with mitochondrial genes. This is a simplified interface for the most common QC analysis in single-cell studies. It automatically identifies mitochondrial genes and computes standard QC metrics. ## What it does 1. **Identifies mitochondrial genes**: Finds genes starting with "MT-" or "mt-" 2. **Adds mito annotation**: Creates boolean column i

(adata: &IMAnnData)

/// If your data uses different naming conventions, use `calculate_qc_metrics` directly
/// with custom gene sets.
pub fn qc_metrics(adata: &IMAnnData) -> anyhow::Result<()> {
    let var_names = adata.var_names();
    let mito_mask: Vec<bool> = var_names
        .iter()
        .map(|name| name.starts_with("MT-") || name.starts_with("mt-"))
        .collect();

    let mut var_df = adata.var().get_data();
    var_df.with_column(Column::new("mito".into(), mito_mask))?;
    adata.var().set_data(var_df)?;

    calculate_qc_metrics(
        adata,
        Some("counts"),
        Some("genes"),
        Some(vec!["mito"]),
        Some(vec![50, 100, 200, 500]),
        None,
        false,
        true,
        true,
    )?;

    Ok(())
}