Rev Language Reference


fnGTR - The General Time-Reversible rate matrix

The GTR rate matrix.

Usage

fnGTR(Simplex exchangeRates, Simplex baseFrequencies)

Arguments

exchangeRates : Simplex (pass by const reference)
The exchangeability rates between states.
baseFrequencies : Simplex (pass by const reference)
The stationary frequencies of the states.

Return Type

Details

The general time-reversible rate matrix: Q(i,j) = exchangeRates(i,j) * baseFrequencies[j] The exchangeRates are symmetric.

Example

er ~ dnDirichlet( v(1,1,1,1,1,1) )
pi ~ dnDirichlet( v(1,1,1,1) )
Q := fnGTR(er,pi)