TO BE ADDED.
This function only works for models with a fixed number of mixture components.
NMixCluster(object, ...)# S3 method for default
NMixCluster(object, ...)
# S3 method for 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, ...)
A data.frame
with three (when HPD
is FALSE
) or five
(when HPD
is TRUE
) columns.
an object of apropriate class.
character string which identifies estimates of the component probabilities to be used for clustering.
when prob
is either “quant.comp.prob” or
“quant.comp.prob_b”, argument pquant
is the
probability of the quantile of the component probabilities to be
used for clustering.
logical value. If TRUE
then only those subjects are
classified for which the lower limit of the pHPD
*100% HPD
credible interval of the component probability exceeds the value of ptrash
.
credible level of the HPD credible interval, see argument HPD
.
an optional threshold for the estimated component
probability (when HPD
is FALSE
) or for the lower limit
of the HPD credible interval (when HPD
is TRUE
) to
classify a subject. No effect when pthresh
is negative.
logical value taken into account when pthresh
is positive. If unclass.na
is TRUE
, unclassified subjects get classification
NA
. If unclass.na
is FALSE
, unclassified
subjects create a separate (last) group.
optional additional arguments.
Arnošt Komárek arnost.komarek@mff.cuni.cz
NMixMCMC
, GLMM_MCMC
.