TN93: The TN93 class

Description

This class implements the Tamura-Nei 93 GTR-submodel.

The rate parameters are the following: "Alpha1", "Alpha2","Beta". Package: Class TN93

Directly known subclasses:

public static class TN93 extends GTR

Usage

TN93(name="Anonymous", rate.params=list(Alpha1 = 1, Alpha2 = 1, Beta = 1),
  base.freqs=c(0.25, 0.25, 0.25, 0.25), ...)

Arguments

name

Object name.

rate.params

Rate parameters.

base.freqs

Base frequency parameters.

...

Not used.

Fields and Methods

Methods:

	`checkConsistency`	-
	`getBaseFreqs`	-
	`getRateParam`	-
	`getRateParamList`	-
	`setBaseFreqs`	-
	`setRateParam`	-
	`setRateParamList`	-
	`summary`	-

Methods inherited from GTR: checkConsistency, getBaseFreqs, getRateParam, getRateParamList, setBaseFreqs, setRateParam, setRateParamList, summary

Methods inherited from UNREST: checkConsistency, summary

Methods inherited from GeneralSubstitution: as.character, checkConsistency, clone, getAlphabet, getEquDist, getEventRate, getEventRateAtSite, getEventsAtSite, getQMatrix, getRate, getRateList, hasUndefinedRate, is, plot, rescaleQMatrix, sampleState, setAlphabet, setEquDist, setQMatrix, setRate, setRateList, summary

Methods inherited from Process: !=, ==, as.character, checkConsistency, clone, getAlphabet, getEventsAtSite, getId, getName, getParameterAtSite, getSiteSpecificParamIds, getSiteSpecificParamList, getWriteProtected, hasSiteSpecificParameter, hasUndefinedRate, is, setAlphabet, setId, setName, setParameterAtSite, setSiteSpecificParamIds, setSiteSpecificParamList, setWriteProtected, summary

Methods inherited from PSRoot: checkConsistency, enableVirtual, getComments, getMethodsList, globalConsistencyCheck, intersect.list, is, is.na, ll, my.all.equal, plot, setComments, setMethodsList, summary, virtualAssignmentForbidden

Methods inherited from Object: $, $<-, [[, [[<-, as.character, attach, attachLocally, clearCache, clearLookupCache, clone, detach, equals, extend, finalize, getEnvironment, getFieldModifier, getFieldModifiers, getFields, getInstantiationTime, getStaticInstance, hasField, hashCode, ll, load, names, objectSize, print, save

References

Tamura, K, and Nei, M (1993) Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees - Molecular Biology and Evolution 10:512-526 http://bit.ly/bNkCqn

Examples

Run this code

# NOT RUN {
	# create substitution process object
       p<-TN93(rate.params=list( "Alpha1"=4,"Alpha2"=3,"Beta"=2),
		base.freqs=c(2,2,1,3)/9
	)
       # get a summary
       summary(p)
	# display a bubble plot
	plot(p)

	# The following code demonstrates how to use
	# the process in a simulation.

	# create a sequence, attach process p
	s<-NucleotideSequence(length=20,processes=list(list(p)))
	# sample states
	sampleStates(s)
	# make the first five positions invariable
	setRateMultipliers(s,p,0,1:5)
	# get rate multipliers
	getRateMultipliers(s,p)
	# create a simulation object
	sim<-PhyloSim(root.seq=s,phylo=rcoal(2))
	# run simulation
	Simulate(sim)
	# print alignment
	sim$alignment
 
# }