see the paper here http://genomebiology.com/content/14/4/R38
it is a R package, may consider to use it.
jMOSAiCS: joint analysis of multiple ChIP-seq datasets
Originally published on my old Blogger blog on July 23, 2013.
Read the original post and its comments.