## draw the coverage profile for a coverate matrix that was previously calculated
# using the cov.matrix function for a BAM file containing ChIP-seq data
# from a H3K36me3 histone modification experiment
data(DF_H3K36me3)
draw.profile(DF_H3K36me3,ylab="coverage",outfile="testProfile.png",main="testProfile",col="red")
## Now, draw two profiles for 2 different coverage matrices previously
# obtained for the same H3K36me3 histone modification experiment and its
# respective 'control' file
data(DF_H3K36me3_control)
#create a list with the two matrices
input_l=list(DF_H3K36me3,DF_H3K36me3_control)
draw.profile(input_l,ylab="coverage",outfile="testProfile.png",main="testProfile")
Run the code above in your browser using DataLab