Some measure of # of selective genes

It’s worth noting that those values are out of date/incomplete.

You can find many examples of genes labeled as LRT >= 100 on the web browser that are lower in that dataset. Sime genes are missing entirely.

What I did is simply re-run the data myself using the R code shared by @pmontgom on the most recent release. There’s an error I havent been able to pin down still where errors/warnings from the try() function in LRT_test() aren’t being properly caught, which allows for some bad skewed T fits on certain genes that are clearly not skewed to return a value instead of NA (HSPA5 and a few RNA polymerase genes). I was able to bandaid this with the “testit” R library and the has_warning() function, but it’s fairly computationally intensive.