The Goldman and Yang (1994) substitution model for codons. More...

#include <Bpp/Phyl/Model/Codon/GY94.h>

Inheritance diagram for bpp::GY94:

Collaboration diagram for bpp::GY94:

Public Member Functions
	GY94 (const GeneticCode gc, FrequenciesSet codonFreqs)

	~GY94 ()

	GY94 (const GY94 &gy94)

GY94 &	operator= (const GY94 &gy94)

GY94 *	clone () const

std::string	getName () const
	Get the name of the model. More...

const SubstitutionModel &	getModel () const

const std::vector< int > &	getAlphabetStates () const

const StateMap &	getStateMap () const

int	getAlphabetStateAsInt (size_t i) const

std::string	getAlphabetStateAsChar (size_t i) const

std::vector< size_t >	getModelStates (int code) const
	Get the state in the model corresponding to a particular state in the alphabet. More...

std::vector< size_t >	getModelStates (const std::string &code) const
	Get the state in the model corresponding to a particular state in the alphabet. More...

virtual double	freq (size_t i) const

virtual double	Qij (size_t i, size_t j) const

virtual double	Pij_t (size_t i, size_t j, double t) const

virtual double	dPij_dt (size_t i, size_t j, double t) const

virtual double	d2Pij_dt2 (size_t i, size_t j, double t) const

virtual const Vdouble &	getFrequencies () const

const Matrix< double > &	getGenerator () const

const Matrix< double > &	getExchangeabilityMatrix () const

double	Sij (size_t i, size_t j) const

const Matrix< double > &	getPij_t (double t) const

const Matrix< double > &	getdPij_dt (double t) const

const Matrix< double > &	getd2Pij_dt2 (double t) const

void	enableEigenDecomposition (bool yn)
	Set if eigenValues and Vectors must be computed. More...

bool	enableEigenDecomposition ()
	Tell if eigenValues and Vectors must be computed. More...

bool	isDiagonalizable () const

bool	isNonSingular () const

const Vdouble &	getEigenValues () const

const Vdouble &	getIEigenValues () const

const Matrix< double > &	getRowLeftEigenVectors () const

const Matrix< double > &	getColumnRightEigenVectors () const

double	getRate () const
	Get the rate. More...

void	setRate (double rate)
	Set the rate of the model (must be positive). More...

void	addRateParameter ()

void	setFreqFromData (const SequenceContainer &data, double pseudoCount=0)
	Set equilibrium frequencies equal to the frequencies estimated from the data. More...

void	setFreq (std::map< int, double > &frequ)
	Set equilibrium frequencies. More...

const Alphabet *	getAlphabet () const

size_t	getNumberOfStates () const
	Get the number of states. More...

double	getInitValue (size_t i, int state) const throw (BadIntException)

const FrequenciesSet *	getFrequenciesSet () const
	If the model owns a FrequenciesSet, returns a pointer to it, otherwise return 0. More...

virtual void	fireParameterChanged (const ParameterList &parameters)
	Tells the model that a parameter value has changed. More...

double	getScale () const
	Get the scalar product of diagonal elements of the generator and the frequencies vector. If the generator is normalized, then scale=1. Otherwise each element must be multiplied by 1/scale. More...

void	setScale (double scale)
	Multiplies the current generator by the given scale. More...

Protected Member Functions
virtual void	updateMatrices ()

Protected Attributes
std::map< std::string, std::string >	mapParNamesFromPmodel_
	Tools to make the link between the Parameters of the object and those of pmixmodel_. More...

ParameterList	lParPmodel_

Private Member Functions
SubstitutionModel &	getModel ()

Private Attributes
GranthamAAChemicalDistance	gacd_

std::auto_ptr< CodonDistanceFrequenciesSubstitutionModel >	pmodel_

Detailed Description

The Goldman and Yang (1994) substitution model for codons.

Author: Laurent Guéguen

This model has one rate of transitions and one rate of transversion. It also allows distinct equilibrium frequencies between codons. A multiplicative factor accounts for the selective restraints at the amino acid level. This factor applies on the distance $d$ between amino acids given by Grantham (1974).

For codons $i=i_1i_2i_3$ and $j=j_1j_2j_3$ , the generator term $Q_{ij} (i \neq j)$ is:

0 if 2 or 3 of the pair $(i_1,j_1) (i_2,j_2) (i_3,j_3) $ are different.

$\mu \pi_j \exp(-d_{aa_i,aa_j}/V)$ if exactly 1 of the pairs $(i_1,j_1) (i_2,j_2) (i_3,j_3) $ is different, and that difference is a transversion.

$\mu \kappa \pi_j \exp(-d_{aa_i,aa_j}/V)$ if exactly 1 of the pairs $(i_1,j_1) (i_2,j_2) (i_3,j_3) $ is different, and that difference is a transition.

$\mu$ is a normalization factor.

This model includes 2 parameters ( $\kappa$ and $V$ ). The codon frequencies $\pi_j$ are either observed or infered.

Reference:

Goldman N. and Yang Z. (1994), Molecular Biology And Evolution 11(5) 725–736.

Definition at line 82 of file GY94.h.

Constructor & Destructor Documentation

◆ GY94() [1/2]

GY94::GY94	(	const GeneticCode *	gc,
		FrequenciesSet *	codonFreqs
	)

Definition at line 49 of file GY94.cpp.

References bpp::AbstractBiblioSubstitutionModel::lParPmodel_, bpp::AbstractBiblioSubstitutionModel::mapParNamesFromPmodel_, pmodel_, and bpp::AbstractBiblioSubstitutionModel::updateMatrices().

Referenced by clone().

◆ ~GY94()

GY94::~GY94 ( )

Definition at line 85 of file GY94.cpp.

◆ GY94() [2/2]

GY94::GY94 ( const GY94 & gy94 )

Definition at line 72 of file GY94.cpp.

Member Function Documentation

◆ addRateParameter()

void AbstractBiblioSubstitutionModel::addRateParameter ( )

virtualinherited

i	the index of the state in the model.
state	An observed state in the sequence/site.

IndexOutOfBoundsException	if array position is out of range.
BadIntException	if states are not allowed in the associated alphabet.

data	The sequences to use.
pseudoCount	A quantity $\psi$ to add to adjust the observed values in order to prevent issues due to missing states on small data set. The corrected frequencies shall be computed as $\pi_i = \frac{n_i+\psi}{\sum_j (f_j+\psi)}$

Public Member Functions

Protected Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ GY94() [1/2]

◆ ~GY94()

◆ GY94() [2/2]

Member Function Documentation

◆ addRateParameter()

◆ clone()

◆ d2Pij_dt2()

◆ dPij_dt()

◆ enableEigenDecomposition() [1/2]

◆ enableEigenDecomposition() [2/2]

◆ fireParameterChanged()

◆ freq()

◆ getAlphabet()

◆ getAlphabetStateAsChar()

◆ getAlphabetStateAsInt()

◆ getAlphabetStates()

◆ getColumnRightEigenVectors()

◆ getd2Pij_dt2()

◆ getdPij_dt()

◆ getEigenValues()

◆ getExchangeabilityMatrix()

◆ getFrequencies()

◆ getFrequenciesSet()

◆ getGenerator()

◆ getIEigenValues()

◆ getInitValue()

◆ getModel() [1/2]

◆ getModel() [2/2]

◆ getModelStates() [1/2]

◆ getModelStates() [2/2]

◆ getName()

◆ getNumberOfStates()

◆ getPij_t()

◆ getRate()

◆ getRowLeftEigenVectors()

◆ getScale()

◆ getStateMap()

◆ isDiagonalizable()

◆ isNonSingular()

◆ operator=()

◆ Pij_t()

◆ Qij()

◆ setFreq()

◆ setFreqFromData()

◆ setRate()

◆ setScale()

◆ Sij()

◆ updateMatrices()

Member Data Documentation

◆ gacd_

◆ lParPmodel_

◆ mapParNamesFromPmodel_

◆ pmodel_