The function read.fstat reads Fstat data files (.dat) and
convert them into a genind object.
Note: read.fstat is meant for DIPLOID DATA ONLY. Haploid
data with the Hierfstat format can be read into R using
read.table or read.csv after removing headers and 'POP'
lines, and then converted using df2genind.
Usage
read.fstat(file,missing=NA,quiet=FALSE)
Arguments
file
a character string giving the path to the file to
convert, with the appropriate extension.
missing
can be NA, 0 or "mean". See details section.
quiet
logical stating whether a conversion message must be
printed (TRUE,default) or not (FALSE).
Value
an object of the class genind
encoding
UTF-8
Details
There are 3 treatments for missing values:
- NA: kept as NA.
- 0: allelic frequencies are set to 0 on all alleles of the concerned
locus. Recommended for a PCA on compositionnal data.
- "mean": missing values are replaced by the mean frequency of the
corresponding allele, computed on the whole set of
individuals. Recommended for a centred PCA.
References
Fstat (version 2.9.3). Software by Jerome
Goudet. http://www2.unil.ch/popgen/softwares/fstat.htm