code: def test_log_binomial_stirling(tol):
      k = torch.arange(200.0)
      n_minus_k = k.unsqueeze(-1)
      n = k + n_minus_k
  
      # Test binomial coefficient choose(n, k).
      expected = (
          (n + 1).lgamma() - (k + 1).lgamma() - (n_minus_k + 1).lgamma()
      )
      actual = log_binomial(n, k, tol=tol)
  
      assert (actual - expected).abs().max() < tol
  
file_path: pyro/tests/ops/test_special.py
imports: import pytest
  import torch
  from scipy.special import iv
  from torch import tensor
  from torch.autograd import grad
  from pyro.ops.special import get_quad_rule, log_beta, log_binomial, log_I1, safe_log
  from tests.common import assert_equal
label_desc: Use image node_id 3 for calling a global function with example usage: test_log_binomial_stirling(tol) without return types