Currently experimental function that creates a list object containing one table and two plots showing the range of plausible values at each ordinal sum score level (corresponds to possible theta values). Means and 95% HDCIs are shown in figures and in table. One plot makes mean adjustments for WL estimate thetas, since the plausible values are based on EAP estimates. This is similar to the bias adjustment made when bootstrapping.
Usage
RIcrel(
data,
draws = 500,
n = 10,
conf_level = 0.95,
estim = "WLE",
theta_range = c(-6, 6)
)