Make a quick first guess where reasonable beta would be
estimate_betas_roughly(Y, model_matrix, offset_matrix, pseudo_count = 1)
a matrix with one column for each coefficient