code: def test_sample_pymc():
      distribs_pymc = [Bernoulli("B", 0.2), Binomial("N", 5, 0.4)]
      size = 3
      pymc = import_module("pymc")
      if not pymc:
          skip("PyMC is not installed. Abort tests for _sample_pymc.")
      else:
          for X in distribs_pymc:
              samps = sample(X, size=size, library="pymc")
              for sam in samps:
                  assert sam in X.pspace.domain.set
          raises(
              NotImplementedError,
              lambda: (sample(Die("D"), library="pymc")),
          )
  
file_path: sympy/sympy/stats/sampling/tests/test_sample_finite_rv.py
imports: from sympy.core.numbers import Rational
  from sympy.core.singleton import S
  from sympy.external import import_module
  from sympy.stats import Binomial, sample, Die, FiniteRV, DiscreteUniform, Bernoulli, BetaBinomial, Hypergeometric, Rademacher
  from sympy.testing.pytest import skip, raises
label_desc: Use image node_id 4 for calling a global function with example usage: test_sample_pymc() without return types