Package: Class ECMrest
Object
~~|
~~+--
PSRoot
~~~~~~~|
~~~~~~~+--
Process
~~~~~~~~~~~~|
~~~~~~~~~~~~+--
GeneralSubstitution
~~~~~~~~~~~~~~~~~|
~~~~~~~~~~~~~~~~~+--
CodonSubst
~~~~~~~~~~~~~~~~~~~~~~|
~~~~~~~~~~~~~~~~~~~~~~+--
ECMrest
Directly known subclasses:
public static class ECMrest extends CodonSubst
ECMrest(equ.dist=NA, ...)
Equilibrium distribution.
Not used.
Methods: No methods defined.
Methods inherited from CodonSubst: buildFromPAML, checkConsistency, newMatrix, setEquDist, 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
Kosiol C., Holmes I., Goldman, N. (2007) An empirical codon model for protein sequence evolution - Mol Biol Evol. 24(7):1464-79 DOI: 10.1093/molbev/msm064 http://bit.ly/1ia8gWm
CodonSubst GeneralSubstitution CodonUNREST
# NOT RUN {
# create substitution model object
p<-ECMrest()
# get object 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<-CodonSequence(length=10,processes=list(list(p)) )
# sample states
sampleStates(s)
# make the first three positions invariable
setRateMultipliers(s,p,0,1:3)
# 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
# }
Run the code above in your browser using DataLab