R: Pairwise bivariate densities: plug-in estimate

NMixPlugDensJoint2 {mixAK}

R Documentation

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, ...)

## Default S3 method:
NMixPlugDensJoint2(x, scale, w, mu, Sigma, ...)

## S3 method for class 'NMixMCMC'
NMixPlugDensJoint2(x, grid, lgrid=50, scaled=FALSE, ...)

## S3 method for class 'GLMM_MCMC'
NMixPlugDensJoint2(x, grid, lgrid=50, scaled=FALSE, ...)

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) for `NMixPlugDensJoint2.default` function.
`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.

Author(s)