plot.norMix: Plotting Methods for 'norMix' Objects

In nor1mix: Normal aka Gaussian 1-d Mixture Models

Plotting Methods for 'norMix' Objects

Description

The plot and lines methods for norMix objects draw the normal mixture density, optionally additonally with a fitted normal density.

Usage

## S3 method for class 'norMix'
plot(x, type = "l", n = 511, xout = NULL, xlim = NULL, ylim,
      xlab = "x", ylab = "f(x)", main = attr(x, "name"), lwd = 1.4,
      p.norm = !p.comp, p.h0 = TRUE, p.comp = FALSE,
      parNorm = list(col = 2, lty = 2, lwd = 0.4),
      parH0 = list(col = 3, lty = 3, lwd = 0.4),
      parComp = list(col= "blue3", lty = 3, lwd = 0.4), ...)

## S3 method for class 'norMix'
lines(x, type = "l", n = 511, xout = NULL,
      lwd = 1.4, p.norm = FALSE, parNorm = list(col = 2, lty = 2, lwd = 0.4),
      ...)

Arguments

x	object of class norMix.
type	character denoting type of plot, see, e.g. lines.
n	number of points to generate if xout is unspecified.
xout	numeric or NULL giving the abscissae at which to draw the density.
xlim	range of x values to use; particularly important if xout is not specified where xlim is passed to dnorMix and gets a smart default if unspecified.
ylim	range of y values to use; by default, if not specified (or containing NA), a smart default is used.
xlab, ylab	labels for the x and y axis with defaults.
main	main title of plot, defaulting to the norMix name.
lwd	line width for plotting with a non-standard default.
p.norm	logical indicating if the normal density with the same mean and variance should be drawn as well.
p.h0	logical indicating if the line y = 0 should be drawn.
p.comp	logical indicating if the Gaussian components should also be drawn individually.
parNorm	graphical parameters for drawing the normal density if p.norm is true.
parH0	graphical parameters for drawing the line y = 0 if p.h0 is true.
parComp	graphical parameters for drawing the single components if p.comp is true.
...	further arguments passed to and from methods.

Author(s)

Martin Maechler

See Also

norMix for the construction and further methods, particularly dnorMix which is used here.

Examples

plot(norMix(m=c(0,3), sigma = c(2,1))) # -> var = c(2^2, 1) = c(4, 1)

plot(MW.nm4, p.norm=FALSE, p.comp = TRUE)
plot(MW.nm4, p.norm=FALSE, p.comp = TRUE, ylim = c(0, 2))# now works
stopifnot(all.equal(c(0,2), par("yaxp")[1:2], tol= 1e-15))

## Further examples in  ?norMix and  ?rnorMix

nor1mix documentation built on May 29, 2024, 8:28 a.m.