code: def test_betabinom_bernoulli():
      # test limiting case that betabinom(1, a, b) = bernoulli(a / (a + b))
      a = 2.3
      b = 0.63
      k = np.arange(2)
      p = betabinom(1, a, b).pmf(k)
      expected = bernoulli(a / (a + b)).pmf(k)
      assert_almost_equal(p, expected)
  
imports: import pytest
  from scipy.stats import betabinom, hypergeom, nhypergeom, bernoulli, boltzmann, skellam, zipf, zipfian, binom, nbinom, nchypergeom_fisher, nchypergeom_wallenius, randint
  import numpy
  from numpy.testing import assert_almost_equal, assert_equal, assert_allclose
  from scipy.special import binom
  import pytest
  from scipy.optimize import root_scalar
  from scipy.integrate import quad
label_desc: Use image node_id 9 for calling a global function with example usage: test_betabinom_bernoulli() without return types