code: def pvalues(self):
          
          self.compute_statistics()
  
          lml_alts = self.alt_lmls()
          lml_null = self.null_lml()
  
          lrs = -2 * lml_null + 2 * asarray(lml_alts)
  
          from scipy.stats import chi2
          chi2 = chi2(df=1)
  
          return chi2.sf(lrs)
docs: Association p-value for candidate markers.