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. seqVCF2GDSfile.in, 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?

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.