It looks like the GSVA enrichment scores are bound between -1 and 1. Is this true? If so, is this ensured in step 3?
Lastly, though I set the scores to be distributed as the difference between the largest positive and negative deviations, the output was not normally distributed (potentially because I have a relatively small number of samples). Thus, I was planning to transform the values using Fisher transformation. Does this seem like an appropriate follow-up approach to make the scores more normally distributed so that I can perform a Two-Way ANOVA test or would this over-transform the scores to point that they are perhaps less meaningful representations of the data? On a side note, this transformation only worked for one of my expression sets. Are there any other transformations you'd recommend for enrichment scores?
I appreciate any advice you can offer.