code: def test_very_high_moments(self):
      running_moments = sample_stats.RunningCentralMoments.from_shape(
          shape=(), moment=np.arange(15) + 1
      )
      for sample in range(5):
          running_moments = running_moments.update(sample)
      moments = self.evaluate(running_moments.moments())
      self.assertAllClose(
          stats.moment(np.arange(5), moment=np.arange(15) + 1),
          moments,
          atol=5e-5,
          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 3 for calling the RunningCentralMomentsTest obj's underlying member method code with example usage: obj.test_very_high_moments() without return types