code: def test_score_test(self):
      res1 = self.res1
      resr = self.resr
      params_restr = np.concatenate([resr.params, [0]])
      r_matrix = np.zeros((1, len(params_restr)))
      r_matrix[0, -1] = 1
      exog_prec_extra = res1.model.exog_precision[:, 1:]
  
      from statsmodels.base._parameter_inference import score_test
  
      sc1 = score_test(
          res1, params_constrained=params_restr, k_constraints=1
      )
      sc2 = score_test(resr, exog_extra=(None, exog_prec_extra))
      assert_allclose(sc2[:2], sc1[:2])
  
      sc1_hc = score_test(
          res1,
          params_constrained=params_restr,
          k_constraints=1,
          r_matrix=r_matrix,
          cov_type="HC0",
      )
      sc2_hc = score_test(
          resr, exog_extra=(None, exog_prec_extra), cov_type="HC0"
      )
      assert_allclose(sc2_hc[:2], sc1_hc[: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 2 for calling the TestBetaIncome obj's underlying member method code with example usage: obj.test_score_test() without return types