diffHic package plotPlaid. Visualize multiple inter-chromosomal and intra-chromosomal interactions in the same heatmap
Entering edit mode
Last seen 5.8 years ago


Is it possible to draw a heatmap with plotPlaid from some (or all) chromosomes together? How?


diffhic GRanges • 643 views
Entering edit mode
Aaron Lun ★ 27k
Last seen 12 hours ago
The city by the bay

No, it's not possible to do that. The plotPlaid function requires one genomic interval on the x- and y-axes. This means you can't specify multiple chromosomes, because that doesn't fit into a single GRanges entry.

Personally, I don't find multi-chromosome plaid plots to be particularly interesting. You get boxes along the diagonal for the strong intra-chromosomal contacts, and nothing much else. There are some exceptions, e.g., yeast data where centromeres of different chromosomes interact with each other, or cancer cell lines with translocations. But in general, most of the inter-chromosomal space is pretty empty and boring to look at.

If you really need it, you could probably mimic it with something like:

all.chrs <- range(param$fragments)
nchrs <- length(all.chrs)
par(mfrow=c(nchrs, nchrs))
for (x in all.chrs) {
    for (y in all.chrs) {
        plotPlaid("myfile.h5", param, x, y) # plot between current chr pair

This will make a plot of nchrs by nchrs subplots, where each subplot represents the plot between two chromosomes. You could probably get it looking reasonably nice by reducing the margins for each plot (using mar= in par) to cut down the whitespace.


Login before adding your answer.

Traffic: 303 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6