Solution for easyRNASeq RPKM value = NA

[weasteam]

I was having problems in generate the rpkm values from easyRNASeq, everything I got is "NA".

The problem is because the easyRNASeq could not calculate the total reads from one bam file. The libSizes generated from the below code is "NA". It looks like the column names of the read table is bam file name, but the name of libSizes is the bam file name with absolute file address. Remove the additional filesDirectory, the problem will be solved.

libSizes <- libSizes[match(colnames(mCounts),names(libSizes))]

To do that, first download the easyRNASeq package from Bioconductor: http://www.bioconductor.org/packages...asyRNASeq.html

wget http://www.bioconductor.org/packages/2.12/bioc/src/contrib/easyRNASeq_1.6.0.tar.gz

Then unzip it

tar zxvf easyRNASeq_1.6.0.tar.gz

Edit the file easyRNASeq-normalize.R under R folder.
Find the following code:

libSizes <- librarySize(obj)
libSizes <- libSizes[match(colnames(mCounts),names(libSizes))]

Modify the above code into below:

libSizes <- librarySize(obj)
#insert these code
libSizeName=names(libSizes)
for (i in 1:length(libSizeName)){
     tmp=unlist(strsplit(libSizeName[i],"/"))
     libSizeName[i]=tmp[length(tmp)]
}
names(libSizes) = libSizeName

libSizes <- libSizes[match(colnames(mCounts),names(libSizes))]

Save the file and reinstall the easyRNASeq package with (u might need root to do so)

R CMD INSTALL easyRNASeq

Rerun the R code you are able to see the normalized value.



==============================================
PHDcell: the best online free reference manager