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

Trying to get LOCO=True to work #440

Open
Ankeet3 opened this issue Sep 26, 2023 · 0 comments
Open

Trying to get LOCO=True to work #440

Ankeet3 opened this issue Sep 26, 2023 · 0 comments

Comments

@Ankeet3
Copy link

Ankeet3 commented Sep 26, 2023

Hi,
So I've been using SAIGE for a while, and while performing the GBAT (gene based association tests), I would always set LOCO=False, because setting LOCO=True would cause some sort of parsing issue where the script would stdout "chunk is empty". For context, my genotypic data is in the PLINK format. Recently, based on some research I read, it seems like when working with a very large sample space which has high chances of having a significant amount of related individuals, having LOCO set to false would lead to a reduction in power. But how am I supposed to have LOCO set to true if it isn't able to parse the input data appropriately? Does LOCO=True only work when the input genotypic data is of VCF/BGEN format? What even is the point of having the option for LOCO if it doesn't even seem to work for the association tests? PS - I have set the chr=1 flag as well

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

1 participant