calculating RPKM and VST by DEseq2
2
0
Entering edit mode
Angel ▴ 40
@angel-7981
Last seen 7.3 years ago
Berlin

friends,

i have a file containing raw counts produced by featurecounts. now i need RPKM and VST...i read DEseq2 vignette but really i could not get what to do...supposing my raw counts file named counts.txt and placed in  /usr/data/nfs6/izadi/microarray/subread-1.4.6-p5-source/ and i installed DEseq2 package...please help me to calculate the RPKM and VST.

for example in the below i don't know what i should replace in the code

library("pasilla")

library("Biobase")

data("pasillaGenes")

countData <- counts(pasillaGenes)

colData <- pData(pasillaGenes)[,c("condition","type")] 

i am not familiar enough with R but i need RPKM and VST files...

thank you

DEseq2 • 5.8k views
ADD COMMENT
1
Entering edit mode
@mikelove
Last seen 12 hours ago
United States

We don't offer any functionality for helping you to calculate the feature size (the gene length), if you are importing a matrix of counts. So you would have to calculate the exonic basepairs for each gene on your own, using import() from rtracklayer and functions from the GRanges package. Or you could use software like RSEM or Cufflinks to calculate the FPKM estimate outside of R.

(Although if you do the fragment counting with the summarizeOverlaps() function, this information is included as rowRanges and we have a function fpkm() which will do the calculation for you.)

You can calculate VST values by first creating a DESeqDataSet and then using the varianceStabilizingTransformation() function on that object.

So you should first follow the vignette or workflow[1] example on creating a DESeqDataSet from matrix input:

http://www.bioconductor.org/help/workflows/rnaseqGene/

ADD COMMENT
0
Entering edit mode

thank you Michael i think i have to more study in R because everything is harsh for me to understand...

ADD REPLY
0
Entering edit mode
@andrewjskelton73-7074
Last seen 4 months ago
United Kingdom

hmmmm, different user, but very similar question.... 

I've provided an answer here

ADD COMMENT

Login before adding your answer.

Traffic: 464 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6