.. DO NOT EDIT. .. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY. .. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE: .. "auto_examples/plot_running_median.py" .. LINE NUMBERS ARE GIVEN BELOW. .. only:: html .. note:: :class: sphx-glr-download-link-note :ref:`Go to the end ` to download the full example code. .. rst-class:: sphx-glr-example-title .. _sphx_glr_auto_examples_plot_running_median.py: Running median example ======================== Plot running median on a data set .. GENERATED FROM PYTHON SOURCE LINES 7-33 .. image-sg:: /auto_examples/images/sphx_glr_plot_running_median_001.png :alt: plot running median :srcset: /auto_examples/images/sphx_glr_plot_running_median_001.png :class: sphx-glr-single-img .. rst-class:: sphx-glr-script-out .. code-block:: none (10.0, 50.0) | .. code-block:: Python from pylab import * from sequana.running_median import RunningMedian N = 1000 X = linspace(0, N - 1, N) # Create some interesting data with SNP and longer over # represented section. data = 20 + randn(N) + sin(X * 2 * pi / 1000.0 * 5) data[300:350] += 10 data[500:505] += 100 data[700] = 1000 plot(X, data, "k", label="data") rm = RunningMedian(data, 101) plot(X, rm.run(), "r", label="median W=201") from sequana.stats import moving_average as ma plot(X[100:-100], ma(data, 201), "g", label="mean W=201") grid() legend() ylim([10, 50]) .. rst-class:: sphx-glr-timing **Total running time of the script:** (0 minutes 0.049 seconds) .. _sphx_glr_download_auto_examples_plot_running_median.py: .. only:: html .. container:: sphx-glr-footer sphx-glr-footer-example .. container:: sphx-glr-download sphx-glr-download-jupyter :download:`Download Jupyter notebook: plot_running_median.ipynb ` .. container:: sphx-glr-download sphx-glr-download-python :download:`Download Python source code: plot_running_median.py ` .. container:: sphx-glr-download sphx-glr-download-zip :download:`Download zipped: plot_running_median.zip ` .. only:: html .. rst-class:: sphx-glr-signature `Gallery generated by Sphinx-Gallery `_