Fit response values using a binomial/logit model
fit_logistic(y, X)
y
numeric vector
X
numeric matrix (covariates)