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

Hi,

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

diffhic GRanges • 643 views
0
Entering edit mode
Aaron Lun ★ 27k
@alun
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.