code: def test_nhypergeom_rvs_shape():
      # Check that when given a size with more dimensions than the
      # dimensions of the broadcast parameters, rvs returns an array
      # with the correct shape.
      x = nhypergeom.rvs(22, [7, 8, 9], [[12], [13]], size=(5, 1, 2, 3))
      assert x.shape == (5, 1, 2, 3)
  
imports: import pytest
  import itertools
  from scipy.stats import betabinom, betanbinom, 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, suppress_warnings
  from scipy.special import binom
  from scipy.optimize import root_scalar
  from scipy.integrate import quad
label_desc: Use image node_id 7 for calling a global function with example usage: test_nhypergeom_rvs_shape() without return types