powered by
Convert genotype calls in the form of witten out alleles (e.g. 'A/T') to dosages (0, 1, 2) based on provided risk alleles from a PGS.
convert.alleles.to.pgs.dosage(called.alleles, risk.alleles)
A vector of dosages corresponding to each genotype in called.alleles. Hemizygous genotypes (one allele e.g. 'A') are counted as 1.
A vector of genotypes in allelic notation separated by a slash or pipe.
A vector of risk alleles from a polygenic score corresponding to each genotype (by locus) in called.alleles.
called.alleles <- c('A/A', 'A/T', 'T/T'); risk.alleles <- c('T', 'T', 'T'); convert.alleles.to.pgs.dosage(called.alleles, risk.alleles);
Run the code above in your browser using DataLab