'max' function in bigWigMerge

[KAWAJI Hideya]

Dear the developers in UCSC,

I needed “max” function over a series of bigWig files, and I’d like to
propose ‘-max’ option in utils/bigWigMerge. I gave a try to implement
this as seen below.

https://gist.github.com/hkawaji/6d0913fdc4c3e2db6427

I’d be happy if you find this idea as useful and include (or re-implement) it
in future release of the source tree.

best wishes,
Kawaji H

(I am re-sending this in case that my previous e-mail is rejected due to attachment).


--- message of the ‘patched’ version of bigWigMerge ---

$ ~/bin/x86_64/bigWigMerge
bigWigMerge v2 - Merge together multiple bigWigs into a single output bedGraph.
You'll have to run bedGraphToBigWig to make the output bigWig.
The signal values are just added together to merge them
usage:
bigWigMerge in1.bw in2.bw .. inN.bw out.bedGraph
options:
-threshold=0.N - don't output values at or below this threshold. Default is 0.0
-adjust=0.N - add adjustment to each value
-clip=NNN.N - values higher than this are clipped to this value
-inList - input file are lists of file names of bigWigs
-max - calculate max values, instead of sum (added feature)

[Steve Heitner]

Hello, Kawaji.

Thank you for letting us know of your modification to bigWigMerge. We believe this is a valuable option and one of our developers is currently looking into incorporating this into a future release. We'll let you know when this becomes available.

Please contact us again at gen...@soe.ucsc.edu if you have any further questions. All messages sent to that address are archived on a publicly-accessible Google Groups forum. If your question includes sensitive data, you may send it instead to genom...@soe.ucsc.edu.

---
Steve Heitner
UCSC Genome Bioinformatics Group

--

[Steve Heitner]

Hello, Kawaji.

I just wanted to let you know that the new version of bigWigMerge with the max function went out with our latest software release yesterday.

[KAWAJI Hideya]

Dear Steve and the development team,

Thanks a lot for letting me know. I’m glad my contribution is incorporated
and I am able to use the latest tool for that calculation.

kind regrads,
Kawaji H