Markov chain
This article may be too long to read and navigate comfortably. The readable prose size is 74 kilobytes. (February 2017)

A Markov chain is a stochastic model describing a sequence of possible events in which the probability of each event depends only on the state attained in the previous event.^{[1]}^{[2]}^{[3]} In continuoustime, it is known as a Markov process. It is named after the Russian mathematician Andrey Markov.
Markov chains have many applications as statistical models of realworld processes,^{[1]}^{[4]}^{[5]}^{[6]} such as studying cruise control systems in motor vehicles, queues or lines of customers arriving at an airport, exchange rates of currencies, and population growths of certain animal species.^{[7]}
Markov processes are the basis for general stochastic simulation methods known as Markov chain Monte Carlo, which are used for simulating sampling from complex probability distributions, and have found extensive application in Bayesian statistics.^{[7]}^{[8]}^{[9]}
The adjective Markovian is used to describe something that is related to a Markov process.^{[1]}^{[10]}
Contents
 1 Introduction
 2 History
 3 Examples
 4 Formal definition

5 Properties
 5.1 Reducibility
 5.2 Periodicity
 5.3 Transience and recurrence
 5.4 Ergodicity
 5.5 Steadystate analysis and limiting distributions
 5.6 Reversible Markov chain
 5.7 Markovian representations
 5.8 Transient evolution
 5.9 Transient behaviour
 5.10 Stationary distribution
 5.11 Hitting times
 5.12 Time reversal
 5.13 Embedded Markov chain
 5.14
 6 Special types of Markov chains

7 Applications
 7.1 Physics
 7.2 Chemistry
 7.3 Testing
 7.4 Solar irradiance variability
 7.5 Speech recognition
 7.6 Information and computer science
 7.7 Queueing theory
 7.8 Internet applications
 7.9 Statistics
 7.10 Economics and finance
 7.11 Social sciences
 7.12 Mathematical biology
 7.13 Genetics
 7.14 Games
 7.15 Music
 7.16 Baseball
 7.17 Markov text generators
 7.18 Bioinformatics
 7.19 Probabilistic forecasting
 8 See also
 9 Notes
 10 References
 11 External links
Introduction[edit]
Definition[edit]
A Markov process is a stochastic process that satisfies the Markov property^{[1]} (sometimes characterized as "memorylessness"). In simpler terms, a Markov process is a process for which one can make for its future based solely on its present state just as well as one could knowing the process's full history.^{[11]} In other words, conditional on the present state of the system, its future and past states are independent.
A Markov chain is a type of Markov process that has either a discrete state space or a discrete index set (often representing time), but the precise definition of a Markov chain varies.^{[12]} For example, it is common to define a Markov chain as a Markov process in either discrete or continuous time with a countable state space (thus regardless of the nature of time),^{[13]}^{[14]}^{[15]}^{[16]} but it is also common to define a Markov chain as having discrete time in either countable or continuous state space (thus regardless of the state space).^{[12]}
Types of Markov chains[edit]
The system's state space and time parameter index need to be specified. The following table gives an overview of the different instances of Markov processes for different levels of state space generality and for discrete time v. continuous time:
Countable state space  Continuous or general state space  

Discretetime  (discretetime) Markov chain on a countable or finite state space  Markov chain on a measurable state space (for example, Harris chain) 
Continuoustime  Continuoustime Markov process or Markov jump process  Any continuous stochastic process with the Markov property (for example, the Wiener process) 
Note that there is no definitive agreement in the literature on the use of some of the terms that signify special cases of Markov processes. Usually the term "Markov chain" is reserved for a process with a discrete set of times, that is, a discretetime Markov chain (DTMC),^{[1]}^{[17]}^{[17]} but a few authors use the term "Markov process" to refer to a continuoustime Markov chain (CTMC) without explicit mention.^{[18]}^{[19]}^{[20]} In addition, there are other extensions of Markov processes that are referred to as such but do not necessarily fall within any of these four categories (see Markov model). Moreover, the time index need not necessarily be realvalued; like with the state space, there are conceivable processes that move through index sets with other mathematical constructs. Notice that the general state space continuoustime Markov chain is general to such a degree that it has no designated term.
While the time parameter is usually discrete, the state space of a Markov chain does not have any generally agreedon restrictions: the term may refer to a process on an arbitrary state space.^{[21]} However, many applications of Markov chains employ finite or countably infinite state spaces, which have a more straightforward statistical analysis. Besides timeindex and statespace parameters, there are many other variations, extensions and generalizations (see Variations). For simplicity, most of this article concentrates on the discretetime, discrete statespace case, unless mentioned otherwise.
Transitions[edit]
The changes of state of the system are called transitions.^{[1]} 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 discretetime random process involves a system which is in a certain state at each step, with the state changing randomly between steps.^{[1]} 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.^{[22]} 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.^{[22]} However, the statistical properties of the system's future can be predicted.^{[22]} In many applications, it is these statistical properties that are important.
History[edit]
Markov studied Markov processes in the early 20th century, publishing his first paper on the topic in 1906.^{[23]}^{[24]}^{[25]}^{[26]} Markov processes in continuous time were discovered long before Andrey Markov's work in the early 20th century^{[1]} in the form of the Poisson process.^{[27]}^{[28]}^{[29]} Markov was interested in studying an extension of independent random sequences, motivated by a disagreement with Pavel Nekrasov who claimed independence was necessary for the weak law of large numbers to hold.^{[1]}^{[30]} In his first paper on Markov chains, published in 1906, Markov showed that under certain conditions the average outcomes of the Markov chain would converge to a fixed vector of values, so proving a weak law of large numbers without the independence assumption,^{[1]}^{[24]}^{[25]}^{[26]} which had been commonly regarded as a requirement for such mathematical laws to hold.^{[26]} Markov later used Markov chains to study the distribution of vowels in Eugene Onegin, written by Alexander Pushkin, and proved a central limit theorem for such chains.^{[1]}^{[24]}
In 1912 Henri Poincaré studied Markov chains on finite groups with an aim to study card shuffling. Other early uses of Markov chains include a diffusion model, introduced by Paul and Tatyana Ehrenfest in 1907, and a branching process, introduced by Francis Galton and Henry William Watson in 1873, preceding the work of Markov.^{[24]}^{[25]} After the work of Galton and Watson, it was later revealed that their branching process had been independently discovered and studied around three decades earlier by IrénéeJules Bienaymé.^{[31]} Starting in 1928, Maurice Fréchet became interested in Markov chains, eventually resulting in him publishing in 1938 a detailed study on Markov chains.^{[24]}^{[32]}
Andrei Kolmogorov developed in a 1931 paper a large part of the early theory of continuoustime Markov processes.^{[33]}^{[34]} Kolmogorov was partly inspired by Louis Bachelier's 1900 work on fluctuations in the stock market as well as Norbert Wiener's work on Einstein's model of Brownian movement.^{[33]}^{[35]} He introduced and studied a particular set of Markov processes known as diffusion processes, where he derived a set of differential equations describing the processes.^{[33]}^{[36]} Independent of Kolmogorov's work, Sydney Chapman derived in a 1928 paper an equation, now called the Chapman–Kolmogorov equation, in a less mathematically rigorous way than Kolmogorov, while studying Brownian movement.^{[37]} The differential equations are now called the Kolmogorov equations^{[38]} or the Kolmogorov–Chapman equations.^{[39]} Other mathematicians who contributed significantly to the foundations of Markov processes include William Feller, starting in 1930s, and then later Eugene Dynkin, starting in the 1950s.^{[34]}
Examples[edit]
Random walks based on integers and the gambler's ruin problem are examples of Markov processes.^{[40]}^{[41]} Some variations of these processes were studied hundreds of years earlier in the context of independent variables.^{[42]}^{[43]}^{[44]} Two important examples of Markov processes are the Wiener process, also known as the Brownian motion process, and the Poisson process,^{[27]} which are considered the most important and central stochastic processes in the theory of stochastic processes.^{[45]}^{[46]}^{[47]} These two processes are Markov processes in continuous time, while random walks on the integers and the gambler's ruin problem are examples of Markov processes in discrete time.^{[40]}^{[41]}
A famous Markov chain is the socalled "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, tomorrow it will eat grapes with probability 4/10 or cheese with probability 6/10. It will not eat lettuce again tomorrow.
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 any other time 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 nontrivially on the current state.
A nonMarkov example[edit]
Suppose that there is a coin purse containing five quarters (each worth 25¢), five dimes (each worth 10¢), and five nickels (each worth 5¢), and one by one, coins are randomly drawn from the purse and are set on a table. If represents the total value of the coins set on the table after n draws, with , then the sequence is not a Markov process.
To see why this is the case, suppose that in the first six draws, all five nickels and a quarter are drawn. Thus . If we know not just , but the earlier values as well, then we can determine which coins have been drawn, and we know that the next coin will not be a nickel; so we can determine that with probability 1. But if we do not know the earlier values, then based only on the value we might guess that we had drawn four dimes and two nickels, in which case it would certainly be possible to draw another nickel next. Thus, our guesses about are impacted by our knowledge of values prior to .
However, it is possible to model this scenario as a Markov process. Instead of defining to represent the total value of the coins on the table, we could define to represent the count of the various coin types on the table. For instance, could be defined to represent the state where there is one quarter, zero dimes, and five nickels on the table after 6 onebyone draws. This new model would be represented by 216 possible states (that is, 6x6x6 states, since each of the three coin types could have zero to five coins on the table by the end of the 6 draws). Suppose that the first draw results in state . The probability of achieving now depends on ; for example, the state is not possible. After the second draw, the third draw depends on which coins have so far been drawn, but no longer only on the coins that were drawn for the first state (since probabilistically important information has since been added to the scenario). In this way, the likelihood of the state depends exclusively on the outcome of the state.
Formal definition[edit]
Discretetime Markov chain[edit]
A discretetime Markov chain is a sequence of random variables X_{1}, X_{2}, X_{3}, ... with the Markov property, namely that the probability of moving to the next state depends only on the present state and not on the previous states:
 if both conditional probabilities are well defined, that is, if
The possible values of X_{i} form a countable set S called the state space of the chain.
Markov chains are often described by a sequence of directed graphs, where the edges of graph n are labeled by the probabilities of going from one state at time n to the other states at time n + 1, The same information is represented by the transition matrix from time n to time n + 1. However, Markov chains are frequently assumed to be timehomogeneous (see variations below), in which case the graph and matrix are independent of n and are thus not presented as sequences.
These descriptions highlight the structure of the Markov chain that is independent of the initial distribution When timehomogeneous, the chain can be interpreted as a state machine assigning a probability of hopping from each vertex or state to an adjacent one. The probability of the machine's state can be analyzed as the statistical behavior of the machine with an element of the state space as input, or as the behavior of the machine with the initial distribution of states as input, where is the Iverson bracket.
The fact that some sequences of states might have zero probability of occurring corresponds to a graph with multiple connected components, where we omit edges that would carry a zero transition probability. For example, if a has a nonzero probability of going to b, but a and x lie in different connected components of the graph, then is defined, while is not.
Variations[edit]
 Timehomogeneous Markov chains (or stationary Markov chains) are processes where
 for all n. The probability of the transition is independent of n.
 A Markov chain with memory (or a Markov chain of order m)
 where m is finite, is a process satisfying
 In other words, the future state depends on the past m states. It is possible to construct a chain from which has the 'classical' Markov property by taking as state space the ordered mtuples of X values, ie. .
Continuoustime Markov chain[edit]
A continuoustime Markov chain (X_{t})_{t ≥ 0} is defined by a finite or countable state space S, a transition rate matrix Q with dimensions equal to that of the state space and initial probability distribution defined on the state space. For i ≠ j, the elements q_{ij} are nonnegative and describe the rate of the process transitions from state i to state j. The elements q_{ii} are chosen such that each row of the transition rate matrix sums to zero, while the rowsums of a probability transition matrix in a (discrete) Markov chain are all equal to one.
There are three equivalent definitions of the process.^{[48]}
Infinitesimal definition[edit]
Let be the random variable describing the state of the process at time t, and assume the process is in a state i at time t. Then, knowing , is independent of previous values , and as h → 0 for all j and for all t,
 ,
where is the Kronecker delta, using the littleo notation. The can be seen as measuring how quickly the transition from i to j happens.
Jump chain/holding time definition[edit]
Define a discretetime Markov chain Y_{n} to describe the nth jump of the process and variables S_{1}, S_{2}, S_{3}, ... to describe holding times in each of the states where S_{i} follows the exponential distribution with rate parameter −q_{YiYi}.
Transition probability definition[edit]
For any value n = 0, 1, 2, 3, ... and times indexed up to this value of n: t_{0}, t_{1}, t_{2}, ... and all states recorded at these times i_{0}, i_{1}, i_{2}, i_{3}, ... it holds that
where p_{ij} is the solution of the forward equation (a firstorder differential equation)
with initial condition P(0) is the identity matrix.
Finite state space[edit]
If the state space is finite, the transition probability distribution can be represented by a matrix, called the transition matrix, with the (i, j)th element of P equal to
Since each row of P sums to one and all elements are nonnegative, P is a right stochastic matrix.
Stationary distribution relation to eigenvectors and simplices[edit]
A stationary distribution π is a (row) vector, whose entries are nonnegative and sum to 1, is unchanged by the operation of transition matrix P on it and so is defined by
By comparing this definition with that of an eigenvector we see that the two concepts are related and that
is a normalized () 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 a stationary state that is the limit of the sequence of distributions for some initial distribution.
The values of a stationary distribution 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 we see that the dot product of π with a vector whose components are all 1 is unity and that π lies on a simplex.
Timehomogeneous Markov chain with a finite state space[edit]
This section includes a list of references, related reading or external links, but its sources remain unclear because it lacks inline citations. (February 2012) (Learn how and when to remove this template message)

If the Markov chain is timehomogeneous, then the transition matrix P is the same after each step, so the kstep transition probability can be computed as the kth power of the transition matrix, P^{k}.
If the Markov chain is irreducible and aperiodic, then there is a unique stationary distribution π.^{[49]} Additionally, in this case P^{k} converges to a rankone matrix in which each row is the stationary distribution π:
where 1 is the column vector with all entries equal to 1. This is stated by the Perron–Frobenius theorem. If, by whatever means, 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 does not exist while the stationary distribution does, as shown by this example:
(This example illustrates a periodic Markov chain.)
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
It is always true that
Subtracting Q from both sides and factoring then yields
where I_{n} is the identity matrix of size n, and 0_{n,n} 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 substitutes each of its elements by one, and on the other one substitutes the corresponding element (the one in the same column) in the vector 0, and next leftmultiplies 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 rightmost column replaced with all 1's. If [f(P − I_{n})]^{−1} exists then^{[50]}^{[51]}
 Explain: The original matrix equation is equivalent to a 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 rightmost column of (PIn)” have been replaced by the n stochastic ones.
One thing to notice is that if P has an element P_{i,i} on its main diagonal that is equal to 1 and the ith row or column is otherwise filled with 0's, then that row or column will remain unchanged in all of the subsequent powers P^{k}. Hence, the ith 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[edit]
This section includes a list of references, related reading or external links, but its sources remain unclear because it lacks inline citations. (February 2012) (Learn how and when to remove this template message)

As stated earlier, from the equation (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 nondiagonalizable, that is, defective matrices, one may start with the Jordan normal form of P and proceed with a bit more involved set of arguments in a similar way.^{[52]})
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, that is, Σ = diag(λ_{1},λ_{2},λ_{3},...,λ_{n}). Then by eigendecomposition
Let the eigenvalues be enumerated such that:
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_{i} be the ith column of U matrix, that is, u_{i} is the left eigenvector of P corresponding to λ_{i}. Also let x be a length n row vector that represents a valid probability distribution; since the eigenvectors u_{i} span we can write
If we multiply x with P from right and continue this operation with the results, in the end we get the stationary distribution π. In other words, π = u_{i} ← xPP...P = xP^{k} as k → ∞. That means
Since π = u_{1}, π^{(k)} approaches to π as k → ∞ with a speed in the order of λ_{2}/λ_{1} exponentially. This follows because hence λ_{2}/λ_{1} is the dominant term. Random noise in the state distribution π can also speed up this convergence to the stationary distribution.^{[53]}
General state space[edit]
For an overview of Markov chains on a general state space, see Markov chains on a measurable state space.
Harris chains[edit]
Many results for Markov chains with finite state space can be generalized to chains with uncountable state space through Harris chains. 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
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 chains 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[edit]
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 (probabilistic cellular automata). See for instance Interaction of Markov Processes^{[54]} or^{[55]}
Properties[edit]
Reducibility[edit]
A Markov chain is said to be irreducible if it is possible to get to any state from any state.^{[1]} The following explains this definition more formally.
A state j is said to be accessible from a state i (written i → j) if a system started in state i has a nonzero probability of transitioning into state j at some point. Formally, state j is accessible from state i if there exists an integer n_{ij} ≥ 0 such that
This integer is allowed to be different for each pair of states, hence the subscripts in n_{ij}. Allowing n to be zero means that every state is accessible from itself by definition. The accessibility relation is reflexive and transitive, but not necessarily symmetric.
A state i is said to communicate with state j (written i ↔ j) if both i → j and j → i. A communicating class is a maximal set of states C such that every pair of states in C communicates with each other. Communication is an equivalence relation, and communicating classes are the equivalence classes of this relation.
A communicating class is closed if the probability of leaving the class is zero, namely if i is in C but j is not, then j is not accessible from i. The set of communicating classes forms a directed, acyclic graph by inheriting the arrows from the original state space. A communicating class is closed if and only if it has no outgoing arrows in this graph.
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.^{[56]} A state is final if and only if its communicating class is closed.
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[edit]
A state i has period k if any return to state i must occur in multiples of k time steps. Formally, the period of a state is defined as
(where "gcd" is the greatest common divisor) provided that this set is not empty. Otherwise the period is not defined. 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. 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.
Transience and recurrence[edit]
A state i is said to be transient if, given that we start in state i, there is a nonzero probability that we will never return to i. Formally, let the random variable T_{i} be the first return time to state i (the "hitting time"):
The number
is the probability that we return to state i for the first time after n steps. Therefore, state i is transient if
State i is recurrent (or persistent) if it is not transient. Recurrent states are guaranteed (with probability 1) to have a finite hitting time. Recurrence and transience are class properties, that is, they either hold or do not hold equally for all members of a communicating class.
Mean recurrence time[edit]
Even if the hitting time is finite with probability 1, it need not have a finite expectation. The mean recurrence time at state i is the expected return time M_{i}:
State i is positive recurrent (or nonnull persistent) if M_{i} is finite; otherwise, state i is null recurrent (or null persistent).
Expected number of visits[edit]
It can be shown that a state i is recurrent if and only if the expected number of visits to this state is infinite:
Absorbing states[edit]
A state i is called absorbing if it is impossible to leave this state. Therefore, the state i is absorbing if and only if
If every state can reach an absorbing state, then the Markov chain is an absorbing Markov chain.^{[1]}
Ergodicity[edit]
A state i is said to be 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 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. More generally, a Markov chain is ergodic if there is a number N such that any state can be reached from any other state in any number of steps greater than or equal to a number N. In case of a fully connected transition matrix, where all transitions have a nonzero probability, this condition is fulfilled with N = 1.
A Markov chain with more than one state and just one outgoing transition per state is either not irreducible or not aperiodic, hence cannot be ergodic.
Steadystate analysis and limiting distributions[edit]
If the Markov chain is a timehomogeneous Markov chain,^{[57]} so that the process is described by a single, timeindependent matrix , then the vector is called a stationary distribution (or invariant measure) if it satisfies
An irreducible chain has a positive stationary distribution (a stationary distribution such that ) if and only if all of its states are positive recurrent.^{[58]} In that case, π is unique and is related to the expected return time:
where is the normalizing constant. Further, if the positive recurrent chain is both irreducible and aperiodic, it is said to have a limiting distribution; for any i and j,
(There is no assumption on the starting distribution; the chain converges to the stationary distribution regardless of where it begins. Such 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 closed communicating class in the chain; each one will have its own unique stationary distribution . 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 's). However, if a state j is aperiodic, then and for any other state i, let f_{ij} be the probability that the chain ever visits state j if it starts at i,
If a state i is periodic with period k > 1 then the limit does not exist, although the limit does exist for every integer r.
Steadystate analysis and the timeinhomogeneous Markov chain[edit]
A Markov chain need not necessarily be timehomogeneous to have an equilibrium distribution. If there is a probability distribution over states such that
for every state j and every time n then is an equilibrium distribution of the Markov chain. Such can occur in 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.
Reversible Markov chain[edit]
A Markov chain is said to be reversible if there is a probability distribution π over its states such that
for all times n and all states i and j. This condition is known as the detailed balance condition (some books call it the local balance equation).
Considering a fixed arbitrary time n and using the shorthand
the detailed balance equation can be written more compactly as
The single timestep from n to n + 1 can be thought of as each person i having π_{i} dollars initially and paying each person j a fraction p_{ij} of it. The detailed balance condition states that upon each payment, the other person pays exactly the same amount of money back.^{[59]} Clearly the total amount of money π each person has remains the same after the timestep, since every dollar spent is balanced by a corresponding dollar received. This can be shown more formally by the equality
which essentially states that the total amount of money person j receives (including from himself) during the timestep equals the amount of money he pays others, which equals all the money he initially had because it was assumed that all money is spent (that is, p_{ji} sums to 1 over i). The assumption is a technical one, because the money not really used is simply thought of as being paid from person j to himself (that is, p_{jj} is not necessarily zero).
As n was arbitrary, this reasoning holds for any n, and therefore for reversible Markov chains π is always a steadystate distribution of Pr(X_{n+1} = j  X_{n} = i) for every n.
If the Markov chain begins in the steadystate distribution, that is, if , then for all and the detailed balance equation can be written as
The left and righthand 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 (MCMC) approaches because the detailed balance equation for a desired distribution π necessarily implies that the Markov chain has been constructed so that π is a steadystate distribution. Even with timeinhomogeneous 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 steadystate distribution of the Markov chain.
Closest reversible Markov chain[edit]
For any timehomogeneous Markov chain given by a transition matrix , any norm on which is induced by a scalar product, and any probability vector , there exists a unique transition matrix which is reversible according to and which is closest to according to the norm The matrix can be computed by solving a quadraticconvex optimization problem.^{[60]} A GNU licensed Matlab script that computes the nearest reversible Markov chain can be found here.
For example, consider the following Markov chain:
This Markov chain is not reversible. According to the Frobenius norm the closest reversible Markov chain according to can be computed as
If we choose the probability vector as , then the closest reversible Markov chain according to the Frobenius norm is approximately given by
Markovian representations[edit]
In some cases, apparently nonMarkovian processes may still have Markovian representations, constructed by expanding the concept of the 'current' and 'future' states. For example, let X be a nonMarkovian process. Then define a process Y, such that each state of Y represents a timeinterval of states of X. Mathematically, this takes the form:
If Y has the Markov property, then it is a Markovian representation of X.
An example of a nonMarkovian process with a Markovian representation is an autoregressive time series of order greater than one.^{[61]}
Transient evolution[edit]
The probability of going from state i to state j in n time steps is
and the singlestep transition is
For a timehomogeneous Markov chain:
and
The nstep transition probabilities satisfy the Chapman–Kolmogorov equation, that for any k such that 0 < k < n,
where S is the state space of the Markov chain.
The marginal distribution Pr(X_{n} = x) is the distribution over states at time n. The initial distribution is Pr(X_{0} = x). The evolution of the process through one time step is described by
(The superscript (n) is an index, and not an exponent.
Transient behaviour[edit]
Write P(t) for the matrix with entries p_{ij} = P(X_{t} = j  X_{0} = i). Then the matrix P(t) satisfies the forward equation, a firstorder differential equation
where the prime denotes differentiation with respect to t. The solution to this equation is given by a matrix exponential
In a simple case such as a CTMC on the state space {1,2}. The general Q matrix for such a process is the following 2 × 2 matrix with α,β > 0
The above relation for forward matrix can be solved explicitly in this case to give
However, direct solutions are complicated to compute for larger matrices. The fact that Q is the generator for a semigroup of matrices
is used.
Stationary distribution[edit]
The stationary distribution for an irreducible recurrent CTMC is the probability distribution to which the process converges for large values of t. Observe that for the twostate process considered earlier with P(t) given by
as t → ∞ the distribution tends to
Observe that each row has the same distribution as this does not depend on starting state. The row vector π may be found by solving^{[62]}
with the additional constraint that
Example 1[edit]
The image to the right describes a continuoustime Markov chain with statespace {Bull market, Bear market, Stagnant market} and transition rate matrix
The stationary distribution of this chain can be found by solving π (QI) = 0, where I is the identity matrix, subject to the constraint that elements must sum to 1 to obtain
Example 2[edit]
The image to the right describes a discretetime Markov chain with statespace {1,2,3,4,5,6,7,8,9}. The player controls PacMan through a maze, eating pacdots. Meanwhile, he is being hunted by ghosts. For convenience, the maze shall be a small 3x3grid and the monsters move randomly in horizontal and vertical directions. A secret passageway between states 2 and 8 can be used in both directions. Entries with probability zero are removed in the following transition matrix:
This Markov chain is irreducible, because the ghosts can fly from every state to every state in a finite amount of time. Due to the secret passageway, the Markov chain is also aperiodic, because the monsters can move from any state to any state both in an even and in an uneven number of state transitions. Therefore, a unique stationary distribution exists and can be found by solving π Q = 0 subject to the constraint that elements must sum to 1. The solution of this linear equation subject to the constraint is The central state and the border states 2 and 8 of the adjacent secret passageway are visited most and the corner states are visited least.
Hitting times[edit]
The hitting time is the time, starting in a given set of states until the chain arrives in a given state or set of states. The distribution of such a time period has a phase type distribution. The simplest such distribution is that of a single exponentially distributed transition.
Expected hitting times[edit]
For a subset of states A ⊆ S, the vector k^{A} of hitting times (where element represents the expected value, starting in state i that the chain enters one of the states in the set A) is the minimal nonnegative solution to^{[62]}
Time reversal[edit]
For a CTMC X_{t}, the timereversed process is defined to be . By Kelly's lemma this process has the same stationary distribution as the forward process.
A chain is said to be reversible if the reversed process is the same as the forward process. Kolmogorov's criterion states that the necessary and sufficient condition for a process to be reversible is that the product of transition rates around a closed loop must be the same in both directions.
Embedded Markov chain[edit]
[edit]
One method of finding the stationary probability distribution, π, of an ergodic continuoustime Markov chain, Q, is by first finding its embedded Markov chain (EMC). Strictly speaking, the EMC is a regular discretetime Markov chain, sometimes referred to as a jump process. Each element of the onestep transition probability matrix of the EMC, S, is denoted by s_{ij}, and represents the conditional probability of transitioning from state i into state j. These conditional probabilities may be found by
From this, S may be written as
where I is the identity matrix and diag(Q) is the diagonal matrix formed by selecting the main diagonal from the matrix Q and setting all other elements to zero.
To find the stationary probability distribution vector, we must next find such that
with being a row vector, such that all elements in are greater than 0 and = 1. From this, π may be found as
(S may be periodic, even if Q is not. Once π is found, it must be normalized to a unit vector.)
Another discretetime process that may be derived from a continuoustime Markov chain is a δskeleton—the (discretetime) Markov chain formed by observing X(t) at intervals of δ units of time. The random variables X(0), X(δ), X(2δ), ... give the sequence of states visited by the δskeleton.
Special types of Markov chains[edit]
Markov model[edit]
Markov models are used to model changing systems. There are 4 main types of models, that generalize Markov chains depending on whether every sequential state is observable or not, and whether the system is to be adjusted on the basis of observations made:
System state is fully observable  System state is partially observable  

System is autonomous  Markov chain  Hidden Markov model 
System is controlled  Markov decision process  Partially observable Markov decision process 
Bernoulli scheme[edit]
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.
Applications[edit]
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[edit]
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 timeinvariant, and that no relevant history need be considered which is not already included in the state description.^{[63]}^{[64]} For example, a thermodynamic state operates under a probability distribution that is difficult or expensive to acquire. Therefore, Markov Chain Monte Carlo method can be used to draw samples randomly from a blackbox to approximate the probability distribution of attributes over a range of objects.^{[65]}
The paths, in the path integral formulation of quantum mechanics, are Markov chains.^{[66]}
Markov chains are used in lattice QCD simulations.^{[67]}
Chemistry[edit]
A reaction network is a chemical system involving multiple reactions and chemical species. The simplest stochastic models of such networks treat the system as a continuous time Markov chain with the state being the number of molecules of each species and with reactions modeled as possible transitions of the chain.^{[68]} Markov chains and continuoustime Markov processes are useful in chemistry when physical systems closely approximate the Markov property. For example, imagine a large number n of molecules in solution in state A, each of which can undergo a chemical reaction to state B with a certain average rate. Perhaps the molecule is an enzyme, and the states refer to how it is folded. The state of any single enzyme follows a Markov chain, and since the molecules are essentially independent of each other, the number of molecules in state A or B at a time is n times the probability a given molecule is in that state.
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 MichaelisMenten is fairly straightforward, far more complicated reaction networks can also be modeled with Markov chains.^{[69]}
An algorithm based on a Markov chain was also used to focus the fragmentbased growth of chemicals in silico towards a desired class of compounds such as drugs or natural products.^{[70]} As a molecule is grown, a fragment is selected from the nascent molecule as the "current" state. It is not aware of its past (that is, 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.^{[71]}
Also, the growth (and composition) of copolymers 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 (for example, whether monomers tend to add in alternating fashion or in long runs of the same monomer). Due to steric effects, secondorder 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.^{[72]}
Testing[edit]
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 statebased 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.
Solar irradiance variability[edit]
Solar irradiance variability assessments are useful for solar power applications. Solar irradiance variability at any location over time is mainly a consequence of the deterministic variability of the sun's path across the sky dome and the variability in cloudiness. The variability of accessible solar irradiance on Earth's surface has been modeled using Markov chains,^{[73]}^{[74]}^{[75]}^{[76]} also including modeling the two states of clear and cloudiness as a twostate Markov chain.^{[77]}^{[78]}
Speech recognition[edit]
Hidden Markov models are the basis for most modern automatic speech recognition systems.
Information and computer science[edit]
Markov chains are used throughout information processing. 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 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 (such as in rearrangements detection^{[79]}).
The LZMA lossless data compression algorithm combines Markov chains with LempelZiv compression to achieve very high compression ratios.
Queueing theory[edit]
Markov chains are the basis for the analytical treatment of queues (queueing theory). Agner Krarup Erlang initiated the subject in 1917.^{[80]} This makes them critical for optimizing the performance of telecommunications networks, where messages must often compete for limited resources (such as bandwidth).^{[81]}
Numerous queueing models use continuoustime Markov chains. For example, an M/M/1 queue is a CTMC on the nonnegative integers where upward transitions from i to i + 1 occur at rate λ according to a Poisson process and describe job arrivals, while transitions from i to i – 1 (for i > 1) occur at rate μ (job service times are exponentially distributed) and describe completed services (departures) from the queue.
Internet applications[edit]
The PageRank of a webpage as used by Google is defined by a Markov chain.^{[82]}^{[83]}^{[84]} It is the probability to be at page in the stationary distribution on the following Markov chain on all (known) webpages. If is the number of known webpages, and a page has links to it then it has transition probability for all pages that are linked to and for all pages that are not linked to. The parameter is taken to be about 0.85.^{[85]}
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 secondorder Markov models and can be used to make predictions regarding future navigation and to personalize the web page for an individual user.
Statistics[edit]
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 distributions to be simulated and their parameters found numerically.
Economics and finance[edit]
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.^{[dubious – discuss]}^{[86]} Another was the regimeswitching model of James D. Hamilton (1989), in which a Markov chain is used to model switches between periods high and low GDP growth (or alternatively, economic expansions and recessions).^{[87]} A more recent example is the Markov switching multifractal model of Laurent E. Calvet and Adlai J. Fisher, which builds upon the convenience of earlier regimeswitching models.^{[88]}^{[89]} 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.^{[90]}
Credit rating agencies produce annual tables of the transition probabilities for bonds of different credit ratings.^{[91]}
Social sciences[edit]
Markov chains are generally used in describing pathdependent 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 to model how once a country reaches a specific level of economic development, the configuration of structural factors, such as size of the middle class, 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.^{[92]}
Mathematical biology[edit]
Markov chains also have many applications in biological modelling, particularly population processes, which are useful in modelling processes that are (at least) analogous to biological populations. The Leslie matrix, is one such example used to describe the population dynamics of many species, 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.^{[93]} Yet another example is the state of ion channels in cell membranes.
Markov chains are also used in simulations of brain function, such as the simulation of the mammalian neocortex.^{[94]} Markov chains have also been used to model viral infection of single cells.^{[95]}
Genetics[edit]
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 the diffusion equation method described by Motoo Kimura.^{[96]}
Games[edit]
Markov chains can be used to model many games of chance.^{[1]} The children's games Snakes and Ladders and "Hi Ho! CherryO", 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[edit]
Markov chains are employed in algorithmic music composition, particularly in software such as Csound, Max, and SuperCollider. In a firstorder 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 output note values based on the transition matrix weightings, which could be MIDI note values, frequency (Hz), or any other desirable metric.^{[97]}
Note  A  C♯  E♭ 

A  0.1  0.6  0.3 
C♯  0.25  0.05  0.7 
E♭  0.7  0.3  0 
Notes  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 
A secondorder Markov chain can be introduced by considering the current state and also the previous state, as indicated in the second table. Higher, nthorder chains tend to "group" particular notes together, while 'breaking off' into other patterns and sequences occasionally. These higherorder chains tend to generate results with a sense of phrasal structure, rather than the 'aimless wandering' produced by a firstorder system.^{[98]}
Markov chains can be used structurally, as in Xenakis's Analogique A and B.^{[99]} Markov chains are also used in systems which use a Markov model to react interactively to music input.^{[100]}
Usually musical systems need to enforce specific control constraints on the finitelength sequences they generate, but control constraints are not compatible with Markov models, since they induce longrange dependencies that violate the Markov hypothesis of limited memory. In order to overcome this limitation, a new approach has been proposed.^{[101]}
Baseball[edit]
Markov chain models have been used in advanced baseball analysis since 1960, although their use is still rare. Each halfinning of a baseball game fits the Markov chain state when the number of runners and outs are considered. During any atbat, 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.^{[102]} 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 bunting and base stealing and differences when playing on grass vs. AstroTurf.^{[103]}
Markov text generators[edit]
Markov processes can also be used to generate superficially reallooking text given a sample document. Markov processes are used in a variety of recreational "parody generator" software (see dissociated press, Jeff Harrison,^{[104]} Mark V. Shaney,^{[105]}^{[106]} and Academias Neutronium).
Bioinformatics[edit]
In the bioinformatics field, they can be used to simulate DNA sequences.^{[107]}^{[108]}
Probabilistic forecasting[edit]
Markov chains have been used for forecasting in several areas: for example, price trends,^{[109]} wind power,^{[110]} and solar irradiance.^{[111]} The Markov chain forecasting models utilize a variety of settings, from discretizing the time series,^{[110]} to hidden Markov models combined with wavelets,^{[109]} and the Markov chain mixture distribution model (MCM).^{[111]}
See also[edit]
 Dynamics of Markovian particles
 Markov chain approximation method
 Markov chain geostatistics
 Markov chain mixing time
 Markov decision process
 Markov information source
 Markov random field
 Quantum Markov chain
 SemiMarkov process
 Stochastic cellular automaton
 Telescoping Markov chain
 Variableorder Markov model
Notes[edit]
 ^ ^{a} ^{b} ^{c} ^{d} ^{e} ^{f} ^{g} ^{h} ^{i} ^{j} ^{k} ^{l} ^{m} ^{n} Gagniuc, Paul A. (2017). Markov Chains: From Theory to Implementation and Experimentation. USA, NJ: John Wiley & Sons. pp. 1–235. ISBN 9781119387558.
 ^ "Markov chain  Definition of Markov chain in US English by Oxford Dictionaries". Oxford Dictionaries  English. Retrieved 20171214.
 ^ Definition at Brilliant.org "Brilliant Math and Science Wiki". Retrieved on 12 May 2019
 ^ Samuel Karlin; Howard E. Taylor (2 December 2012). A First Course in Stochastic Processes. Academic Press. p. 47. ISBN 9780080570419. Archived from the original on 23 March 2017.
 ^ Bruce Hajek (12 March 2015). Random Processes for Engineers. Cambridge University Press. ISBN 9781316241240. Archived from the original on 23 March 2017.
 ^ G. Latouche; V. Ramaswami (1 January 1999). Introduction to Matrix Analytic Methods in Stochastic Modeling. SIAM. pp. 4–. ISBN 9780898714258. Archived from the original on 23 March 2017.
 ^ ^{a} ^{b} Sean Meyn; Richard L. Tweedie (2 April 2009). Markov Chains and Stochastic Stability. Cambridge University Press. p. 3. ISBN 9780521731829. Archived from the original on 23 March 2017.
 ^ Reuven Y. Rubinstein; Dirk P. Kroese (20 September 2011). Simulation and the Monte Carlo Method. John Wiley & Sons. p. 225. ISBN 9781118210529. Archived from the original on 23 March 2017.
 ^ Dani Gamerman; Hedibert F. Lopes (10 May 2006). Markov Chain Monte Carlo: Stochastic Simulation for Bayesian Inference, Second Edition. CRC Press. ISBN 9781584885870. Archived from the original on 23 March 2017.
 ^ "Markovian". Oxford English Dictionary (3rd ed.). Oxford University Press. September 2005. (Subscription or UK public library membership required.)
 ^ Øksendal, B. K. (Bernt Karsten), 1945 (2003). Stochastic differential equations : an introduction with applications (6th ed ed.). Berlin: Springer. ISBN 3540047581. OCLC 52203046.CS1 maint: multiple names: authors list (link) CS1 maint: extra text (link)
 ^ ^{a} ^{b} Søren Asmussen (15 May 2003). Applied Probability and Queues. Springer Science & Business Media. p. 7. ISBN 9780387002118. Archived from the original on 23 March 2017.
 ^ Emanuel Parzen (17 June 2015). Stochastic Processes. Courier Dover Publications. p. 188. ISBN 9780486796888. Archived from the original on 20 November 2017.
 ^ Samuel Karlin; Howard E. Taylor (2 December 2012). A First Course in Stochastic Processes. Academic Press. pp. 29 and 30. ISBN 9780080570419. Archived from the original on 23 March 2017.
 ^ John Lamperti (1977). Stochastic processes: a survey of the mathematical theory. SpringerVerlag. pp. 106–121. ISBN 9783540902751. Archived from the original on 20170323.
 ^ Sheldon M. Ross (1996). Stochastic processes. Wiley. pp. 174 and 231. ISBN 9780471120629. Archived from the original on 20170323.
 ^ ^{a} ^{b} Everitt, B.S. (2002) The Cambridge Dictionary of Statistics. CUP. ISBN 052181099X
 ^ Parzen, E. (1962) Stochastic Processes, HoldenDay. ISBN 0816266646 (Table 6.1)
 ^ Dodge, Y. (2003) The Oxford Dictionary of Statistical Terms, OUP. ISBN 0199206139 (entry for "Markov chain")
 ^ Dodge, Y. The Oxford Dictionary of Statistical Terms, OUP. ISBN 0199206139
 ^ Meyn, S. Sean P., and Richard L. Tweedie. (2009) Markov chains and stochastic stability. Cambridge University Press. (Preface, p. iii)
 ^ ^{a} ^{b} ^{c} Gagniuc, Paul A. (2017). Markov Chains: From Theory to Implementation and Experimentation. USA, NJ: John Wiley & Sons. pp. 159–163. ISBN 9781119387558.
 ^ Gagniuc, Paul A. (2017). Markov Chains: From Theory to Implementation and Experimentation. USA, NJ: John Wiley & Sons. pp. 2–8. ISBN 9781119387558.
 ^ ^{a} ^{b} ^{c} ^{d} ^{e} Charles Miller Grinstead; James Laurie Snell (1997). Introduction to Probability. American Mathematical Soc. pp. 464–466. ISBN 9780821807491. Archived from the original on 20170323.
 ^ ^{a} ^{b} ^{c} Pierre Bremaud (9 March 2013). Markov Chains: Gibbs Fields, Monte Carlo Simulation, and Queues. Springer Science & Business Media. p. ix. ISBN 9781475731248. Archived from the original on 23 March 2017.
 ^ ^{a} ^{b} ^{c} Hayes, Brian (2013). "First links in the Markov chain". American Scientist. 101 (2): 92–96. doi:10.1511/2013.101.92.
 ^ ^{a} ^{b} Sheldon M. Ross (1996). Stochastic processes. Wiley. pp. 235 and 358. ISBN 9780471120629. Archived from the original on 20170323.
 ^ Jarrow, Robert; Protter, Philip (2004). A short history of stochastic integration and mathematical finance: The early years, 1880–1970. A Festschrift for Herman Rubin. Institute of Mathematical Statistics Lecture Notes – Monograph Series. pp. 75–91. CiteSeerX 10.1.1.114.632. doi:10.1214/lnms/1196285381. ISBN 9780940600614. ISSN 07492170.
 ^ Guttorp, Peter; Thorarinsdottir, Thordis L. (2012). "What Happened to Discrete Chaos, the Quenouille Process, and the Sharp Markov Property? Some History of Stochastic Point Processes". International Statistical Review. 80 (2): 253–268. doi:10.1111/j.17515823.2012.00181.x. ISSN 03067734.
 ^ Seneta, E. (1996). "Markov and the Birth of Chain Dependence Theory". International Statistical Review / Revue Internationale de Statistique. 64 (3): 255–257. doi:10.2307/1403785. ISSN 03067734. JSTOR 1403785.
 ^ Seneta, E. (1998). "I.J. Bienaymé [1796–1878]: Criticality, Inequality, and Internationalization". International Statistical Review / Revue Internationale de Statistique. 66 (3): 291–292. doi:10.2307/1403518. ISSN 03067734. JSTOR 1403518.
 ^ Bru B, Hertz S (2001). "Maurice Fréchet". In Heyde CC, Seneta E, Crépel P, Fienberg SE, Gani J (eds.). Statisticians of the Centuries. New York, NY: Springer. pp. 331–334. doi:10.1007/9781461301790_71. ISBN 9780387952833.
 ^ ^{a} ^{b} ^{c} Kendall, D. G.; Batchelor, G. K.; Bingham, N. H.; Hayman, W. K.; Hyland, J. M. E.; Lorentz, G. G.; Moffatt, H. K.; Parry, W.; Razborov, A. A.; Robinson, C. A.; Whittle, P. (1990). "Andrei Nikolaevich Kolmogorov (1903–1987)". Bulletin of the London Mathematical Society. 22 (1): 33. doi:10.1112/blms/22.1.31. ISSN 00246093.
 ^ ^{a} ^{b} Cramer, Harald (1976). "Half a Century with Probability Theory: Some Personal Recollections". The Annals of Probability. 4 (4): 509–546. doi:10.1214/aop/1176996025. ISSN 00911798.
 ^ Marc Barbut; Bernard Locker; Laurent Mazliak (23 August 2016). Paul Lévy and Maurice Fréchet: 50 Years of Correspondence in 107 Letters. Springer London. p. 5. ISBN 9781447172628. Archived from the original on 23 March 2017.
 ^ Valeriy Skorokhod (5 December 2005). Basic Principles and Applications of Probability Theory. Springer Science & Business Media. p. 146. ISBN 9783540263128. Archived from the original on 23 March 2017.
 ^ Bernstein, Jeremy (2005). "Bachelier". American Journal of Physics. 73 (5): 395–398. Bibcode:2005AmJPh..73..395B. doi:10.1119/1.1848117. ISSN 00029505.
 ^ William J. Anderson (6 December 2012). ContinuousTime Markov Chains: An ApplicationsOriented Approach. Springer Science & Business Media. p. vii. ISBN 9781461230380. Archived from the original on 23 March 2017.
 ^ Kendall, D. G.; Batchelor, G. K.; Bingham, N. H.; Hayman, W. K.; Hyland, J. M. E.; Lorentz, G. G.; Moffatt, H. K.; Parry, W.; Razborov, A. A.; Robinson, C. A.; Whittle, P. (1990). "Andrei Nikolaevich Kolmogorov (1903–1987)". Bulletin of the London Mathematical Society. 22 (1): 57. doi:10.1112/blms/22.1.31. ISSN 00246093.
 ^ ^{a} ^{b} Ionut Florescu (7 November 2014). Probability and Stochastic Processes. John Wiley & Sons. pp. 373 and 374. ISBN 9781118593202. Archived from the original on 23 March 2017.
 ^ ^{a} ^{b} Samuel Karlin; Howard E. Taylor (2 December 2012). A First Course in Stochastic Processes. Academic Press. p. 49. ISBN 9780080570419. Archived from the original on 23 March 2017.
 ^ Gagniuc, Paul A. (2017). Markov Chains: From Theory to Implementation and Experimentation. USA, NJ: John Wiley & Sons. pp. 1–2. ISBN 9781119387558.
 ^ Weiss, George H. (2006). "Random Walks". Encyclopedia of Statistical Sciences. p. 1. doi:10.1002/0471667196.ess2180.pub2. ISBN 9780471667193.
 ^ Michael F. Shlesinger (1985). The Wonderful world of stochastics: a tribute to Elliott W. Montroll. NorthHolland. pp. 8–10. ISBN 9780444869371. Archived from the original on 20170323.
 ^ Emanuel Parzen (17 June 2015). Stochastic Processes. Courier Dover Publications. p. 7 and 8. ISBN 9780486796888. Archived from the original on 20 November 2017.
 ^ Joseph L. Doob (1990). Stochastipoic processes. Wiley. p. 46 and 47. Archived from the original on 20171120.
 ^ Donald L. Snyder; Michael I. Miller (6 December 2012). Random Point Processes in Time and Space. Springer Science & Business Media. p. 32. ISBN 9781461231660. Archived from the original on 20 November 2017.
 ^ Norris, J. R. (1997). "Continuoustime Markov chains I". Markov Chains. pp. 60–107. doi:10.1017/CBO9780511810633.004. ISBN 9780511810633.
 ^ Serfozo, Richard (2009). "Basics of Applied Stochastic Processes". Probability and Its Applications. doi:10.1007/9783540893325. ISSN 14317028.
 ^ Chapter 11 "Markov Chains". Retrieved from"Archived copy" (PDF). Archived (PDF) from the original on 20170215. Retrieved 20170602.CS1 maint: archived copy as title (link)
 ^ Serfozo, Richard (2009). "Basics of Applied Stochastic Processes". Probability and Its Applications. doi:10.1007/9783540893325. ISSN 14317028.

^ 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. "Archived copy". CiteSeerX 10.1.1.28.6191. Cite journal requires
journal=
(help)  ^ Franzke, Brandon; Kosko, Bart (1 October 2011). "Noise can speed convergence in Markov chains". Physical Review E. 84 (4): 041112. Bibcode:2011PhRvE..84d1112F. doi:10.1103/PhysRevE.84.041112. PMID 22181092.
 ^ Spitzer, Frank (1970). "Interaction of Markov Processes". Advances in Mathematics. 5 (2): 246–290. doi:10.1016/00018708(70)900344.
 ^ R. L. Dobrushin; V. I. Kri︠u︡kov; A. L. Toom (1978). Stochastic Cellular Systems: Ergodicity, Memory, Morphogenesis. ISBN 9780719022067. Archived from the original on 20170205. Retrieved 20160304.
 ^ Asher Levin, David (2009). Markov chains and mixing times. p. 16. ISBN 9780821847398. Archived from the original on 20170205. Retrieved 20160304.
 ^ Wan, Lijie; Lou, Wenjie; Abner, Erin; Kryscio, Richard J. (2016). "A comparison of timehomogeneous Markov chain and Markov process multistate models". Communications in Statistics: Case Studies, Data Analysis and Applications. 2 (3–4): 92–100. doi:10.1080/23737484.2017.1361366. ISSN 23737484. PMC 5803756. PMID 29430521.
 ^ Serfozo, Richard (2009), "Basics of Applied Stochastic Processes", Probability and its Applications: 35, doi:10.1007/9783540893325, ISBN 9783540893318, MR 2484222, archived from the original on 20150319
 ^ Richard Durrett (19 May 2012). Essentials of Stochastic Processes. Springer Science & Business Media. p. 37. ISBN 9781461436157. Archived from the original on 6 February 2017.
 ^ : A. Nielsen and M. Weber, "Computing the nearest reversible Markov chain". Numerical Linear Algebra with Applications, 22(3):483–499, 2015.

^ Doblinger, G., 1998. Smoothing of Noise AR Signals Using an Adaptive Kalman Filter. In EUSIPCO 98. pp. 781–784. Available at: "Archived copy". CiteSeerX 10.1.1.251.3078. Cite journal requires
journal=
(help) [Accessed January 15, 2015].  ^ ^{a} ^{b} Norris, J. R. (1997). "Continuoustime Markov chains II". Markov Chains. pp. 108–127. doi:10.1017/CBO9780511810633.005. ISBN 9780511810633.
 ^ Thermodynamics and Statistical Mechanics. Richard Fitzpatrick. Retrieved from "Archived copy" (PDF). Archived (PDF) from the original on 20161130. Retrieved 20170602.CS1 maint: archived copy as title (link)
 ^ van Ravenzwaaij, Don; Cassey, Pete; Brown, Scott D. (20160311). "A simple introduction to Markov Chain Monte–Carlo sampling". Psychonomic Bulletin & Review. 25 (1): 143–154. doi:10.3758/s1342301610158. ISSN 10699384.
 ^ van Ravenzwaaij, Don; Cassey, Pete; Brown, Scott D. (20160311). "A simple introduction to Markov Chain Monte–Carlo sampling". Psychonomic Bulletin & Review. 25 (1): 143–154. doi:10.3758/s1342301610158. ISSN 10699384.
 ^ Ryder, Lewis H. (1985). Quantum field theory. Cambridge [Cambridgeshire]: Cambridge University Press. p. 160. ISBN 9780521338592. OCLC 10533049.
 ^ Gattringer, Christof; Lang, Christian B (2010). Quantum Chromodynamics on the Lattice. Lecture Notes in Physics. 788. SpringerVerlag Berlin Heidelberg. doi:10.1007/9783642018503. ISBN 9783642018497. Archived from the original on 20170801.
 ^ Anderson, David F.; Kurtz, Thomas G. (2011), "Continuous Time Markov Chain Models for Chemical Reaction Networks", Design and Analysis of Biomolecular Circuits, Springer New York, pp. 3–42, ISBN 9781441967657, retrieved 20191116
 ^ Du, Chao; Kou, S. C. (September 2012). "Correlation analysis of enzymatic reaction of a single protein molecule". The Annals of Applied Statistics. 6 (3): 950–976. doi:10.1214/12aoas541. ISSN 19326157.
 ^ Kutchukian, Peter; Lou, David; Shakhnovich, Eugene (2009). "FOG: Fragment Optimized Growth Algorithm for the de Novo Generation of Molecules occupying Druglike Chemical". Journal of Chemical Information and Modeling. 49 (7): 1630–1642. doi:10.1021/ci9000458. PMID 19527020.
 ^ Kutchukian, Peter S.; Lou, David; Shakhnovich, Eugene I. (20090615). "FOG: Fragment Optimized Growth Algorithm for the de Novo Generation of Molecules Occupying Druglike Chemical Space". Journal of Chemical Information and Modeling. 49 (7): 1630–1642. doi:10.1021/ci9000458. ISSN 15499596.
 ^ Kopp, V. S.; Kaganer, V. M.; Schwarzkopf, J.; Waidick, F.; Remmele, T.; Kwasniewski, A.; Schmidbauer, M. (2011). "Xray diffraction from nonperiodic layered structures with correlations: Analytical calculation and experiment on mixed Aurivillius films". Acta Crystallographica Section A. 68 (Pt 1): 148–155. Bibcode:2012AcCrA..68..148K. doi:10.1107/S0108767311044874. PMID 22186291.
 ^ Aguiar, R. J.; CollaresPereira, M.; Conde, J. P. (1988). "Simple procedure for generating sequences of daily radiation values using a library of Markov transition matrices". Solar Energy. 40: 269–279. doi:10.1016/0038092X(88)900497.
 ^ Ngoko, B. O.; Sugihara, H.; Funaki, T. (2014). "Synthetic generation of high temporal resolution solar radiation data using Markov models". Solar Energy. 103: 160–170. doi:10.1016/j.solener.2014.02.026.
 ^ Bright, J. M.; Smith, C. I.; Taylor, P. G.; Crook, R. (2015). "Stochastic generation of synthetic minutely irradiance time series derived from mean hourly weather observation data". Solar Energy. 115: 229–242. doi:10.1016/j.solener.2015.02.032.
 ^ Munkhammar, J.; Widén, J. (2018). "An Nstate Markovchain mixture distribution model of the clearsky index,". Solar Energy. 173: 487–495. doi:10.1016/j.solener.2018.07.056.
 ^ Morf, H. (1998). "The stochastic twostate solar irradiance model (STSIM)". Solar Energy. 62: 101–112. Bibcode:1998SoEn...62..101M. doi:10.1016/S0038092X(98)000048.
 ^ Munkhammar, J.; Widén, J. (2018). "A Markovchain probability distribution mixture approach to the clearsky index". Solar Energy. 170: 174–183. doi:10.1016/j.solener.2018.05.055.
 ^ Pratas, D; Silva, R; Pinho, A; Ferreira, P (May 18, 2015). "An alignmentfree method to find and visualise rearrangements between pairs of DNA sequences". Scientific Reports (Group Nature). 5 (10203): 10203. Bibcode:2015NatSR...510203P. doi:10.1038/srep10203. PMC 4434998. PMID 25984837.
 ^ O'Connor, John J.; Robertson, Edmund F., "Markov chain", MacTutor History of Mathematics archive, University of St Andrews.
 ^ S. P. Meyn, 2007. Control Techniques for Complex Networks Archived 20150513 at the Wayback Machine, Cambridge University Press, 2007.
 ^ U.S. Patent 6,285,999
 ^ Gupta, Brij; Agrawal, Dharma P.; Yamaguchi, Shingo (16 May 2016). Handbook of Research on Modern Cryptographic Solutions for Computer and Cyber Security. IGI Global. pp. 448–. ISBN 9781522501060. Archived from the original on 23 March 2017.
 ^ Langville, Amy N.; Meyer, Carl D. (2006). "A Reordering for the PageRank Problem" (PDF). SIAM Journal on Scientific Computing. 27 (6): 2112–2113. CiteSeerX 10.1.1.58.8652. doi:10.1137/040607551. ISSN 10648275. Archived (PDF) from the original on 20170921.
 ^ Page, Lawrence and Brin, Sergey and Motwani, Rajeev and Winograd, Terry (1999). The PageRank Citation Ranking: Bringing Order to the Web (Technical report). Archived from the original on 20160303. Retrieved 20160304.CS1 maint: multiple names: authors list (link)
 ^ Prasad, NR; RC Ender; ST Reilly; G Nesgos (1974). "Allocation of resources on a minimized cost basis". 1974 IEEE Conference on Decision and Control Including the 13th Symposium on Adaptive Processes. 13: 402–3. doi:10.1109/CDC.1974.270470.
 ^ Hamilton, James (1989). "A new approach to the economic analysis of nonstationary time series and the business cycle". Econometrica. 57 (2): 357–84. CiteSeerX 10.1.1.397.3582. doi:10.2307/1912559. JSTOR 1912559.
 ^ Calvet, Laurent E.; Fisher, Adlai J. (2001). "Forecasting Multifractal Volatility". Journal of Econometrics. 105 (1): 27–58. doi:10.1016/S03044076(01)000690.
 ^ Calvet, Laurent; Adlai Fisher (2004). "How to Forecast longrun volatility: regimeswitching and the estimation of multifractal processes". Journal of Financial Econometrics. 2: 49–83. CiteSeerX 10.1.1.536.8334. doi:10.1093/jjfinec/nbh003.
 ^ Brennan, Michael; Xiab, Yihong. "Stock Price Volatility and the Equity Premium" (PDF). Department of Finance, the Anderson School of Management, UCLA. Archived from the original (PDF) on 20081228.
 ^ A Markov Chain Example in Credit Risk Modelling Columbia University lectures Archived 20160324 at the Wayback Machine
 ^ Acemoglu, Daron; Georgy Egorov; Konstantin Sonin (2011). "Political model of social evolution". Proceedings of the National Academy of Sciences. 108: 21292–21296. Bibcode:2011PNAS..10821292A. CiteSeerX 10.1.1.225.6090. doi:10.1073/pnas.1019454108. PMC 3271566. PMID 22198760. Archived from the original on 20130415.
 ^ Gibson, Matthew C; Patel, Ankit P.; Perrimon, Norbert; Perrimon, Norbert (2006). "The emergence of geometric order in proliferating metazoan epithelia". Nature. 442 (7106): 1038–1041. Bibcode:2006Natur.442.1038G. doi:10.1038/nature05014. PMID 16900102.
 ^ George, Dileep; Hawkins, Jeff (2009). Friston, Karl J. (ed.). "Towards a Mathematical Theory of Cortical Microcircuits". PLoS Comput Biol. 5 (10): e1000532. Bibcode:2009PLSCB...5E0532G. doi:10.1371/journal.pcbi.1000532. PMC 2749218. PMID 19816557.
 ^ Gupta, Ankur; Rawlings, James B. (April 2014). "Comparison of Parameter Estimation Methods in Stochastic Chemical Kinetic Models: Examples in Systems Biology". AIChE Journal. 60 (4): 1253–1268. doi:10.1002/aic.14409. PMC 4946376. PMID 27429455.
 ^ Watterson, G (1996). "Motoo Kimura's Use of Diffusion Theory in Population Genetics". Theoretical Population Biology. 49 (2): 154–188. doi:10.1006/tpbi.1996.0010. PMID 8813021.
 ^ K McAlpine; E Miranda; S Hoggar (1999). "Making Music with Algorithms: A CaseStudy System". Computer Music Journal. 23 (2): 19–30. doi:10.1162/014892699559733.
 ^ Curtis Roads (ed.) (1996). The Computer Music Tutorial. MIT Press. ISBN 9780262181587.CS1 maint: extra text: authors list (link)
 ^ Xenakis, Iannis; Kanach, Sharon (1992) Formalized Music: Mathematics and Thought in Composition, Pendragon Press. ISBN 1576470792
 ^ Continuator Archived July 13, 2012, at the Wayback Machine
 ^ Pachet, F.; Roy, P.; Barbieri, G. (2011) "FiniteLength Markov Processes with Constraints" Archived 20120414 at the Wayback Machine, Proceedings of the 22nd International Joint Conference on Artificial Intelligence, IJCAI, pages 635–642, Barcelona, Spain, July 2011
 ^ Pankin, Mark D. "MARKOV CHAIN MODELS: THEORETICAL BACKGROUND". Archived from the original on 20071209. Retrieved 20071126.
 ^ Pankin, Mark D. "BASEBALL AS A MARKOV CHAIN". Archived from the original on 20090513. Retrieved 20090424.
 ^ Poet's Corner – Fieralingue Archived December 6, 2010, at the Wayback Machine
 ^ Kenner, Hugh; O'Rourke, Joseph (November 1984). "A Travesty Generator for Micros". BYTE. 9 (12): 129–131, 449–469.
 ^ Hartman, Charles (1996). Virtual Muse: Experiments in Computer Poetry. Hanover, NH: Wesleyan University Press. ISBN 9780819522399.
 ^ Gagniuc, Paul A. (2017). Markov chains: from theory to implementation and experimentation. Hoboken, NJ: John Wiley & Sons. pp. 131–163. ISBN 9781119387572. OCLC 982373850.
 ^ Pratas, Diogo; Bastos, Carlos; Pinho, Armando; Neves, Antonio; Matos, Luis (June 2011). DNA synthetic sequences generation using multiple competing Markov models. Statistical Signal Processing Workshop (SSP), 2011 IEEE. 9 (12). pp. 133–136. doi:10.1109/SSP.2011.5967639.
 ^ ^{a} ^{b} de Souza e Silva, E.G.; Legey, L.F.L.; de Souza e Silva, E.A. (2010). "Forecasting oil price trends using wavelets and hidden Markov models". Energy Economics. 32.
 ^ ^{a} ^{b} Carpinone, A; Giorgio, M; Langella, R.; Testa, A. (2015). "Markov chain modeling for veryshortterm wind power forecasting". Electric Power Systems Research. 122.
 ^ ^{a} ^{b} Munkhammar, J.; van der Meer, D.W.; Widén, J. (2019). "Probabilistic forecasting of highresolution clearsky index timeseries using a Markovchain mixture distribution model". Solar Energy. 184.
References[edit]
 A. A. Markov (1906) "Rasprostranenie zakona bol'shih chisel na velichiny, zavisyaschie drug ot druga". Izvestiya Fizikomatematicheskogo obschestva pri Kazanskom universitete, 2ya seriya, tom 15, pp. 135–156.
 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: Markov, A. A. (2006). Translated by Link, David. "An Example of Statistical Investigation of the Text Eugene Onegin Concerning the Connection of Samples in Chains". Science in Context. 19 (4): 591–600. doi:10.1017/s0269889706001074.
 Leo Breiman (1992) [1968] Probability. Original edition published by AddisonWesley; reprinted by Society for Industrial and Applied Mathematics ISBN 0898712963. (See Chapter 7)
 J. L. Doob (1953) Stochastic Processes. New York: John Wiley and Sons ISBN 0471523690.
 S. P. Meyn and R. L. Tweedie (1993) Markov Chains and Stochastic Stability. London: SpringerVerlag ISBN 0387198326. online: https://web.archive.org/web/20100619010320/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 9780521884419. Appendix contains abridged Meyn & Tweedie. online: https://web.archive.org/web/20100619011046/https://netfiles.uiuc.edu/meyn/www/spm_files/CTCN/CTCN.htmlBooth, Taylor L. (1967). Sequential Machines and Automata Theory (1st ed.). New York, NY: John Wiley and Sons, Inc. Library of Congress Card Catalog Number 6725924. Extensive, wideranging 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 Ztransforms, D transforms in their context.
 Kemeny, John G.; Hazleton Mirkil; J. Laurie Snell; Gerald L. Thompson (1959). Finite Mathematical Structures (1st ed.). Englewood Cliffs, NJ: PrenticeHall, Inc. Library of Congress Card Catalog Number 5912841. Classical text. cf Chapter 6 Finite Markov Chains pp. 384ff.
 John G. Kemeny & J. Laurie Snell (1960) Finite Markov Chains, D. van Nostrand Company ISBN 0442043287
 E. Nummelin. "General irreducible Markov chains and nonnegative operators". Cambridge University Press, 1984, 2004. ISBN 052160494X
 Seneta, E. Nonnegative 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 9780387297651
 Kishor S. Trivedi, Probability and Statistics with Reliability, Queueing, and Computer Science Applications, John Wiley & Sons, Inc. New York, 2002. ISBN 0471333417.
 K. S. Trivedi and R.A.Sahner, SHARPE at the age of twentytwo, 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 examplebased approach using the SHARPE software package, Kluwer Academic Publishers, 1996. ISBN 0792396502.
 G. Bolch, S. Greiner, H. de Meer and K. S. Trivedi, Queueing Networks and Markov Chains, John Wiley, 2nd edition, 2006. ISBN 9780792396505.
External links[edit]
 Introduction to Markov Chains on YouTube
 Hazewinkel, Michiel, ed. (2001) [1994], "Markov chain", Encyclopedia of Mathematics, Springer Science+Business Media B.V. / Kluwer Academic Publishers, ISBN 9781556080104
 Techniques to Understand Computer Simulations: Markov Chain Analysis
 Markov Chains chapter in American Mathematical Society's introductory probability book(pdf)
 A beautiful visual explanation of Markov Chains
 Making Sense and Nonsense of Markov Chains
 Original paper by A.A Markov(1913): An Example of Statistical Investigation of the Text Eugene Onegin Concerning the Connection of Samples in Chains (translated from Russian)