code: def test_score():
      uniq, load, corr, par = _toy()
      fa = Factor(n_factor=2, corr=corr)
  
      def f(par):
          return fa.loglike(par)
  
      par2 = np.r_[
          0.1, 0.2, 0.3, 0.4, 0.3, 0.1, 0.2, -0.2, 0, 0.8, 0.5, 0
      ]
  
      for pt in (par, par2):
          g1 = approx_fprime(pt, f, 1e-8)
          g2 = fa.score(pt)
          assert_allclose(g1, g2, atol=1e-3)
  
file_path: statsmodels/statsmodels/multivariate/tests/test_ml_factor.py
imports: import numpy
  from statsmodels.multivariate.factor import Factor
  from numpy.testing import assert_allclose, assert_equal
  from scipy.optimize import approx_fprime
  import warnings
label_desc: Use image node_id 3 for calling a global function with example usage: test_score() and returns: fa