stationary: Compute the stationary distribution of a homogeneous Markov chain
Description
A homogeneous, finite state Markov chain that is irreducible and aperiodic converges to a unique stationary distribution, here called \(\delta\).
As it is stationary, this distribution satisfies
$$\delta \Gamma = \delta,$$ subject to \(\sum_{j=1}^N \delta_j = 1\),
where \(\Gamma\) is the transition probability matrix.
This function solves the linear system of equations above.
Usage
stationary(Gamma)
Value
either a single stationary distribution of the Markov chain (vector of length N) or a matrix of stationary distributions of dimension c(nTracks,N) with one stationary distribution in each row
Arguments
Gamma
transition probability matrix of dimension c(N,N) or array of such matrices of dimension c(N,N,nTracks) if the stationary distribution should be computed for several matrices at once
See Also
tpm to create a transition probabilty matrix using the multinomial logistic link (softmax)
Other stationary distribution functions:
stationary_cont(),
stationary_p()