Question: PeakNarrow ATAC-Seq analysis
2
gravatar for Spacebio
2.2 years ago by
Spacebio180
Spacebio180 wrote:

Hello biostars,

I am new to this type of analysis. I have IDR.NarrowPeak data and I would like to normalise it to FPKM using R. I read that the packages limma and DESeq are helpful but I still have problems creating the count matrix. Is there a step I am missing?

Any help/suggestion will be appreciated.

narrowpeak atac-seq • 932 views
ADD COMMENTlink modified 4 weeks ago by ATpoint25k • written 2.2 years ago by Spacebio180
1
gravatar for igor
2.2 years ago by
igor8.8k
United States
igor8.8k wrote:

You should try DiffBind, which uses DESeq and edgeR in the background and automates this process: https://bioconductor.org/packages/release/bioc/html/DiffBind.html

ADD COMMENTlink written 2.2 years ago by igor8.8k
1
gravatar for ATpoint
4 weeks ago by
ATpoint25k
Germany
ATpoint25k wrote:

The peak file is the template for the count matrix. You can e.g. use featureCounts to get the counts per peak and sample directly from the BAM files:

## 1. Transform narrowPeak to the featureCount SAF format:
awk 'OFS="\t" {print $1"_"$2+1"_"$3, $1, $2+1, $3, "+"}' in.narrowPeak > out.saf

## 2. Run fC:
featureCounts -a out.saf -F SAF -o countmatrix.txt *.bam

This countmatrix you can normalize with e.g. DESeq2 or edgeR or transformations such as vst and rlog. Stay away from naive methods such as R/FPKM. They are not suitable for inter-sample comparisons.

ADD COMMENTlink modified 4 weeks ago • written 4 weeks ago by ATpoint25k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1795 users visited in the last hour