Skip to content
New issue

Have a question about this project? # for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “#”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? # to your account

computeMatrix on 2 regions with replicates #1318

Open
asmariyaz23 opened this issue Jul 15, 2024 · 0 comments
Open

computeMatrix on 2 regions with replicates #1318

asmariyaz23 opened this issue Jul 15, 2024 · 0 comments

Comments

@asmariyaz23
Copy link

asmariyaz23 commented Jul 15, 2024

Hello,

I am new to deepTools and want to run computeMatrix (twice) to study gene promoter effect during activation on 1 mark (H3K4me3) by comparing 2 cell types:

  1. activated WT vs quiescent WT
  2. activated mdx vs quiescent mdx. I have 3 replicates per sample.

I merged the BED files (promoter regions) for replicates and provide the bigWig files like so:

rep1_cond1="${gt1}_${st1}_${mark}_R1.bigWig"
rep2_cond1="${gt1}_${st1}_${mark}_R2.bigWig"
rep3_cond1="${gt1}_${st1}_${mark}_R3.bigWig"
rep1_cond2="${gt2}_${st2}_${mark}_R1.bigWig"
rep2_cond2="${gt2}_${st2}_${mark}_R2.bigWig"
rep3_cond2="${gt2}_${st2}_${mark}_R3.bigWig" 
computeMatrix reference-point \
    --regionsFileName "${peak_ip}/${gt1}-${st1}_mLb.clN_${mark}_promoters.bed" "${peak_ip}/${gt2}-${st2}_mLb.clN_${mark}_promoters.bed" \
    --scoreFileName ${rep1_cond1} ${rep2_cond1} ${rep3_cond1} ${rep1_cond2} ${rep2_cond2} ${rep3_cond2} \
    --referencePoint TSS \
    --blackListFileName /scratch/asmab/mm10-blacklist.v2.bed \
    --missingDataAsZero \
    --outFileName "${outdir}/${gt1}_${st1}_vs_${gt2}_${st2}_${mark}_computeMatrix" \
    --numberOfProcessors 1

I expected 2 rows representing the 2 bed files provided (which I get), however I get 6 columns (I expected 3 columns). Am I running computeMatrix correctly? Could you please help? I am running deeptools version 3.5.5

Thank you,
Asma

# for free to join this conversation on GitHub. Already have an account? # to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant