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).

thanks

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?

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.

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.

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

https://www.rdocumentation.org/packages/samr/versions/2.0/topics/SAMseq

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