Posterior probabilities/intercepts, and mass point classifications
Description
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.
Usage
1  post(object, level="upper")

Arguments
object 
an object of class 
level 

Value
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 
Author(s)
Jochen Einbeck and John Hinde (2006)
References
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.
See Also
alldist
, allvc
Examples
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
