SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



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

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

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

Hi,

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)

output:

$samples
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!

A
alisonewr is offline   Reply With Quote
Old 05-04-2014, 11:48 PM   #2
Gordon Smyth
Member
 
Location: Melbourne, Australia

Join Date: Apr 2011
Posts: 91
Default

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, 08:48 AM   #3
Gonza
Member
 
Location: Ithaca, NY

Join Date: Mar 2013
Posts: 78
Default

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?
thanks

normalization:

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

Tags
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:27 AM.


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