Hello Sara and Amy,

I am trying to use the tool with one transcriptome that presents ~16.000 genes. I decided to filter out on the basis of total amount of reads, keeping ~1100 genes.

From these, I have multiple reference sets that I wanted to test out to see how the behaviour changes for an specific gene that I know how it covaries with the iron concentration.

To speed up everything I selected only 400 genes as a testing ground to understand what I was doing, and this took some time.

Then I selected these 1100 genes and I ran the program with multiple denominators. What stroke me is that having 1100 genes, the call was way faster. Does this mean that above a certain threshold the computation is entirely different?

Additionally, the final coefficient was quite different from both calculations. In here I am displaying the same gene and the coefficient variation:

Let me know what you think, thanks a ton!