Maybe variational inference is a possibility? You can try searching for black box stochastic gradient variational inference, easy to use and scales well.