Function fit_logistic 

pub fn fit_logistic(
    x: &[f64],
    y: &[f64],
    n: usize,
    p: usize,
    fit_intercept: bool,
) -> Result<GlmModel, LearningError>

Convenience: logistic regression (Bernoulli y ∈ {0,1}).