RNAseq tximport Deseq2
1
0
Entering edit mode
tanyabioinfo ▴ 20
@tanyabioinfo-14091
Last seen 5.2 years ago

 

HI

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 • 1.0k views
ADD COMMENT
0
Entering edit mode

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

ADD REPLY
1
Entering edit mode
@mikelove
Last seen 18 hours ago
United States

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 COMMENT

Login before adding your answer.

Traffic: 1048 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