deepTools
deepTools copied to clipboard
Compute matrix -a -b influence the signal
When I use deepTools computeMatrix 3.1.2 to calculate the different TF binding signal around the peaks,I meet some problem. When we extend the -a -b region from 1kb,2.5kb,5kb,10kb,25kb to 50kb with the bin size also increase, the enrich pattern of different peaks changed.(python 2.7.15).
Code: computeMatrix reference-point -S FOXA1.bigwig -R SE_enhancer.bed TE_enhancer.bed Open_region.bed -o Path_out/FOXA1.mat.gz --referencePoint center -b 1000 -a 1000 --binSize 10 -p 15
(total 2kb region bin 10)
(total 5kb region bin 10)
why these will happen,and we icontinue to increase the region length ,the signal for TE and OR continue to drop.
Oh, that's incredibly strange. Can you post your bigWig files somewhere so we can debug this?
I downloaded bigwigs from Encode database. the Enocde ID for FXOA1 is ENCFF795BHZ.