bootstrapHet_byloc: Calculate confidence intervals on observed and expected heterozygosity in each population by bootstrapping markers

Description

This function uses a Arlequin project file (.arp) as input. It resamples loci with replacement (i.e. bootstrap resampling) and calculates the corresponding confidence intervals on Ho (observed heterozygosity) and He (expected heterozygosity) in each population.

Usage

bootstrapHet_byloc(datafile, ndigit = 3, nrepet = 1000)

Arguments

datafile

Name of the Arlequin file to use as input, with the extension (e.g. 'glyptemys3.arp')

ndigit

Number of digits per allele in input file. Can be 2 or 3 (default is 3)

nrepet

Number of bootstrap replicates to perform. Default is 1000. Can be slow for large datafiles.

Value

Results are printed to a file saved in the working directory, with the suffix 'bootHet_byloc' in the filename. First part of the table reports real Ho and He values for each locus in each population. Bottom part reports bootstrap results for each population: mean calculated Ho and He as well as the lower and upper bounds of their 95 % confidence intervals.

Details

Some .arp files might not be imported correctly depending on which software was used to generate them. If this is the case, start from the Genepop file and convert it first with convert. It is possible to perform bootstrapping on individuals instead of loci. To do so, use the function bootstrapHet instead.

References