def METHOD_NAME():
  # Test statistic calculation
  peak_list = [5]+[0]*5+[50]+[0]*12+[25]+[0]*10+[10]+[0]*12+[5]+[0]*10+[10]+[0]*17
  binned_peaks = score.calculate_binned_counts(peak_list)
  assert binned_peaks==[50,25,10,5,10,5]
  zscore, rotamer_ratio = score.statistic(binned_peaks)
  assert approx_equal(rotamer_ratio, 2.0/3)
  assert approx_equal(zscore, 2.570679)
  new_rotamer_ratio, new_zscore  = score.calc_ratio(peak_list)
  print(zscore)
  print(new_zscore)
  assert approx_equal(new_rotamer_ratio, rotamer_ratio)
  assert approx_equal(new_zscore,zscore)