How to imports DS field from VCF files with SeqArray seqVCF2GDS()
Entering edit mode
ealopera • 0
Last seen 24 months ago

Hello, I have VCF file coming from imputed data that contains three genotype info in three formats GT:DS:GP. Specifically, my dosage (DS) field is composed of CONTINUOUS values from 0 to 2. I use the function to convert this file to GDS., file.out,scenario="imputation")

and later I check the dosage field in the GDS with this command

seqGetData(genofile, "$dosage")

and I notice the dosage field contains only discrete values 0, 1 or 2, but not continuous. Which gives me the feeling that DS filed is not being imported but rather re-calculated. How can I change this?

seqArray seqVCF2GDS • 698 views
Entering edit mode
Last seen 11 days ago
University of Washington

seqGetData(genofile, "$dosage") is indeed calculating dosage based on the "genotype" node. To read from the dosage node, you can use seqGetData(genofile, "annotation/format/DS"). Or you can use SeqVarTools::imputedDosage(), which does the same thing.


Login before adding your answer.

Traffic: 197 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6