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

Error running SAIGE using unimputed data #130

Open
zrayw opened this issue Dec 19, 2023 · 1 comment
Open

Error running SAIGE using unimputed data #130

zrayw opened this issue Dec 19, 2023 · 1 comment

Comments

@zrayw
Copy link

zrayw commented Dec 19, 2023

Dear Developer,

We were trying to run an older SAIGE (0.45) on unimputed data. In SAIGE step 2, We imposed IsDropMissingDosages=TRUE since I don't have any imputed genotypes. We met the following error:
Error in solve.default(XVX) :
system is computationally singular: reciprocal condition number = 6.31709e-20
Calls: SPAGMMATtest ... ScoreTest_NULL_Model -> solve -> solve -> solve.default

Would it be the reason that so many genotypes were discarded? Do you have any suggestions running SAIGE on unimputed data? Btw, we could successfully run SAIGE with sDropMissingDosages=FALSE, but in that case SAIGE would impute missingness with mean values and cause bias.

Thanks so much!!

@mattimpat
Copy link

It seems the problem is caused by SAIGE calculating an equivalent MAF for the MAC, based on the total number of samples. However, this does not work when there are many missing genotypes, since the equivalent MAF threshold can be met, even when the MAC one is not. It can be fixed by checking the MAC threshold for each variant, rather than the equivalent MAF.

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