Coal.Theta.Sim.File.G: Coal.Theta.Sim.File.G: function to simulate a vector (k.vec)...

In radamsRHA/ThetaMater: ThetaMater: Bayesian estimation of the population size parameter theta from genomic data

Coal.Theta.Sim.File.G: function to simulate a vector (k.vec) of mutation counts given a vector of locus lengths (l.vec), vector of sample numbers (n.vec) with gamma rate variation

Description

This function returns k.vec, l.vec, n.vec from simulating under the coalescent model, with the results written to 'out.file'

Usage

Coal.Theta.Sim.File.G(theta, n.vec, l.vec, num.loci, out.file, alpha.param, K)

Arguments

theta	Population size parameter theta
num.loci	Number of loci to simulate
out.file	Outfile to write in alleles format
alpha.param	Shape parameter of the gamma distribution for among-locus rate variation
K	Number of discrete classes to approximate the gamma distribution
n.list	List containing the number of samples for each locus
l.list	List containing the number of samples for each locus

radamsRHA/ThetaMater documentation built on April 25, 2024, 10:11 a.m.