Run the Collapsed Gibbs sampler for the keyATM covariates (Polya-Gamma)
keyATM_fit_covPG(model, resume = FALSE)
model
A initialized model
resume
resume or not