ctlab / fgsea

Fast Gene Set Enrichment Analysis
Other
379 stars 67 forks source link

NA in pval column #97

Closed pstew closed 3 years ago

pstew commented 3 years ago

Hello, I am running fgsea on some RNA-seq data with parameters minSize = 15, maxSize = 500, scoreType = "std", and eps = 1E-50. I am getting some NA's in the pval results column. These gene sets are a little on the large side, but they are under my maxSize. They seem to be the only ones with this issue. Below is the data sorted so that you can see the NA's. Could you help me understand what the issue might be? Thank you.

Screen Shot 2021-07-06 at 1 26 29 PM
assaron commented 3 years ago

Hi, there should be a warning message, describing why you get NAs and probably even a way fix the problem. Likely, your gene-level statisitic has a skewed distribution.

On Tue, Jul 6, 2021, 20:27 Paul Stewart @.***> wrote:

Hello, I am running fgsea on some RNA-seq data with parameters minSize = 15, maxSize = 500, scoreType = "std", and eps = 1E-50. I am getting some NA's in the pval results column. These gene sets are a little on the large side, but they are under my maxSize. They seem to be the only ones with this issue. Below is the data sorted so that you can see the NA's. Could you help me understand what the issue might be? Thank you.

[image: Screen Shot 2021-07-06 at 1 26 29 PM] https://user-images.githubusercontent.com/5620536/124642392-cb5b1180-de5d-11eb-8dc3-d462221816dd.png

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/ctlab/fgsea/issues/97, or unsubscribe https://github.com/notifications/unsubscribe-auth/AAH56WPD6GHVGB7LIYMA56DTWM4INANCNFSM475AHDGQ .

pstew commented 3 years ago

Ahhh, my knitr settings were suppressing the warning messages. I see it now. Thank you.