code: def test_mvn_mvt_1(self):
      a, b = self.a, self.b
      df = self.df
      corr_equal = self.corr_equal
      # result from R, mvtnorm with option
      # algorithm = GenzBretz(maxpts = 100000, abseps = 0.000001, releps = 0)
      #     or higher
      probmvt_R = 0.60414  # report, ed error approx. 7.5e-06
      probmvn_R = 0.673970  # reported error approx. 6.4e-07
      assert_almost_equal(
          probmvt_R, mvstdtprob(a, b, corr_equal, df), 4
      )
      assert_almost_equal(
          probmvn_R, mvstdnormcdf(a, b, corr_equal, abseps=1e-5), 4
      )
  
      mvn_high = mvstdnormcdf(
          a, b, corr_equal, abseps=1e-8, maxpts=10000000
      )
      assert_almost_equal(probmvn_R, mvn_high, 5)
  
file_path: statsmodels/statsmodels/sandbox/distributions/tests/test_multivariate.py
imports: import numpy
  from numpy.testing import assert_almost_equal, assert_allclose
  from statsmodels.sandbox.distributions.multivariate import mvstdtprob, mvstdnormcdf
  from statsmodels.sandbox.distributions.mv_normal import MVT, MVNormal
label_desc: Use image node_id 2 for calling the Test_MVN_MVT_prob obj's underlying member method code with example usage: obj.test_mvn_mvt_1() without return types