mixAK (version 2.2)

NMixPlugDensJoint2: Pairwise bivariate densities: plug-in estimate

Description

This function serves as an inference tool for the MCMC output obtained using the function NMixMCMC. It computes marginal (pairwise bivariate) plug-in densities obtained by using posterior summary statistics (e.g., posterior means) of mixture weights, means and variances.

Usage

NMixPlugDensJoint2(x, ...)
## S3 method for class 'default':
NMixPlugDensJoint2(x, scale, w, mu, Sigma, \dots)
## S3 method for class 'NMixMCMC':
NMixPlugDensJoint2(x, grid, lgrid=50, scaled=FALSE, \dots)
## S3 method for class 'GLMM_MCMC':
NMixPlugDensJoint2(x, grid, lgrid=50, scaled=FALSE, \dots)

Arguments

x

an object of class NMixMCMC for NMixPlugDensJoint2.NMixMCMC function.

An object of class GLMM_MCMC for NMixPlugDensJoint2.GLMM_MCMC function. A list with the grid values (see below)

scale

a two component list giving the shift and the scale. If not given, shift is equal to zero and scale is equal to one.

w

a numeric vector with posterior summary statistics for the mixture weights. The length of this vector determines the number of mixture components.

mu

a matrix with posterior summary statistics for mixture means in rows. That is, mu has $K$ rows and $p$ columns, where $K$ denotes the number of mixture components and $p$ is dimension of the mixture distribution.

Sigma

a list with posterior summary statistics for for mixture covariance matrices.

grid

a list with the grid values for each margin in which the density should be evaluated.

If grid is not specified, it is created automatically using the information from the posterior summary statistics stored in x.

lgrid

a length of the grid used to create the grid if that is not specified.

scaled

if TRUE, the density of shifted and scaled data is summarized. The shift and scale vector are taken from the scale component of the object x.

...

optional additional arguments.

Value

• An object of class NMixPlugDensJoint2 which has the following components:
• x
  a list with the grid values for each margin. The components of the list are named x1, ...or take names from grid argument.
• dens
  a list with the computed densities for each pair of margins. The components of the list are named 1-2, 1-3, ..., i.e., dens[[1]]$=$dens[["1-2"]] is the pairwise predictive density for margins 1 and 2, etc. Each component of the list is a matrix in such a form that it can be directly passed together with the proper components of x to the plotting functions like contour or image.
• There is also a plot method implemented for the resulting object.

See Also

plot.NMixPlugDensJoint2, NMixMCMC, GLMM_MCMC, NMixPredDensJoint2.