I agree that the profile or metagene plot is useful and widely applied.
Creating such a plot consists of several steps, including:
- quantifying and accumulating signal around sites of interest
- normalization
- visualization
The qProfile() function in QuasR was designed to do the
quantification/
Bsmooth/bsseq (and QuasR, as well) are great for manipulating and analyzing
BSseq data, but there's something in this email that I've never really
found a great solution for, and that is so-called "profile" or "metagene"
plots (the latter is a confusing bit of terminology, given that it is
re-used
