simulate_sequence_clusters: Simulate sequences from a mixture of transition systems

In simulab: Simulates a Wide Variety of Data Shapes, Sizes and Distributions

Simulate sequences from a mixture of transition systems

Description

Simulate sequences from a mixture of transition systems

Usage

simulate_sequence_clusters(
  n,
  transitions,
  chain_length,
  proportions = NULL,
  initial = NULL,
  labels = NULL,
  states = NULL,
  seed = NULL,
  batch = NULL
)

Arguments

n	Number of sequences.
transitions	List of transition matrices, or a tidy data frame with columns cluster, from, to and probability.
chain_length	Sequence length.
proportions	Cluster proportions. When NULL, clusters are equally likely.
initial	Optional common initial probabilities, shared by every cluster. When NULL, every sequence starts in the first state.
labels	Optional sequence-cluster labels.
states	Optional state labels. When the matrices carry no dimnames and states is NULL, states are labeled by their row index.
seed	Optional random seed.
batch	Optional single positive whole number. When given, the simulator runs batch times and returns a plain list of batch results, each exactly what the same call without batch returns. With a seed, every dataset gets its own seed drawn from seed, so the whole batch is reproducible; without one the datasets are consecutive draws from the session's random-number stream. The default NULL returns a single result.

Value

A long-form simulab_sim base data.frame with columns id, period, state, and the true sequence_cluster, one row per sequence position. A transitions table (sequence_cluster, from, to, probability) is available through as.data.frame().

Examples

# One transition matrix per latent cluster.
result <- simulate_sequence_clusters(
  n = 60,
  transitions = list(
    matrix(c(0.8, 0.2, 0.3, 0.7), nrow = 2, byrow = TRUE),
    matrix(c(0.3, 0.7, 0.6, 0.4), nrow = 2, byrow = TRUE)
  ),
  chain_length = 12,
  seed = 1
)
head(result)
components(result)

simulab documentation built on Oct. 15, 2026, 5:08 p.m.