Question: calculating logfold change from variance stabilizing transformation
gravatar for lirongrossmann
13 months ago by
lirongrossmann20 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 13 months ago by Michael Love20k • written 13 months ago by lirongrossmann20
gravatar for Michael Love
13 months ago by
Michael Love20k
United States
Michael Love20k 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 13 months ago by Michael Love20k

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 13 months ago by lirongrossmann20

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 13 months ago by Michael Love20k

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 13 months ago by lirongrossmann20

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 13 months ago by Michael Love20k
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: 227 users visited in the last hour