code: def test_confint_poisson_2indep():
      # test other confint methods
  
      count1, exposure1, count2, exposure2 = 60, 51477.5, 30, 54308.7
  
      ci = confint_poisson_2indep(
          count1,
          exposure1,
          count2,
          exposure2,
          method="mover",
          compare="ratio",
          alpha=0.1,
          method_mover="jeff",
      )
      ci1 = (1.4667, 3.0608)  # LTW 2014
      assert_allclose(ci, ci1, atol=0.05)
      # moverci(60, 51477.5, 30, 54308.7, type = "jeff", distrib="poi",
      # contrast="RR", cc=0, level=0.9)
      #   Lower Estimate    Upper
      # est 1.466768 2.104135 3.058634
      # ratesci uses different center for mover-jeffreys
      ci1 = (1.466768, 3.058634)
      assert_allclose(ci, ci1, rtol=0.001)
  
      ci = confint_poisson_2indep(
          count1,
          exposure1,
          count2,
          exposure2,
          method="mover",
          compare="ratio",
          alpha=0.1,
          method_mover="score",
      )
      ci1 = (1.4611, 3.0424)  # LTW 2014
      assert_allclose(ci, ci1, atol=0.05)
  
      ci = confint_poisson_2indep(
          count1,
          exposure1,
          count2,
          exposure2,
          method="waldcc",
          compare="ratio",
          alpha=0.1,
      )
      ci1 = (1.4523, 3.0154)  # LTW 2014
      assert_allclose(ci, ci1, atol=0.0005)
  
      ci = confint_poisson_2indep(
          count1,
          exposure1,
          count2,
          exposure2,
          method="score",
          compare="ratio",
          alpha=0.05,
      )
      ci1 = (1.365962, 3.259306)
      assert_allclose(ci, ci1, atol=5e-6)
  
      # with diff
      # use larger rates for diff
      exposure1 /= 1000
      exposure2 /= 1000
      ci = confint_poisson_2indep(
          count1,
          exposure1,
          count2,
          exposure2,
          method="mover",
          compare="diff",
          alpha=0.05,
          method_mover="jeff",
      )
      # moverci(60, 51477.5/1000, 30, 54308.7/1000, type = "jeff", distrib="poi",
      # contrast="RD", cc=0, level=0.95)
      ci1 = (0.2629322, 0.9786493)
      assert_allclose(ci, ci1, atol=0.005)
      ci = confint_poisson_2indep(
          count1,
          exposure1,
          count2,
          exposure2,
          method="score",
          compare="diff",
          alpha=0.05,
      )
      # scoreci(60, 51477.5/1000, 30, 54308.7/1000, distrib="poi", contrast="RD",
      # skew=FALSE, cc=0, level=0.95)
      ci1 = (0.265796, 0.989192)
      assert_allclose(ci, ci1, atol=5e-6)
  
      # example in Li et al 2011, scaled
      ci = confint_poisson_2indep(
          count2,
          exposure2,
          count1,
          exposure1,
          method="mover",
          compare="diff",
          alpha=0.1,
          method_mover="jeff",
      )
      # moverci(30, 54308.7/1000, 60, 51477.5/1000, type = "jeff", distrib="poi",
      # contrast="RD", cc=0, level=0.90)
      ci1 = (-0.9183272231752, -0.3188611692202)
      assert_allclose(ci, ci1, atol=0.005)
      ci1 = (-0.9195, -0.3193)  # from Li et al 2011
      assert_allclose(ci, ci1, atol=0.005)
  
      ci = confint_poisson_2indep(
          count2,
          exposure2,
          count1,
          exposure1,
          method="mover",
          compare="diff",
          alpha=0.1,
          method_mover="jeff",
      )
      # scoreci(30, 54308.7/1000, 60, 51477.5/1000, distrib="poi", contrast="RD",
      # skew=FALSE, cc=0, level=0.90)
      ci1 = (-0.9232, -0.3188)  # from Li et al 2011
      assert_allclose(ci, ci1, atol=0.006)
  
imports: import pytest
  import warnings
  import numpy
  from numpy import arange
  from numpy.testing import assert_allclose, assert_equal
  from scipy import stats
  import statsmodels.stats.rates
  from statsmodels.stats.rates import confint_poisson, tolerance_int_poisson, confint_quantile_poisson, etest_poisson_2indep, confint_poisson_2indep, nonequivalence_poisson_2indep, power_poisson_ratio_2indep, power_equivalence_poisson_2indep, power_poisson_diff_2indep, power_equivalence_neginb_2indep, power_negbin_ratio_2indep, method_names_poisson_1samp, method_names_poisson_2indep
label_desc: Use image node_id 10 for calling a global function with example usage: test_confint_poisson_2indep() without return types