Hello all,
I'm trying to compare coverage within specific introns across different RNA-seq experiments, and I was hoping to do so through bedgraph files. Would deeptools bamCoverage normalization be adequate for enabling comparisons between bedgraphs? I've been doing so using the normalizeTo1x flag.
•
written
by
albert.dayuan.yu •
1 answer
You could use the normalization factors produced by e.g. edgeR to scale the raw counts. For this you would need to calculate norm. factors in R using a standard count matrix e.g. with the exon counts (you could use featureCounts for creation of the matrix directly from the BAM files) and then use edgeR for normalization followed by providing these factors to bamCoverage. See the link below for example code:
Log in to answer this question.
-
Comparing P-Values & FDR adjusted p-values between RNA-Seq experiments?
written by Saran •
-
Technical single-cell question
written by karen •
-
Comparison of ChIP-seq tracks from two cell lines using FC over control (ENCODE data)
written by erica.fary •
-
Normalizing RNA-seq read counts within a single sample - gut check
written by plberry •
-
Comparison of ratio of expression between two genes across samples
written by guillepalou4 •
-
BAM file and Bedgraph file don't match after bedtools genomecov?
written by albert.dayuan.yu •
-
Best-suited unit of expression to perform between-sample comparison
written by leaodel
-
Comparison of ChIPseq enrichment in specific regions
written by Roman Hillje •
-
Within-sample comparison of TMM normalized values?
written by Ekarl2
-
Microarray Across Experiments
written by Greg Clark •
