code: def test_random_higher_rank_samples(self):
      rng = test_util.test_np_rng()
      x_orig = rng.rand(100, 10)
      x = tf.convert_to_tensor(x_orig, dtype=tf.float32)
      running_moments = sample_stats.RunningCentralMoments.from_shape(
          shape=(10,), moment=np.arange(5) + 1
      )
      running_moments = consume(running_moments, x)
      moments = self.evaluate(running_moments.moments())
      self.assertAllClose(
          stats.moment(x_orig, moment=[1, 2, 3, 4, 5]),
          moments,
          rtol=1e-6,
      )
  
imports: from absl.testing import parameterized
  import numpy
  import scipy.stats
  import tensorflow.compat.v2
  from tensorflow_probability.python.experimental.stats import sample_stats
  from tensorflow_probability.python.internal import test_util
  from tensorflow_probability.python.mcmc import diagnostic
label_desc: Use image node_id 5 for calling the RunningCentralMomentsTest obj's underlying member method code with example usage: obj.test_random_higher_rank_samples() without return types