I need it as a likelihood distribution. I did try something like you suggest:
with pm.Model() as m:
...
shift = pm......
x = pm.Lognormal('y_obs' ... observed = data + shift)
But this leads to:
Cannot resolve inputs for ['y_obs']
When sampling the PP.