Likelihood and Approximate Likelihood Analyses of Genetic Structure in a Linear Habitat: Performance and Robustness to Model Mis-Specification Export

@article{citeulike:2215659, abstract = {We evaluate the performance of maximum likelihood (ML) analysis of allele frequency data in a linear array of populations. The parameters are a mutation rate and either the dispersal rate in a stepping stone model or a dispersal rate and a scale parameter in a geometric dispersal model. An approximate procedure known as maximum product of approximate conditional (PAC) likelihood is found to perform as well as ML. Mis-specification biases may occur because the importance sampling algorithm is formally defined in term of mutation and migration rates scaled by the total size of the population, and this size may differ widely in the statistical model and in reality. As could be expected, ML generally performs well when the statistical model is correctly specified. Otherwise, mutation rate estimates are much closer to mutation probability scaled by number of demes in the statistical model than scaled by number of demes in reality when mutation probability is high and dispersal is most limited. This mis-specification bias actually has practical benefits. However, opposite results are found in opposite conditions. Migration rate estimates show roughly similar trends, but they may not always be easily interpreted as low-bias estimates of dispersal rate under any scaling. Estimation of the dispersal scale parameter is also affected by mis-specification of the number of demes, and the different biases compensate each other in such a way that good estimation of the so-called neighborhood size (or more precisely the product of population density and mean-squared parentoffspring dispersal distance) is achieved. Results congruent with these findings are found in an application to a damselfly data set. 10.1093/molbev/msm206}, author = {Rousset, Francois and Leblois, Raphael}, citeulike-article-id = {2215659}, citeulike-linkout-0 = {http://dx.doi.org/10.1093/molbev/msm206}, citeulike-linkout-1 = {http://mbe.oxfordjournals.org/cgi/content/abstract/24/12/2730}, citeulike-linkout-2 = {http://view.ncbi.nlm.nih.gov/pubmed/17893401}, citeulike-linkout-3 = {http://www.hubmed.org/display.cgi?uids=17893401}, day = {1}, doi = {10.1093/molbev/msm206}, journal = {Mol Biol Evol}, keywords = {geographic, modeling, simulation, to-get}, month = {December}, number = {12}, pages = {2730--2745}, posted-at = {2008-01-10 18:21:16}, priority = {2}, title = {Likelihood and Approximate Likelihood Analyses of Genetic Structure in a Linear Habitat: Performance and Robustness to Model Mis-Specification}, url = {http://dx.doi.org/10.1093/molbev/msm206}, volume = {24}, year = {2007} }

TY - JOUR ID - citeulike:2215659 L3 - citeulike-article-id:2215659 N2 - We evaluate the performance of maximum likelihood (ML) analysis of allele frequency data in a linear array of populations. The parameters are a mutation rate and either the dispersal rate in a stepping stone model or a dispersal rate and a scale parameter in a geometric dispersal model. An approximate procedure known as maximum product of approximate conditional (PAC) likelihood is found to perform as well as ML. Mis-specification biases may occur because the importance sampling algorithm is formally defined in term of mutation and migration rates scaled by the total size of the population, and this size may differ widely in the statistical model and in reality. As could be expected, ML generally performs well when the statistical model is correctly specified. Otherwise, mutation rate estimates are much closer to mutation probability scaled by number of demes in the statistical model than scaled by number of demes in reality when mutation probability is high and dispersal is most limited. This mis-specification bias actually has practical benefits. However, opposite results are found in opposite conditions. Migration rate estimates show roughly similar trends, but they may not always be easily interpreted as low-bias estimates of dispersal rate under any scaling. Estimation of the dispersal scale parameter is also affected by mis-specification of the number of demes, and the different biases compensate each other in such a way that good estimation of the so-called neighborhood size (or more precisely the product of population density and mean-squared parentoffspring dispersal distance) is achieved. Results congruent with these findings are found in an application to a damselfly data set. 10.1093/molbev/msm206 IS - 12 JF - Mol Biol Evol EP - 2745 TI - Likelihood and Approximate Likelihood Analyses of Genetic Structure in a Linear Habitat: Performance and Robustness to Model Mis-Specification VL - 24 SP - 2730 KW - geographic KW - modeling KW - simulation KW - to-get AU - Rousset, Francois AU - Leblois, Raphael PY - 2007/12/01/ UR - http://dx.doi.org/10.1093/molbev/msm206 ER -