Simulate a Markov chain — simulate_markov_chain • fHMM

This function simulates a Markov chain.

Usage

simulate_markov_chain(
  Gamma,
  T,
  delta = Gamma2delta(Gamma),
  seed = NULL,
  total_length = T
)

Arguments

Gamma: A tpm (transition probability matrix).
T: The length of the Markov chain.
delta: A probability vector, the initial distribution. If not specified, delta is set to the stationary distribution vector.
seed: Set a seed.
total_length: An integer, the total length of the output vector. Must be greater or equal than T.

Value

A numeric vector of length T with states.

Examples

Gamma <- matrix(c(0.5, 0.3, 0.5, 0.7), 2, 2)
T <- 10
fHMM:::simulate_markov_chain(Gamma = Gamma, T = T)
#>  [1] 1 1 2 2 2 2 1 1 1 1