Something to keep in mind is that ADVI with the default option (i.e., mean-field approximation) will underestimated the uncertainty of the parameters (much narrower posterior).
You can have a look at this code of using different sampling on a mixed (hierarchical) model: https://github.com/junpenglao/GLMM-in-Python/blob/master/pymc3_different_sampling.py
2 Likes