Question: PeakNarrow ATAC-Seq analysis
2
gravatar for Spacebio
3.0 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 • 1.1k views
ADD COMMENTlink modified 11 months ago by ATpoint39k • written 3.0 years ago by Spacebio180
1
gravatar for igor
3.0 years ago by
igor11k
United States
igor11k 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 3.0 years ago by igor11k
1
gravatar for ATpoint
11 months ago by
ATpoint39k
Germany
ATpoint39k 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 11 months ago • written 11 months ago by ATpoint39k
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: 1210 users visited in the last hour