Rlog performed on tpm or raw transcript abundances?
2
0
Entering edit mode
@94e8e072
Last seen 18 months ago
United States

Background: Following DESeq2 tutorials, I'm importing Kallisto abundance.tsv files using tximport, creating a DESeq2 dataset, and rlog transforming the data for downstream figure generation. Question: What column of Kallisto abundance.tsv is transformed using "rlog(dds)"? In other words, is DESeq2 "rlog" acting on Transcript per million (tpm) or raw abundance values (est_counts) from the Kallisto abundance.tsv files?

txi<- tximport(files, type="kallisto", tx2gene = tx2gene)
dds <-DESeqDataSetFromTximport(txi, colData=samples, design = ~Batch + Condition)
dds <- DESeq(dds)
rld <- rlog(dds, blind=FALSE)

#not relevant to the question, simply showing an example of using rlog transformed data in figures
plotPCA(rld, intgroup=c("Condition", "Batch"))
DESeq2 rlog • 1.0k views
ADD COMMENT
2
Entering edit mode

Abundance. If you look at the counts in dds, they should match abundance, not TPM (unless there is a slight shift due to the transcript-length offset value)

ADD REPLY

Login before adding your answer.

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