Difficulties in using .dist and pm.Potential for inference

Ah yeah I had a typo, you need to give it a name, like: log_norm = pm.Lognormal('log_norm', mu=μlog, sigma=σlog)

Then you should be using the variable log_norm in the scan, not pm.Lognormal.