ip <- generate_ip(model = sample(c("GPCM", "2PL"), 20, TRUE))
n_examinee <- 1000
true_theta <- rnorm(n_examinee)
observed_theta <- true_theta + runif(n_examinee, -.5, .5)
theta_cs <- c(-1, 0, 1.5)
raw_cs <- round(rsss(ip = ip, scale_score = theta_cs))
true_raw_score <- round(rsss(ip = ip, scale_score = true_theta))
observed_raw_score <- round(rsss(ip = ip, scale_score = observed_theta))
# Theta scores
classification_agreement_index(true_score = true_theta,
estimated_score = observed_theta,
cut_scores = theta_cs)
# Summed scores
classification_agreement_index(true_score = true_raw_score,
estimated_score = observed_raw_score,
cut_scores = raw_cs)
# Add labels
classification_agreement_index(true_score = true_theta,
estimated_score = observed_theta,
cut_scores = theta_cs,
cat_labels = c("Unsatisfactory", "Basic",
"Mastery", "Advanced"))
Run the code above in your browser using DataLab