hdrcde (version 3.5.0)

plot.cde: Plots conditional densities

Description

Produces stacked density plots or highest density region plots for a univariate density conditional on one covariate.

Usage

# S3 method for cde
plot(
  x,
  firstvar = 1,
  mfrow = n2mfrow(dim(x$z)[firstvar]),
  plot.fn = "stacked",
  x.name,
  margin = NULL,
  ...
)

Value

If plot.fn=="stacked" and there is only one conditioning variable, the function returns the output from graphics::persp(). If plot.fn=="hdr" and there is only one conditioning variable, the function returns the output from hdr.cde(). When there is more than one conditioning variable, nothing is returned.

Arguments

x

Output from cde().

firstvar

If there is more than one conditioning variable, firstvar specifies which variable to fix first.

mfrow

If there is more than one conditioning variable, mfrow is passed to graphics::par() before plotting.

plot.fn

Specifies which plotting function to use: "stacked" results in stacked conditional densities and "hdr" results in highest density regions.

x.name

Name of x (conditioning) variable for use on x-axis.

margin

Marginal density of conditioning variable. If present, only conditional densities corresponding to non-negligible marginal densities will be plotted.

...

Additional arguments to plot.

Author

Rob J Hyndman

References

Hyndman, R.J., Bashtannyk, D.M. and Grunwald, G.K. (1996) "Estimating and visualizing conditional densities". Journal of Computational and Graphical Statistics, 5, 315-336.

See Also

hdr.cde(), cde(), hdr()

Examples

faithful.cde <- cde(faithful$waiting, faithful$eruptions,
  x.name = "Waiting time", y.name = "Duration time"
)
plot(faithful.cde)
plot(faithful.cde, plot.fn = "hdr")