Interpreting the flock algorithm from a statistical perspective.

Autor: Anderson EC; Fisheries Ecology Division, Southwest Fisheries Science Center, National Marine Fisheries Service, NOAA, 110 Shaffer Road, Santa Cruz, CA, 95060, USA., Barry PD; School of Fisheries and Ocean Science, University of Alaska, Fairbanks, 17101 Point Lena Loop Road, Juneau, AK, 99801, USA.
Jazyk: angličtina
Zdroj: Molecular ecology resources [Mol Ecol Resour] 2015 Sep; Vol. 15 (5), pp. 1020-30. Date of Electronic Publication: 2015 May 06.
DOI: 10.1111/1755-0998.12417
Abstrakt: We show that the algorithm in the program flock (Duchesne & Turgeon 2009) can be interpreted as an estimation procedure based on a model essentially identical to the structure (Pritchard et al. 2000) model with no admixture and without correlated allele frequency priors. Rather than using MCMC, the flock algorithm searches for the maximum a posteriori estimate of this structure model via a simulated annealing algorithm with a rapid cooling schedule (namely, the exponent on the objective function →∞). We demonstrate the similarities between the two programs in a two-step approach. First, to enable rapid batch processing of many simulated data sets, we modified the source code of structure to use the flock algorithm, producing the program flockture. With simulated data, we confirmed that results obtained with flock and flockture are very similar (though flockture is some 200 times faster). Second, we simulated multiple large data sets under varying levels of population differentiation for both microsatellite and SNP genotypes. We analysed them with flockture and structure and assessed each program on its ability to cluster individuals to their correct subpopulation. We show that flockture yields results similar to structure albeit with greater variability from run to run. flockture did perform better than structure when genotypes were composed of SNPs and differentiation was moderate (FST= 0.022-0.032). When differentiation was low, structure outperformed flockture for both marker types. On large data sets like those we simulated, it appears that flock's reliance on inference rules regarding its 'plateau record' is not helpful. Interpreting flock's algorithm as a special case of the model in structure should aid in understanding the program's output and behaviour.
(Published 2015. This article is a U.S. Government work and is in the public domain in the USA.)
Databáze: MEDLINE
Externí odkaz: