Skip to content
New issue

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

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

Already on GitHub? # to your account

Error in snpgdsVCF2GDS......Genotype code is out of range #95

Open
gaochengPRC opened this issue Jul 4, 2022 · 0 comments
Open

Error in snpgdsVCF2GDS......Genotype code is out of range #95

gaochengPRC opened this issue Jul 4, 2022 · 0 comments

Comments

@gaochengPRC
Copy link

Hi,

When I run the code below,

library("SNPRelate")
snpgdsVCF2GDS("xxxx.vcf", "ccm.gds", method="biallelic.only")
genofile <- openfn.gds("ccm.gds")
ccm_pca<-snpgdsPCA(genofile)

I get the following error.


Loading required package: gdsfmt
SNPRelate -- supported by Streaming SIMD Extensions 2 (SSE2)
Start file conversion from VCF to SNP GDS ...
Method: extracting biallelic SNPs
Number of samples: xxx
Parsing "xxxx.vcf" ...
Error in snpgdsVCF2GDS("xxxx.vcf", :
Genotype code is out of range "2".
FILE: xxxx.vcf
LINE: 46, COLUMN: 20, 2
Execution halted


I feel very confused. The (imputed) genotypes are all in the range [0, 2], why does it say that genotype code is out of range?

Thank you in advance!

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

No branches or pull requests

1 participant