x <- freqtab(KBneat$x, scales = list(0:36, 0:12))
y <- freqtab(KBneat$y, scales = list(0:36, 0:12))
set.seed(111213)
eq1 <- equate(x, y, type = "linear", method = "tuck",
boot = TRUE, reps = 10)
eq2 <- equate(x, y, type = "equip", method = "freq",
boot = TRUE, smooth = "log", reps = 10,
showWarnings = FALSE)
eq3 <- equate(x, y, type = "linear", method = "braun",
boot = TRUE, reps = 10, showWarnings = FALSE)
plot(eq1, eq2, eq3)
plot(eq1, eq2, eq3, out = "se", legendplace = "top",
addident = FALSE)
Run the code above in your browser using DataLab