scaledTPM for DESeq2 using tximport
1
0
Entering edit mode
@ffcc7fa5
Last seen 11 months ago
India

Hello!

I am performing a transcript level differential expression analysis. For the same, I have quantified my samples using salmon and am using tximport as an offset before DESeq2. I wanted to know whether passing scaledTPM values is recommended to DESeqDataSetFromTximport? I am attaching my code below. Any help is appreciated. Thanks!!


txi <- tximport(files, type= "salmon", txOut=TRUE, countsFromAbundance="scaledTPM")

sampleTable <- data.frame(condition = factor(samples$condition))
rownames(sampleTable) <- colnames(txi$counts)
dds <- DESeqDataSetFromTximport(txi, sampleTable, ~condition)
DESeq2 tximport • 638 views
ADD COMMENT
2
Entering edit mode
@mikelove
Last seen 1 hour ago
United States

Either is fine (countsFromAbundance (cFA) = "no" or the other options which scale TPM to act like counts):

https://bioconductor.org/packages/release/bioc/vignettes/tximport/inst/doc/tximport.html#Downstream_DGE_in_Bioconductor

I personally recommend cFA = "no" when using DESeq2 or edgeR.

ADD COMMENT

Login before adding your answer.

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