Function rank_gene_groups

src/memory/processing/diffexp/mod.rs:189–240

(
    adata: &IMAnnData,
    groupby: &str,
    reference: Option<&str>,
    groups: Option<&[&str]>,
    key_added: Option<&str>,
    method: Option<TestMethod>,
    n_genes: Option<usize>,
    corre

/// ```
#[allow(clippy::too_many_arguments)]
pub fn rank_gene_groups(
    adata: &IMAnnData,
    groupby: &str,
    reference: Option<&str>,
    groups: Option<&[&str]>,
    key_added: Option<&str>,
    method: Option<TestMethod>,
    n_genes: Option<usize>,
    correction_method: CorrectionMethod,
    compute_logfoldchanges: Option<bool>,
    pseudocount: Option<f64>,
) -> anyhow::Result<()> {
    let method = method.unwrap_or(TestMethod::TTest(TTestType::Welch));
    let key = key_added.unwrap_or("").to_string();
    let compute_lfc = compute_logfoldchanges.unwrap_or(true);
    let pseudocount = pseudocount.unwrap_or(1.0);
    let n_genes = n_genes.unwrap_or(adata.n_vars());

    let computation = compute_rank_gene_groups(
        adata,
        groupby,
        reference,
        groups,
        method,
        n_genes,
        correction_method,
        compute_lfc,
        pseudocount,
    )?;

    let RankGeneGroupsComputed {
        groups: groups_to_test,
        results,
        reference,
    } = computation;

    store_results(
        adata,
        &key,
        &groups_to_test,
        results.scores,
        results.pvals,
        results.pvals_adj,
        results.logfoldchanges,
        results.gene_names,
        method,
        groupby,
        reference,
    )?;

    Ok(())
}

#[allow(clippy::too_many_arguments)]
pub fn rank_gene_groups_dataframe(