How to create Bayesian data fusion in python with pymc3?

The arguments to pm.sample have to do with helping the MCMC sampler converge to the true posterior. Are your models converging? What are the r_hat values for each variable you estimate? Are there divergences in the sampler? Asking these questions will guide you to an answer about sampler settings.