Non-equilibrium allele frequency spectra via spectral methods
S. Lukić, J. Hey, and K. Chen. Theoretical Population Biology, 79 (4):
203 - 219(2011)
Abstract
A major challenge in the analysis of population genomics data consists of isolating signatures of natural selection from background noise caused by random drift and gene flow. Analyses of massive amounts of data from many related populations require high-performance algorithms to determine the likelihood of different demographic scenarios that could have shaped the observed neutral single nucleotide polymorphism (SNP) allele frequency spectrum. In many areas of applied mathematics, Fourier Transforms and Spectral Methods are firmly established tools to analyze spectra of signals and model their dynamics as solutions of certain Partial Differential Equations (PDEs). When spectral methods are applicable, they have excellent error properties and are the fastest possible in high dimension; see Press et al. (2007). In this paper we present an explicit numerical solution, using spectral methods, to the forward Kolmogorov equations for a Wright–Fisher process with migration of K populations, influx of mutations, and multiple population splitting events.
Description
ScienceDirect.com - Theoretical Population Biology - Non-equilibrium allele frequency spectra via spectral methods
%0 Journal Article
%1 luki2011nonequilibrium
%A Lukić, Sergio
%A Hey, Jody
%A Chen, Kevin
%D 2011
%J Theoretical Population Biology
%K diffusion_approximation population_structure selection site_frequency_spectrum special_functions spectral_theory
%N 4
%P 203 - 219
%T Non-equilibrium allele frequency spectra via spectral methods
%U http://www.sciencedirect.com/science/article/pii/S0040580911000165
%V 79
%X A major challenge in the analysis of population genomics data consists of isolating signatures of natural selection from background noise caused by random drift and gene flow. Analyses of massive amounts of data from many related populations require high-performance algorithms to determine the likelihood of different demographic scenarios that could have shaped the observed neutral single nucleotide polymorphism (SNP) allele frequency spectrum. In many areas of applied mathematics, Fourier Transforms and Spectral Methods are firmly established tools to analyze spectra of signals and model their dynamics as solutions of certain Partial Differential Equations (PDEs). When spectral methods are applicable, they have excellent error properties and are the fastest possible in high dimension; see Press et al. (2007). In this paper we present an explicit numerical solution, using spectral methods, to the forward Kolmogorov equations for a Wright–Fisher process with migration of K populations, influx of mutations, and multiple population splitting events.
@article{luki2011nonequilibrium,
abstract = {A major challenge in the analysis of population genomics data consists of isolating signatures of natural selection from background noise caused by random drift and gene flow. Analyses of massive amounts of data from many related populations require high-performance algorithms to determine the likelihood of different demographic scenarios that could have shaped the observed neutral single nucleotide polymorphism (SNP) allele frequency spectrum. In many areas of applied mathematics, Fourier Transforms and Spectral Methods are firmly established tools to analyze spectra of signals and model their dynamics as solutions of certain Partial Differential Equations (PDEs). When spectral methods are applicable, they have excellent error properties and are the fastest possible in high dimension; see Press et al. (2007). In this paper we present an explicit numerical solution, using spectral methods, to the forward Kolmogorov equations for a Wright–Fisher process with migration of K populations, influx of mutations, and multiple population splitting events.},
added-at = {2013-02-26T19:41:52.000+0100},
author = {Lukić, Sergio and Hey, Jody and Chen, Kevin},
biburl = {https://www.bibsonomy.org/bibtex/2dd6c748d78accf02e9e203206476b712/peter.ralph},
description = {ScienceDirect.com - Theoretical Population Biology - Non-equilibrium allele frequency spectra via spectral methods},
interhash = {c0401fabe2de931dfb0de6b17cb40ff1},
intrahash = {dd6c748d78accf02e9e203206476b712},
journal = {Theoretical Population Biology},
keywords = {diffusion_approximation population_structure selection site_frequency_spectrum special_functions spectral_theory},
number = 4,
pages = {203 - 219},
timestamp = {2013-02-26T19:41:52.000+0100},
title = {Non-equilibrium allele frequency spectra via spectral methods},
url = {http://www.sciencedirect.com/science/article/pii/S0040580911000165},
volume = 79,
year = 2011
}