code: def _(dist: NegativeBinomialDistribution):
      return pymc.NegativeBinomial(
          "X",
          mu=float((dist.p * dist.r) / (1 - dist.p)),
          alpha=float(dist.r),
      )
  
imports: from functools import singledispatch
  from sympy.external import import_module
  from sympy.stats.crv_types import BetaDistribution, CauchyDistribution, ChiSquaredDistribution, ExponentialDistribution, GammaDistribution, LogNormalDistribution, NormalDistribution, ParetoDistribution, UniformDistribution, GaussianInverseDistribution
  from sympy.stats.drv_types import PoissonDistribution, GeometricDistribution, NegativeBinomialDistribution
  from sympy.stats.frv_types import BinomialDistribution, BernoulliDistribution
label_desc: Use image node_id 13 for calling a global function with example usage: _(dist) and returns: pymc