code: def a_st70_upp(stat, nobs):
      nobsinv = 1.0 / nobs
      stat_modified = stat - 0.7 * nobsinv + 0.9 * nobsinv**2
      stat_modified *= 1 + 1.23 * nobsinv
      pval = 1.273 * np.exp(-2 * stat_modified / 2.0 * np.pi**2)
      digits = np.sum(stat > np.array([0.11, 0.11, 0.452]))
      # repeat low to get {0,2,3}
      return stat_modified, pval, digits
  
imports: from statsmodels.compat.python import lmap
  import numpy
  from scipy.stats import distributions
  from statsmodels.tools.decorators import cache_readonly
  from scipy.special import kolmogorov
label_desc: Use image node_id 8 for calling a global function with example usage: a_st70_upp(stat, nobs) and returns: stat_modified,  pval,  digits