Run the Collapsed Gibbs sampler for the keyATM covariates (Dir-Multi)
keyATM_fit_cov(model, resume = FALSE)
model
A initialized model
resume
resume or not