Question: How to run the rpkm function of edgeR
gravatar for Gary
2.7 years ago by
Gary0 wrote:


I would like to use edgeR to estimate the RPKM values.  I have (1) read counts files estimated by HTSeq-count, and (2) a transcript length file. I know how to estimate CPM in edgeR, using below command lines. However, I don’t know how to estimate RPKM values based on the files I have. Could you show me how to run the command lines in edgeR? Many thanks.




Transcript count files

Transcript length file


> targets=readTargets()

> targets

             files group  description

1 T8_Early_R_T.txt Early T8_Early_R_T

2  T8_Late_R_T.txt  Late  T8_Late_R_T

3 T9_Early_R_T.txt Early T9_Early_R_T

4  T9_Late_R_T.txt  Late  T9_Late_R_T

> d=readDGE(targets,skip=0,comment.char="!")

> keep=rowSums(cpm(d)>1)>=2

> d=d[keep,,keep.lib.sizes=FALSE]

> d=calcNormFactors(d,method="TMM")

> dim(d)

[1] 8202    4

> detags=rownames(d)

> detagsfile=cpm(d)[detags,]

> write.table(detagsfile,file="CPM_MSC_8202ene.txt")

> rpkm<-rpkm(d)

Error in rpkm.DGEList(d) : Gene lengths not found

> rpkm<-rpkm(d,gene.lengths)

Error in rpkm.DGEList(d, gene.lengths) : object 'gene.lengths' not found

ADD COMMENTlink modified 2.7 years ago by Aaron Lun21k • written 2.7 years ago by Gary0
gravatar for Aaron Lun
2.7 years ago by
Aaron Lun21k
Cambridge, United Kingdom
Aaron Lun21k wrote:

Well, it would be a good start to load the gene lengths into memory. <- read.table("hg19refGene20160206_TranscriptLength.txt", header=TRUE)

You then need to match the gene names in d to the ordering of

m <- match(rownames(d),$Transcript)
gene.lengths <-$TranscriptLength[m]

... and then you should be able to run rpkm on that.

ADD COMMENTlink modified 2.7 years ago • written 2.7 years ago by Aaron Lun21k

Hi Aaron,

Thank you so much.

ADD REPLYlink written 2.7 years ago by Gary0
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.2.0
Traffic: 189 users visited in the last hour