|
|
| Line 1: |
Line 1: |
| {{About|Markov chains in discrete time|Markov chains in continuous time|continuous-time Markov chain}}
| | I'm Celeste and I live with my husband and our 3 children in Premer, in the NSW south part. My hobbies are Inline Skating, Coloring and Worldbuilding.<br><br>Stop by my site ... [http://quintomyap.blogspot.com/2010/03/930am-in-morning-2-igo-seniors-playing.html Fifa 15 coin generator] |
| [[File:Markovkate 01.svg|thumb|right|A simple two-state Markov chain]]
| |
| A '''Markov chain''' ('''discrete-time Markov chain''' or '''DTMC'''<ref>{{cite book | last = Norris | first = James R. | authorlink = James R. Norris | title = Markov chains | publisher = Cambridge University Press | year = 1998 | url = http://www.statslab.cam.ac.uk/~james/Markov/}}</ref>) named after [[Andrey Markov]], is a mathematical system that undergoes transitions from one state to another on a [[state space]]. It is a [[stochastic process|random process]] usually characterized as [[Memorylessness|memoryless]]: the next state depends only on the current state and not on the sequence of events that preceded it. This specific kind of "memorylessness" is called the [[Markov process#Markov property|Markov property]]. Markov chains have many applications as [[statistical modeling|statistical models]] of real-world processes.
| |
| | |
| ==Introduction==
| |
| [[File:AAMarkov.jpg|thumb|220px|right|Russian mathematician [[Andrey Markov]], the namesake.]]
| |
| A Markov chain is a [[stochastic process]] with the [[Markov property]]. The term "Markov chain" refers to the sequence of random variables such a process moves through, with the Markov property defining [[serial dependence]] only between adjacent periods (as in a "chain"). It can thus be used for describing systems that follow a chain of linked events, where what happens next depends only on the current state of the system.
| |
| | |
| In literature different [[Markov process]]es are designated as "Markov chains". Usually however, the term is reserved for a process with a discrete set of times (i.e. a [[discrete-time]] Markov chain (DTMC))<ref>Everitt,B.S. (2002) ''The Cambridge Dictionary of Statistics''. CUP. ISBN 0-521-81099-X</ref> although some authors use the same terminology to refer to a [[continuous-time Markov chain]] without explicit mention.<ref>Parzen, E. (1962) ''Stochastic Processes'', Holden-Day. ISBN 0-8162-6664-6 (Table 6.1))</ref><ref>Dodge, Y. (2003) ''The Oxford Dictionary of Statistical Terms'', OUP. ISBN 0-19-920613-9 (entry for "Markov chain")</ref> While the time parameter definition is mostly agreed upon to mean discrete-time, the Markov chain [[state space]] does not have an established definition: the term may refer to a process on an arbitrarily general state space.<ref>Meyn, S. Sean P., and Richard L. Tweedie. (2009) ''Markov chains and stochastic stability''. Cambridge University Press. (Preface, p. iii)</ref> However, many uses of Markov chains employ finite or [[countable]] (or [[Discrete space|discrete]] on the [[real line]]) state spaces, which have a more straightforward statistical analysis. Because there are—besides time index and state space parameters—many other variations, extensions and generalisations (see [[#Variations|Variations]]), the remainder of this article concentrates on the simplest discrete-time, discrete state-space case, unless mentioned otherwise.
| |
| | |
| The changes of state of the system are called transitions, and the probabilities associated with various state changes are called transition probabilities. The process is characterized by a state space, a [[transition matrix]] describing the probabilities of particular transitions, and an initial state (or initial distribution) across the state space. By convention, we assume all possible states and transitions have been included in the definition of the process, so there is always a next state, and the process does not terminate.
| |
| | |
| A discrete-time random process involves a system which is in a certain state at each step, with the state changing randomly between steps. The steps are often thought of as moments in time, but they can equally well refer to physical distance or any other discrete measurement. Formally, the steps are the [[integers]] or [[natural numbers]], and the random process is a mapping of these to states. The Markov property states that the [[conditional probability distribution]] for the system at the next step (and in fact at all future steps) depends only on the current state of the system, and not additionally on the state of the system at previous steps.
| |
| | |
| Since the system changes randomly, it is generally impossible to predict with certainty the state of a Markov chain at a given point in the future. However, the statistical properties of the system's future can be predicted. In many applications, it is these statistical properties that are important.
| |
| | |
| A famous Markov chain is the so-called "drunkard's walk", a [[random walk]] on the [[number line]] where, at each step, the position may change by +1 or −1 with equal probability. From any position there are two possible transitions, to the next or previous integer. The transition probabilities depend only on the current position, not on the manner in which the position was reached. For example, the transition probabilities from 5 to 4 and 5 to 6 are both 0.5, and all other transition probabilities from 5 are 0. These probabilities are independent of whether the system was previously in 4 or 6.
| |
| | |
| Another example is the dietary habits of a creature who eats only grapes, cheese, or lettuce, and whose dietary habits conform to the following rules:
| |
| | |
| * It eats exactly once a day.
| |
| * If it ate cheese today, tomorrow it will eat lettuce or grapes with equal probability.
| |
| * If it ate grapes today, tomorrow it will eat grapes with probability 1/10, cheese with probability 4/10 and lettuce with probability 5/10.
| |
| * If it ate lettuce today, it will not eat lettuce again tomorrow but will eat grapes with probability 4/10 or cheese with probability 6/10.
| |
| | |
| This creature's eating habits can be modeled with a Markov chain since its choice tomorrow depends solely on what it ate today, not what it ate yesterday or even farther in the past. One statistical property that could be calculated is the expected percentage, over a long period, of the days on which the creature will eat grapes.
| |
| | |
| A series of independent events (for example, a series of coin flips) satisfies the formal definition of a Markov chain. However, the theory is usually applied only when the probability distribution of the next step depends non-trivially on the current state.
| |
| | |
| Many other [[examples of Markov chains]] exist.
| |
| | |
| ==Formal definition==
| |
| A Markov chain is a sequence of [[random variable]]s ''X''<sub>1</sub>, ''X''<sub>2</sub>, ''X''<sub>3</sub>, ... with the Markov property, namely that, given the present state, the future and past states are independent. Formally,
| |
| | |
| :<math>\Pr(X_{n+1}=x|X_1=x_1, X_2=x_2, \ldots, X_n=x_n) = \Pr(X_{n+1}=x|X_n=x_n)</math>, if both sides of the equation are well defined.
| |
| | |
| The possible values of ''X''<sub>''i''</sub> form a [[countable set]] ''S'' called the '''state space''' of the chain.
| |
| | |
| Markov chains are often described by a [[directed graph]], where the edges are labeled by the probabilities of going from one state to the other states.
| |
| | |
| ===Variations===
| |
| * [[Continuous-time Markov process]]es have a continuous index.
| |
| * {{anchor|homogeneous}}'''Time-homogeneous Markov chains''' (or '''stationary Markov chains''') are processes where
| |
| | |
| ::<math>\Pr(X_{n+1}=x|X_n=y) = \Pr(X_n=x|X_{n-1}=y)\,</math>
| |
| | |
| : for all ''n''. The probability of the transition is independent of ''n''.
| |
| * A '''Markov chain of order ''m'' ''' (or a Markov chain with memory ''m''), where ''m'' is finite, is a process satisfying
| |
| | |
| ::<math>
| |
| \begin{align}
| |
| {} &\Pr(X_n=x_n|X_{n-1}=x_{n-1}, X_{n-2}=x_{n-2}, \dots , X_1=x_1) \\
| |
| = &\Pr(X_n=x_n|X_{n-1}=x_{n-1}, X_{n-2}=x_{n-2}, \dots, X_{n-m}=x_{n-m})
| |
| \text{ for }n > m
| |
| \end{align}
| |
| </math>
| |
| | |
| : In other words, the future state depends on the past ''m'' states. It is possible to construct a chain (''Y<sub>n</sub>'') from (''X<sub>n</sub>'') which has the 'classical' Markov property by taking as state space the ordered ''m''-tuples of ''X'' values, ie. ''Y<sub>n</sub>'' = (''X<sub>n</sub>'', ''X''<sub>''n''−1</sub>, ..., ''X''<sub>''n''−''m''+1</sub>).
| |
| | |
| * An [[additive Markov chain]] of order ''m'' is determined by an additive conditional probability,
| |
| | |
| ::<math>\Pr(X_n=x_n|X_{n-1}=x_{n-1}, X_{n-2}=x_{n-2}, \dots, X_{n-m}=x_{n-m}) = \sum_{r=1}^{m} f(x_n,x_{n-r},r) .</math>
| |
| | |
| The value ''f''(''x<sub>n</sub>,x<sub>n-r</sub>,r'') is the additive contribution of the variable ''x<sub>n-r</sub>'' to the conditional probability.<ref>Usatenko, O. V.; Apostolov, S. S.; Mayzelis, Z. A.; Melnik, S. S. (2010) ''Random finite-valued dynamical systems: additive Markov chain approach''. Cambridge Scientific Publisher. ISBN 978-1-904868-74-3</ref>{{Clarify|reason=why is this form useful|date=June 2009}}
| |
| | |
| ===Example===
| |
| {{Main|Examples of Markov chains}}
| |
| [[File:Finance Markov chain example state space.svg|thumb|400px|right|<!--Directed graph representation of a Markov chain for the state of financial markets (note : numbers are made-up).-->]]
| |
| A [[state diagram]] for a simple example is shown in the figure on the right, using a directed graph to picture the [[state transition]]s. The states represent whether a hypothetical stock market is exhibiting a [[Market trend#Bull market|bull market]], [[Market trend#Bear market|bear market]], or stagnant market trend during a given week. According to the figure, a bull week is followed by another bull week 90% of the time, a bear week 7.5% of the time, and a stagnant week the other 2.5% of the time. Labelling the state space {1 = bull, 2 = bear, 3 = stagnant} the [[transition matrix]] for this example is
| |
| :<math>P = \begin{bmatrix}
| |
| 0.9 & 0.075 & 0.025 \\
| |
| 0.15 & 0.8 & 0.05 \\
| |
| 0.25 & 0.25 & 0.5
| |
| \end{bmatrix}.</math>
| |
| The distribution over states can be written as a [[stochastic row vector]] ''x'' with the relation ''x''<sup>(''n'' + 1)</sup> = ''x''<sup>(''n'')</sup>''P''. So if at time ''n'' the system is in state 2 (bear), then three time periods later, at time ''n'' + 3 the distribution is
| |
| :<math>\begin{align}
| |
| x^{(n+3)} &= x^{(n+2)} P = \left(x^{(n+1)} P\right) P \\\\
| |
| &= x^{(n+1)} P^2 = \left( x^{(n)} P^2 \right) P\\
| |
| &= x^{(n)} P^3 \\
| |
| &= \begin{bmatrix} 0 & 1 & 0 \end{bmatrix} \begin{bmatrix}
| |
| 0.9 & 0.075 & 0.025 \\
| |
| 0.15 & 0.8 & 0.05 \\
| |
| 0.25 & 0.25 & 0.5
| |
| \end{bmatrix}^3 \\
| |
| &= \begin{bmatrix} 0 & 1 & 0 \end{bmatrix} \begin{bmatrix}
| |
| 0.7745 & 0.17875 & 0.04675 \\
| |
| 0.3575 & 0.56825 & 0.07425 \\
| |
| 0.4675 & 0.37125 & 0.16125 \\
| |
| \end{bmatrix} \\
| |
| & = \begin{bmatrix} 0.3575 & 0.56825 & 0.07425 \end{bmatrix}.
| |
| \end{align}</math>
| |
| Using the transition matrix it is possible to calculate, for example, the long-term fraction of weeks during which the market is stagnant, or the average number of weeks it will take to go from a stagnant to a bull market. Using the transition probabilities, the steady-state probabilities indicate that 62.5% of weeks will be in a bull market, 31.25% of weeks will be in a bear market and 6.25% of weeks will be stagnant, since:
| |
| | |
| <math>\lim_{N\to \infty } \, P^N=
| |
| \begin{bmatrix}
| |
| 0.625 & 0.3125 & 0.0625 \\
| |
| 0.625 & 0.3125 & 0.0625 \\
| |
| 0.625 & 0.3125 & 0.0625 \\
| |
| \end{bmatrix}</math>
| |
| | |
| A thorough development and many examples can be found in the on-line monograph
| |
| Meyn & Tweedie 2005.<ref name=MCSS>S. P. Meyn and R.L. Tweedie, 2005. [https://netfiles.uiuc.edu/meyn/www/spm_files/book.html Markov Chains and Stochastic Stability].
| |
| Second edition to appear, Cambridge University Press, 2008.</ref>
| |
| | |
| The appendix of Meyn 2007,<ref name=CTCN>S. P. Meyn, 2007. [http://decision.csl.uiuc.edu/~meyn/pages/CTCN/CTCN.html Control Techniques for Complex Networks], Cambridge University Press, 2007.</ref> also available on-line, contains an abridged Meyn & Tweedie.
| |
| | |
| A [[finite state machine]] can be used as a representation of a Markov chain. Assuming a sequence of [[Independent and identically distributed random variables|independent and identically distributed]] input signals (for example, symbols from a binary alphabet chosen by coin tosses), if the machine is in state ''y'' at time ''n'', then the probability that it moves to state ''x'' at time ''n'' + 1 depends only on the current state.
| |
| | |
| ===System Dynamics approach to Markov Chain===
| |
| | |
| [[File:Markov system represented by Stock & Flow.jpg|thumb|Stock and flow diagram of the Markov chain example for stock market]]
| |
| [[File:SD Output BOT.jpg|thumb|Transient behavior of the Markov chain example for stock market]]
| |
| | |
| The Markov chain can be easily represented by the [[stock and flow]] diagram of [[System dynamics]]. One of the benefits of SD approach to Markov chain is that the transition probability can be dynamic (i.e., non-constant). The output from SD shows that the stock market will reach the steady state after 10 weeks.
| |
| | |
| The equations for the stock market model can also be represented by systems of ordinary differential equations. This dynamic system is a form of '''dX/dt=Ax''', thus the eigenvalues and eigenvectors of '''A''' will show the variation modes of the stock market. :
| |
| | |
| <math>
| |
| \begin{align}
| |
| \frac{d(Bear\;to\;Bull)}{dt} & = P(Bear\;to\;Bull) * Bear \\
| |
| \frac{d(Bull\;to\;Bear)}{dt} & = P(Bull\;to\;Bear) * Bull \\
| |
| \frac{d(Bear\;to\;Stagnant)}{dt} & = P(Bear\;to\;Stagnant) * Bear \\
| |
| \frac{d(Stagnant\;to\;Bear)}{dt} & = P(Stagnant\;to\;Bear) * Stagnant \\
| |
| \frac{d(Bull\;to\;Stagnant)}{dt} & = P(Bull\;to\;Stagnant) * Bull \\
| |
| \frac{d(Stagnant\;to\;Bull)}{dt} & = P(Stagnant\;to\;Bull) * Stagnant
| |
| \end{align}
| |
| </math>
| |
| | |
| with initial values for the three states
| |
| | |
| <math>
| |
| \begin{align}
| |
| Bull(t=0) & = 1 \\
| |
| Bear(t=0) & = 0 \\
| |
| Stagnant(t=0) & = 0
| |
| \end{align}
| |
| </math>
| |
| | |
| ==Transient evolution==
| |
| | |
| The probability of going from state ''i'' to state ''j'' in ''n'' time steps is
| |
| | |
| :<math>p_{ij}^{(n)} = \Pr(X_n=j\mid X_0=i) \,</math>
| |
| | |
| and the single-step transition is
| |
| | |
| :<math>p_{ij} = \Pr(X_1=j\mid X_0=i). \,</math>
| |
| | |
| For a time-homogeneous Markov chain:
| |
| | |
| :<math>p_{ij}^{(n)} = \Pr(X_{k+n}=j \mid X_{k}=i) \,</math>
| |
| | |
| and
| |
| | |
| :<math>p_{ij} = \Pr(X_{k+1}=j \mid X_k=i). \,</math>
| |
| | |
| The ''n''-step transition probabilities satisfy the [[Chapman–Kolmogorov equation]], that for any ''k'' such that 0 < ''k'' < ''n'',
| |
| | |
| :<math>p_{ij}^{(n)} = \sum_{r \in S} p_{ir}^{(k)} p_{rj}^{(n-k)}</math>
| |
| | |
| where ''S'' is the state space of the Markov chain.
| |
| | |
| The [[marginal distribution]] Pr(''X''<sub>''n''</sub> = ''x'') is the distribution over states at time ''n''. The initial distribution is Pr(''X''<sub>0</sub> = ''x''). The evolution of the process through one time step is described by
| |
| | |
| : <math> \Pr(X_{n}=j) = \sum_{r \in S} p_{rj} \Pr(X_{n-1}=r) = \sum_{r \in S} p_{rj}^{(n)} \Pr(X_0=r).</math>
| |
| | |
| Note: The superscript (''n'') is an [[Indexed family|index]] and not an [[exponent]].
| |
| | |
| ==Properties==
| |
| | |
| ===Reducibility===
| |
| | |
| A state ''j'' is said to be '''accessible''' from a state ''i'' (written ''i'' → ''j'') if a system started in state ''i'' has a non-zero probability of transitioning into state ''j'' at some point. Formally, state ''j'' is accessible from state ''i'' if there exists an integer ''n<sub>ij</sub>'' ≥ 0 such that
| |
| | |
| : <math> \Pr(X_{n}=j | X_0=i) = p_{ij}^{(n_{ij})} > 0.\, </math>
| |
| | |
| This integer is allowed to be different for each pair of states, hence the subscripts in n<sub>ij</sub>.
| |
| Allowing ''n'' to be zero means that every state is defined to be accessible from itself.
| |
| | |
| A state ''i'' is said to '''communicate''' with state ''j'' (written ''i'' ↔ ''j'') if both ''i'' → ''j'' and ''j'' → ''i''. A set of states ''C'' is a '''communicating class''' if every pair of states in ''C'' communicates with each other, and no state in ''C'' communicates with any state not in ''C''. It can be shown that communication in this sense is an [[equivalence relation]] and thus that communicating classes are the [[equivalence classes]] of this relation. A communicating class is '''closed''' if the probability of leaving the class is zero, namely that if ''i'' is in ''C'' but ''j'' is not, then ''j'' is not accessible from ''i''.
| |
| | |
| A state ''i'' is said to be '''essential''' or '''final''' if for all ''j'' such that ''i'' → ''j'' it is also true that ''j'' → ''i''. A state ''i'' is '''inessential''' if it is not essential.<ref>{{cite book|last=Asher Levin|first=David|title=Markov chains and mixing times|pages=16|url=http://books.google.co.uk/books?id=6Cg5Nq5sSv4C&lpg=PA16&dq=essential%20class%20markov&pg=PA16#v=onepage|isbn=978-0-8218-4739-8|year=2009}}</ref>
| |
| | |
| A Markov chain is said to be '''irreducible''' if its state space is a single communicating class; in other words, if it is possible to get to any state from any state.
| |
| | |
| ===Periodicity===
| |
| A state ''i'' has '''period''' ''k'' if any return to state ''i'' must occur in multiples of ''k'' time steps. Formally, the [[Periodic function|period]] of a state is defined as
| |
| | |
| : <math> k = \operatorname{gcd}\{ n: \Pr(X_n = i | X_0 = i) > 0\}</math>
| |
| | |
| (where "gcd" is the [[greatest common divisor]]). Note that even though a state has period ''k'', it may not be possible to reach the state in ''k'' steps. For example, suppose it is possible to return to the state in {6, 8, 10, 12, ...} time steps; ''k'' would be 2, even though 2 does not appear in this list.
| |
| | |
| If ''k'' = 1, then the state is said to be '''aperiodic''': returns to state ''i'' can occur at irregular times. In other words, a state ''i'' is '''aperiodic''' if there exists ''n'' such that for all ''n' ≥ n'',
| |
| | |
| : <math> \Pr(X_{n'} = i | X_0 = i) > 0.</math>
| |
| | |
| Otherwise (''k'' > 1), the state is said to be '''periodic with period ''k'''''. A Markov chain is aperiodic if every state is aperiodic. An irreducible markov chain only needs one aperiodic state to imply all states are aperiodic.
| |
| | |
| Every state of a [[bipartite graph]] has an even '''period'''.
| |
| | |
| ===Recurrence===
| |
| A state ''i'' is said to be '''transient''' if, given that we start in state ''i'', there is a non-zero probability that we will never return to ''i''. Formally, let the [[random variable]] ''T<sub>i</sub>'' be the first return time to state ''i'' (the "hitting time"):
| |
| | |
| : <math> T_i = \inf \{ n\ge1: X_n = i | X_0 = i\}.</math>
| |
| | |
| The number
| |
| | |
| : <math> f_{ii}^{(n)} = \Pr(T_i = n)</math>
| |
| | |
| is the probability that we return to state ''i'' for the first time after ''n'' steps.
| |
| Therefore, state ''i'' is transient if
| |
| | |
| : <math> \Pr(T_i < {\infty}) = \sum_{n=1}^{\infty} f_{ii}^{(n)} < 1. </math>
| |
| | |
| State ''i'' is '''recurrent''' (or '''persistent''') if it is not transient.
| |
| Recurrent states have finite hitting time with probability ''1''.
| |
| | |
| ====Mean recurrence time====
| |
| Even if the hitting time is finite with probability ''1'', it need not have a finite [[Expected value|expectation]]<!-- Should provide example later for this and reference it here -->.
| |
| The '''mean recurrence time''' at state ''i'' is the expected return time ''M<sub>i</sub>'':
| |
| | |
| : <math> M_i = E[T_i]=\sum_{n=1}^{\infty} n\cdot f_{ii}^{(n)}.\, </math>
| |
| | |
| State ''i'' is '''positive recurrent''' (or '''non-null persistent''') if ''M<sub>i</sub>'' is finite; otherwise, state ''i'' is '''null recurrent''' (or '''null persistent''').
| |
| | |
| ====Expected number of visits====
| |
| It can be shown that a state ''i'' is recurrent [[if and only if]] the expected number of visits to this state is infinite, i.e.,
| |
| | |
| : <math>\sum_{n=0}^{\infty} p_{ii}^{(n)} = \infty.</math>
| |
| | |
| ====Absorbing states====
| |
| A state ''i'' is called '''absorbing''' if it is impossible to leave this state. Therefore, the state ''i'' is absorbing if and only if
| |
| | |
| : <math> p_{ii} = 1\text{ and }p_{ij} = 0\text{ for }i \not= j.</math>
| |
| | |
| If every state can reach an absorbing state, then the Markov chain is an [[absorbing Markov chain]].
| |
| | |
| ===Ergodicity===
| |
| A state ''i'' is said to be '''[[ergodic theory|ergodic]]''' if it is aperiodic and positive recurrent. In other words, a state ''i'' is ergodic if it is recurrent, has a period of ''1'' and it has finite mean recurrence time. If all states in an irreducible Markov chain are ergodic, then the chain is said to be ergodic.
| |
| | |
| It can be shown that a finite state irreducible Markov chain is ergodic if it has an aperiodic state. A model has the ergodic property if there's a finite number ''N'' such that any state can be reached from any other state in exactly ''N'' steps. In case of a fully connected transition matrix where all transitions have a non-zero probability, this condition is fulfilled with ''N''=1. A model with more than one state and just one out-going transition per state cannot be ergodic.
| |
| | |
| ===Steady-state analysis and limiting distributions===
| |
| If the Markov chain is a [[time-homogeneous Markov chain]], so that the process is described by a single, [[time-independent matrix]] <math>p_{ij}</math>, then the vector <math>\boldsymbol{\pi}</math> is called a '''stationary distribution''' (or '''[[invariant measure]]''') if <math>\forall j \in S</math> it satisfies
| |
| | |
| : <math>0\leq\pi_j\leq1.</math>
| |
| | |
| : <math>\sum_{j \in S}\pi_j = 1.</math>
| |
| | |
| : <math>\pi_j = \sum_{i \in S} \pi_i p_{ij}.</math>
| |
| | |
| An irreducible chain has a stationary distribution if and only if all of its states are positive recurrent. In that case, ''π'' is unique and is related to the expected return time:
| |
| | |
| : <math>\pi_j = \frac{C}{M_j}\,,</math>
| |
| | |
| where <math>C</math> is the normalizing constant. Further, if the chain is both irreducible and aperiodic, then for any ''i'' and ''j'',
| |
| | |
| : <math>\lim_{n \rarr \infty} p_{ij}^{(n)} = \frac{C}{M_j}.</math>
| |
| | |
| Note that there is no assumption on the starting distribution; the chain converges to the stationary distribution regardless of where it begins. Such ''<math>\pi</math>'' is called the '''equilibrium distribution''' of the chain.
| |
| | |
| If a chain has more than one closed communicating class, its stationary distributions will not be unique (consider any [[#Reducibility|closed communicating class]] <math>C_i</math> in the chain; each one will have its own unique stationary distribution <math>\pi_i</math>. Extending these distributions to the overall chain, setting all values to zero outside the communication class, yields that the set of invariant measures of the original chain is the set of all convex combinations of the <math>\pi_i</math>'s). However, if a state ''j'' is aperiodic, then
| |
| | |
| : <math>\lim_{n \rarr \infty} p_{jj}^{(n)} = \frac{C}{M_j}</math>
| |
| | |
| and for any other state ''i'', let ''f<sub>ij</sub>'' be the probability that the chain ever visits state ''j'' if it starts at ''i'',
| |
| | |
| : <math>\lim_{n \rarr \infty} p_{ij}^{(n)} = C \frac{f_{ij}}{M_j}.</math>
| |
| | |
| If a state ''i'' is periodic with period ''k'' > 1 then the limit
| |
| | |
| : <math>\lim_{n \rarr \infty} p_{ii}^{(n)}</math>
| |
| | |
| does not exist, although the limit
| |
| | |
| : <math>\lim_{n \rarr \infty} p_{ii}^{(kn+r)}</math>
| |
| | |
| does exist for every integer ''r''.
| |
| | |
| ====Steady-state analysis and the time-inhomogeneous Markov chain====
| |
| A Markov chain need not necessarily be time-homogeneous to have an equilibrium distribution. If there is a probability distribution over states <math>\boldsymbol{\pi}</math> such that
| |
| : <math>\pi_j = \sum_{i \in S} \pi_i \, \Pr(X_{n+1} = j \mid X_{n} = i)</math>
| |
| for every state ''j'' and every time ''n'' then <math>\boldsymbol{\pi}</math> is an equilibrium distribution of the Markov chain. Such can occur in [[Markov chain Monte Carlo|Markov chain Monte Carlo (MCMC)]] methods in situations where a number of different transition matrices are used, because each is efficient for a particular kind of mixing, but each matrix respects a shared equilibrium distribution.
| |
| | |
| ==Finite state space==
| |
| If the state space is [[finite set|finite]], the transition probability distribution can be represented by a [[matrix (mathematics)|matrix]], called the '''transition matrix''', with the (''i'', ''j'')th [[Element (mathematics)|element]] of '''P''' equal to
| |
| | |
| :<math>p_{ij} = \Pr(X_{n+1}=j\mid X_n=i). \,</math>
| |
| | |
| Since each row of '''P''' sums to one and all elements are non-negative, '''P''' is a [[right stochastic matrix]].
| |
| | |
| ===Stationary distribution relation to eigenvectors and simplexes===
| |
| The stationary distribution '''π''' is a (row) vector, whose entries are non-negative and sum to 1, is unchanged by the operation of transition matrix '''P''' on it and so is defined by
| |
| | |
| :<math> \pi\mathbf{P} = \pi.\,</math>
| |
| | |
| By comparing this definition with that of an [[eigenvector]] we see that the two concepts are related and that
| |
| | |
| :<math>\pi=\frac{e}{\sum_i{e_i}}</math>
| |
| | |
| is a normalized (<math>\textstyle \sum_i \pi_i=1</math>) multiple of a left eigenvector '''e''' of the transition matrix '''P''' with an [[eigenvalue]] of 1. If there is more than one unit eigenvector then a weighted sum of the corresponding stationary states is also a stationary state. But for a Markov chain one is usually more interested in the stationary state that is the limit of the sequence distributions for some initial distribution.
| |
| | |
| The values of stationary distribution <math> \textstyle \pi_i </math> are associated with the state space of '''P''' and its eigenvectors have their relative proportions preserved. Since the components of '''π''' are positive and the constraint that their sum is unity can be rewritten as <math>\textstyle \sum_i 1 \cdot \pi_i=1</math> we see that the dot product of '''π''' with a vector whose components are all 1 is unity and that '''π''' lies on a [[Simplex#The standard simplex|simplex]].
| |
| | |
| ===Time-homogeneous Markov chain with a finite state space===
| |
| {{no footnotes|section|date=February 2012}}
| |
| If the Markov chain is time-homogeneous, then the transition matrix '''P''' is the same after each step, so the ''k''-step transition probability can be computed as the ''k''-th power of the transition matrix, '''P'''<sup>''k''</sup>.
| |
| | |
| If the Markov chain is irreducible and aperiodic, then there is a unique stationary distribution '''π'''. Additionally, in this case '''P'''<sup>''k''</sup> converges to a rank-one matrix in which each row is the stationary distribution '''π''', that is,
| |
| | |
| :<math>\lim_{k\rightarrow\infty}\mathbf{P}^k=\mathbf{1}\pi</math>
| |
| | |
| where '''1''' is the column vector with all entries equal to 1. This is stated by the [[Perron–Frobenius theorem]]. If, by whatever means, <math>\scriptstyle \lim_{k\to\infty}\mathbf{P}^k</math> is found, then the stationary distribution of the Markov chain in question can be easily determined for any starting distribution, as will be explained below.
| |
| | |
| For some stochastic matrices '''P''', the limit <math>\scriptstyle \lim\limits_{k\to\infty}\mathbf{P}^k</math> does not exist while the stationary distribution does, as shown by this example:
| |
| | |
| : <math>\mathbf P=\begin{pmatrix} 0& 1\\ 1& 0 \end{pmatrix} \qquad \mathbf P^{2k}=I \qquad \mathbf P^{2k+1}=\mathbf P</math>
| |
| | |
| : <math>\begin{pmatrix}\frac{1}{2}&\frac{1}{2}\end{pmatrix}\begin{pmatrix} 0& 1\\ 1& 0 \end{pmatrix}=\begin{pmatrix}\frac{1}{2}&\frac{1}{2}\end{pmatrix}</math>
| |
| | |
| Because there are a number of different special cases to consider, the process of finding this limit if it exists can be a lengthy task. However, there are many techniques that can assist in finding this limit. Let '''P''' be an ''n''×''n'' matrix, and define <math>\scriptstyle \mathbf{Q} = \lim_{k\to\infty}\mathbf{P}^k.</math>
| |
| | |
| It is always true that
| |
| | |
| :<math>\mathbf{QP} = \mathbf{Q}.</math>
| |
| | |
| Subtracting '''Q''' from both sides and factoring then yields
| |
| | |
| :<math>\mathbf{Q}(\mathbf{P} - \mathbf{I}_{n}) = \mathbf{0}_{n,n} ,</math>
| |
| | |
| where '''I'''<sub>''n''</sub> is the [[identity matrix]] of size ''n'', and '''0'''<sub>''n'',''n''</sub> is the [[zero matrix]] of size ''n''×''n''. Multiplying together stochastic matrices always yields another stochastic matrix, so '''Q''' must be a [[stochastic matrix]] (see the definition above). It is sometimes sufficient to use the matrix equation above and the fact that '''Q''' is a stochastic matrix to solve for '''Q'''. Including the fact that the sum of each the rows in '''P''' is 1, there are ''n+1'' equations for determining ''n'' unknowns, so it is computationally easier if on the one hand one selects one row in '''Q''' and substitute each of its elements by one, and on the other one substitute the corresponding element (the one in the same column) in the vector '''0''', and next left-multiply this latter vector by the inverse of transformed former matrix to find '''Q'''.
| |
| | |
| Here is one method for doing so: first, define the function ''f''('''A''') to return the matrix '''A''' with its right-most column replaced with all 1's. If [''f''('''P''' − '''I'''<sub>n</sub>)]<sup>−1</sup> exists then{{Citation needed|date=March 2009}}
| |
| | |
| :<math>\mathbf{Q}=f(\mathbf{0}_{n,n})[f(\mathbf{P}-\mathbf{I}_n)]^{-1}.</math>
| |
| :Explain: The original matrix equation is equivalent to a [[System of linear equations|system of n×n linear equations]] in n×n variables. And there are n more linear equations from the fact that Q is a right [[stochastic matrix]] whose each row sums to 1. So it needs any n×n independent linear equations of the (n×n+n) equations to solve for the n×n variables. In this example, the n equations from “Q multiplied by the right-most column of (P-In)” have been replaced by the n stochastic ones.
| |
| One thing to notice is that if '''P''' has an element '''P'''<sub>''i'',''i''</sub> on its main diagonal that is equal to 1 and the ''i''th row or column is otherwise filled with 0's, then that row or column will remain unchanged in all of the subsequent powers '''P'''<sup>''k''</sup>. Hence, the ''i''th row or column of '''Q''' will have the 1 and the 0's in the same positions as in '''P'''.
| |
| | |
| ===Convergence speed to the stationary distribution===
| |
| {{No footnotes|section|date=February 2012}}
| |
| As stated earlier, from the equation <math> \mathbf{\pi} = \mathbf{\pi P} </math>, (if exists) the stationary (or steady state) distribution '''π''' is a left eigenvector of row [[stochastic matrix]] '''P'''. Then assuming that '''P''' is diagonalizable or equivalently that '''P''' has n linearly independent eigenvectors, speed of convergence is elaborated as follows. For non-diagonalizable matrices, one may start with "Jordan Canonical Form" (''almost'' diagonal form) of '''P''' and proceed with a bit more involved set of arguments in a similar way.<ref>Florian Schmitt and Franz Rothlauf, "On the Mean of the Second Largest Eigenvalue on the Convergence Rate of Genetic Algorithms", Working Paper 1/2001, Working Papers in Information Systems, 2001. http://citeseerx.ist.psu.edu/viewdoc/summary?doi=10.1.1.28.6191</ref>
| |
| | |
| Let '''U''' be the matrix of eigenvectors (each normalized to having an L2 norm equal to 1) where each column is a left eigenvector of '''P''' and let '''Σ''' be the diagonal matrix of left eigenvalues of '''P''', i.e. '''Σ''' = diag(λ<sub>1</sub>,λ<sub>2</sub>,λ<sub>3</sub>,...,λ<sub>n</sub>). Then by [[eigendecomposition]]
| |
| | |
| :<math> \mathbf{P} = \mathbf{U\Sigma U}^{-1} .</math>
| |
| | |
| Let the eigenvalues be enumerated such that 1=|λ<sub>1</sub>|>|λ<sub>2</sub>|≥|λ<sub>3</sub>|≥...≥|λ<sub>n</sub>|. Since '''P''' is a row stochastic matrix, its largest left eigenvalue is 1. If there is a unique stationary distribution, then the largest eigenvalue and the corresponding eigenvector is unique too (because there is no other '''π''' which solves the stationary distribution equation above). Let '''u'''<sub>''i''</sub> be the ''i''th column of '''U''' matrix, i.e. '''u'''<sub>''i''</sub> is the left eigenvector of '''P''' corresponding to λ<sub>''i''</sub>. Also let '''x''' be an arbitrary length n row vector in the span of the eigenvectors '''u'''<sub>''i''</sub>, that is
| |
| | |
| :<math> \mathbf{x}^T = \sum_{i=1}^n a_i \mathbf{u}_i </math>
| |
| | |
| for some set of a<sub>''i''</sub>∈ℝ. If we start multiplying '''P''' with '''x''' from left and continue this operation with the results, in the end we get the stationary distribution '''π'''. In other words '''π''' = '''u'''<sub>''i''</sub> ← '''xPPP'''...'''P''' = '''xP'''<sup>''k''</sup> as ''k'' goes to infinity. That means
| |
| | |
| :<math> \mathbf{\pi}^{(k)} = \mathbf{x}(\mathbf{U\Sigma U}^{-1})(\mathbf{U\Sigma U}^{-1})...(\mathbf{U\Sigma U}^{-1}) </math>
| |
| :<math> = \mathbf{xU\Sigma}^k \mathbf{U}^{-1} </math>
| |
| since '''UU'''<sup>''-1''</sup> = '''I''' the identity matrix and power of a diagonal matrix is also a diagonal matrix where each entry is taken to that power.
| |
| :<math> = (a_1\mathbf{u}_1^T + a_2\mathbf{u}_2^T + ... + a_n\mathbf{u}_n^T)\mathbf{U\Sigma}^k\mathbf{U}^{-1} ,</math>
| |
| :<math> = a_1\lambda_1^k\mathbf{u}_1 + a_2\lambda_2^k\mathbf{u}_2 + ... + a_n\lambda_n^k\mathbf{u}_n ,</math>
| |
| since the eigenvectors are [[orthonormal]]. Then<ref>Gene H. Golub, Charles F. Van Loan, "Matrix computations", Third Edition, The Johns Hopkins University Press, Baltimore and London, 1996.</ref>
| |
| :<math> = \lambda_1^k\left\{a_1\mathbf{u}_1 + a_2\left(\frac{\lambda_2}{\lambda_1}\right)^k\mathbf{u}_2 + a_3\left(\frac{\lambda_3}{\lambda_1}\right)^k\mathbf{u}_3 + ... + a_n\left(\frac{\lambda_n}{\lambda_1}\right)^k\mathbf{u}_n\right\} .</math>
| |
| | |
| Since '''π''' = '''u'''<sub>1</sub>, '''π'''<sup>(k)</sup> approaches to '''π''' as ''k'' goes to infinity with a speed in the order of λ<sub>2</sub>/λ<sub>1</sub> exponentially. This follows because |λ<sub>2</sub>|≥|λ<sub>3</sub>|≥...≥|λ<sub>n</sub>|, hence λ<sub>2</sub>/λ<sub>1</sub> is the dominant term. Random noise in the state distribution '''π''' can also speed up this convergence to the stationary distribution.<ref>{{cite journal|last=Franzke|first=Brandon|coauthors=Kosko, Bart|title=Noise can speed convergence in Markov chains|journal=Physical Review E|date=1 October 2011|volume=84|issue=4|doi=10.1103/PhysRevE.84.041112}}</ref>
| |
| | |
| ==Reversible Markov chain{{Anchor|detailed balance}}==
| |
| | |
| A Markov chain is said to be '''reversible''' if there is a probability distribution over states, '''π''', such that
| |
| :<math>\pi_i \Pr(X_{n+1} = j \mid X_{n} = i) = \pi_j \Pr(X_{n+1} = i \mid X_{n} = j)</math>
| |
| for all times ''n'' and all states ''i'' and ''j''.
| |
| This condition is also known as the '''[[detailed balance]]''' condition (some books refer the local balance equation).
| |
| With a time-homogeneous Markov chain, Pr(''X''<sub>n+1</sub> = ''j'' | ''X''<sub>n</sub> = ''i'') does not change with time ''n'' and it can be written more simply as <math>p_{ij}</math>. In this case, the detailed balance equation can be written more compactly as
| |
| :<math>\pi_i p_{ij} = \pi_j p_{ji}\,.</math>
| |
| | |
| Summing the original equation over ''i'' gives
| |
| | |
| :<math>\begin{align}\sum_i \pi_i \Pr(X_{n+1} = j \mid X_{n} = i) &= \sum_i \pi_j \Pr(X_{n+1} = i \mid X_{n} = j) \\ &= \pi_j \sum_i \Pr(X_{n+1} = i \mid X_{n} = j) = \pi_j\,,\end{align}</math>
| |
| so, for reversible Markov chains, '''π''' is always a steady-state distribution of Pr(''X''<sub>n+1</sub> = ''j'' | ''X''<sub>n</sub> = ''i'') for every ''n''.
| |
| | |
| If the Markov chain begins in the steady-state distribution, ''i.e.'', if Pr(''X''<sub>0</sub> = ''i'') = π<sub>''i''</sub>, then Pr(''X''<sub>''n''</sub> = ''i'') = π<sub>''i''</sub> for all ''n'' and the detailed balance equation can be written as
| |
| :<math>\Pr(X_{n} = i, X_{n+1} = j) = \Pr(X_{n+1} = i, X_{n} = j)\,.</math>
| |
| The left- and right-hand sides of this last equation are identical except for a reversing of the time indices ''n'' and ''n'' + 1.
| |
| | |
| [[Kolmogorov's criterion]] gives a necessary and sufficient condition for a Markov chain to be reversible directly from the transition matrix probabilities. The criterion requires that the products of probabilities around every closed loop are the same in both directions around the loop.
| |
| | |
| Reversible Markov chains are common in [[Markov chain Monte Carlo|Markov chain Monte Carlo (MCMC)]] approaches because the detailed balance equation for a desired distribution '''π''' necessarily implies that the Markov chain has been constructed so that '''π''' is a steady-state distribution. Even with time-inhomogeneous Markov chains, where multiple transition matrices are used, if each such transition matrix exhibits detailed balance with the desired '''π''' distribution, this necessarily implies that '''π''' is a steady-state distribution of the Markov chain.
| |
| | |
| ==Bernoulli scheme==
| |
| | |
| A [[Bernoulli scheme]] is a special case of a Markov chain where the transition probability matrix has identical rows, which means that the next state is even independent of the current state (in addition to being independent of the past states). A Bernoulli scheme with only two possible states is known as a [[Bernoulli process]].
| |
| | |
| ==General state space==
| |
| | |
| Many results for Markov chains with finite state space can be generalized to chains with uncountable state space through [[Harris chain]]s. The main idea is to see if there is a point in the state space that the chain hits with probability one. Generally, it is not true for continuous state space, however, we can define sets ''A'' and ''B'' along with a positive number ''ε'' and a probability
| |
| measure ''ρ'', such that
| |
| | |
| # <math>\text{If }\tau_A = \inf\{n\geq 0: X_n \in A\},\text{ then } P_z(\tau_A<\infty)>0\text{ for all }z.</math>
| |
| # <math>\text{If }x \in A\text{ and }C\subset B,\text{ then } p(x, C)\geq \varepsilon \rho(C).</math>
| |
| | |
| Then we could collapse the sets into an auxiliary point ''α'', and a recurrent [[Harris chain]] can be modified to contain ''α''. Lastly, the collection of [[Harris chain]]s is a comfortable level of generality, which is broad enough to contain a large number of interesting examples, yet restrictive enough to allow for a rich theory.
| |
| | |
| The use of Markov chains in [[Markov chain Monte Carlo]] methods covers cases where the process follows a continuous state space.
| |
| | |
| ===Locally interacting Markov chains===
| |
| | |
| Considering a collection of Markov chains whose evolution takes in account the state of other Markov chains, is related to the notion
| |
| of [[locally interacting Markov chains]]. This corresponds to the situation when the state space has a (Cartesian-) product form.
| |
| See [[interacting particle system]] and [[stochastic cellular automata]].
| |
| See for instance ''Interaction of Markov Processes''<ref>{{cite journal|last=Spitzer|first=Frank|title=Interaction of Markov Processes|journal=Advances in Mathematics|year=1970|volume=5|pages=246–290|doi=10.1016/0001-8708(70)90034-4|issue=2}}</ref>
| |
| or<ref>{{cite book|title=Stochastic Cellular Systems: Ergodicity, Memory, Morphogenesis|author= R. L. Dobrushin, V. I. Kri︠u︡kov, A. L. Toom|year=1978|url=http://books.google.com/?id=0Wa7AAAAIAAJ&pg=PA181&lpg=PA181&dq=locally+interacting+markov+chains+toom+Dobrushin#v=onepage&q=locally%20interacting%20markov%20chains%20toom%20Dobrushin&f=false|isbn=9780719022067}}</ref>
| |
| | |
| ==Applications==
| |
| {{Unreferenced section|date=May 2012}}
| |
| Research has reported the application and usefulness of Markov chains in a wide range of topics such as physics, chemistry, medicine, music, game theory and sports.
| |
| | |
| ===Physics===
| |
| Markovian systems appear extensively in [[thermodynamics]] and [[statistical mechanics]], whenever probabilities are used to represent unknown or unmodelled details of the system, if it can be assumed that the dynamics are time-invariant, and that no relevant history need be considered which is not already included in the state description.{{Citation needed|date=March 2012}}
| |
| | |
| ===Chemistry===
| |
| [[File:Simple mechanism.svg|thumb|left|200px|[[Michaelis-Menten kinetics]]. The enzyme (E) binds a substrate (S) and produces a product (P). Each reaction is a state transition in a Markov chain.]]
| |
| Chemistry is often a place where Markov chains and continuous-time Markov processes are especially useful because these simple physical systems tend to satisfy the Markov property quite well. The classical model of enzyme activity, [[Michaelis-Menten kinetics]], can be viewed as a Markov chain, where at each time step the reaction proceeds in some direction. While Michaelis-Menten is fairly straightforward, far more complicated reaction networks can also be modeled with Markov chains.
| |
| | |
| An algorithm based on a Markov chain was also used to focus the fragment-based growth of chemicals in silico towards a desired class of compounds such as drugs or natural products.<ref>{{cite journal|title=FOG: Fragment Optimized Growth Algorithm for the de Novo Generation of Molecules occupying Druglike Chemical | last=Kutchukian | first=Peter | coauthors=Lou, David; Shakhnovich, Eugene |journal=Journal of Chemical Information and Modeling | year=2009 |volume=49 | pages=1630–1642|doi=10.1021/ci9000458|pmid=19527020|issue=7 }}</ref> As a molecule is grown, a fragment is selected from the nascent molecule as the "current" state. It is not aware of its past (i.e., it is not aware of what is already bonded to it). It then transitions to the next state when a fragment is attached to it. The transition probabilities are trained on databases of authentic classes of compounds.
| |
| | |
| Also, the growth (and composition) of [[copolymer]]s may be modeled using Markov chains. Based on the reactivity ratios of the monomers that make up the growing polymer chain, the chain's composition may be calculated (e.g., whether monomers tend to add in alternating fashion or in long runs of the same monomer). Due to [[steric effects]], second-order Markov effects may also play a role in the growth of some polymer chains.
| |
| | |
| Similarly, it has been suggested that the crystallization and growth of some epitaxial [[superlattice]] oxide materials can be accurately described by Markov chains.<ref>{{cite doi|10.1107/S0108767311044874}}</ref>
| |
| | |
| ===Testing===
| |
| Several theorists have proposed the idea of the Markov chain statistical test (MCST), a method of conjoining Markov chains to form a "[[Markov blanket]]", arranging these chains in several recursive layers ("wafering") and producing more efficient test sets—samples—as a replacement for exhaustive testing. MCSTs also have uses in temporal state-based networks; Chilukuri et al.'s paper entitled "Temporal Uncertainty Reasoning Networks for Evidence Fusion with Applications to Object Detection and Tracking" (ScienceDirect) gives a background and case study for applying MCSTs to a wider range of applications.
| |
| | |
| ===Information sciences===
| |
| Markov chains are used throughout information processing. [[Claude Shannon|Claude Shannon's]] famous 1948 paper ''[[A mathematical theory of communication]]'', which in a single step created the field of [[information theory]], opens by introducing the concept of [[information entropy|entropy]] through Markov modeling of the English language. Such idealized models can capture many of the statistical regularities of systems. Even without describing the full structure of the system perfectly, such signal models can make possible very effective [[data compression]] through [[entropy encoding]] techniques such as [[arithmetic coding]]. They also allow effective [[state estimation]] and [[pattern recognition]]. Markov chains also play an important role in [[reinforcement learning]].
| |
| | |
| Markov chains are also the basis for [[Hidden Markov Models]], which are an important tool in such diverse fields as telephone networks (which use the [[Viterbi algorithm]] for error correction), speech recognition and [[bioinformatics]].
| |
| | |
| ===Queueing theory===
| |
| Markov chains are the basis for the analytical treatment of queues ([[queueing theory]]). [[Agner Krarup Erlang]] initiated the subject in 1917.<ref name="MacTutor Biography|id=Erlang">{{MacTutor Biography|id=Erlang}}</ref> This makes them critical for optimizing the performance of telecommunications networks, where messages must often compete for limited resources (such as bandwidth).<ref name=CTCN/>
| |
| | |
| ===Internet applications===
| |
| The [[PageRank]] of a webpage as used by [[Google]] is defined by a Markov chain.<ref>{{US patent|6285999}}</ref> It is the probability to be at page <math>i</math> in the stationary distribution on the following Markov chain on all (known) webpages. If <math>N</math> is the number of known webpages, and a page <math>i</math> has <math>k_i</math> links to it then it has transition probability <math>\frac{\alpha}{k_i} + \frac{1-\alpha}{N}</math> for all pages that are linked to and <math>\frac{1-\alpha}{N}</math> for all pages that are not linked to. The parameter <math>\alpha</math> is taken to be about 0.85.<ref name="pagerank">{{cite techreport |
| |
| author = Page, Lawrence and Brin, Sergey and Motwani, Rajeev and Winograd, Terry | title = The PageRank Citation Ranking: Bringing Order to the Web | url = http://citeseerx.ist.psu.edu/viewdoc/summary?doi=10.1.1.31.1768 | year = 1999}}</ref>
| |
| | |
| Markov models have also been used to analyze web navigation behavior of users. A user's web link transition on a particular website can be modeled using first- or second-order Markov models and can be used to make predictions regarding future navigation and to personalize the web page for an individual user.
| |
| | |
| ===Statistics===
| |
| Markov chain methods have also become very important for generating sequences of random numbers to accurately reflect very complicated desired probability distributions, via a process called [[Markov chain Monte Carlo]] (MCMC). In recent years this has revolutionized the practicability of [[Bayesian inference]] methods, allowing a wide range of [[posterior distribution]]s to be simulated and their parameters found numerically.
| |
| | |
| ===Economics and finance===
| |
| Markov chains are used in Finance and Economics to model a variety of different phenomena, including asset prices and market crashes. The first financial model to use a Markov chain was from Prasad ''et al.'' in 1974.<ref>{{cite journal|title=Allocation of resources on a minimized cost basis | last=Prasad | first=NR | coauthors=RC Ender,ST Reilly,G Nesgos |journal=1974 IEEE Conference on Decision and Control including the 13th Symposium on Adaptive Processes | year=1974 |volume=13 | pages=402–3 |doi=10.1109/CDC.1974.270470 |url= http://ieeexplore.ieee.org/Xplore/defdeny.jsp?url=http://ieeexplore.ieee.org/stamp/stamp.jsp%3Ftp%3D%26arnumber%3D4045263%26userType%3Dinst&denyReason=-133&arnumber=4045263&productsMatched=null&userType=inst }}</ref> Another was the regime-switching model of [[James D. Hamilton]] (1989), in which a Markov chain is used to model switches between periods of high volatility and low volatility of asset returns.<ref>{{cite journal|title=A new approach to the economic analysis of nonstationary time series and the business cycle | last=Hamilton | first=James | coauthors= |journal=Econometrica | year=1989 |volume=57 | pages=357–84|doi=10.2307/1912559|jstor=1912559|issue=2|publisher=Econometrica, Vol. 57, No. 2 }}</ref> A more recent example is the [[Markov Switching Multifractal]] asset pricing model, which builds upon the convenience of earlier regime-switching models.<ref>{{cite journal|title=How to Forecast long-run volatility: regime-switching and the estimation of multifractal processes | last=Calvet | first=Laurent | coauthors=Adlai Fisher |journal=Journal of Financial Econometrics | year=2004 |volume=2 | pages=49–83|doi=10.1093/jjfinec/nbh003 }}</ref> It uses an arbitrarily large Markov chain to drive the level of volatility of asset returns.
| |
| | |
| Dynamic macroeconomics heavily uses Markov chains. An example is using Markov chains to exogenously model prices of equity (stock) in a [[general equilibrium]] setting.<ref>{{cite web |last=Brennan |first=Michael |first2=Yihong |last2=Xiab |title=Stock Price Volatility and the Equity Premium |work=Department of Finance, the Anderson School of Management, UCLA |url=http://bbs.cenet.org.cn/uploadImages/200352118122167693.pdf }}</ref>
| |
| | |
| ===Social sciences===
| |
| Markov chains are generally used in describing [[path-dependent]] arguments, where current structural configurations condition future outcomes. An example is the reformulation of the idea, originally due to [[Karl Marx]]'s [[Das Kapital]], tying [[economic development]] to the rise of [[capitalism]]. In current research, it is common to use a Markov chain <!-- this is actually a Markov perefect equilibria, not simply a Markov chain, I'll try to remember get back to this ~~~~ --> to model how once a country reaches a specific level of economic development, the configuration of structural factors, such as size of the commercial [[bourgeoisie]], the ratio of urban to rural residence, the rate of [[political]] mobilization, etc., will generate a higher probability of transitioning from [[authoritarian]] to [[democratic]] regime.<ref>{{cite journal | last = Acemoglu | first = Daron | coauthors = Georgy Egorov, Konstantin Sonin| title = Political model of social evolution | journal = Proceedings of the National Academy of Sciences | year = 2011| volume = 108 | pages = 21292–21296 | url = http://www.pnas.org/content/108/suppl.4/21292.short | doi = 10.1073/pnas.1019454108 }}</ref>
| |
| | |
| ===Mathematical biology===
| |
| | |
| Markov chains also have many applications in biological modelling, particularly [[population process]]es, which are useful in modelling processes that are (at least) analogous to biological populations. The [[Leslie matrix]] is one such example, though some of its entries
| |
| are not probabilities (they may be greater than 1). Another example is the modeling of cell shape in dividing sheets of [[epithelial cells]].<ref>{{cite journal | last=Gibson | first=Matthew C |first2=Ankit P. |last2=Patel |first3=Norbert |last3=Perrimon |year=2006 |title=The emergence of geometric order in proliferating metazoan epithelia |journal=Nature |volume=442 |pages=1038–1041 |doi=10.1038/nature05014 | last4=Perrimon | first4=Norbert | issue=7106 | pmid=16900102}}</ref> Yet another example is the state of [[ion channel]]s in cell membranes.
| |
| | |
| Markov chains are also used in simulations of brain function, such as the simulation of the mammalian neocortex.<ref>{{cite journal |last=George |first=Dileep |first2=Jeff |last2=Hawkins |year=2009 |title=Towards a Mathematical Theory of Cortical Micro-circuits |journal=PLoS Comput Biol |volume=5 |issue=10 |pages=e1000532 |doi=10.1371/journal.pcbi.1000532 |editor1-last=Friston |editor1-first=Karl J. |pmid=19816557 |pmc=2749218 }}</ref>
| |
| | |
| ===Genetics===
| |
| Markov chains have been used in [[population genetics]] in order to describe the change in gene frequencies in small populations affected by [[genetic drift]], for example in [[diffusion equation]] method described by [[Motoo Kimura]].{{citation needed|date=April 2013}} Another application is Using a VOM Model for Reconstructing Potential Coding Regions in EST Sequences <ref>{{cite paper |url= http://www.eng.tau.ac.il/~bengal/VOM_EST.pdf|title= Using specific Markov models for Reconstructing Potential Coding Regions in EST Sequences, |author= Shmilovici A. and Ben-Gal I.|publisher= Journal of Computational Statistics, vol. 22, no. 1, 49 - 69.|year=2007}}</ref> [http://www.eng.tau.ac.il/~bengal/VOM_EST.pdf]
| |
| | |
| ===Games===
| |
| Markov chains can be used to model many games of chance. The children's games [[Snakes and Ladders]] and "[[Hi Ho! Cherry-O]]", for example, are represented exactly by Markov chains. At each turn, the player starts in a given state (on a given square) and from there has fixed odds of moving to certain other states (squares).
| |
| | |
| ===Music===
| |
| Markov chains are employed in [[algorithmic composition|algorithmic music composition]], particularly in [[software]] programs such as [[CSound]], [[Max (software)|Max]] or [[SuperCollider]]. In a first-order chain, the states of the system become note or pitch values, and a [[probability vector]] for each note is constructed, completing a transition probability matrix (see below). An algorithm is constructed to produce and output note values based on the transition matrix weightings, which could be [[MIDI]] note values, frequency ([[Hertz|Hz]]), or any other desirable metric.<ref>{{cite journal |title=Making Music with Algorithms: A Case-Study System |author=K McAlpine, E Miranda, S Hoggar |url=http://www.mitpressjournals.org/doi/abs/10.1162/014892699559733 |journal=Computer Music Journal |issue=2 |year=1999 |volume=23 |doi=10.1162/014892699559733 |pages=19 }}</ref>
| |
| | |
| {| class="wikitable" style="float: left"
| |
| |+ 1st-order matrix
| |
| ! Note !! A !! C♯ !! E♭
| |
| |-
| |
| ! A
| |
| | 0.1 || 0.6 || 0.3
| |
| |-
| |
| ! C♯
| |
| | 0.25 || 0.05 || 0.7
| |
| |-
| |
| ! E♭
| |
| | 0.7 || 0.3 || 0
| |
| |}
| |
| | |
| {| class="wikitable" style="float: left"
| |
| |+ 2nd-order matrix
| |
| ! Note !! A !! D !! G
| |
| |-
| |
| ! AA
| |
| | 0.18 || 0.6 || 0.22
| |
| |-
| |
| ! AD
| |
| | 0.5 || 0.5 || 0
| |
| |-
| |
| ! AG
| |
| | 0.15 || 0.75 || 0.1
| |
| |-
| |
| ! DD
| |
| | 0 || 0 || 1
| |
| |-
| |
| ! DA
| |
| | 0.25 || 0 || 0.75
| |
| |-
| |
| ! DG
| |
| | 0.9 || 0.1 || 0
| |
| |-
| |
| ! GG
| |
| | 0.4 || 0.4 || 0.2
| |
| |-
| |
| ! GA
| |
| | 0.5 || 0.25 || 0.25
| |
| |-
| |
| ! GD
| |
| | 1 || 0 || 0
| |
| |-
| |
| |}
| |
| {{Clear}}
| |
| | |
| A second-order Markov chain can be introduced by considering the current state ''and'' also the previous state, as indicated in the second table. Higher, ''n''th-order chains tend to "group" particular notes together, while 'breaking off' into other patterns and sequences occasionally. These higher-order chains tend to generate results with a sense of [[Phrase (music)|phrasal]] structure, rather than the 'aimless wandering' produced by a first-order system.<ref name=Roads>{{cite book|author=Curtis Roads (ed.)|title=The Computer Music Tutorial| year=1996|publisher=MIT Press|isbn= 0-262-18158-4 }}</ref>
| |
| | |
| Markov chains can be used structurally, as in Xenakis's Analogique A and B.<ref>Xenakis, Iannis; Kanach, Sharon (1992) ''Formalized Music: Mathematics and Thought in Composition'', Pendragon Press. ISBN 1576470792</ref> Markov chains are also used in systems which use a Markov model to react interactively to music input.<ref>[http://www.csl.sony.fr/~pachet/ Continuator]{{dead link|date=April 2013}}</ref>
| |
| | |
| Usually musical systems need to enforce specific control constraints on the finite-length sequences they generate, but control constraints are not compatible with Markov models, since they induce long-range dependencies that violate the Markov hypothesis of limited memory. In order to overcome this limitation, a new approach has been proposed.<ref>Pachet, F.; Roy, P.; Barbieri, G. (2011) [http://www.csl.sony.fr/downloads/papers/2011/pachet-11b.pdf "Finite-Length Markov Processes with Constraints"], ''Proceedings of the 22nd International Joint Conference on Artificial Intelligence'', IJCAI, pages 635-642,Barcelona, Spain, July 2011</ref>
| |
| | |
| ===Baseball===
| |
| Markov chain models have been used in advanced baseball analysis since 1960, although their use is still rare. Each half-inning of a baseball game fits the Markov chain state when the number of runners and outs are considered. During any at-bat, there are 24 possible combinations of number of outs and position of the runners. Mark Pankin shows that Markov chain models can be used to evaluate runs created for both individual players as well as a team.<ref>{{cite web| last=Pankin |first=Mark D. |title=MARKOV CHAIN MODELS: THEORETICAL BACKGROUND| url=http://www.pankin.com/markov/theory.htm |accessdate = 2007-11-26}}
| |
| </ref>
| |
| He also discusses various kinds of strategies and play conditions: how Markov chain models have been used to analyze statistics for game situations such as [[bunt (baseball)|bunting]] and [[base stealing]] and differences when playing on grass vs. [[astroturf]].<ref>{{cite web| last=Pankin |first=Mark D. |title=BASEBALL AS A MARKOV CHAIN| url=http://www.pankin.com/markov/intro.htm |accessdate = 2009-04-24}}
| |
| </ref>
| |
| | |
| ===Markov text generators===
| |
| Markov processes can also be used to generate superficially real-looking text given a sample document: they are used in a variety of recreational "[[parody generator]]" software (see [[dissociated press]], Jeff Harrison,<ref>[http://www.fieralingue.it/modules.php?name=Content&pa=list_pages_categories&cid=111 Poet's Corner - Fieralingue]</ref> [[Mark V Shaney]]<ref name="Travesty">
| |
| {{Cite journal
| |
| | last1 = Kenner
| |
| | first1 = Hugh
| |
| | last2 = O'Rourke
| |
| | first2 = Joseph | authorlink2 = Joseph O'Rourke (professor)
| |
| | title = A Travesty Generator for Micros
| |
| | date = November 1984
| |
| | journal = BYTE
| |
| | volume = 9
| |
| | issue = 12
| |
| | pages = 129–131, 449–469
| |
| | postscript = <!--None-->}}
| |
| </ref><ref name="Hartman">
| |
| {{Cite book
| |
| | last = Hartman
| |
| | first = Charles
| |
| | title = Virtual Muse: Experiments in Computer Poetry
| |
| | place= Hanover, NH
| |
| | publisher = Wesleyan University Press
| |
| | year = 1996
| |
| | isbn = 0-8195-2239-2
| |
| | postscript = <!--None-->}}
| |
| </ref>
| |
| ).
| |
| | |
| These processes are also used by [[E-mail spam|spammers]] to inject real-looking hidden paragraphs into [[unsolicited email]] and [[spam in blogs|post comments]] in an attempt to get these messages past [[spam filter]]s.
| |
| | |
| ==Fitting==
| |
| | |
| When fitting a Markov chain to data, situations where parameters poorly describe the situation may highlight interesting trends.<ref>{{cite jstor|2680818}}</ref><ref>{{cite paper |url= http://www.eng.tau.ac.il/~bengal/VOM_EST.pdf|title= Using a VOM Model for Reconstructing Potential Coding Regions in EST Sequences, |author= Shmilovici A. and Ben-Gal I.|publisher= Journal of Computational Statistics, vol. 22, no. 1, 49 - 69.|year=2007}}</ref> [http://www.eng.tau.ac.il/~bengal/VOM_EST.pdf]
| |
| | |
| ==History==
| |
| [[Andrey Markov]] produced the first results (1906) for these processes, purely theoretically.
| |
| A generalization to countably infinite state spaces was given by [[Andrey Nikolaevich Kolmogorov|Kolmogorov]] (1936).
| |
| Markov chains are related to [[Brownian motion]] and the [[ergodic hypothesis]], two topics in physics which were important in the early years of the twentieth century. However, Markov first pursued this in 1906 as part of his argument against Pavel Nekrasov, in particular to make the case that the [[law of large numbers]] can be extended to dependent events.<ref name="Hayes 2013 92–97">{{cite journal |last=Hayes |first=Brian |year=2013 |month=March–April |title=First Links in the Markov Chain |journal=American Scientist |volume=101 |pages=92–97}}</ref> In 1913, he applied his findings to the first 20,000 letters of Pushkin's ''[[Eugene Onegin]]''.<ref name="Hayes 2013 92–97"/> By 1917, more practical application of his work was made by [[Agner Krarup Erlang|Erlang]] to obtain formulas for call loss and waiting time in telephone networks.<ref name="MacTutor Biography|id=Erlang"/>
| |
| | |
| Seneta provides an account of Markov's motivations and the theory's early development.<ref>{{cite journal |last=Seneta |first=E. |year=1996 |title=Markov and the Birth of Chain Dependence Theory |journal=International Statistical Review |volume=64 |issue=3 |pages=255–263 |jstor=1403785 |doi=10.2307/1403785 }}</ref> The term "chain" was used by Markov (1906) to suggest a sequence of pairwise dependent variables.<ref>{{cite book |last=Upton |first=G. |last2=Cook |first2=I. |year=2008 |title=Oxford Dictionary of Statistics |location= |publisher=OUP |isbn=978-0-19-954145-4 }}</ref>
| |
| | |
| ==See also==
| |
| {{refbegin|4}}
| |
| * [[Hidden Markov model]]
| |
| * [[Markov blanket]]
| |
| * [[Markov chain geostatistics]]
| |
| * [[Markov chain mixing time]]
| |
| * [[Markov chain Monte Carlo]]
| |
| * [[Markov decision process]]
| |
| * [[Markov information source]]
| |
| * [[Markov network]]
| |
| * [[Markov process]]
| |
| * [[Quantum Markov chain]]
| |
| * [[Semi-Markov process]]
| |
| * [[Telescoping Markov chain]]
| |
| * [[Variable-order Markov model]]
| |
| {{refend}}
| |
| | |
| ==Notes==
| |
| {{Reflist|30em}}
| |
| | |
| ==References==
| |
| {{Refbegin}}
| |
| * A.A. Markov. "Rasprostranenie zakona bol'shih chisel na velichiny, zavisyaschie drug ot druga". ''Izvestiya Fiziko-matematicheskogo obschestva pri Kazanskom universitete'', 2-ya seriya, tom 15, pp. 135–156, 1906.
| |
| * A.A. Markov. "Extension of the limit theorems of probability theory to a sum of variables connected in a chain". reprinted in Appendix B of: R. Howard. ''Dynamic Probabilistic Systems, volume 1: Markov Chains''. John Wiley and Sons, 1971.
| |
| * Classical Text in Translation: A. A. Markov, An Example of Statistical Investigation of the Text Eugene Onegin Concerning the Connection of Samples in Chains, trans. David Link. Science in Context 19.4 (2006): 591–600. Online: http://journals.cambridge.org/production/action/cjoGetFulltext?fulltextid=637500
| |
| * Leo Breiman. ''Probability''. Original edition published by Addison-Wesley, 1968; reprinted by Society for Industrial and Applied Mathematics, 1992. ISBN 0-89871-296-3. ''(See Chapter 7.)''
| |
| * J.L. Doob. ''Stochastic Processes''. New York: John Wiley and Sons, 1953. ISBN 0-471-52369-0.
| |
| * S. P. Meyn and R. L. Tweedie. ''Markov Chains and Stochastic Stability''. London: Springer-Verlag, 1993. ISBN 0-387-19832-6. online: https://netfiles.uiuc.edu/meyn/www/spm_files/book.html . Second edition to appear, Cambridge University Press, 2009.
| |
| * S. P. Meyn. ''Control Techniques for Complex Networks''. Cambridge University Press, 2007. ISBN 978-0-521-88441-9. Appendix contains abridged Meyn & Tweedie. online: https://netfiles.uiuc.edu/meyn/www/spm_files/CTCN/CTCN.html
| |
| * {{Cite book |last=Booth| first=Taylor L. |coauthors=| title=Sequential Machines and Automata Theory |edition=1st |publisher=John Wiley and Sons, Inc. |location=New York |year=1967| id=Library of Congress Card Catalog Number 67-25924}} Extensive, wide-ranging book meant for specialists, written for both theoretical computer scientists as well as electrical engineers. With detailed explanations of state minimization techniques, FSMs, Turing machines, Markov processes, and undecidability. Excellent treatment of Markov processes pp. 449ff. Discusses Z-transforms, D transforms in their context.
| |
| * {{Cite book |last=Kemeny| first=John G. |coauthors=Hazleton Mirkil, J. Laurie Snell, Gerald L. Thompson |title=Finite Mathematical Structures| edition=1st |publisher=Prentice-Hall, Inc. |location=Englewood Cliffs, N.J. |year=1959| id=Library of Congress Card Catalog Number 59-12841}} Classical text. cf Chapter 6 ''Finite Markov Chains'' pp. 384ff.
| |
| * E. Nummelin. "General irreducible Markov chains and non-negative operators". Cambridge University Press, 1984, 2004. ISBN 0-521-60494-X
| |
| * Seneta, E. ''Non-negative matrices and Markov chains''. 2nd rev. ed., 1981, XVI, 288 p., Softcover Springer Series in Statistics. (Originally published by Allen & Unwin Ltd., London, 1973) ISBN 978-0-387-29765-1
| |
| * Kishor S. Trivedi, ''Probability and Statistics with Reliability, Queueing, and Computer Science Applications'', John Wiley & Sons, Inc. New York, 2002. ISBN 0-471-33341-7.
| |
| * K.S.Trivedi and R.A.Sahner, ''SHARPE at the age of twenty-two'', vol. 36, no. 4, pp.-52-57, ACM SIGMETRICS Performance Evaluation Review, 2009.
| |
| * R.A.Sahner, K.S.Trivedi and A. Puliafito, ''Performance and reliability analysis of computer systems: an example-based approach using the SHARPE software package'', Kluwer Academic Publishers, 1996. ISBN 0-7923-9650-2.
| |
| * G.Bolch, S.Greiner, H.de Meer and K.S.Trivedi, ''Queueing Networks and Markov Chains'', John Wiley, 2nd edition, 2006. ISBN 978-0-7923-9650-5.
| |
| {{Refend}}
| |
| | |
| ==External links==
| |
| {{refbegin}}
| |
| * {{YouTube|o-jdJxXL_W4|Introduction to Markov Chains}}
| |
| * {{springer|title=Markov chain|id=p/m062350}}
| |
| * [http://jasss.soc.surrey.ac.uk/12/1/6.html Techniques to Understand Computer Simulations: Markov Chain Analysis]
| |
| * [http://www.dartmouth.edu/~chance/teaching_aids/books_articles/probability_book/Chapter11.pdf Markov Chains chapter in American Mathematical Society's introductory probability book](pdf)
| |
| *[http://www.math.rutgers.edu/courses/338/coursenotes/chapter5.pdf Chapter 5: Markov Chain Models]
| |
| {{refend}}
| |
| | |
| {{Stochastic processes}}
| |
| | |
| {{DEFAULTSORT:Markov Chain}}
| |
| [[Category:Markov processes]]
| |
| [[Category:Markov models]]
| |
| [[Category:Graph theory]]
| |