# NOT RUN {
data(hiv)
npobj <- npmle(hiv, family = gaussian, maxiter = 4)
### Compute unit-specific posterior means
pmean <- npmixapply(npobj, function(x) { x })
### Compute post. prob that \theta_i < .1
pp <- npmixapply(npobj, function(x) { x < .1})
# }
Run the code above in your browser using DataLab