            d['subject'] = subjects
            d['outcome'] = outcomes
            d['preds'] = preds

            for f_id in range(0, num_top_features):
                feat = features_best[f_id]
                curr_beta = betas[:, order[f_id]]
                curr_shap = shap_values[:, order[f_id]]
                d[f"{feat}_beta"] = curr_beta
                d[f"{feat}_shap"] = curr_shap