Clustering based on the MCMC output of the mixture model
TO BE ADDED.
This function only works for models with a fixed number of mixture components.
NMixCluster(object, ...) ## Default S3 method: NMixCluster(object, ...) ## S3 method for class 'GLMM_MCMC' NMixCluster(object, prob = c("poster.comp.prob", "quant.comp.prob", "poster.comp.prob_b", "quant.comp.prob_b", "poster.comp.prob_u"), pquant = 0.5, HPD = FALSE, pHPD = 0.95, pthresh = -1, unclass.na = FALSE, ...)
object |
an object of apropriate class. |
prob |
character string which identifies estimates of the component probabilities to be used for clustering. |
pquant |
when |
HPD |
logical value. If |
pHPD |
credible level of the HPD credible interval, see argument |
pthresh |
an optional threshold for the estimated component
probability (when |
unclass.na |
logical value taken into account when |
... |
optional additional arguments. |
A data.frame
with three (when HPD
is FALSE
) or five
(when HPD
is TRUE
) columns.
Arnošt Komárek arnost.komarek[AT]mff.cuni.cz
## TO BE ADDED.
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.