Sampling posterior predictive from shifted Lognormal

Questions
6

I am not convinced that it is necessary to have a new class for your use case - shifting does not change the logp so you might consider doing below in your model instead:

with pm.Model() as m:
    ...
    x = pm.Lognormal(...)
    x_shift = x + shift