Function to generate a sequence of states from homogeneous Markov chains.

Provided any markovchain object, it returns a sequence of states coming from the underlying stationary distribution.

markovchainSequence(
  n,
  markovchain,
  t0 = sample(markovchain@states, 1),
  include.t0 = FALSE,
  useRCpp = TRUE
)

Arguments

n

Sample size

markovchain

markovchain object

t0

The initial state

include.t0

Specify if the initial state shall be used

useRCpp

Boolean. Should RCpp fast implementation being used? Default is yes.

Value

A Character Vector

Details

A sequence of size n is sampled.

References

A First Course in Probability (8th Edition), Sheldon Ross, Prentice Hall 2010

See also

markovchainFit

Author

Giorgio Spedicato

Examples

# define the markovchain object
statesNames <- c("a", "b", "c")
mcB <- new("markovchain", states = statesNames, 
   transitionMatrix = matrix(c(0.2, 0.5, 0.3, 0, 0.2, 0.8, 0.1, 0.8, 0.1), 
   nrow = 3, byrow = TRUE, dimnames = list(statesNames, statesNames)))

# show the sequence
outs <- markovchainSequence(n = 100, markovchain = mcB, t0 = "a")