Thanks for the reply!
Running main.py here should perform the calculation, it only needs pymc3 installed.
The likelihood part is interesting, and I haven’t thought about it either. That might be the reason why I had trouble with NUTS in the first place.
Will implement the changes you recommended! Thanks again