340 likes | 355 Views
This chapter explores Markov Chains and their properties, using nucleotide evolution as an example. Topics include types of point mutations, Kimura's 2-parameter model, and Markov Chain properties.
E N D
A First Course in Stochastic Processes Chapter Two: Markov Chains
X2= X2=2 X3=1 X4=3 X1=1
X1 X2 X3 X4 X5 etc
G C T Example Two: Nucleotide evolution A
α β β β β α Types of point mutation A G Purine Transitions Transversions T C Pyramidine Transitions
A G C T A G P = C T Kimura’s 2 parameter model (K2P)
G C C G A C T G A T G C G C T G T A C A G T A T C A T
G C C G A C T G A T G C A G T A T C A T The Markov Property
Markov Chain properties accessible aperiodic communicate recurrent irreducible transient
A G C T A G P = C T Accessible 0
A G C T A G P = C T Accessible A (and G) are no longer accessible from C (or T). 0 0 0 0
A G C T A G P = C T Accessible But C (and T) are still accessible from A (or G). 0 0 0 0
A G C T A G P = C T Communicate Reciprocal accessibility
A G C T A G P = C T Irreducible All elements communicate
A G C T A P1 0 0 0 G P = = 0 P2 0 0 C 0 0 T 0 0 A G C T P1 = P2 = A C G T Non-irreducible
Repercussions of communication • Reflexivity • Symmetry • Transitivity
P = Periodicity
Periodicity • The period d(i) of an element i is defined as the greatest common divisor of the numbers of the generations in which the element is visited. • Most Markov Chains that we deal with do not exhibit periodicity. • A Markov Chain is aperiodicif d(i) = 1 for all i.
recurrent transient Recurrence
More on Recurrence • and i is recurrent then j is recurrent • In a one-dimensional symmetric random walk the origin is recurrent • In a two-dimensional symmetric random walk the origin is recurrent • In a three-dimensional symmetric random walk the origin is transient
Markov Chain properties accessible aperiodic communicate recurrent irreducible transient
Markov Chains Examples
X1=1 X1 X2 X3 X4 X5 etc
Haldane (1927) branching process model of fixation probability 2 3 4 4 4 4 2
Haldane (1927) branching process model of fixation probability
Haldane (1927) branching process model of fixation probability Pi,j = coefficient of sj in the above generating function
Haldane (1927) branching process model of fixation probability Probability of fixation = 2s
Markov Chain properties accessible aperiodic communicate recurrent irreducible transient