Run the Collapsed Gibbs sampler for weighted LDA with covariates
keyATM_fit_LDAcov(model, resume = FALSE)
model
A initialized model
resume
resume or not