Vector of allele counts: each element either 0,1,2 or NA.
p
Matrix of allele frequencies. Each row corresponds with a SNP. Number of rows must equal length of data. Each column is a population
nchains
Number of Markov Chains to perform gibbs sampling
iterations
Number of samples for each Markov Chain
alpha
1 minus the width of credible interval taken around the posterior mode.
prior
Concentration parameter for the dirichlet distribution. 1 is uninformative, small values suggest first-generation admixture (uninformative as to which source populations involved), and higher values suggest low parental divergence.
burn
Burn- in period. The number of MCMC samples to discard.
SampSizes
Vector of the same length as the number of source populations. Each element is the number of individuals used to calculate allele frequencies in that population. Parameterises the prior on allele frequencies
Details
BEAPFrOG requires jags and rjags to be installed. http://mcmc-jags.sourceforge.net/
Credible intervals are centered around the mode.
Values for parameter vectors m1 and m2 can be exchanged to give identical likelihoods. This may give rise to bimodal posterior distributions, particularly with first-generation admixture, and the resulting credible intervals are not useful. Therefore, for all MCMC samples, we redefine m1 as the admixture proportions for the parent with admixture from population 1 less than 0.5, and vice-versa for m2.
Value
A list including elements
P1est
A vector: Posterior mode for admixture proportions in parent 'A'.
P2
A vector: Posterior mode for admixture proportions in parent 'B'.
P1i
A matrix of two columns and number of rows equal to number of populations: upper and lower credible intervals of width defined by the argument alpha, for admixture proportions in parent 'A'.
P2i
A matrix of two columns and number of rows equal to number of populations: upper and lower credible intervals of width defined by the argument alpha, for admixture proportions in parent 'B'.
Monitor
All stored MCMC samples. Use plot() on this object for visualisation of posterior distribution and MCMC trace
Author(s)
Daniel Crouch & Michael Weale, Department of Medical and Molecular Genetics, King's College London