|
bpp-phyl3 3.0.0
|
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 SubstitutionModelInterface > | getSubstitutionModel () const |
| Get the substitution model associated to the mutation process. More... | |
Protected Attributes | |
| std::shared_ptr< const SubstitutionModelInterface > | model_ |
| 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... | |
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.
|
inline |
Definition at line 298 of file MutationProcess.h.
|
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.
| initialState | The state before beginning evolution. |
| time | The time during which evolution must occur. |
Implements bpp::MutationProcess.
Definition at line 70 of file MutationProcess.cpp.
References bpp::MutationPath::addEvent(), getTimeBeforeNextMutationEvent(), model_, and mutate().
Referenced by bpp::NonHomogeneousSequenceSimulator::dEvolveInternal(), and bpp::SimpleSubstitutionProcessSiteSimulator::evolveInternal().
|
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.
| initialState | The state before beginning evolution. |
| finalState | The state after evolution. |
| time | The time during which evolution must occur. |
Implements bpp::MutationProcess.
Definition at line 89 of file MutationProcess.cpp.
References bpp::MutationPath::addEvent(), bpp::MutationPath::clear(), bpp::exp(), getTimeBeforeNextMutationEvent(), bpp::RandomTools::giveRandomNumberBetweenZeroAndEntry(), bpp::log(), model_, and mutate().
|
virtual |
Simulation a character evolution during a specified time according to the given substitution model and send the final state.
| initialState | The state before beginning evolution. |
| time | The time during which evolution must occur. |
Implements bpp::MutationProcess.
Reimplemented in bpp::SimpleMutationProcess.
Definition at line 55 of file MutationProcess.cpp.
References getTimeBeforeNextMutationEvent(), and mutate().
|
inlinevirtual |
Get the substitution model associated to the mutation process.
Implements bpp::MutationProcess.
Definition at line 309 of file MutationProcess.h.
References model_.
|
virtual |
Get the time before next mutation event.
| state | The actual state of the chain; |
Implements bpp::MutationProcess.
Definition at line 48 of file MutationProcess.cpp.
References model_, and bpp::RandomTools::randExponential().
Referenced by detailedEvolve(), and evolve().
|
virtual |
Mutate a character in state i.
| state | The current state of the character. |
Implements bpp::MutationProcess.
Definition at line 15 of file MutationProcess.cpp.
References bpp::RandomTools::giveRandomNumberBetweenZeroAndEntry(), repartition_, size_, and bpp::TextTools::toString().
Referenced by detailedEvolve(), and evolve().
|
virtual |
Mutate a character in state i n times.
| state | The current state of the character. |
| n | The number of mutations to perform. |
Implements bpp::MutationProcess.
Definition at line 28 of file MutationProcess.cpp.
References bpp::RandomTools::giveRandomNumberBetweenZeroAndEntry(), repartition_, and size_.
|
protected |
The substitution model to use:
Definition at line 282 of file MutationProcess.h.
Referenced by detailedEvolve(), bpp::SimpleMutationProcess::evolve(), getSubstitutionModel(), and getTimeBeforeNextMutationEvent().
|
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 mutate(), bpp::SelfMutationProcess::SelfMutationProcess(), and bpp::SimpleMutationProcess::SimpleMutationProcess().
|
protected |
The number of states allowed for the character to mutate.
Definition at line 287 of file MutationProcess.h.
Referenced by bpp::SimpleMutationProcess::evolve(), mutate(), bpp::SelfMutationProcess::SelfMutationProcess(), and bpp::SimpleMutationProcess::SimpleMutationProcess().