robinweide / GENOVA

GENome Organisation Visual Analytics
GNU General Public License v3.0
68 stars 15 forks source link

Issues related to "chip" in function "hic_matrixplot" #326

Open ccmeth opened 1 year ago

ccmeth commented 1 year ago

Dear developers,

Thanks for providing a powerful tool for us to deal with Hi-C data. Recently, I used hic_matrixplot to exhibit Hi-C contact matrices as well as H3K4me3/H3K9me3 signals, however, when I used "chip = list(h3k4me3,h3k9me3)", the two tracks seemed to overlap together.

Code: h3k4me3 <- "~/methylhic/DSS_results/chromHMM/GSM1000124_mm9_wgEncodeLicrHistoneEse14H3k04me3ME0129olaStdSig.bigWig" h3k9me3 <- "~/methylhic/DSS_results/chromHMM/GSM1003751_mm9_wgEncodeSydhHistEse14H3k09me3StdSig.bigWig"

hic_matrixplot(exp1 = wt_tettko[[1]],exp2 = wt_tettko[[2]],chrom = 'chr7',start = 29e6,end = 30e6,colour_lim = c(0, 150),chip = list(h3k9me3),chip.yMax = 25) image

hic_matrixplot(exp1 = wt_tettko[[1]],exp2 = wt_tettko[[2]],chrom = 'chr7',start = 29e6,end = 30e6,colour_lim = c(0, 150),chip = list(h3k4me3),chip.yMax = 25) image

hic_matrixplot(exp1 = wt_tettko[[1]],exp2 = wt_tettko[[2]],chrom = 'chr7',start = 29e6,end = 30e6,colour_lim = c(0, 150),chip = list(h3k9me3,h3k4me3),chip.yMax = 25) image