Function fit_penalized_binomial_multi 

pub fn fit_penalized_binomial_multi(
    inputs: BinomialMultiFitInputs<'_>,
) -> Result<BinomialMultiFitOutputs, EstimationError>

Fit K independent penalized binomial-logit GLMs sharing the design X and penalty S. See the module docs for the optimization problem.