plot.cint: Plot confidence intervals for a Gaussian scalar mixture...

In mixreg: Functions to Fit Mixtures of Regressions

Description

Plots confidence and prediction intervals, corresponding to each component of a Gaussian scalar mixture model.

Usage

## S3 method for class 'cint'
plot(x, cints=TRUE, pints=TRUE,rugged=TRUE,
                     type=c("both","upper","lower"),...)

Arguments

x	An object of class "cint", specifying the confidence and prediction intervals to be plotted, as produced by the undocumented function cint(). Such an object would usually be produced (indirectly) from a call to cband(). If the class "mixreg" object, to which cband() is applied, specifies a fit to a model involving no predictors, then confidence and prediction bands make no sense, and the object is handed over to cint() to produce confidence and prediction intervals.
cints	Logical scalar; should the confidence intervals be plotted?
pints	Logical scalar; should the prediction intervals be plotted?
rugged	Logical scalar. Should a “rug” be added to the plot, displaying the values of the variable to which the model was fitted? See rug(). If rugged is TRUE then a rug is added at the right hand margin of the plot.
type	Character string specifying whether the intervals plotted should be two-sided ("both") or if only the upper ("upper" or lower ("lower") bounds should be plotted. May be abbreviated (e.g. to "b", "u" or "l").
...	Optional extra arguments for plot() and title() These may include xlim, ylim, xlab, ylab, col, pch, and main. Note that col is replicated to have length K where K is the number of components in the model. If main is not supplied then a “sensible” default main title is created. If no main title is desired, specify main="".

Details

This function is a method for plot. Note that a simple plot of the fit may be produced by calling plot(object) where object is an object of class "mixreg" and the x argument of this function was produced by applying cband() to this object.

Value

None. This function is called for its side effect of producing a plot.

Author(s)

Rolf Turner [e-mail]

See Also

cband(), plot.cband(), plot.mixreg(), plot.mixresid(), qqMix(), residuals.mixreg()

Examples

    fit  <- mixreg(plntsInf~1,ncomp=2,seed=42,covMat=TRUE,data=aphids)
    cis  <- cband(fit) # *Could* call cint(fit) directly 
    plot(cis)
    plot(cis,pintss=FALSE)
    plot(cis,pintss=FALSE,type="u")

mixreg documentation built on Oct. 14, 2021, 9:12 a.m.