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


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?


ADD COMMENTlink modified 6.2 years ago by Devon Ryan95k • written 6.2 years ago by AW350

Any questions relating to edgeR may be best answered via the 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 modified 6 months ago by RamRS27k • written 6.2 years ago by Jason20
gravatar for Devon Ryan
6.2 years ago by
Devon Ryan95k
Freiburg, Germany
Devon Ryan95k 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 modified 6 months ago by RamRS27k • written 6.2 years ago by Devon Ryan95k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 590 users visited in the last hour