simulate.dmm: Simulate a sequence under a drifting Markov model

View source: R/dmm.R

simulate.dmmR Documentation

Simulate a sequence under a drifting Markov model

Description

Simulate a sequence under a k-th order DMM.

Usage

## S3 method for class 'dmm'
simulate(x, output_file = NULL, model_size = NULL, ncpu = 2)

Arguments

x

An object of class dmm

output_file

(Optional) File containing the simulated sequence (e.g, "C:/.../SIM.txt")

model_size

Size of the model

ncpu

Default=2. Represents the number of cores used to parallelized computation. If ncpu=-1, then it uses all available cores.

Value

the vector of simulated sequence

Author(s)

Annthomy Gilles, Alexandre Seiller

References

\insertRef

BaVe2018drimmR \insertRefVer08drimmR

See Also

fitdmm, getTransitionMatrix, getStationaryLaw

Examples

data(lambda, package = "drimmR")
dmm <- fitdmm(lambda, 1, 1, c('a','c','g','t'), init.estim = "freq", fit.method="sum")
simulate(dmm, model_size=100)

drimmR documentation built on Oct. 9, 2026, 5:08 p.m.