# \donttest{
J <- 5
K <- 4
pi <- rep(1 / K, K)
mu <- matrix(0, nrow = K, ncol = K)
diag(mu) <- 5
sigma <- matrix(sqrt(2) / sqrt(pi), nrow = K, ncol = K)
sim_data <- data.frame(item = rep(1:2, each = 5), rater = rep(1:5, 2))
sim_result <- simulate_hier_dawid_skene_model(pi, mu, sigma, sim_data)
sim_result$sim
sim_result$theta
# }
Run the code above in your browser using DataLab