Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Got a nan in the p-value results #273

Open
voichek opened this issue Apr 11, 2023 · 1 comment
Open

Got a nan in the p-value results #273

voichek opened this issue Apr 11, 2023 · 1 comment

Comments

@voichek
Copy link

voichek commented Apr 11, 2023

Describe the bug
When I run a GWAS I got one SNP that have a nan value in the three columns logl_H1, l_mle, and p_lrt

chr     rs      ps      n_miss  allele1 allele0 af      logl_H1 l_mle   p_lrt
1       1_5861128       5861128 62      C       T       0.161   nan     nan     nan

Gemma version - 0.98.5 (2021-08-25)
Command line run:
gemma-0.98.5-linux-static-AMD64 -bfile table -lmm 2 -k table.kinship -snps snps.txt -maf 0.05 -miss 0.5

To Reproduce
I have made a small example to reproduce the error, it has all the files needed to run it (attached).
The exact run command can be found in the log file of running the gemma, as well as all other exact parameters.
20230411_gemma_bug_report.zip

I couldn't find a previous issue like this in the GitHub issues page

Best wishes,
Yoav

@yefanglee
Copy link

Hi, did you address it? I met with same problem

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants