rmeove accidental peak-utils edit

anopheles-genomic-surveillance · Jan 22, 2024 · e413751 · e413751
1 parent 5f16823
commit e413751
Showing 1 changed file with 0 additions and 10 deletions.
diff --git a/workflow/notebooks/peak-utils.ipynb b/workflow/notebooks/peak-utils.ipynb
@@ -326,16 +326,6 @@
     "    fit_gcenter = peak_result.params['center'].value\n",
     "\n",
     "    peak_in_scan_interval = ((gcenter - scan_interval) < fit_gcenter < (gcenter + scan_interval))\n",
-    "\n",
-    "    # Signals are defined by fitting an exponential peak to the data moving along the genome in windows. \n",
-    "    # this is done both to flat 'null' lines, and H12 GWSS scans\n",
-    "    # the different between the AIC of the null and H12 model, provides a summary statistic (delta_i) of how strong the signal is\n",
-    "    # Overlapping signals are de-duplicated, with the highest delta_i kept\n",
-    "    # the focus is the decay parameter\n",
-    "    \n",
-    "    focus_right = center + .25 * (2**(-skew_param) * decay_param) \n",
-    "    focus_left = center - .25 * (2**skew_param * decay_param) \n",
-    "    \n",
     "    \n",
     "    fit_params = peak_result.params\n",
     "    fit_skew = fit_params['skew'].value\n",