code: def _ll_geom(y, X, beta):
      """Geometric regression"""
      ll = _ll_nbp(y, X, beta, alph=1, Q=0)
      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 4 for calling a global function with example usage: _ll_geom(y, X, beta) and returns: ll