How to run plotPlaid in diffHic Package using output from hic Pro
1
0
Entering edit mode
Andrew • 0
@andrew-11520
Last seen 5.0 years ago

Hi Aaron,

Is there any way to run plotPlain in diffHic package using output from HiC-Pro? saying using a precomputed HiC matrix from HiC-Pro, and a bed file with bin coordinates. Thanks.

plotPlaid diffHic • 619 views
0
Entering edit mode
Aaron Lun ★ 27k
@alun
Last seen 2 hours ago
The city by the bay

The easiest approach is to take the BAM files produced by HiC-Pro, use preparePairs to generate HDF5 files, and then use those files in plotPlaid. Trying to do so from a precomputed matrix is somewhat confusing as you'll have to pretend that your bins are fragments in order to use plotPlaid.

Alternatively, I guess you could use plotDI for plotting instead. This requires that you generate an InteractionSet from the HiC-Pro contact matrix, and then pretend the normalised intensities are log-fold changes. It should be trivial to construct a ContactMatrix object from the matrix and coordinates; you can then use the deflate method to easily convert this to an InteractionSet.