Rev Language Reference
fnGTR - The General Time-Reversible rate matrix
The GTR rate matrix.
Usage
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)
Edit this page