This script performs an AMOVA based on the genetic distance matrix from stamppNeisD() [package StAMPP] using the amova() function from the package PEGAS for exploring within and between population variation. For detailed information use their help pages: ?pegas::amova, ?StAMPP::stamppAmova. Be aware due to a conflict of the amova functions from various packages I had to 'hack' StAMPP::stamppAmova to avoid a namespace conflict.
Arguments
- x
Name of the genlight containing the SNP genotypes, with population information [required].
- distance
Distance matrix between individuals (if not provided NeisD from StAMPP::stamppNeisD is calculated) [default NULL].
- permutations
Number of permutations to perform for hypothesis testing [default 100]. Please note should be set to 1000 for analysis.
- verbose
Verbosity: 0, silent or fatal errors; 1, begin and end; 2, progress log ; 3, progress and results summary; 5, full report [default 2, unless specified using gl.set.verbosity].
Value
An object of class 'amova' which is a list with a table of sums of square deviations (SSD), mean square deviations (MSD), and the number of degrees of freedom, and a vector of variance components.
Author
Bernd Gruber (bugs? Post to https://groups.google.com/d/forum/dartr)