pub fn select_model(
y: &Array2<f64>,
designlist: &[Array2<f64>],
criterion: SelectCriterion,
df_prior: f64,
s2_prior: Option<f64>,
s2_true: Option<&[f64]>,
) -> SelectModelResultExpand 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").