Density or uncertainty surface for bivariate mixtures
Plots a density or uncertainty surface given bivariate data and parameters of a MVN mixture model for the data.
surfacePlot(data, parameters,
what = c("density", "uncertainty"),
type = c("contour", "hdr", "image", "persp"),
transformation = c("none", "log", "sqrt"),
grid = 200, nlevels = 11, levels = NULL,
prob = c(0.25, 0.5, 0.75),
col = gray(0.5),
col.palette = function(...) hcl.colors(..., "blues", rev = TRUE),
hdr.palette = blue2grey.colors,
xlim = NULL, ylim = NULL, xlab = NULL, ylab = NULL,
main = FALSE, scale = FALSE, swapAxes = FALSE,
verbose = FALSE, ...)data |
A matrix, or data frame of bivariate observations. Categorical variables are not allowed. If a matrix or data frame, rows correspond to observations and columns correspond to variables. |
parameters |
A named list giving the parameters of an MCLUST model, used to produce superimposing ellipses on the plot. The relevant components are as follows:
|
what |
Choose from one of the following options: |
type |
Choose from one of the following three options: |
transformation |
Choose from one of the following three options: |
grid |
The number of grid points (evenly spaced on each axis).
The mixture density and uncertainty is computed at
|
nlevels |
The number of levels to use for a contour plot.
Default: |
levels |
A vector of levels at which to draw the lines in a contour plot. |
prob |
A vector of probability levels for computing HDR.
Only used if |
col |
A string specifying the colour to be used for |
col.palette |
A function which defines a palette of colours to be used for
|
hdr.palette |
A function which defines a palette of colours to be used for
|
xlim, ylim |
Optional argument specifying bounds for the ordinate, abscissa of the plot. This may be useful for when comparing plots. |
xlab, ylab |
Optional argument specifying labels for the x-axis and y-axis. |
main |
A logical variable or |
scale |
A logical variable indicating whether or not the two dimensions should be plotted on the same scale, and thus preserve the shape of the distribution. The default is not to scale. |
swapAxes |
A logical variable indicating whether or not the axes should be swapped for the plot. |
verbose |
A logical variable telling whether or not to print an indication that the function is in the process of computing values at the grid points, which typically takes some time to complete. |
... |
Other graphics parameters. |
For an image plot, a color scheme may need to be selected on the display device in order to view the plot.
A plots showing (a transformation of) the density or uncertainty for the given mixture model and data.
The function also returns an invisible list with components x,
y, and z in which x and y are the values used to
define the grid and z is the transformed density or uncertainty at the
grid points.
C. Fraley and A. E. Raftery (2002). Model-based clustering, discriminant analysis, and density estimation. Journal of the American Statistical Association 97:611-631.
C. Fraley, A. E. Raftery, T. B. Murphy and L. Scrucca (2012). mclust Version 4 for R: Normal Mixture Modeling for Model-Based Clustering, Classification, and Density Estimation. Technical Report No. 597, Department of Statistics, University of Washington.
faithfulModel <- Mclust(faithful)
surfacePlot(faithful, parameters = faithfulModel$parameters,
type = "contour", what = "density", transformation = "none",
drawlabels = FALSE)
surfacePlot(faithful, parameters = faithfulModel$parameters,
type = "persp", what = "density", transformation = "log")
surfacePlot(faithful, parameters = faithfulModel$parameters,
type = "contour", what = "uncertainty", transformation = "log")Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.