MACS to Signac
Hey guys--
Once MACS is completed running, how would I feed the peaks/bed file into Signac for each sample up to merging and finally clustering? I believe the vignette only demonstrates if one were to add *fragment files and I need the peaks file from MACS for a stronger signal.
• 1,735 views
•
link
updated
by
ATpoint
•
written
by
bioinformatics.girl •
•
written
by
bioinformatics.girl •
1 answer
If you have externally defined peaks, you can use
library(rtracklayer)
peaks = import('/path/to/peaks.bed')
FeatureMatrix(
fragments=CreateFragmentObject('/path/to/fragments.tsv'),
features=peaks
)
This will generate a peak quantification matrix, which you can push into a Seurat object for downstream analysis.
Log in to answer this question.
More posts like this
-
Joint RNA and ATAC analysis: 10x multiomic - Fragment Path
written by cgp09741 •
-
Zero peaks detected using findMotifsGenome.pl
written by ichmagbaeume •
-
Get MACs2 peaks associated with specific genes
written by ahmed.saadawi •
-
eCLIP peak calling
written by GiusiG •
-
Ideas to split peaks into subpeaks - beyond PeakSplitter
written by Sergio Martínez Cuesta
-
Calculate enrichment values of Chip-seq peaks for hierarchical clustering
written by florian.noack •
-
ChIP-seq peak calling & visualized in IGV
written by izzy.yichao.cai
-
How to find overlapping regions among three bed files [Solved]
written by Bioradical •
-
Paired-end ChIP-Seq Peak Calling
written by simon.pearce •
-
Combined Tss Plots For Chip Seq Peaks
written by aditi.qamra
Please go through your previous questions and upvote helpful answers, mark the one that solved the issue as accepted. You opened an extensive number of questions and barely or did not do any upvoting st all. That's not how an open community works.