bpp-phyl3  3.0.0
bpp::AbstractMutationProcess Class Reference

Partial implementation of the MutationProcess interface. More...

#include <Bpp/Phyl/Simulation/MutationProcess.h>

+ Inheritance diagram for bpp::AbstractMutationProcess:
+ Collaboration diagram for bpp::AbstractMutationProcess:

Public Member Functions

 AbstractMutationProcess (std::shared_ptr< const SubstitutionModelInterface > model)
 
size_t mutate (size_t state) const
 Mutate a character in state i. More...
 
size_t mutate (size_t state, unsigned int n) const
 Mutate a character in state i n times. More...
 
double getTimeBeforeNextMutationEvent (size_t state) const
 Get the time before next mutation event. More...
 
size_t evolve (size_t initialState, double time) const
 Simulation a character evolution during a specified time according to the given substitution model and send the final state. More...
 
MutationPath detailedEvolve (size_t initialState, double time) const
 Simulation a character evolution during a specified time according to the given substitution model and send the total path with all intermediate states and times between mutation events. More...
 
MutationPath detailedEvolve (size_t initialState, size_t finalState, double time) const
 Simulation a character evolution during a specified time according to the given substitution model and send the total path with all intermediate states and times between mutation events, conditional to the final state.q. More...
 
std::shared_ptr< const SubstitutionModelInterfacegetSubstitutionModel () const
 Get the substitution model associated to the mutation process. More...
 

Protected Attributes

std::shared_ptr< const SubstitutionModelInterfacemodel_
 The substitution model to use: More...
 
size_t size_
 The number of states allowed for the character to mutate. More...
 
VVdouble repartition_
 The repartition function for states probabilities. More...
 

Detailed Description

Partial implementation of the MutationProcess interface.

This class provides an implementation of the MutationProcess interface. It assumes that there are size_ states allowed for the character of interest, and that the distribution of probabilities are in repartition_. As a matter of facts, probabilities must be cumulative, so that repartition_ contains values of the repartition function. The mutate function hence draws a random number between 0 and 1 and gives the corresponding character using the bijection of the repartition function.

All derived classes must initialize the repartition_ and size_ fields.

Definition at line 275 of file MutationProcess.h.

Constructor & Destructor Documentation

◆ AbstractMutationProcess()

bpp::AbstractMutationProcess::AbstractMutationProcess ( std::shared_ptr< const SubstitutionModelInterface model)
inline

Definition at line 298 of file MutationProcess.h.

Member Function Documentation

◆ detailedEvolve() [1/2]

MutationPath AbstractMutationProcess::detailedEvolve ( size_t  initialState,
double  time 
) const
virtual

Simulation a character evolution during a specified time according to the given substitution model and send the total path with all intermediate states and times between mutation events.

Parameters
initialStateThe state before beginning evolution.
timeThe time during which evolution must occur.
Returns
The resulting mutation path.

Implements bpp::MutationProcess.

Definition at line 70 of file MutationProcess.cpp.

References bpp::MutationPath::addEvent().

Referenced by bpp::NonHomogeneousSequenceSimulator::dEvolveInternal(), and bpp::SimpleSubstitutionProcessSiteSimulator::evolveInternal().

◆ detailedEvolve() [2/2]

MutationPath AbstractMutationProcess::detailedEvolve ( size_t  initialState,
size_t  finalState,
double  time 
) const
virtual

Simulation a character evolution during a specified time according to the given substitution model and send the total path with all intermediate states and times between mutation events, conditional to the final state.q.

Parameters
initialStateThe state before beginning evolution.
finalStateThe state after evolution.
timeThe time during which evolution must occur.
Returns
The resulting mutation path.

Implements bpp::MutationProcess.

Definition at line 89 of file MutationProcess.cpp.

References bpp::MutationPath::addEvent(), bpp::MutationPath::clear(), bpp::exp(), bpp::RandomTools::giveRandomNumberBetweenZeroAndEntry(), and bpp::log().

◆ evolve()

size_t AbstractMutationProcess::evolve ( size_t  initialState,
double  time 
) const
virtual

Simulation a character evolution during a specified time according to the given substitution model and send the final state.

Parameters
initialStateThe state before beginning evolution.
timeThe time during which evolution must occur.
Returns
The resulting state after evolution is completed.

Implements bpp::MutationProcess.

Reimplemented in bpp::SimpleMutationProcess.

Definition at line 55 of file MutationProcess.cpp.

◆ getSubstitutionModel()

std::shared_ptr<const SubstitutionModelInterface> bpp::AbstractMutationProcess::getSubstitutionModel ( ) const
inlinevirtual

Get the substitution model associated to the mutation process.

Returns
The SubstitutionModel associated to this instance.

Implements bpp::MutationProcess.

Definition at line 309 of file MutationProcess.h.

References model_.

◆ getTimeBeforeNextMutationEvent()

double AbstractMutationProcess::getTimeBeforeNextMutationEvent ( size_t  state) const
virtual

Get the time before next mutation event.

Parameters
stateThe actual state of the chain;
Returns
A random time before next mutation event.

Implements bpp::MutationProcess.

Definition at line 48 of file MutationProcess.cpp.

References bpp::RandomTools::randExponential().

◆ mutate() [1/2]

size_t AbstractMutationProcess::mutate ( size_t  state) const
virtual

Mutate a character in state i.

Parameters
stateThe current state of the character.

Implements bpp::MutationProcess.

Definition at line 15 of file MutationProcess.cpp.

References bpp::RandomTools::giveRandomNumberBetweenZeroAndEntry(), and bpp::TextTools::toString().

◆ mutate() [2/2]

size_t AbstractMutationProcess::mutate ( size_t  state,
unsigned int  n 
) const
virtual

Mutate a character in state i n times.

Parameters
stateThe current state of the character.
nThe number of mutations to perform.

Implements bpp::MutationProcess.

Definition at line 28 of file MutationProcess.cpp.

References bpp::RandomTools::giveRandomNumberBetweenZeroAndEntry().

Member Data Documentation

◆ model_

std::shared_ptr<const SubstitutionModelInterface> bpp::AbstractMutationProcess::model_
protected

The substitution model to use:

Definition at line 282 of file MutationProcess.h.

Referenced by bpp::SimpleMutationProcess::evolve(), and getSubstitutionModel().

◆ repartition_

VVdouble bpp::AbstractMutationProcess::repartition_
protected

The repartition function for states probabilities.

repartition_[i][j] = probability that, being in state i at time t, we'll be in state <= j at time t+1.

Definition at line 295 of file MutationProcess.h.

Referenced by bpp::SelfMutationProcess::SelfMutationProcess(), and bpp::SimpleMutationProcess::SimpleMutationProcess().

◆ size_

size_t bpp::AbstractMutationProcess::size_
protected

The number of states allowed for the character to mutate.

Definition at line 287 of file MutationProcess.h.

Referenced by bpp::SimpleMutationProcess::evolve(), bpp::SelfMutationProcess::SelfMutationProcess(), and bpp::SimpleMutationProcess::SimpleMutationProcess().


The documentation for this class was generated from the following files: