simulateHMM: Simulate a hidden Markov model path (HMM)

In plbaldoni/mixNBHMM: Differential Peak Calling via Hidden Markov Model With Mixture of Negative Binomial Distributions

Description

This function simulates a HMM path of size 'M' given a matrix of transition probabilities 'gamma'. By default, the path always begins at the first state of the HMM.

Usage

simulateHMM(gamma, M)

Arguments

gamma	a square matrix of transition probabiltiies such that all(rowSums(P)==1) is TRUE
M	length of the HMM path

Value

A vector of integers of length 'M' representing states of the HMM

Author(s)

Pedro L. Baldoni, [e-mail]

References

https://github.com/plbaldoni/mixNBHMM

Examples

rle(simulateHMM(gamma = matrix(c(0.95,0.05,0.05,0.95),2,2),M = 1e4))

plbaldoni/mixNBHMM documentation built on Dec. 24, 2019, 1:31 p.m.