Skip to main content

select_model

Function select_model 

Source
pub fn select_model(
    y: &Array2<f64>,
    designlist: &[Array2<f64>],
    criterion: SelectCriterion,
    df_prior: f64,
    s2_prior: Option<f64>,
    s2_true: Option<&[f64]>,
) -> SelectModelResult
Expand description

Score designlist against expression matrix y (genes × arrays).

df_prior / s2_prior supply a prior for the residual variance (AIC/BIC only; s2_prior is required when df_prior > 0). s2_true (length 1 or ngenes) is the known variance required by Mallows’ Cp. Panics on any NaN in y, mirroring limma’s stop("NAs not allowed").