code: def _ll_nb1(y, X, beta, alph):
      """Negative Binomial regression (type 1 likelihood)"""
      ll = _ll_nbp(y, X, beta, alph, Q=1)
      return ll
  
imports: from urllib.request import urlopen
  import numpy
  from numpy.testing import assert_almost_equal
  from scipy.special import digamma
  from scipy.stats import nbinom
  import pandas
  import patsy
  from statsmodels.base.model import GenericLikelihoodModel
  from statsmodels.base.model import GenericLikelihoodModelResults
label_desc: Use image node_id 2 for calling a global function with example usage: _ll_nb1(y, X, beta, alph) and returns: ll