Question: calculating logfold change from variance stabilizing transformation
gravatar for lirongrossmann
7 months ago by
lirongrossmann10 wrote:

Hi All,

I am interested in calculating the logfold change from the transformed expression file obtained by using variance stabilizing transformation function. 

I used the transformed data to normalize my count data and compare two groups within the data using a different statistical test than the one in Deseq, thus I cannot use the logfold change given by the differential analysis of Deseq. 

Is there a straightforward conversion between the logfold change and the values obtained by variance stabilizing transformation?  (I tried taking the exp of the difference between the averages of the two groups, but it did not work).


ADD COMMENTlink modified 7 months ago by Michael Love18k • written 7 months ago by lirongrossmann10
gravatar for Michael Love
7 months ago by
Michael Love18k
United States
Michael Love18k wrote:

Transformed data is on the log2 scale for large counts. 

Note that we recommend use of DESeq, results, and lfcShrink for the optimal estimation of log fold changes.

ADD COMMENTlink written 7 months ago by Michael Love18k

My groups are small (9 vs 9 samples) and for some reason when I tried a different statistical test to compare them (Wilcoxon rank sum test, to be specific), I got a more uniform expression within each group for the genes that were significant. When I used Deseq to perform the statistical test, the genes that were found to be differentially expression had a large variance within each groups. 

I would really like to use Deseq for my statistical analysis as well, is there a way to overcome this problem?

ADD REPLYlink written 7 months ago by lirongrossmann10

I’d recommend using a package specifically designed for the data, as it will have known operating characteristics, whereas custom solutions are just less known quantities. So if you find that you prefer working with ranks for this dataset, I’d recommend using SAMseq.

ADD REPLYlink written 7 months ago by Michael Love18k

Thank you so much! 

I see that SamSeq requires normalized data to run. From your experience, is it “safe” to give SamSeq the output of variance stabilizing transformation from Deseq2 as the input for the analysis.

ADD REPLYlink written 7 months ago by lirongrossmann10

You should provide SAMseq with the raw counts as 'x':

It uses its own resampling strategy in lieu of parametric normalization. It works well.

ADD REPLYlink written 7 months ago by Michael Love18k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 106 users visited in the last hour