Question: RNAseq tximport Deseq2
gravatar for tanyabioinfo
21 months ago by
tanyabioinfo20 wrote:



I am using the tximport module to process the output of Salmon. Then I am using the vst(dds) to normalise the data. I understand that the assay for dds is count. How can I process the abundance(TPM) using the vst so that my TPM data is normalised. 

My R code looks:

txi <- tximport(files, type="salmon", tx2gene=tx2gene, ignoreTxVersion=TRUE,dropInfReps=TRUE)
sampleTable <- data.frame(condition =samples$condition,time=factor(samples$time))
rownames(sampleTable) <- colnames(txi$counts)
dds <- DESeqDataSetFromTximport(txi, sampleTable, ~condition+time)

vsd <- vst(dds) 

deseq2 tximport • 385 views
ADD COMMENTlink modified 21 months ago by Michael Love25k • written 21 months ago by tanyabioinfo20

I missed this post because the tags were not separated, so didn’t get an email alert.

ADD REPLYlink written 21 months ago by Michael Love25k
Answer: RNAseq tximport Deseq2
gravatar for Michael Love
21 months ago by
Michael Love25k
United States
Michael Love25k wrote:

The VST is only designed for counts. We recommend using variance stabilized, transformed counts for downstream applications like visualization, distances, clustering or classification/prediction methods that benefit from homoskedastic data.

TPM is useful for comparing abundance across transcripts or genes. It should be proportional to the true abundance of RNA molecules.

ADD COMMENTlink written 21 months ago by Michael Love25k
Please log in to add an answer.


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