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

test_that("sim_response returns binary matrix with expected dimensions", {
  Q <- sim_qmatrix(5, 3)
  alpha <- sim_students(10, 3)
  Y <- sim_response(Q, alpha, slip = rep(0.1, 5), guess = rep(0.2, 5))
  expect_true(is.matrix(Y))
  expect_equal(dim(Y), c(10, 5))
  expect_true(all(Y %in% c(0, 1)))
})