pairs.profile {ellipse}R Documentation

Profile pairs

Description

This function produces pairwise plots of profile traces, profile sketches, and ellipse approximations to confidence intervals.

Usage

## S3 method for class 'profile'
pairs(x, labels = c(names(x), "Profile tau"), panel = lines, invert = TRUE, 
    plot.tau = TRUE, plot.trace = TRUE, plot.sketch = TRUE, 
	plot.ellipse = FALSE, level = 0.95, ...)

Arguments

x

An object of class profile, generally the result of the profile() function.

labels

The labels to use for each variable. These default to the variable names.

panel

The function to use to draw the sketch in each panel.

invert

Whether to swap the axes so things look better.

plot.tau

Whether to do the profile tau (profile t) plots.

plot.trace

Whether to do the profile trace plots.

plot.sketch

Whether to do the profile sketch plots.

plot.ellipse

Whether to do the ellipse approximations.

level

The nominal confidence level for the profile sketches and ellipses.

...

Other plotting parameters.

Details

This function implements the plots used in Bates and Watts (1988) for nonlinear regression diagnostics.

Unfortunately, the MASS package also defines a pairs.profile method. The ellipse::pairs generic is supplied to allow users to choose to use the version in this package.

If x is a profile object, then ellipse::pairs(x) is guaranteed to call the method from this package. If you'd rather use the MASS method, then two steps are needed: you need to be sure that MASS is loaded, then call graphics::pairs(x). (If MASS is not loaded, you'll get the default method from the graphics package, which will trigger an error since it doesn't know about profile objects.)

If both ellipse and MASS are loaded, then pairs(x) is slightly ambiguous: the method that is called depends on the ordering of ellipse and graphics (not MASS!) in the search list. In almost all cases ellipse will precede graphics, so the ellipse method will usually be the default.

If x is any object other than a profile object, ellipse::pairs(x) will pass the call on to graphics::pairs(x), and the standard method will be called.

Side Effects

Produces a plot on the current device for each pair of variables in the profile object.

References

Bates and Watts (1988), Nonlinear Regression Analysis \& its Applications.

See Also

profile, ellipse.profile, ellipse.nls

Examples

 # Plot everything for the Puromycin data
 data(Puromycin)
 Purboth <- nls(formula = rate ~ ((Vm + delV * (state == "treated"))
   * conc)/(K + conc), data = Puromycin,
   start = list(Vm = 160, delV = 40, K = 0.05))
 Pur.prof <- profile(Purboth)
 ellipse::pairs(Pur.prof, plot.ellipse = TRUE)
 
 # Show the corresponding plot from MASS:
 if (requireNamespace("MASS"))
   graphics::pairs(Pur.prof)

[Package ellipse version 0.4.2 Index]