Dear Bioconductors, 

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). Does transforming the values 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 no longer meaningful representations of the data?

I'd appreciate any advice you can offer, thank you.