HiCExplorer
HiCExplorer copied to clipboard
What are the differences of the pearson matrix generated by hicPCA and hicTransform?
Hello, For A / B compartments analysis, I plot the A / B compartments based on the pearson matrix. The pearson matrix was generated by two tools (hicPCA and hicTransform) separately. I found that visualizaion of A / B compartments is a little different.
a.) the commandlines for hicPCA:
hicPCA -m Mus.100Kb.corrected.h5
--numberOfEigenvectors 2 --method lieberman
--outputFileName Mus.100Kb.pca1.bw Mus.100Kb.pca2.bw --format bigwig
--pearsonMatrix Mus.100Kb.hicPCA.pearson.h5
--obsexpMatrix Mus.100Kb.hicPCA.obsexp.h5
b.) and then plot A / B compartments based on the pearson matrix generated by hicPCA :
hicPlotMatrix -m Mus.100Kb.hicPCA.pearson.h5
--outFileName Mus.100Kb.hicPCA.png --perChr --bigwig Mus.100Kb.pca1.bw
the checkerboard pattern of A / B compartments is abvious (using one chromosome as the case):

c.) the commandlines for hicTransform:
hicTransform -m Mus.100Kb.corrected.h5 --perChromosome
--outFileName Mus.100Kb.hicTransform.pearson.h5 --method pearson
Of note, accroding to #580, i add the prameter --perChromosome but no changes for final plot. In addition, --method all in this document (https://hicexplorer.readthedocs.io/en/latest/content/example_usage.html#a-b-compartment-analysis) is not valid options in my used version of hicexplorer (v3.6). So I use --method pearson instead.
d.) and then plot A / B compartments based on the pearson matrix generated by hicTransform:
hicPlotMatrix -m Mus.100Kb.hicTransform.pearson.h5
--outFileName Mus.100Kb.hicTransform.png --perChr --bigwig Mus.100Kb.pca1.bw
However, the checkerboard pattern of A / B compartments is not very abvious compared the plot above (using the same chromosome as the case):

e.) And the matrix files Mus.100Kb.hicPCA.pearson.h5 and Mus.100Kb.hicTransform.pearson.h5 are not the same files. According to my understanding, the peason matix was claculated by hicPCA and hicTransform (--perChromosome) with equivalent method (https://github.com/deeptools/HiCExplorer/issues/580) per chromosom. If so, why the plots of A / B compartments are a little different? Which pearson matrix is recommended for A / B compartments analysis?