Implementig new distributions

How can I implement the Generalized gamma distribution and the Gompertz distribution using TransformedDistribution (if it’s possible) ?

Probably not possible, would imply a bijective transform from 2 to 3 param family

Not sure that a transformation into this exists.

You can just write log_prob functions manually, for example adapting scipy.stats code