recover_marginals doesn't work with implicit model from modelcontext #610
Description
with marginal_model:
recover_marginals(idata) # missing argument idata (because it assumes idata is model)