How to compute DS (diploid dosage) from HDS (haploid dosage) in VCF?

Question

[Nadia]

Hello, I have a VCF file with GT:HDS (haploid dosage) values, but I need to generate GT:DS (diploid dosage), since this is the most common format in vcf files and used commonly for association analyses.

Is there any way, command or plugin to compute or convert the GT:DS field from GT:HDS?

The VCF FORMAT column is GT:HDS. A sample line looks like: chr22 10510303 ... GT:HDS 0|0:0,0 0|0:0.001,0 ...

VCF Header Snippet:

FORMAT=<ID=HDS,Number=2,Type=Float,Description="Allele Probs for two haplotypes=Haploid Alternate Allele Dosage">

Thank you in advance for your guidance.

[Jeremy Leipzig]

OK this is really a Minimac4 / Michigan Imputation Server concept, not a general VCF thing

Answer 1

[LChart]

for biallelic sites, the trivial way is DS = HDS[0] + HDS[1]