code: def test_confidence_prob_func(data, pars):
      """Test conf_interval with alternate prob_func."""
      minimizer = lmfit.Minimizer(residual, pars, fcn_args=data)
      out = minimizer.minimize(method="leastsq")
      called = 0
  
      def my_f_compare(best_fit, new_fit):
          nonlocal called
          called += 1
          nfree = best_fit.nfree
          nfix = best_fit.nvarys - new_fit.nvarys
          dchi = new_fit.chisqr / best_fit.chisqr - 1.0
          return f.cdf(dchi * nfree / nfix, nfix, nfree)
  
      lmfit.conf_interval(
          minimizer, out, sigmas=[1], prob_func=my_f_compare
      )
      assert called > 10
  
imports: import copy
  import numpy
  from numpy.testing import assert_allclose
  import pytest
  from scipy.stats import f
  import lmfit
label_desc: Use image node_id 15 for calling a global function with example usage: test_confidence_prob_func(data, pars) and returns: f