Documentation for SLiM function calcWattersonsTheta, which is a method of
the SLiM class SLiMBuiltin.
Note that the R function is a stub, it does not do anything in R (except bring
up this documentation). It will only do
anything useful when used inside a slim_block function further
nested in a slim_script
function call, where it will be translated into valid SLiM code as part of a
full SLiM script.
calcWattersonsTheta(genomes, muts, start, end)An object of type Genome object. See details for description.
An object of type null or Mutation object. The default value is
NULL. See details for description.
An object of type null or integer. Must be of length 1 (a
singleton). The default value is NULL. See details for description.
An object of type null or integer. Must be of length 1 (a singleton).
The default value is NULL. See details for description.
An object of type float. Return will be of length 1 (a singleton)
Documentation for this function can be found in the official SLiM manual: page 753.
Calculates Watterson's theta (a metric of genetic diversity comparable to heterozygosity) for a vector of genomes, based upon the mutations in the genomes. Often genomes will be all of the genomes in a subpopulation, or in the entire population, but any genome vector may be used. By default, with muts=NULL, the calculation is based upon all mutations in the simulation; the calculation can instead be based upon a subset of mutations, such as mutations of a specific mutation type, by passing the desired vector of mutations for muts. The calculation can be narrowed to apply to only a window - a subrange of the full chromosome - by passing the interval bounds [start, end] for the desired window. In this case, the vector of mutations used for the calculation will be subset to include only mutations within the specified window. The default behavior, with start and end of NULL, provides the genome-wide Watterson's theta. The implementation of calcWattersonsTheta(), viewable with functionSource(), treats every mutation as independent in the heterozygosity calculations. One could regard this choice as embodying an infinite-sites interpretation of the segregating mutations, as with calcHeterozygosity(). In most biologically realistic models, such genetic states will be quite rare, and so the impact of this assumption will be negligible; however, in some models this distinction may be important. See calcPairHeterozygosity() for further discussion.
This is documentation for a function in the SLiM software, and has been reproduced from the official manual, which can be found here: http://benhaller.com/slim/SLiM_Manual.pdf. This documentation is Copyright © 2016-2020 Philipp Messer. All rights reserved. More information about SLiM can be found on the official website: https://messerlab.org/slim/
Other SLiMBuiltin:
SB,
calcFST(),
calcHeterozygosity(),
calcInbreedingLoad(),
calcPairHeterozygosity(),
calcVA(),
codonsToAminoAcids(),
mm16To256(),
mmJukesCantor(),
mmKimura(),
nucleotideCounts(),
nucleotideFrequencies(),
nucleotidesToCodons(),
summarizeIndividuals(),
treeSeqMetadata()
## This just brings up the documentation:
calcWattersonsTheta()