Repository 'pandas_rolling_window'
hg clone https://toolshed.g2.bx.psu.edu/repos/bgruening/pandas_rolling_window

Changeset 0:a06f7b5c4dc7 (2019-05-20)
Commit message:
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/pandas_rolling_window commit bdbedf42854d16bb00c396045007d4baece0a869
added:
pandas_rolling.xml
test-data/1 (Kopie).bedgraph
test-data/1.bedgraph
test-data/1_boxcar.bedgraph
test-data/1_hamming.bedgraph
test-data/2_boxcar.bedgraph
test-data/2_hamming.bedgraph
b
diff -r 000000000000 -r a06f7b5c4dc7 pandas_rolling.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/pandas_rolling.xml Mon May 20 08:42:54 2019 -0400
[
b'@@ -0,0 +1,175 @@\n+<tool id="pandas_rolling_window" name="Rolling window" version="0.1">\n+    <description>over a dataframe (e.g. for data smoothing)</description>\n+    <requirements>\n+        <requirement type="package" version="1.16.3">numpy</requirement>\n+        <requirement type="package" version="1.2.1">scipy</requirement>\n+        <requirement type="package" version="0.24.2">pandas</requirement>\n+    </requirements>\n+    <command>\n+<![CDATA[\n+\n+    cat \'$pandas_script\' &&\n+    python \'$pandas_script\'\n+\n+]]>\n+    </command>\n+    <configfiles>\n+        <configfile name="pandas_script"><![CDATA[\n+import argparse\n+import sys\n+\n+import pandas as pd\n+\n+kwargs = dict()\n+window_type = \'$smooth_function.smooth_function_opts_selector\'\n+\n+#if $smooth_function.smooth_function_opts_selector == \'gaussian\':\n+kwargs.update({\'std\': $smooth_function.gaussian_std})\n+#elif $smooth_function.smooth_function_opts_selector == \'general_gaussian\':\n+kwargs = ({\'power\': $smooth_function.ggaussian_power, \'width\': $smooth_function.ggaussian_width})\n+#elif $smooth_function.smooth_function_opts_selector == \'kaiser\':\n+kwargs.update({\'beta\': $smooth_function.kaiser_beta})\n+#elif $smooth_function.smooth_function_opts_selector == \'slepian\':\n+kwargs.update({\'width\': $smooth_function.slepian_width})\n+#end if\n+\n+df = pd.read_csv(\'${infile}\', sep=\'\\t\', index_col=None, header=None, dtype={\'strand\': object} )\n+\n+#if $group_column:\n+df[\'aggregate\'] = df.groupby( int($group_column)-1, sort=False )[int($value_column)-1].rolling(${window_len}, win_type=window_type, center=$centering ).${statistics}(**kwargs).reset_index(drop=True)\n+#else:\n+df[\'aggregate\'] = df[int($value_column)-1].rolling(${window_len}, win_type=window_type, center=$centering ).${statistics}(**kwargs).reset_index(drop=True)\n+#end if\n+\n+df.to_csv(\'${outfile}\', index=False, header=False, sep=\'\\t\', na_rep=\'0\', float_format=\'%.2f\')\n+    ]]> </configfile>\n+    </configfiles>\n+    <inputs>\n+        <param name="infile" type="data" format="tabular,bed.interval" label="Select input file in tabular or BED format"/>\n+        <param name="group_column" type="data_column" data_ref="infile" optional="true" label="Optional column to group"\n+            help="For example if you have a chromosome column you probably want to group each chromosome before you apply any function." />\n+\n+        <param name="value_column" type="data_column" data_ref="infile" label="Column with the value of interest"\n+            help="" />\n+\n+        <conditional name="smooth_function">\n+            <param name="smooth_function_opts_selector" type="select" label="Provide a window type"\n+                help="For more information please see https://en.wikipedia.org/wiki/Window_function">\n+                <option value="boxcar" selected="True">Boxcar or Dirichlet, all points are evenly weighted</option>\n+                <option value="triang">triang</option>\n+                <option value="blackman">blackman</option>\n+                <option value="hamming">hamming</option>\n+                <option value="bartlett">bartlett</option>\n+                <option value="parzen">parzen</option>\n+                <option value="bohman">bohman</option>\n+                <option value="blackmanharris">blackmanharris</option>\n+                <option value="nuttall">nuttall</option>\n+                <option value="barthann">barthann</option>\n+                <!--option value="kaiser">kaiser</option>\n+                <option value="gaussian">gaussian</option>\n+                <option value="general_gaussian">general gaussian</option>\n+                <option value="slepian">slepian</option-->\n+            </param>\n+            <when value="boxcar" />\n+            <when value="triang" />\n+            <when value="blackman" />\n+            <when value="hamming" />\n+            <when value="bartlett" />\n+            <when value="parzen" />\n+            <when value="bohman" />\n+            <when value="blackmanharris" />\n+            <when value="nuttal'..b'     <when value="general_gaussian">\n+                <param name="ggaussian_power" type="integer" value="2" min=\'1\' label="power" />\n+                <param name="ggaussian_width" type="integer" value="2" min="1" label="width" />\n+            </when>\n+            <when value="slepian">\n+                <param name="slepian_width" type="integer" value="2" min="1" label="width" />\n+            </when>\n+        </conditional>\n+\n+        <param name="statistics" type="select" label="Provide a statistical function">\n+            <option value="count">Number of non-null observations (count)</option>\n+            <option value="sum">Sum of values (sum)</option>\n+            <option value="mean" selected="true">Mean of values (mean)</option>\n+            <option value="median">Arithmetic median of values (median)</option>\n+            <option value="min">Minimum (min)</option>\n+            <option value="max">max (max)</option>\n+            <option value="std">Bessel-corrected sample standard deviation (std)</option>\n+            <option value="var">Unbiased variance (var)</option>\n+            <option value="skew">Sample skewness (3rd moment)</option>\n+            <option value="kurt">Sample kurtosis (4th moment)</option>\n+            <option value="quantil">Sample quantile (value at %)</option>\n+            <option value="cov">Unbiased covariance (binary) (cov)</option>\n+            <option value="corr">Correlation (corr)</option>\n+        </param>\n+\n+        <param name="centering" type="boolean" truevalue="True" falsevalue="False" label="center smoothed values"\n+            help="By default the labels are set to the right edge of the window. Here you can change that to the center." />\n+        <!-- Options for all formats.-->\n+        <param name="window_len" type="integer" value="3" min="2" label="Window length"/>\n+    </inputs>\n+    <outputs>\n+        <data name="outfile" format_source="infile" />\n+    </outputs>\n+    <tests>\n+        <test>\n+            <param name="infile" value="1.bedgraph"/>\n+            <param name="group_column" value="1"/>\n+            <param name="value_column" value="5"/>\n+            <conditional name="smooth_function">\n+                <param name="smooth_function_opts_selector" value="boxcar"/>\n+            </conditional>\n+            <param name="window_len" value="3"/>\n+            <output name="outfile" value="1_boxcar.bedgraph"/>\n+        </test>\n+        <test>\n+            <!-- None test -->\n+            <param name="infile" value="1.bedgraph"/>\n+            <param name="value_column" value="5"/>\n+            <conditional name="smooth_function">\n+                <param name="smooth_function_opts_selector" value="boxcar"/>\n+            </conditional>\n+            <param name="window_len" value="3"/>\n+            <output name="outfile" value="2_boxcar.bedgraph"/>\n+        </test>\n+        <test>\n+            <param name="infile" value="1.bedgraph"/>\n+            <param name="group_column" value="1"/>\n+            <param name="value_column" value="5"/>\n+            <conditional name="smooth_function">\n+                <param name="smooth_function_opts_selector" value="hamming"/>\n+            </conditional>\n+            <param name="window_len" value="3"/>\n+            <param name="statistics" value="sum"/>\n+            <output name="outfile" value="1_hamming.bedgraph"/>\n+        </test>\n+        <test>\n+            <param name="infile" value="1.bedgraph"/>\n+            <param name="value_column" value="5"/>\n+            <conditional name="smooth_function">\n+                <param name="smooth_function_opts_selector" value="hamming"/>\n+            </conditional>\n+            <param name="window_len" value="3"/>\n+            <param name="statistics" value="sum"/>\n+            <output name="outfile" value="2_hamming.bedgraph"/>\n+        </test>\n+    </tests>\n+    <help>\n+<![CDATA[\n+\n+**What it does**\n+\n+Provides rolling window calculations, e.g. for smoothing values.\n+\n+\n+]]>\n+    </help>\n+</tool>\n'
b
diff -r 000000000000 -r a06f7b5c4dc7 test-data/1 (Kopie).bedgraph
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/1 (Kopie).bedgraph Mon May 20 08:42:54 2019 -0400
b
@@ -0,0 +1,100 @@
+chr10 3000180 3000182 100 5 0 5.0
+chr10 3000185 3000187 100 5 0 5.0
+chr10 3000270 3000272 73 11 4 15.0
+chr10 3000294 3000296 100 14 0 14.0
+chr10 3000353 3000355 81 9 2 11.0
+chr10 3000488 3000490 66 4 2 6.0
+chr10 3000657 3000659 66 4 2 6.0
+chr10 3000709 3000711 80 4 1 5.0
+chr10 3000718 3000720 100 5 0 5.0
+chr10 3000923 3000925 80 4 1 5.0
+chr10 3001027 3001029 88 8 1 9.0
+chr10 3001029 3001031 100 9 0 9.0
+chr10 3001055 3001057 100 8 0 8.0
+chr10 3001188 3001190 75 3 1 4.0
+chr10 3001586 3001588 83 5 1 6.0
+chr10 3001747 3001749 100 7 0 7.0
+chr10 3001898 3001900 37 3 5 8.0
+chr10 3001908 3001910 55 5 4 9.0
+chr10 3001923 3001925 100 5 0 5.0
+chr10 3001928 3001930 100 5 0 5.0
+chr10 3001989 3001991 69 9 4 13.0
+chr10 3001996 3001998 84 11 2 13.0
+chr10 3002030 3002032 92 12 1 13.0
+chr10 3002077 3002079 77 7 2 9.0
+chr10 3002089 3002091 100 9 0 9.0
+chr10 3002093 3002095 44 4 5 9.0
+chr10 3002116 3002118 100 7 0 7.0
+chr10 3002147 3002149 80 8 2 10.0
+chr10 3002159 3002161 50 4 4 8.0
+chr10 3002261 3002263 100 16 0 16.0
+chr10 3002320 3002322 92 12 1 13.0
+chr10 3002514 3002516 100 4 0 4.0
+chr10 3002526 3002528 83 5 1 6.0
+chr10 3002539 3002541 100 8 0 8.0
+chr10 3002560 3002562 88 8 1 9.0
+chr10 3002594 3002596 100 8 0 8.0
+chr10 3002607 3002609 88 8 1 9.0
+chr10 3002628 3002630 100 12 0 12.0
+chr10 3002641 3002643 100 9 0 9.0
+chr10 3002904 3002906 100 9 0 9.0
+chr10 3002932 3002934 83 5 1 6.0
+chr10 3003427 3003429 61 8 5 13.0
+chr10 3003435 3003437 84 11 2 13.0
+chr10 3003465 3003467 90 10 1 11.0
+chr10 3003494 3003496 100 9 0 9.0
+chr10 3003584 3003586 60 6 4 10.0
+chr10 3003924 3003926 100 4 0 4.0
+chr10 3003976 3003978 100 12 0 12.0
+chr10 3004035 3004037 71 10 4 14.0
+chr10 3004640 3004642 90 9 1 10.0
+chr10 3005170 3005172 28 2 5 7.0
+chr10 3005441 3005443 21 3 11 14.0
+chr10 3005994 3005996 50 7 7 14.0
+chr10 3008524 3008526 75 9 3 12.0
+chr10 3008594 3008596 90 9 1 10.0
+chr10 3008737 3008739 50 5 5 10.0
+chr10 3008827 3008829 80 12 3 15.0
+chr10 3008937 3008939 84 11 2 13.0
+chr10 3008985 3008987 66 6 3 9.0
+chr10 3009020 3009022 66 10 5 15.0
+chr10 3009103 3009105 50 8 8 16.0
+chr10 3009118 3009120 83 15 3 18.0
+chr10 3009207 3009209 66 6 3 9.0
+chr10 3009577 3009579 80 4 1 5.0
+chr10 3009627 3009629 100 4 0 4.0
+chr10 3010000 3010002 100 8 0 8.0
+chr10 3010399 3010401 58 10 7 17.0
+chr10 3010447 3010449 94 16 1 17.0
+chr10 3011056 3011058 85 6 1 7.0
+chr10 3011170 3011172 68 11 5 16.0
+chr10 3011280 3011282 66 10 5 15.0
+chr10 3011403 3011405 64 22 12 34.0
+chr10 3011468 3011470 92 12 1 13.0
+chr10 3011556 3011558 100 13 0 13.0
+chr10 3011580 3011582 68 11 5 16.0
+chr10 3011870 3011872 90 18 2 20.0
+chr10 3011999 3012001 81 18 4 22.0
+chr10 3012040 3012042 93 14 1 15.0
+chr10 3012053 3012055 66 10 5 15.0
+chr10 3012065 3012067 80 12 3 15.0
+chr10 3012214 3012216 61 11 7 18.0
+chr10 3012269 3012271 50 7 7 14.0
+chr10 3012319 3012321 33 5 10 15.0
+chr10 3012389 3012391 77 14 4 18.0
+chr10 3012415 3012417 77 14 4 18.0
+chr10 3012560 3012562 94 17 1 18.0
+chr10 3012609 3012611 85 12 2 14.0
+chr10 3012619 3012621 85 12 2 14.0
+chr10 3012741 3012743 56 13 10 23.0
+chr10 3012786 3012788 78 15 4 19.0
+chr10 3012790 3012792 63 12 7 19.0
+chr10 3012792 3012794 57 11 8 19.0
+chr10 3012794 3012796 94 17 1 18.0
+chr10 3012796 3012798 55 10 8 18.0
+chr10 3012798 3012800 66 12 6 18.0
+chr10 3012818 3012820 35 7 13 20.0
+chr10 3012824 3012826 16 4 20 24.0
+chr10 3012858 3012860 23 6 20 26.0
+chr10 3012940 3012942 20 2 8 10.0
+chr10 3012977 3012979 7 1 13 14.0
b
diff -r 000000000000 -r a06f7b5c4dc7 test-data/1.bedgraph
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/1.bedgraph Mon May 20 08:42:54 2019 -0400
b
@@ -0,0 +1,100 @@
+chr10 3000180 3000182 100 5 0 5.0
+chr10 3000185 3000187 100 5 0 5.0
+chr10 3000270 3000272 73 11 4 15.0
+chr10 3000294 3000296 100 14 0 14.0
+chr10 3000353 3000355 81 9 2 11.0
+chr10 3000488 3000490 66 4 2 6.0
+chr10 3000657 3000659 66 4 2 6.0
+chr10 3000709 3000711 80 4 1 5.0
+chr10 3000718 3000720 100 5 0 5.0
+chr10 3000923 3000925 80 4 1 5.0
+chr10 3001027 3001029 88 8 1 9.0
+chr10 3001029 3001031 100 9 0 9.0
+chr10 3001055 3001057 100 8 0 8.0
+chr10 3001188 3001190 75 3 1 4.0
+chr10 3001586 3001588 83 5 1 6.0
+chr10 3001747 3001749 100 7 0 7.0
+chr10 3001898 3001900 37 3 5 8.0
+chr10 3001908 3001910 55 5 4 9.0
+chr10 3001923 3001925 100 5 0 5.0
+chr10 3001928 3001930 100 5 0 5.0
+chr10 3001989 3001991 69 9 4 13.0
+chr10 3001996 3001998 84 11 2 13.0
+chr10 3002030 3002032 92 12 1 13.0
+chr10 3002077 3002079 77 7 2 9.0
+chr10 3002089 3002091 100 9 0 9.0
+chr10 3002093 3002095 44 4 5 9.0
+chr10 3002116 3002118 100 7 0 7.0
+chr10 3002147 3002149 80 8 2 10.0
+chr10 3002159 3002161 50 4 4 8.0
+chr10 3002261 3002263 100 16 0 16.0
+chr10 3002320 3002322 92 12 1 13.0
+chr10 3002514 3002516 100 4 0 4.0
+chr10 3002526 3002528 83 5 1 6.0
+chr10 3002539 3002541 100 8 0 8.0
+chr10 3002560 3002562 88 8 1 9.0
+chr10 3002594 3002596 100 8 0 8.0
+chr10 3002607 3002609 88 8 1 9.0
+chr10 3002628 3002630 100 12 0 12.0
+chr10 3002641 3002643 100 9 0 9.0
+chr10 3002904 3002906 100 9 0 9.0
+chr10 3002932 3002934 83 5 1 6.0
+chr10 3003427 3003429 61 8 5 13.0
+chr10 3003435 3003437 84 11 2 13.0
+chr10 3003465 3003467 90 10 1 11.0
+chr10 3003494 3003496 100 9 0 9.0
+chr10 3003584 3003586 60 6 4 10.0
+chr10 3003924 3003926 100 4 0 4.0
+chr10 3003976 3003978 100 12 0 12.0
+chr10 3004035 3004037 71 10 4 14.0
+chr10 3004640 3004642 90 9 1 10.0
+chr10 3005170 3005172 28 2 5 7.0
+chr10 3005441 3005443 21 3 11 14.0
+chr10 3005994 3005996 50 7 7 14.0
+chr10 3008524 3008526 75 9 3 12.0
+chr10 3008594 3008596 90 9 1 10.0
+chr10 3008737 3008739 50 5 5 10.0
+chr10 3008827 3008829 80 12 3 15.0
+chr10 3008937 3008939 84 11 2 13.0
+chr10 3008985 3008987 66 6 3 9.0
+chr10 3009020 3009022 66 10 5 15.0
+chr10 3009103 3009105 50 8 8 16.0
+chr10 3009118 3009120 83 15 3 18.0
+chr10 3009207 3009209 66 6 3 9.0
+chr10 3009577 3009579 80 4 1 5.0
+chr10 3009627 3009629 100 4 0 4.0
+chr10 3010000 3010002 100 8 0 8.0
+chr10 3010399 3010401 58 10 7 17.0
+chr10 3010447 3010449 94 16 1 17.0
+chr10 3011056 3011058 85 6 1 7.0
+chr10 3011170 3011172 68 11 5 16.0
+chr10 3011280 3011282 66 10 5 15.0
+chr10 3011403 3011405 64 22 12 34.0
+chr10 3011468 3011470 92 12 1 13.0
+chr10 3011556 3011558 100 13 0 13.0
+chr10 3011580 3011582 68 11 5 16.0
+chr10 3011870 3011872 90 18 2 20.0
+chr10 3011999 3012001 81 18 4 22.0
+chr10 3012040 3012042 93 14 1 15.0
+chr10 3012053 3012055 66 10 5 15.0
+chr10 3012065 3012067 80 12 3 15.0
+chr10 3012214 3012216 61 11 7 18.0
+chr10 3012269 3012271 50 7 7 14.0
+chr10 3012319 3012321 33 5 10 15.0
+chr10 3012389 3012391 77 14 4 18.0
+chr10 3012415 3012417 77 14 4 18.0
+chr10 3012560 3012562 94 17 1 18.0
+chr10 3012609 3012611 85 12 2 14.0
+chr10 3012619 3012621 85 12 2 14.0
+chr10 3012741 3012743 56 13 10 23.0
+chr10 3012786 3012788 78 15 4 19.0
+chr10 3012790 3012792 63 12 7 19.0
+chr10 3012792 3012794 57 11 8 19.0
+chr10 3012794 3012796 94 17 1 18.0
+chr10 3012796 3012798 55 10 8 18.0
+chr10 3012798 3012800 66 12 6 18.0
+chr10 3012818 3012820 35 7 13 20.0
+chr10 3012824 3012826 16 4 20 24.0
+chr10 3012858 3012860 23 6 20 26.0
+chr10 3012940 3012942 20 2 8 10.0
+chr10 3012977 3012979 7 1 13 14.0
b
diff -r 000000000000 -r a06f7b5c4dc7 test-data/1_boxcar.bedgraph
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/1_boxcar.bedgraph Mon May 20 08:42:54 2019 -0400
b
@@ -0,0 +1,100 @@
+chr10 3000180 3000182 100 5 0 5.00 0
+chr10 3000185 3000187 100 5 0 5.00 0
+chr10 3000270 3000272 73 11 4 15.00 7.00
+chr10 3000294 3000296 100 14 0 14.00 10.00
+chr10 3000353 3000355 81 9 2 11.00 11.33
+chr10 3000488 3000490 66 4 2 6.00 9.00
+chr10 3000657 3000659 66 4 2 6.00 5.67
+chr10 3000709 3000711 80 4 1 5.00 4.00
+chr10 3000718 3000720 100 5 0 5.00 4.33
+chr10 3000923 3000925 80 4 1 5.00 4.33
+chr10 3001027 3001029 88 8 1 9.00 5.67
+chr10 3001029 3001031 100 9 0 9.00 7.00
+chr10 3001055 3001057 100 8 0 8.00 8.33
+chr10 3001188 3001190 75 3 1 4.00 6.67
+chr10 3001586 3001588 83 5 1 6.00 5.33
+chr10 3001747 3001749 100 7 0 7.00 5.00
+chr10 3001898 3001900 37 3 5 8.00 5.00
+chr10 3001908 3001910 55 5 4 9.00 5.00
+chr10 3001923 3001925 100 5 0 5.00 4.33
+chr10 3001928 3001930 100 5 0 5.00 5.00
+chr10 3001989 3001991 69 9 4 13.00 6.33
+chr10 3001996 3001998 84 11 2 13.00 8.33
+chr10 3002030 3002032 92 12 1 13.00 10.67
+chr10 3002077 3002079 77 7 2 9.00 10.00
+chr10 3002089 3002091 100 9 0 9.00 9.33
+chr10 3002093 3002095 44 4 5 9.00 6.67
+chr10 3002116 3002118 100 7 0 7.00 6.67
+chr10 3002147 3002149 80 8 2 10.00 6.33
+chr10 3002159 3002161 50 4 4 8.00 6.33
+chr10 3002261 3002263 100 16 0 16.00 9.33
+chr10 3002320 3002322 92 12 1 13.00 10.67
+chr10 3002514 3002516 100 4 0 4.00 10.67
+chr10 3002526 3002528 83 5 1 6.00 7.00
+chr10 3002539 3002541 100 8 0 8.00 5.67
+chr10 3002560 3002562 88 8 1 9.00 7.00
+chr10 3002594 3002596 100 8 0 8.00 8.00
+chr10 3002607 3002609 88 8 1 9.00 8.00
+chr10 3002628 3002630 100 12 0 12.00 9.33
+chr10 3002641 3002643 100 9 0 9.00 9.67
+chr10 3002904 3002906 100 9 0 9.00 10.00
+chr10 3002932 3002934 83 5 1 6.00 7.67
+chr10 3003427 3003429 61 8 5 13.00 7.33
+chr10 3003435 3003437 84 11 2 13.00 8.00
+chr10 3003465 3003467 90 10 1 11.00 9.67
+chr10 3003494 3003496 100 9 0 9.00 10.00
+chr10 3003584 3003586 60 6 4 10.00 8.33
+chr10 3003924 3003926 100 4 0 4.00 6.33
+chr10 3003976 3003978 100 12 0 12.00 7.33
+chr10 3004035 3004037 71 10 4 14.00 8.67
+chr10 3004640 3004642 90 9 1 10.00 10.33
+chr10 3005170 3005172 28 2 5 7.00 7.00
+chr10 3005441 3005443 21 3 11 14.00 4.67
+chr10 3005994 3005996 50 7 7 14.00 4.00
+chr10 3008524 3008526 75 9 3 12.00 6.33
+chr10 3008594 3008596 90 9 1 10.00 8.33
+chr10 3008737 3008739 50 5 5 10.00 7.67
+chr10 3008827 3008829 80 12 3 15.00 8.67
+chr10 3008937 3008939 84 11 2 13.00 9.33
+chr10 3008985 3008987 66 6 3 9.00 9.67
+chr10 3009020 3009022 66 10 5 15.00 9.00
+chr10 3009103 3009105 50 8 8 16.00 8.00
+chr10 3009118 3009120 83 15 3 18.00 11.00
+chr10 3009207 3009209 66 6 3 9.00 9.67
+chr10 3009577 3009579 80 4 1 5.00 8.33
+chr10 3009627 3009629 100 4 0 4.00 4.67
+chr10 3010000 3010002 100 8 0 8.00 5.33
+chr10 3010399 3010401 58 10 7 17.00 7.33
+chr10 3010447 3010449 94 16 1 17.00 11.33
+chr10 3011056 3011058 85 6 1 7.00 10.67
+chr10 3011170 3011172 68 11 5 16.00 11.00
+chr10 3011280 3011282 66 10 5 15.00 9.00
+chr10 3011403 3011405 64 22 12 34.00 14.33
+chr10 3011468 3011470 92 12 1 13.00 14.67
+chr10 3011556 3011558 100 13 0 13.00 15.67
+chr10 3011580 3011582 68 11 5 16.00 12.00
+chr10 3011870 3011872 90 18 2 20.00 14.00
+chr10 3011999 3012001 81 18 4 22.00 15.67
+chr10 3012040 3012042 93 14 1 15.00 16.67
+chr10 3012053 3012055 66 10 5 15.00 14.00
+chr10 3012065 3012067 80 12 3 15.00 12.00
+chr10 3012214 3012216 61 11 7 18.00 11.00
+chr10 3012269 3012271 50 7 7 14.00 10.00
+chr10 3012319 3012321 33 5 10 15.00 7.67
+chr10 3012389 3012391 77 14 4 18.00 8.67
+chr10 3012415 3012417 77 14 4 18.00 11.00
+chr10 3012560 3012562 94 17 1 18.00 15.00
+chr10 3012609 3012611 85 12 2 14.00 14.33
+chr10 3012619 3012621 85 12 2 14.00 13.67
+chr10 3012741 3012743 56 13 10 23.00 12.33
+chr10 3012786 3012788 78 15 4 19.00 13.33
+chr10 3012790 3012792 63 12 7 19.00 13.33
+chr10 3012792 3012794 57 11 8 19.00 12.67
+chr10 3012794 3012796 94 17 1 18.00 13.33
+chr10 3012796 3012798 55 10 8 18.00 12.67
+chr10 3012798 3012800 66 12 6 18.00 13.00
+chr10 3012818 3012820 35 7 13 20.00 9.67
+chr10 3012824 3012826 16 4 20 24.00 7.67
+chr10 3012858 3012860 23 6 20 26.00 5.67
+chr10 3012940 3012942 20 2 8 10.00 4.00
+chr10 3012977 3012979 7 1 13 14.00 3.00
b
diff -r 000000000000 -r a06f7b5c4dc7 test-data/1_hamming.bedgraph
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/1_hamming.bedgraph Mon May 20 08:42:54 2019 -0400
b
@@ -0,0 +1,100 @@
+chr10 3000180 3000182 100 5 0 5.00 0
+chr10 3000185 3000187 100 5 0 5.00 0
+chr10 3000270 3000272 73 11 4 15.00 21.00
+chr10 3000294 3000296 100 14 0 14.00 30.00
+chr10 3000353 3000355 81 9 2 11.00 34.00
+chr10 3000488 3000490 66 4 2 6.00 27.00
+chr10 3000657 3000659 66 4 2 6.00 17.00
+chr10 3000709 3000711 80 4 1 5.00 12.00
+chr10 3000718 3000720 100 5 0 5.00 13.00
+chr10 3000923 3000925 80 4 1 5.00 13.00
+chr10 3001027 3001029 88 8 1 9.00 17.00
+chr10 3001029 3001031 100 9 0 9.00 21.00
+chr10 3001055 3001057 100 8 0 8.00 25.00
+chr10 3001188 3001190 75 3 1 4.00 20.00
+chr10 3001586 3001588 83 5 1 6.00 16.00
+chr10 3001747 3001749 100 7 0 7.00 15.00
+chr10 3001898 3001900 37 3 5 8.00 15.00
+chr10 3001908 3001910 55 5 4 9.00 15.00
+chr10 3001923 3001925 100 5 0 5.00 13.00
+chr10 3001928 3001930 100 5 0 5.00 15.00
+chr10 3001989 3001991 69 9 4 13.00 19.00
+chr10 3001996 3001998 84 11 2 13.00 25.00
+chr10 3002030 3002032 92 12 1 13.00 32.00
+chr10 3002077 3002079 77 7 2 9.00 30.00
+chr10 3002089 3002091 100 9 0 9.00 28.00
+chr10 3002093 3002095 44 4 5 9.00 20.00
+chr10 3002116 3002118 100 7 0 7.00 20.00
+chr10 3002147 3002149 80 8 2 10.00 19.00
+chr10 3002159 3002161 50 4 4 8.00 19.00
+chr10 3002261 3002263 100 16 0 16.00 28.00
+chr10 3002320 3002322 92 12 1 13.00 32.00
+chr10 3002514 3002516 100 4 0 4.00 32.00
+chr10 3002526 3002528 83 5 1 6.00 21.00
+chr10 3002539 3002541 100 8 0 8.00 17.00
+chr10 3002560 3002562 88 8 1 9.00 21.00
+chr10 3002594 3002596 100 8 0 8.00 24.00
+chr10 3002607 3002609 88 8 1 9.00 24.00
+chr10 3002628 3002630 100 12 0 12.00 28.00
+chr10 3002641 3002643 100 9 0 9.00 29.00
+chr10 3002904 3002906 100 9 0 9.00 30.00
+chr10 3002932 3002934 83 5 1 6.00 23.00
+chr10 3003427 3003429 61 8 5 13.00 22.00
+chr10 3003435 3003437 84 11 2 13.00 24.00
+chr10 3003465 3003467 90 10 1 11.00 29.00
+chr10 3003494 3003496 100 9 0 9.00 30.00
+chr10 3003584 3003586 60 6 4 10.00 25.00
+chr10 3003924 3003926 100 4 0 4.00 19.00
+chr10 3003976 3003978 100 12 0 12.00 22.00
+chr10 3004035 3004037 71 10 4 14.00 26.00
+chr10 3004640 3004642 90 9 1 10.00 31.00
+chr10 3005170 3005172 28 2 5 7.00 21.00
+chr10 3005441 3005443 21 3 11 14.00 14.00
+chr10 3005994 3005996 50 7 7 14.00 12.00
+chr10 3008524 3008526 75 9 3 12.00 19.00
+chr10 3008594 3008596 90 9 1 10.00 25.00
+chr10 3008737 3008739 50 5 5 10.00 23.00
+chr10 3008827 3008829 80 12 3 15.00 26.00
+chr10 3008937 3008939 84 11 2 13.00 28.00
+chr10 3008985 3008987 66 6 3 9.00 29.00
+chr10 3009020 3009022 66 10 5 15.00 27.00
+chr10 3009103 3009105 50 8 8 16.00 24.00
+chr10 3009118 3009120 83 15 3 18.00 33.00
+chr10 3009207 3009209 66 6 3 9.00 29.00
+chr10 3009577 3009579 80 4 1 5.00 25.00
+chr10 3009627 3009629 100 4 0 4.00 14.00
+chr10 3010000 3010002 100 8 0 8.00 16.00
+chr10 3010399 3010401 58 10 7 17.00 22.00
+chr10 3010447 3010449 94 16 1 17.00 34.00
+chr10 3011056 3011058 85 6 1 7.00 32.00
+chr10 3011170 3011172 68 11 5 16.00 33.00
+chr10 3011280 3011282 66 10 5 15.00 27.00
+chr10 3011403 3011405 64 22 12 34.00 43.00
+chr10 3011468 3011470 92 12 1 13.00 44.00
+chr10 3011556 3011558 100 13 0 13.00 47.00
+chr10 3011580 3011582 68 11 5 16.00 36.00
+chr10 3011870 3011872 90 18 2 20.00 42.00
+chr10 3011999 3012001 81 18 4 22.00 47.00
+chr10 3012040 3012042 93 14 1 15.00 50.00
+chr10 3012053 3012055 66 10 5 15.00 42.00
+chr10 3012065 3012067 80 12 3 15.00 36.00
+chr10 3012214 3012216 61 11 7 18.00 33.00
+chr10 3012269 3012271 50 7 7 14.00 30.00
+chr10 3012319 3012321 33 5 10 15.00 23.00
+chr10 3012389 3012391 77 14 4 18.00 26.00
+chr10 3012415 3012417 77 14 4 18.00 33.00
+chr10 3012560 3012562 94 17 1 18.00 45.00
+chr10 3012609 3012611 85 12 2 14.00 43.00
+chr10 3012619 3012621 85 12 2 14.00 41.00
+chr10 3012741 3012743 56 13 10 23.00 37.00
+chr10 3012786 3012788 78 15 4 19.00 40.00
+chr10 3012790 3012792 63 12 7 19.00 40.00
+chr10 3012792 3012794 57 11 8 19.00 38.00
+chr10 3012794 3012796 94 17 1 18.00 40.00
+chr10 3012796 3012798 55 10 8 18.00 38.00
+chr10 3012798 3012800 66 12 6 18.00 39.00
+chr10 3012818 3012820 35 7 13 20.00 29.00
+chr10 3012824 3012826 16 4 20 24.00 23.00
+chr10 3012858 3012860 23 6 20 26.00 17.00
+chr10 3012940 3012942 20 2 8 10.00 12.00
+chr10 3012977 3012979 7 1 13 14.00 9.00
b
diff -r 000000000000 -r a06f7b5c4dc7 test-data/2_boxcar.bedgraph
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/2_boxcar.bedgraph Mon May 20 08:42:54 2019 -0400
b
@@ -0,0 +1,100 @@
+chr10 3000180 3000182 100 5 0 5.00 0
+chr10 3000185 3000187 100 5 0 5.00 0
+chr10 3000270 3000272 73 11 4 15.00 7.00
+chr10 3000294 3000296 100 14 0 14.00 10.00
+chr10 3000353 3000355 81 9 2 11.00 11.33
+chr10 3000488 3000490 66 4 2 6.00 9.00
+chr10 3000657 3000659 66 4 2 6.00 5.67
+chr10 3000709 3000711 80 4 1 5.00 4.00
+chr10 3000718 3000720 100 5 0 5.00 4.33
+chr10 3000923 3000925 80 4 1 5.00 4.33
+chr10 3001027 3001029 88 8 1 9.00 5.67
+chr10 3001029 3001031 100 9 0 9.00 7.00
+chr10 3001055 3001057 100 8 0 8.00 8.33
+chr10 3001188 3001190 75 3 1 4.00 6.67
+chr10 3001586 3001588 83 5 1 6.00 5.33
+chr10 3001747 3001749 100 7 0 7.00 5.00
+chr10 3001898 3001900 37 3 5 8.00 5.00
+chr10 3001908 3001910 55 5 4 9.00 5.00
+chr10 3001923 3001925 100 5 0 5.00 4.33
+chr10 3001928 3001930 100 5 0 5.00 5.00
+chr10 3001989 3001991 69 9 4 13.00 6.33
+chr10 3001996 3001998 84 11 2 13.00 8.33
+chr10 3002030 3002032 92 12 1 13.00 10.67
+chr10 3002077 3002079 77 7 2 9.00 10.00
+chr10 3002089 3002091 100 9 0 9.00 9.33
+chr10 3002093 3002095 44 4 5 9.00 6.67
+chr10 3002116 3002118 100 7 0 7.00 6.67
+chr10 3002147 3002149 80 8 2 10.00 6.33
+chr10 3002159 3002161 50 4 4 8.00 6.33
+chr10 3002261 3002263 100 16 0 16.00 9.33
+chr10 3002320 3002322 92 12 1 13.00 10.67
+chr10 3002514 3002516 100 4 0 4.00 10.67
+chr10 3002526 3002528 83 5 1 6.00 7.00
+chr10 3002539 3002541 100 8 0 8.00 5.67
+chr10 3002560 3002562 88 8 1 9.00 7.00
+chr10 3002594 3002596 100 8 0 8.00 8.00
+chr10 3002607 3002609 88 8 1 9.00 8.00
+chr10 3002628 3002630 100 12 0 12.00 9.33
+chr10 3002641 3002643 100 9 0 9.00 9.67
+chr10 3002904 3002906 100 9 0 9.00 10.00
+chr10 3002932 3002934 83 5 1 6.00 7.67
+chr10 3003427 3003429 61 8 5 13.00 7.33
+chr10 3003435 3003437 84 11 2 13.00 8.00
+chr10 3003465 3003467 90 10 1 11.00 9.67
+chr10 3003494 3003496 100 9 0 9.00 10.00
+chr10 3003584 3003586 60 6 4 10.00 8.33
+chr10 3003924 3003926 100 4 0 4.00 6.33
+chr10 3003976 3003978 100 12 0 12.00 7.33
+chr10 3004035 3004037 71 10 4 14.00 8.67
+chr10 3004640 3004642 90 9 1 10.00 10.33
+chr10 3005170 3005172 28 2 5 7.00 7.00
+chr10 3005441 3005443 21 3 11 14.00 4.67
+chr10 3005994 3005996 50 7 7 14.00 4.00
+chr10 3008524 3008526 75 9 3 12.00 6.33
+chr10 3008594 3008596 90 9 1 10.00 8.33
+chr10 3008737 3008739 50 5 5 10.00 7.67
+chr10 3008827 3008829 80 12 3 15.00 8.67
+chr10 3008937 3008939 84 11 2 13.00 9.33
+chr10 3008985 3008987 66 6 3 9.00 9.67
+chr10 3009020 3009022 66 10 5 15.00 9.00
+chr10 3009103 3009105 50 8 8 16.00 8.00
+chr10 3009118 3009120 83 15 3 18.00 11.00
+chr10 3009207 3009209 66 6 3 9.00 9.67
+chr10 3009577 3009579 80 4 1 5.00 8.33
+chr10 3009627 3009629 100 4 0 4.00 4.67
+chr10 3010000 3010002 100 8 0 8.00 5.33
+chr10 3010399 3010401 58 10 7 17.00 7.33
+chr10 3010447 3010449 94 16 1 17.00 11.33
+chr10 3011056 3011058 85 6 1 7.00 10.67
+chr10 3011170 3011172 68 11 5 16.00 11.00
+chr10 3011280 3011282 66 10 5 15.00 9.00
+chr10 3011403 3011405 64 22 12 34.00 14.33
+chr10 3011468 3011470 92 12 1 13.00 14.67
+chr10 3011556 3011558 100 13 0 13.00 15.67
+chr10 3011580 3011582 68 11 5 16.00 12.00
+chr10 3011870 3011872 90 18 2 20.00 14.00
+chr10 3011999 3012001 81 18 4 22.00 15.67
+chr10 3012040 3012042 93 14 1 15.00 16.67
+chr10 3012053 3012055 66 10 5 15.00 14.00
+chr10 3012065 3012067 80 12 3 15.00 12.00
+chr10 3012214 3012216 61 11 7 18.00 11.00
+chr10 3012269 3012271 50 7 7 14.00 10.00
+chr10 3012319 3012321 33 5 10 15.00 7.67
+chr10 3012389 3012391 77 14 4 18.00 8.67
+chr10 3012415 3012417 77 14 4 18.00 11.00
+chr10 3012560 3012562 94 17 1 18.00 15.00
+chr10 3012609 3012611 85 12 2 14.00 14.33
+chr10 3012619 3012621 85 12 2 14.00 13.67
+chr10 3012741 3012743 56 13 10 23.00 12.33
+chr10 3012786 3012788 78 15 4 19.00 13.33
+chr10 3012790 3012792 63 12 7 19.00 13.33
+chr10 3012792 3012794 57 11 8 19.00 12.67
+chr10 3012794 3012796 94 17 1 18.00 13.33
+chr10 3012796 3012798 55 10 8 18.00 12.67
+chr10 3012798 3012800 66 12 6 18.00 13.00
+chr10 3012818 3012820 35 7 13 20.00 9.67
+chr10 3012824 3012826 16 4 20 24.00 7.67
+chr10 3012858 3012860 23 6 20 26.00 5.67
+chr10 3012940 3012942 20 2 8 10.00 4.00
+chr10 3012977 3012979 7 1 13 14.00 3.00
b
diff -r 000000000000 -r a06f7b5c4dc7 test-data/2_hamming.bedgraph
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/2_hamming.bedgraph Mon May 20 08:42:54 2019 -0400
b
@@ -0,0 +1,100 @@
+chr10 3000180 3000182 100 5 0 5.00 0
+chr10 3000185 3000187 100 5 0 5.00 0
+chr10 3000270 3000272 73 11 4 15.00 6.28
+chr10 3000294 3000296 100 14 0 14.00 12.52
+chr10 3000353 3000355 81 9 2 11.00 15.60
+chr10 3000488 3000490 66 4 2 6.00 10.44
+chr10 3000657 3000659 66 4 2 6.00 5.04
+chr10 3000709 3000711 80 4 1 5.00 4.64
+chr10 3000718 3000720 100 5 0 5.00 4.72
+chr10 3000923 3000925 80 4 1 5.00 5.64
+chr10 3001027 3001029 88 8 1 9.00 5.04
+chr10 3001029 3001031 100 9 0 9.00 9.04
+chr10 3001055 3001057 100 8 0 8.00 10.28
+chr10 3001188 3001190 75 3 1 4.00 8.96
+chr10 3001586 3001588 83 5 1 6.00 4.04
+chr10 3001747 3001749 100 7 0 7.00 5.80
+chr10 3001898 3001900 37 3 5 8.00 7.64
+chr10 3001908 3001910 55 5 4 9.00 3.96
+chr10 3001923 3001925 100 5 0 5.00 5.64
+chr10 3001928 3001930 100 5 0 5.00 5.80
+chr10 3001989 3001991 69 9 4 13.00 6.12
+chr10 3001996 3001998 84 11 2 13.00 10.28
+chr10 3002030 3002032 92 12 1 13.00 12.68
+chr10 3002077 3002079 77 7 2 9.00 13.44
+chr10 3002089 3002091 100 9 0 9.00 8.68
+chr10 3002093 3002095 44 4 5 9.00 9.88
+chr10 3002116 3002118 100 7 0 7.00 5.28
+chr10 3002147 3002149 80 8 2 10.00 7.96
+chr10 3002159 3002161 50 4 4 8.00 8.88
+chr10 3002261 3002263 100 16 0 16.00 5.92
+chr10 3002320 3002322 92 12 1 13.00 17.28
+chr10 3002514 3002516 100 4 0 4.00 13.60
+chr10 3002526 3002528 83 5 1 6.00 5.36
+chr10 3002539 3002541 100 8 0 8.00 5.96
+chr10 3002560 3002562 88 8 1 9.00 9.04
+chr10 3002594 3002596 100 8 0 8.00 9.28
+chr10 3002607 3002609 88 8 1 9.00 9.28
+chr10 3002628 3002630 100 12 0 12.00 9.60
+chr10 3002641 3002643 100 9 0 9.00 13.36
+chr10 3002904 3002906 100 9 0 9.00 10.68
+chr10 3002932 3002934 83 5 1 6.00 10.12
+chr10 3003427 3003429 61 8 5 13.00 6.36
+chr10 3003435 3003437 84 11 2 13.00 9.28
+chr10 3003465 3003467 90 10 1 11.00 12.44
+chr10 3003494 3003496 100 9 0 9.00 11.60
+chr10 3003584 3003586 60 6 4 10.00 10.28
+chr10 3003924 3003926 100 4 0 4.00 7.04
+chr10 3003976 3003978 100 12 0 12.00 5.44
+chr10 3004035 3004037 71 10 4 14.00 13.12
+chr10 3004640 3004642 90 9 1 10.00 11.68
+chr10 3005170 3005172 28 2 5 7.00 9.96
+chr10 3005441 3005443 21 3 11 14.00 2.96
+chr10 3005994 3005996 50 7 7 14.00 3.72
+chr10 3008524 3008526 75 9 3 12.00 7.96
+chr10 3008594 3008596 90 9 1 10.00 10.28
+chr10 3008737 3008739 50 5 5 10.00 10.12
+chr10 3008827 3008829 80 12 3 15.00 6.68
+chr10 3008937 3008939 84 11 2 13.00 13.28
+chr10 3008985 3008987 66 6 3 9.00 12.44
+chr10 3009020 3009022 66 10 5 15.00 7.68
+chr10 3009103 3009105 50 8 8 16.00 11.12
+chr10 3009118 3009120 83 15 3 18.00 10.00
+chr10 3009207 3009209 66 6 3 9.00 16.12
+chr10 3009577 3009579 80 4 1 5.00 7.52
+chr10 3009627 3009629 100 4 0 4.00 4.80
+chr10 3010000 3010002 100 8 0 8.00 4.96
+chr10 3010399 3010401 58 10 7 17.00 9.12
+chr10 3010447 3010449 94 16 1 17.00 11.92
+chr10 3011056 3011058 85 6 1 7.00 17.28
+chr10 3011170 3011172 68 11 5 16.00 8.16
+chr10 3011280 3011282 66 10 5 15.00 12.28
+chr10 3011403 3011405 64 22 12 34.00 12.64
+chr10 3011468 3011470 92 12 1 13.00 23.76
+chr10 3011556 3011558 100 13 0 13.00 14.80
+chr10 3011580 3011582 68 11 5 16.00 14.84
+chr10 3011870 3011872 90 18 2 20.00 13.48
+chr10 3011999 3012001 81 18 4 22.00 20.32
+chr10 3012040 3012042 93 14 1 15.00 20.56
+chr10 3012053 3012055 66 10 5 15.00 16.24
+chr10 3012065 3012067 80 12 3 15.00 12.08
+chr10 3012214 3012216 61 11 7 18.00 13.68
+chr10 3012269 3012271 50 7 7 14.00 12.52
+chr10 3012319 3012321 33 5 10 15.00 8.28
+chr10 3012389 3012391 77 14 4 18.00 6.68
+chr10 3012415 3012417 77 14 4 18.00 15.52
+chr10 3012560 3012562 94 17 1 18.00 16.48
+chr10 3012609 3012611 85 12 2 14.00 19.08
+chr10 3012619 3012621 85 12 2 14.00 14.32
+chr10 3012741 3012743 56 13 10 23.00 14.00
+chr10 3012786 3012788 78 15 4 19.00 15.16
+chr10 3012790 3012792 63 12 7 19.00 17.00
+chr10 3012792 3012794 57 11 8 19.00 14.08
+chr10 3012794 3012796 94 17 1 18.00 13.32
+chr10 3012796 3012798 55 10 8 18.00 18.68
+chr10 3012798 3012800 66 12 6 18.00 12.32
+chr10 3012818 3012820 35 7 13 20.00 13.36
+chr10 3012824 3012826 16 4 20 24.00 8.28
+chr10 3012858 3012860 23 6 20 26.00 5.04
+chr10 3012940 3012942 20 2 8 10.00 6.48
+chr10 3012977 3012979 7 1 13 14.00 2.56