R: Plotting Methods for 'norMix' Objects

plot.norMix {nor1mix} | R Documentation

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

[Package nor1mix version 1.2-3 Index]