Takes an object of class glmmNPML
or glmmGQ
and displays the
posterior probabilites w_ik as well as the posterior intercepts
(Sofroniou et. al, 2006). Further it classfies the observations to mass points
according to their posterior probability. The level on which the information
in all three cases is displayed can be chosen by the user via the 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.
1  post(object, level="upper")

object 
an object of class 
level 

A list of the following four items:
prob 
posterior probabilities (identical to 
int 
posterior intercepts (identical to 
classif 
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 
level 
either 
Jochen Einbeck and John Hinde (2006)
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.
alldist
, allvc
1 2 3 4 5  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

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.
All documentation is copyright its authors; we didn't write any of that.