## draw the heatmap 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.heatmap(DF_H3K36me3,outfile="testHeatmap.png")
## Now, draw two heatmaps 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.heatmap(input_l,outfile="testHeatmap.png")
Run the code above in your browser using DataLab