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

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

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

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(theta, n.vec, l.vec, num.loci, out.file)

Arguments

theta	Population size parameter theta
num.loci	Number of loci to simulate
out.file	Outfile to write in alleles format
n.list	List containing the number of samples for each locus
l.list	List containing the number of samples for each locus

Examples

example.data <- Coal.Theta.Sim.File(theta = 0.017, n.vec = c(5,7), l.vec = c(1,1), num.loci = 100, out.file = '~/Desktop/example.alleles')

Coal.Theta.Sim.File Returns:
k.vec = # of mutations for each unique data pattern
n.vec = # of samples for each unique data pattern
l.vec = # of sites for each unique data pattern
c.vec = # of counts of each unique data pattern

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