Fix the description of the number of raters in the anesthesia data. The documentation had erroneously stated there were four anesthetists, not five.
Update the stan code for compatibility with rstan v2.26.0 (@andrjohns)
Updated the Stan implementation, priors, and initialisation points of the hierarchical Dawid-Skene model, leading to much more reliable convergence.
Added the ability to visualise the theta parameter with uncertainty.
Added row names to the output of
class_probabilities().
Added the ability to specify the column names of long format data
passed to rater().
Added simulate_dawid_skene_model() and
simulate_hier_dawid_skene_model() to simulate data from the
Dawid-Skene and hierarchical Dawid-Skene models.
Re-export loo_compare().
Allowed the theta parameter to be extracted from the hierarchical Dawid-Skene model.
Add waic() function for model comparison
Silence warnings with the latest ggplot2 version
Fix validation bug in posterior_predict()
summary() now works with the class conditional and
hierarchical Dawid-Skene models.
All functions applied to fitted class conditional Dawid-Skene models will automatically convert the relevant parameters of the model into a full theta parameter equivalent to the Dawid-Skene model. This is designed to allow easier comparison of the class conditional model with the full Dawid-Skene model.
Plotting via plot() of the rater_fit
object has been changed in several ways. plot.rater_fit
now:
prob, which (called
rater_index) and new item_index arguments in
the plot generic.Add the ability to only plot a subset of items when plotting the
class probabilities. This can be controlled by the new
item_index argument to plot()
Added the function wide_to_long() to convert wide
data to long data.
Add the option data_format = "wide" to
rater() to allow wide data to be passed into
rater() directly.
Added the get_stanfit() function to extract the
underlying stanfit object from a rater fit object.
Added an implementation of the posterior_predict
generic from {rstantools} allowing simulation from the posterior
predictive distribution of fitted standard, and class conditional,
Dawid-Skene models. (The hierarchical Dawid-Skene model is not yet
supported).
Added an implementation of the prior_summary generic
from {rstantools} for rater_fit objects.
Add the loo.rater_fit method to allow the
calculation of loo, a modern Bayesian model comparison metric, for rater
models. loo values can be compared using the excellent {loo}
package.
Added the loo.rater_fit method to allow the
calculation of loo, a modern Bayesian model comparison metric, for rater
models. loo values can be compared using the excellent {loo}
package.
Rater specific prior parameters can now be used in the
Dawid-Skene model for both grouped and long data. In practice this means
that it is now possible to pass a J * K * K array for beta
into dawid_skene() which encodes a K * K prior parameter
for each of the J raters’ error matrices. For backwards compatibility
and ease of use it is still possible to pass a single matrix for
beta which will still be interpreted as the prior parameter
for all the of the raters’ error matrices.
The plot produced for the pi parameter has been changed. The new plot represents the uncertainty in the point estimates when MCMC has been used to fit the model.
Prior parameters for the Dawid-Skene and class conditional Dawid-Skene models have been altered slightly to improve convergence of optimization when the number of classes is small.
summary.mcmc_fit now displays the number of
remaining parameters correctly.
Added the as_mcmc.list() function to convert MCMC
fits to {coda} mcmc.list objects.