DeSeq Strategies for run to run Replication
1
0
Entering edit mode
Giorgio • 0
@b0635dc4
Last seen 4 weeks ago
United States

Hi all,

This is more a general (philosophical?) question:

1. Say I have a dataset analyzed with DeSeq default, and a vst normalized dataset is obtained.

2. The VST is then used to generate a final predictive ML model (binary) with only few genes from the VST dataset.

3. Now I have a few samples that I re-rerun end-to-end with the exact same pipeline as above.

4. After importing in DeSeq what normalization would you think to use in order to obtain a close VST from the original run?

5. The goal is to predict correctly the new repeated samples.

I know there are millions of variables in play, but was curious to see what the folks would answer.

Thank you in advance

DESeq2 Normalization Clustering replicate • 94 views
0
Entering edit mode
@mikelove
Last seen 8 hours ago
United States

You can use the same VST on the new samples. We have some help to do this in ?varianceStabilizingTransformation:

The variance stabilizing transformation from a previous dataset can be "frozen" and reapplied to new samples. The frozen VST is accomplished by saving the dispersion function accessible with dispersionFunction, assigning this to the DESeqDataSet with the new samples, and running varianceStabilizingTransformation with 'blind' set to FALSE. Then the dispersion function from the previous dataset will be used to transform the new sample(s).