View source: R/PLMIXfunctions.R
gibbsPLMIX | R Documentation |
Perform Gibbs sampling simulation for a Bayesian mixture of Plackett-Luce models fitted to partial orderings.
gibbsPLMIX(
pi_inv,
K,
G,
init = list(z = NULL, p = NULL),
n_iter = 1000,
n_burn = 500,
hyper = list(shape0 = matrix(1, nrow = G, ncol = K), rate0 = rep(0.001, G), alpha0 =
rep(1, G)),
centered_start = FALSE
)
pi_inv |
An object of class |
K |
Number of possible items. |
G |
Number of mixture components. |
init |
List of named objects with initialization values: |
n_iter |
Total number of MCMC iterations. |
n_burn |
Number of initial burn-in drawings removed from the returned MCMC sample. |
hyper |
List of named objects with hyperparameter values for the conjugate prior specification: |
centered_start |
Logical: whether a random start whose support parameters and weights should be centered around the observed relative frequency that each item has been ranked top. Default is |
The size L
of the final MCMC sample is equal to n_iter
-n_burn
.
A list of S3 class gsPLMIX
with named elements:
W |
Numeric |
P |
Numeric |
log_lik |
Numeric vector of |
deviance |
Numeric vector of |
objective |
Numeric vector of |
call |
The matched call. |
Cristina Mollica and Luca Tardella
Mollica, C. and Tardella, L. (2017). Bayesian Plackett-Luce mixture models for partially ranked data. Psychometrika, 82(2), pages 442–458, ISSN: 0033-3123, <doi:10.1007/s11336-016-9530-0>.
data(d_carconf)
GIBBS <- gibbsPLMIX(pi_inv=d_carconf, K=ncol(d_carconf), G=3, n_iter=30, n_burn=10)
str(GIBBS)
GIBBS$P
GIBBS$W
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.