code: def test_pearson(self):
      p_x = np.array([0.0, 0.9, 1.8, 2.6, 3.3, 4.4, 5.2, 6.1, 6.5, 7.4])
      p_y = np.array([5.9, 5.4, 4.4, 4.6, 3.5, 3.7, 2.8, 2.8, 2.4, 1.5])
      p_sx = np.array(
          [0.03, 0.03, 0.04, 0.035, 0.07, 0.11, 0.13, 0.22, 0.74, 1.0]
      )
      p_sy = np.array(
          [1.0, 0.74, 0.5, 0.35, 0.22, 0.22, 0.12, 0.12, 0.1, 0.04]
      )
  
      p_dat = RealData(p_x, p_y, sx=p_sx, sy=p_sy)
  
      # Reverse the data to test invariance of results
      pr_dat = RealData(p_y, p_x, sx=p_sy, sy=p_sx)
  
      p_mod = Model(self.pearson_fcn, meta=dict(name="Uni-linear Fit"))
  
      p_odr = ODR(p_dat, p_mod, beta0=[1.0, 1.0])
      pr_odr = ODR(pr_dat, p_mod, beta0=[1.0, 1.0])
  
      out = p_odr.run()
      assert_array_almost_equal(
          out.beta,
          np.array([5.4767400299231674, -0.4796082367610305]),
      )
      assert_array_almost_equal(
          out.sd_beta,
          np.array([0.3590121690702467, 0.0706291186037444]),
      )
      assert_array_almost_equal(
          out.cov_beta,
          np.array(
              [
                  [0.0854275622946333, -0.0161807025443155],
                  [-0.0161807025443155, 0.003306337993922],
              ]
          ),
      )
  
      rout = pr_odr.run()
      assert_array_almost_equal(
          rout.beta,
          np.array([11.4192022410781231, -2.0850374506165474]),
      )
      assert_array_almost_equal(
          rout.sd_beta,
          np.array([0.9820231665657161, 0.3070515616198911]),
      )
      assert_array_almost_equal(
          rout.cov_beta,
          np.array(
              [
                  [0.6391799462548782, -0.1955657291119177],
                  [-0.1955657291119177, 0.0624888159223392],
              ]
          ),
      )
  
imports: import tempfile
  import shutil
  import os
  import numpy
  from numpy import pi
  from numpy.testing import assert_array_almost_equal, assert_equal, assert_warns
  import pytest
  from pytest import raises
  from scipy.odr import Data, Model, ODR, RealData, OdrStop, OdrWarning, multilinear, exponential, unilinear, quadratic, polynomial
label_desc: Use image node_id 13 for calling the TestODR obj's underlying member method code with example usage: obj.test_pearson() without return types