Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
edgeR VS deseq normalization narges Bioinformatics 0 07-04-2013 06:27 AM
RPKM and RNAseq data normalization drghosh1139 Introductions 3 07-24-2012 07:31 AM
EdgeR Normalization Question R_user RNA Sequencing 1 05-07-2012 11:14 PM
RPKM normalization between samples bt7 Bioinformatics 1 08-11-2011 11:30 AM
TMM nomalization in RNA-seq minghui Bioinformatics 1 05-14-2010 03:01 AM

Thread Tools
Old 04-29-2014, 09:49 AM   #1
Junior Member
Location: London

Join Date: Nov 2012
Posts: 1
Default RNA-seq normalization: How to use TMM and rpkm() in EdgeR??


I have some RNA-seq samples that I want to normalize and then output RPKM expression, but I am unsure how to do this.

This is my pipeline so far:

1. Normalise raw read counts with TMM in edgeR

expr <- DGEList(counts=data, group=conditions)
expr <- calcNormFactors(expr)


group lib.size norm.factors
Sample1 F 19770521 1.0462660
Sample2 F 17970679 0.8794805
Sample3 F 19184265 1.0573665

QUESTION: How do I get normalized raw read counts from this? Do I multiply the read counts by the norm.factors?

QUESTION: Ultimately, I want to end up with RPKM values for each gene in each sample. I know I can use the rpkm() function below in edgeR

expr_norm <- rpkm(expr, log=FALSE,gene.length=vector)

but is expr the output from calcNormFactors or something else?

Thanks for your help!

alisonewr is offline   Reply With Quote
Old 05-05-2014, 12:48 AM   #2
Gordon Smyth
Location: Melbourne, Australia

Join Date: Apr 2011
Posts: 91

Yes, you simply use

rpkm(expr, log=FALSE,gene.length=vector)

where expr is the output from calcNormFactors(). edgeR is an object orientated package where the DGEList data object assemblies all the required information.

There is no need for you to compute normalized raw counts (in fact I would argue there isn't any such thing). edgeR uses the normalization factors automatically.
Gordon Smyth is offline   Reply With Quote
Old 10-27-2014, 09:48 AM   #3
Location: Ithaca, NY

Join Date: Mar 2013
Posts: 78

Hello, does it really matter which normalization you use with edgeR? I guess I do not understand really well the differences.
I am using TMM now but I have no 'solid' argument why i chose that instead of upperquartile or RLE.

Any thoughts?


#y <- calcNormFactors(y, method="TMM")
#y <- calcNormFactors(y, method="upperquartile")
#y <- calcNormFactors(y, method="RLE")
Gonza is offline   Reply With Quote

edger, normalisation, rna-seq, tmm

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 07:47 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO