My use case is the insurance industry modeling claims/loss. PyTorch has a loss function for Poisson likelihood, but not Gamma or Poisson-Gamma distributions or Tweedie (where you can specify p as a hyperparameter). Does anyone know of any publicly available code out there implementing this?