code: def test_percentile_of_score(score, axis):
      shape = 10, 20, 30
      np.random.seed(0)
      x = np.random.rand(*shape)
      p = _bootstrap._percentile_of_score(x, score, axis=-1)
  
      def vectorized_pos(a, score, axis):
          return np.apply_along_axis(
              stats.percentileofscore, axis, a, score
          )
  
      p2 = vectorized_pos(x, score, axis=-1) / 100
  
      assert_allclose(p, p2, 1e-15)
  
imports: import numpy
  import pytest
  from scipy.stats import bootstrap
  from numpy.testing import assert_allclose, assert_equal
  from scipy import stats
  from scipy.stats import _bootstrap
  from scipy._lib._util import rng_integers
label_desc: Use image node_id 13 for calling a global function with example usage: test_percentile_of_score(score, axis) and returns: np