Skip to contents

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

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

Value

If plot.fn=="stacked" and there is only one conditioning variable, the function returns the output from 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.

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

Author

Rob J Hyndman

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")