plot.LogConcDEAD: Plot a log-concave maximum likelihood estimator

Description Usage Arguments Details Author(s) See Also

View source: R/plot.LogConcDEAD.R


plot method for class "LogConcDEAD". Plots of various types are available for 1- and 2-d data. For dimension greater than 1, plots of axis-aligned marginal density estimates are available.


## S3 method for class 'LogConcDEAD'
plot(x, uselog=FALSE, type="ic", addp=TRUE,
  drawlabels=TRUE, gridlen=400, g, marg, g.marg, main, xlab, ylab, ...)



Object of class "LogConcDEAD" (typically output from mlelcd)


Scalar logical: should the plot be on the log scale?


Plot type: "p" perspective, "c" contour, "i" image, ic image and contour, r using rgl (the best!)


Scalar logical: should the data points be plotted? (as black dots on the surface for d >= 2; as circles for d=1)


Scalar logical: should labels be added to contour lines? (only relevant for types "ic" and "c")


Integer scalar indicating the number of points at which the maximum likelihood estimator is evaluated in each dimension


(optional) a matrix of density estimate values (the result of a call to interplcd). If many plots of a single dataset are required, it may be quicker to compute the grid using interplcd(x) and pass the result to plot


If non-NULL, this scalar integer determines which marginal should be plotted (should be between 1 and d)


If g is non-NULL, can contain a vector of marginal density estimate values (the output of interpmarglcd). If many plots of a single dataset are required, it may be quicker to compute the marginal values to compute marginal values using interpmarglcd and pass the result to plot




x-axis label


y-axis label


Other arguments to be passed to the generic plot method


The density estimate is evaluated on a grid of points using the interplcd function. If several plots are required, this may be computed separately and passed to plot using the g argument.

For two dimensional data, the default plot type is "ic", corresponding to image and contour plots. These may be obtained separately using plot type "i" or "c" respectively. Where available, the use of plot type "r" is recommended. This uses the rgl package to produce a 3-d plot that may be rotated by the user. The option "p" produces perspective plots.

For data of dimension at least 2, axis-aligned marginals may be plotted by setting the marg argument. This integrates the estimated density over the remaining dimensions. If several plots are required, the estimate may be computed using the function interpmarglcd and passed using the argument g.marg.

Where relevant, the colors were obtained from the function heat_hcl in the package colorspace. Thanks to Achim Zeileis for this suggestion.

For examples, see mlelcd.


Madeleine Cule

Robert B. Gramacy

Richard Samworth

Yining Chen

See Also

mlelcd, interplcd, interpmarglcd, heat_hcl

LogConcDEAD documentation built on April 13, 2018, 9:04 a.m.