How to do MCMC sampling with Multiple likelihood functions?

The model code is valid and multiple likelihood will be accounted for - I dont see np.vstack in your code, where and how do you use it? The error seems to raise from there.