uu.seUppsala universitets publikationer
Ändra sökning
RefereraExporteraLänk till posten
Permanent länk

Direktlänk
Referera
Referensformat
  • apa
  • ieee
  • modern-language-association
  • vancouver
  • Annat format
Fler format
Språk
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Annat språk
Fler språk
Utmatningsformat
  • html
  • text
  • asciidoc
  • rtf
Full likelihood inference from the site frequency spectrum based on the optimal tree resolution
Uppsala universitet, Teknisk-naturvetenskapliga vetenskapsområdet, Matematisk-datavetenskapliga sektionen, Matematiska institutionen, Tillämpad matematik och statistik.
Ecole Polytech, CNRS, CMAP, Palaiseau, France.
2018 (Engelska)Ingår i: Theoretical Population Biology, ISSN 0040-5809, E-ISSN 1096-0325, Vol. 124, s. 1-15Artikel i tidskrift (Refereegranskat) Published
Abstract [en]

We develop a novel importance sampler to compute the full likelihood function of a demographic or structural scenario given the site frequency spectrum (SFS) at a locus free of intra-locus recombination. This sampler, instead of representing the hidden genealogy of a sample of individuals by a labelled binary tree, uses the minimal level of information about such a tree that is needed for the likelihood of the SFS and thus takes advantage of the huge reduction in the size of the state space that needs to be integrated. We assume that the population may have demographically changed and may be non-panmictically structured, as reflected by the branch lengths and the topology of the genealogical tree of the sample, respectively. We also assume that mutations conform to the infinitely-many-sites model. We achieve this by a controlled Markov process that generates 'particles' in the hidden space of SFS histories which are always compatible with the observed SFS. To produce the particles, we use Aldous' Beta-splitting model for a one parameter family of prior distributions over genealogical topologies or shapes (including that of the Kingman coalescent) and allow the branch lengths or epoch times to have a parametric family of priors specified by a model of demography (including exponential growth and bottleneck models). Assuming independence across unlinked loci, we can estimate the likelihood of a population scenario based on a large collection of independent SFS by an importance sampling scheme, using the (unconditional) distribution of the genealogies under this scenario when the latter is available. When it is not available, we instead compute the joint likelihood of the tree balance parameter beta assuming that the tree topology follows Aldous' Beta splitting model, and of the demographic scenario determining the distribution of the inter-coalescence times or epoch times in the genealogy of a sample, in order to at least distinguish different equivalence classes of population scenarios leading to different tree balances and epoch times. Simulation studies are conducted to demonstrate the capabilities of the approach with publicly available code.

Ort, förlag, år, upplaga, sidor
ACADEMIC PRESS INC ELSEVIER SCIENCE , 2018. Vol. 124, s. 1-15
Nyckelord [en]
Importance sampler, Semi-parametric estimation, Optimal tree resolution, Controlled Markov process on hidden genealogical trees
Nationell ämneskategori
Sannolikhetsteori och statistik
Identifikatorer
URN: urn:nbn:se:uu:diva-373021DOI: 10.1016/j.tpb.2018.07.002ISI: 000453111000001PubMedID: 30048667OAI: oai:DiVA.org:uu-373021DiVA, id: diva2:1277345
Tillgänglig från: 2019-01-10 Skapad: 2019-01-10 Senast uppdaterad: 2019-01-10Bibliografiskt granskad

Open Access i DiVA

Fulltext saknas i DiVA

Övriga länkar

Förlagets fulltextPubMed

Personposter BETA

Sainudiin, Raazesh

Sök vidare i DiVA

Av författaren/redaktören
Sainudiin, Raazesh
Av organisationen
Tillämpad matematik och statistik
I samma tidskrift
Theoretical Population Biology
Sannolikhetsteori och statistik

Sök vidare utanför DiVA

GoogleGoogle Scholar

doi
pubmed
urn-nbn

Altmetricpoäng

doi
pubmed
urn-nbn
Totalt: 3 träffar
RefereraExporteraLänk till posten
Permanent länk

Direktlänk
Referera
Referensformat
  • apa
  • ieee
  • modern-language-association
  • vancouver
  • Annat format
Fler format
Språk
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Annat språk
Fler språk
Utmatningsformat
  • html
  • text
  • asciidoc
  • rtf