npmixapply | R Documentation |
Using a nonparametric estimate of the mixing distribution, computes a posterior quantity of interest for each unit.
npmixapply(object, FUN, ...)
object |
an object of class "npmix" |
FUN |
a user provided function |
... |
optional arguments to FUN |
object
is an object of class "npmix" containing
a nonparametric estimate of the mixing distribution F
in the following two-level sampling model
X_i|θ_i ~ p(x|θ_i,η_i)
and θ_i ~ F for i = 1,...,n.
Using npmixapply(object, f)
, then returns the
posterior expectation of f:
E[ f(θ_i) | X_i, η_i] , for i = 1,...,n.
a vector with length equal to n
Nicholas Henderson
npmle
## 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}) ## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.