
Takes an object of class glmmNPML
or glmmGQ
and displays the
posterior probabilites level
argument ("upper"
or "lower"
). The actual information in both cases is
identical, the latter is just an expanded version of the former. In case of
simple overdispersion models, the level
argument is not relevant.
post(object, level="upper")
an object of class glmmNPML
or glmmGQ
.
"upper"
or "lower"
.
A list of the following four items:
posterior probabilities (identical to object$post.prob
in case of "lower"
and for one-level models).
posterior intercepts (identical to object$post.int
in case of "lower"
and for one-level models).
a numerical vector containing the class numbers (the order of the classes corresponds to the
order of the mass points given in the output of alldist
or allvc
).
either "lower"
, "upper"
, or "none"
(for one-level models).
Sofroniou, N., Einbeck, J., and Hinde, J. (2006). Analyzing Irish suicide rates with mixture models. Proceedings of the 21st International Workshop on Statistical Modelling in Galway, Ireland, 2006.
# NOT RUN {
data(galaxies, package="MASS")
gal <- as.data.frame(galaxies)
post(alldist(galaxies/1000~1, random=~1, data=gal, k=5))$classif
# classifies the 82 galaxies to one of the five mass points
# }
Run the code above in your browser using DataLab