Rarefac

This programme (RAREFAC) is written in Turbo-Pascal and relates to the following paper :[br/] Petit, El Mousadik & Pons. 1998. Identifying populations for conservation on the basis of genetic markers. Conservation Biology 12, 844-855. It is a simplified version of the programme CONTRIB providing just the option of rarefaction and the partitioning of allelic richness within and among populations. This is contrasted with equivalent measurements for H, the expected heterozygosity.[br/] It can be used in conjunction with the program haplodiv based on the paper by Pons & Petit 1995, TAG 90, 462-470, which will provide standard errors for the diversity and differentiation parameters.[br/] The input file is a text file, where the first line indicates the number of haplotypes (limited to 100), then the number of populations (limited to 50), and finally the rarefaction size (it should not be larger than the smallest population sample size.[br/] Then follows the data for each population (line), with the number of each haplotype in each population (don't use relative frequencies). Example:[br/] 18 4 10[br/] 1 0 1 0 0 0 1 1 ...(18 columns)[br/] 0 1 2 1 1 0 13 0 ...[br/] 0 0 8 0 0 3 6 0 ...[br/] 1 0 9 0 0 3 7 1[br/] Results can be seen in the output file. General measures are given first: Within population diversity (HS), total diversity (HT), and GST are given, followed by similar measures based on allelic richness. Then you get the results for each population: H, its standard error, and allelic richness after rarefaction (minus one: a monomorphic population has one allele which has to be removed to compute the partitioning of diversity within and among populations).[br/] The program is written for an haploid gene but may be used for nuclear genes, assuming Hardy-Weinberg equilibrium. How to proceed when there are several loci? Do not take the mean across Gst or Rst. They are ratios, so you should take the mean of the numerator and the mean of the denominator separately. Then compute the ratio of the two means.