test_that("multiplication works", {
  expect_equal(2 * 2, 4)
})

test_that("perturb_qmatrix flips entries as expected", {
  Q <- matrix(1, nrow = 5, ncol = 3)
  Q2 <- perturb_qmatrix(Q, flip_prob = 1)
  expect_false(identical(Q, Q2))
  expect_true(all(Q2 %in% c(0, 1)))
})