code: def test_issue_6682():
      # Reference value from R:
      # options(digits=16)
      # print(pnbinom(250, 50, 32/63, lower.tail=FALSE))
      assert_allclose(
          nbinom.sf(250, 50, 32.0 / 63.0), 1.460458510976452e-35
      )
  
imports: import pytest
  from scipy.stats import betabinom, hypergeom, nhypergeom, bernoulli, boltzmann, skellam, zipf, zipfian, binom, nbinom, nchypergeom_fisher, nchypergeom_wallenius, randint
  import numpy
  from numpy.testing import assert_almost_equal, assert_equal, assert_allclose
  from scipy.special import binom
  import pytest
  from scipy.optimize import root_scalar
  from scipy.integrate import quad
label_desc: Use image node_id 18 for calling a global function with example usage: test_issue_6682() without return types