Learn R Programming

hierfstat (version 0.5-11)

fstat2dos: Converts a hierfstat genetic data frame to dosage data

Description

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

Usage

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

Value

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

Run this code
# NOT RUN {
dat<-sim.genot(nbal=5,nbloc=10)
dos<-fstat2dos(dat[,-1])
dim(dos) 
wc(dat)
fst.dosage(dos,pop=dat[,1])

# }

Run the code above in your browser using DataLab