© Copyright 2000, 2001, Donald Simon & Bret Larget, Department of Mathematics and Computer Science, Duquesne University.
We parameterize the instantaneous rate matrix as:
+- -+ | -() K pi_G pi_C pi_T | | | | K pi_A -() pi_C pi_T | theta | | | pi_A pi_G -() K gamma pi_T | | | | pi_A pi_G K gamma pi_C -() | +- -+
where K = 4*ttr/(gamma+1), and we implicitly assume the order A,G,C,T for bases.
There are seven parameters, six of which are free. The model is reversible with stationary distribution given by pi_A, pi_G, pi_C, pi_T which are constrained to sum to 1. The parameter theta controls the overall mutation rate. The transition/transversion ratio is ttr and gamma is the final parameter which affects the ratio of transition/transversion rates among purines and pyrimidines.
Under HKY85 and F84, there is an additional parameter kappa, and ttr and gamma are functions of kappa and the pi's.
When the value kappa is input to F84, it is different than the transition/transversion parameter `T' that PHYLIP takes as input.
Back to the table of contents.