code: def test_oim(self):
      # estimate with default oim, cov_type nonrobust
      res1 = self.res1.model.fit()
      res2 = self.res2
  
      k_mean = 4
      # R betareg uses numerical derivatives from bfgs, has lower precision
      p, se, zv, pv = res2.table_mean_oim.T
      assert_allclose(res1.params[:k_mean], p, rtol=1e-6)
      assert_allclose(res1.bse[:k_mean], se, rtol=1e-5)
      assert_allclose(res1.tvalues[:k_mean], zv, rtol=1e-5)
      assert_allclose(res1.pvalues[:k_mean], pv, atol=1e-6, rtol=1e-5)
  
      p, se, zv, pv = res2.table_precision_oim.T
      assert_allclose(res1.params[k_mean:], p, rtol=1e-6)
      assert_allclose(res1.bse[k_mean:], se, rtol=1e-3)
      assert_allclose(res1.tvalues[k_mean:], zv, rtol=1e-3)
      assert_allclose(res1.pvalues[k_mean:], pv, rtol=1e-2)
  
imports: import io
  import os
  import pytest
  import numpy
  from numpy.testing import assert_allclose, assert_equal
  import pandas
  import patsy
  from statsmodels.api import families
  from statsmodels.tools.sm_exceptions import ValueWarning
  from statsmodels.othermod.betareg import BetaModel
  from .results import results_betareg
label_desc: Use image node_id 4 for calling the TestBetaMeth obj's underlying member method code with example usage: obj.test_oim() without return types