Skip to main content

classify_tests_f

Function classify_tests_f 

Source
pub fn classify_tests_f(
    tstat: &Array2<f64>,
    cor_matrix: Option<&Array2<f64>>,
    df: &[f64],
    p_value: f64,
) -> Array2<i8>
Expand description

Classify each contrast per gene as 1 (up), -1 (down) or 0 via the nested step-down F-test. df is the denominator degrees of freedom, recycled when length 1 (f64::INFINITY selects the chi-squared limit).

Genes with any NaN t-statistic are returned all-zero: limma sets them NA, and this i8 port maps NA → 0 to match crate::decidetests.