DeSeq2 rlog or varianceStabilizingTransformation using ouput file
1
0
Entering edit mode
7.2 years ago
BM ▴ 70

I have a DeSeq2 output file generated by an external RNA-seq. So the output is the usual has the usual columns: gene identifiers, baseMean, log2FoldChange, lfcSE stat, pvalue, padj. However, I want to perform rlog or varianceStabilizingTransformation on my data. Is there a way to do this transformation just using the output file above? When I go back and perform the DeSeq2 analysis on my counts, there are differences between my output and the one given to me by the company. So I would like to use the output file.

RNA-Seq DeSeq2 Transformation rlog vsd • 3.7k views
ADD COMMENT
0
Entering edit mode

It all depends on the model that was used by the external collaborator and the design they have implemented while performing the DE analysis. You can simply ask them how they did and implement the same on your counts and see how are the end results comparable. Or ask them to perform the way you are performing the analysis to the counts (if that is possible). The motivation of the models depends on various factors and the design of the experimental matrix in the DESEQ2 design model. You cannot use an external collaborator output to trace back the counts data from DESeq2 since am guessing your counts data might not be exact integer and also depends on what rounding was applied before feeding the counts to DESeq2. Hope I was clear

ADD REPLY
0
Entering edit mode
7.2 years ago

No, that's not possible. The output with differentially expressed genes doesn't contain the information needed to get per sample normalized expression levels.

ADD COMMENT

Login before adding your answer.

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