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: No valid variant remaining #368

Open
ryshi06 opened this issue Oct 23, 2024 · 0 comments
Open

Error: No valid variant remaining #368

ryshi06 opened this issue Oct 23, 2024 · 0 comments

Comments

@ryshi06
Copy link

ryshi06 commented Oct 23, 2024

Hello,

I am trying to run the program with a publicly available GWAS summary statistics and my PLINK binary file set with the following command:

Rscript $PRSice/PRSice.R --dir . \
    --prsice $PRSice/PRSice_linux \
    --base ${DIR_NAME}/${BASE_NAME}_updated.txt \
    --snp variant_alternate_id2 --chr chromosome --bp base_pair_location --A1 effect_allele --A2 other_allele --stat odds_ratio --pvalue p_value \
    --extract PRSice.valid \
    --target $TARGET_PLINK \
    --pheno $PHENO --pheno-col z_SUVr \
    --clump-kb 250 \
    --clump-r2 0.1 \
    --clump-p 5e-08 \
    --thread 1 \
    --beta \
    --binary-target F

But I received the error message saying:

Reading 100.00%
21101114 variant(s) observed in base file, with: 
21101114 variant(s) excluded based on user input 
0 total variant(s) included from base file 

Error: No valid variant remaining 

I checked the overlap variants of the two file and returns 17,081,910. I also checked the column format as required. The variant_alternate_id2 is in chr:pos:ref:alt format, and the PLINK file also follows this format.

Can you please guide me on how to solve this? Thank you!

Best,
Ruyu

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