fstat2dos function

Converts a hierfstat genetic data frame to dosage data

Converts a hierfstat genetic data frame to dosage. For each allele at each locus, allelic dosage (number of copies of the allele) is reported. The column name is the allele identifier

fstat2dos(dat,diploid=TRUE)

Arguments

• dat: data frame with genetic data without the first column (population identifier)
• diploid: whether the data set is from a diploid organism

Returns

a matrix with $\sum_l n_l^a$ columns (where $n_l^a$ is the number of alleles at locus l), as many rows as individuals, and containing the number of copies (dosage) of the corresponding allele

Examples

## Not run:

dat<-sim.genot(nbal=5,nbloc=10)
dos<-fstat2dos(dat[,-1])
dim(dos) 
wc(dat)
fst.dosage(dos,pop=dat[,1])
## End(Not run)