biall2dos function

Converts bi-allelic SNPs from hierfstat format to dosage format

Converts bi-allelic SNPs from hierfstat format to dosage format

Converts bi-allelic SNPs hierfstat format to dosage format, the number of alternate allele copies at a locus for an individual, i.e. 11 -> 0; 12 or 21 >1 and 22 ->2

biall2dos(dat,diploid=TRUE)

Arguments

  • dat: a hierfstat data frame without the first column (the population identifier), individuals in rows, columns with individual genotypes encoded as 11, 12, 21 and 22
  • diploid: whether the data set is from a diploid organism

Returns

a matrix containing allelic dosages

Examples

## Not run: biall2dos(sim.genot(nbal=2,nbloc=10)[,-1]) # a 10 column matrix ## End(Not run)