Question: EdgeR TMM normalisation
4
gravatar for AW
4.3 years ago by
AW350
United Kingdom
AW350 wrote:

Hi,

I have some RNA-seq samples that I want to normalize and then output RPKM expression, and I will use the following commands from EdgeR.

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

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

Id be very grateful if you could answer these questions.

Q1. When creating the expr <- DGEList(counts=data, group=conditions), what effect does specifying groups have one the TMM normalisation? How does TMM use this information and how would the results differ if you did specify groups versus not?

Q2. The expression data I am using was obtained from mapping reads onto denovo contigs assembled with Trinity. I then chose the most highly expressed contig from each cluster as the "best isoform" and then summed expression across all the contigs in the cluster as the expression value for that cluster. Therefore I do not have one obvious gene length to use. Should I use the longest contig from the cluster?

Thanks!

ADD COMMENTlink modified 4.3 years ago by Devon Ryan82k • written 4.3 years ago by AW350

Any questions relating to edgeR may be best answered via the mailing list.

http://www.bioconductor.org/help/mailing-list/

If I were you I'd search the mailing list first to see if someone else has asked a similar question. If not make a post. The developer of edgeR and other R bioconductor packages are pretty active on the mailing lists. 

ADD REPLYlink written 4.3 years ago by jmillo46860
1
gravatar for Devon Ryan
4.3 years ago by
Devon Ryan82k
Freiburg, Germany
Devon Ryan82k wrote:

calcNormFactors() doesn't use the group information. In fact, if you leave it out when making the DGEList object then every sample is given its own group. In other words, it makes no difference. The group information is used by other functions when you specify a design.

 

I've never used trinity, so I can't give any guidance there (perhaps make that its own question and give it an appropriate title).

ADD COMMENTlink written 4.3 years ago by Devon Ryan82k
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: 1736 users visited in the last hour