plot.nlreg.profiles function

Use plot() on a profile.nlreg' and all.profiles.nlreg' object

These are methods for the function plot for objects inheriting from class "profile.nlreg" or "all.profiles.nlreg" .

## S3 method for class 'nlreg.profile'
plot(x, alpha = 0.05, add.leg = FALSE, stats = c("sk", "fr"), 
     cex = 0.7, cex.lab = 1, cex.axis = 1, cex.main = 1, lwd1 = 1, 
     lwd2 = 2, lty1 = "solid", lty2 = "solid", cl1 = "blue", 
     cl2 = "red", col = "black", ylim = c(-3,3), ...)
## S3 method for class 'all.nlreg.profiles'
plot(x, nframe, alpha = 0.05, stats = c("sk", "fr"), cex = 0.7, 
     cex.lab = 1, cex.axis = 1, cex.main = 1, lwd1 = 1, lwd2 = 2, 
     lty1 = "solid", lty2 = "solid", cl1 = "blue", cl2 = "red", 
     col = "black", ylim = c(-3,3), ...)

Arguments

• x: an object of class profile.nlreg or all.profiles.nlreg such as generated by a call to profile.nlreg.

• nframe: the number of frames into which to split the graphics device; only if x is an all.profiles.nlreg object.

• alpha: numeric vector with the levels used to read off confidence intervals; the default is 5% which corresponds to a confidence level of $1-\code{alpha}=0.95$.

• stats: character value indicating which higher order statistics to plot. Admissible values are "sk" for Skovgaard's (1996)

  proposal and "fr" for Fraser, Reid and Wu's (1999)

  approach. The default is "sk".

• add.leg: logical value indicating whether a legend should be added to the plot; only if x is a profile.nlreg object. The default is FALSE.

• cex, cex.lab, cex.axis, cex.main: the character expansions relative to the standard size of the device to be used for printing text, labels, axes and main title. See par for details.

• lwd1, lwd2: the line widths used to compare different curves in the same plot; default is lwd2 = 2 for higher order solutions and lwd1 = 1 for first order solutions.

• lty1, lty2: line types used to compare different curves in the same plot; default is "solid" for all statistics.

• cl1, cl2, col: colors used to compare different curves in the same plot; default is cl2 = "red" for higher order solutions, and cl1 = "blue" for the remaining first order statistics. The default color of the plot is col = "black".

• ylim: a numerical vector with two elements defining the $y$-axis range; only if x is a profile.nlreg object.

• ...: additional graphics parameters.

Returns

No value is returned.

Side Effects

A plot is produced on the current graphics device.

Details

The function defaults to:

plot.nlreg.profile(x = stop("nothing to plot"), alpha = 0.05, add.leg = FALSE, 
              stats = c("sk", "fr"), cex = 0.7, cex.lab = 1, cex.axis = 1, 
              cex.main = 1, lwd1 = 1, lwd2 = 2, lty1 = "solid", lty2 = "solid", 
              cl1 = "blue", cl2 = "red", col = "black", ylim = c(-3,3), \dots)

plot.all.nlreg.profiles(x = stop("nothing to plot"), nframe, alpha = 0.05, 
                   stats = c("sk", "fr"), cex = 0.7, cex.lab = 1, cex.axis = 1, 
                   cex.main = 1, lwd1 = 1, lwd2 = 2, lty1 = "solid", lty2 = "solid", 
                   cl1 = "blue", cl2 = "red", col = "black", ylim = c(-3,3), \dots)

References

Fraser, D.A.S., Reid, N. and Wu, J. (1999). A simple general formula for tail probabilities for frequentist and Bayesian inference. Biometrika, 86 , 249--264.

Skovgaard, I. (1996) An explicit large-deviation approximation to one-parameter tests. Bernoulli, 2 , 145--165.

See Also

profile.nlreg, nlreg.profile.objects, plot

Examples

## Not run:

data(metsulfuron)
metsulfuron.nl <- 
    nlreg( formula = log(area) ~ log( b1+(b2-b1) / (1+(dose/b4)^b3) ), 
           weights = ~ ( 1+dose^exp(g) )^2, data = metsulfuron, 
           start = c(b1 = 138, b2 = 2470, b3 = 2, b4 = 0.07, g = log(0.3)),
           hoa = TRUE)
##
metsulfuron.prof <- profile( metsulfuron.nl, offset = g, trace = TRUE )
plot( metsulfuron.prof, lwd2 = 2 )
##
metsulfuron.prof <- profile( metsulfuron.nl, trace = TRUE )
plot( metsulfuron.prof, lwd2 = 2, nframe = c(2,3) )
## End(Not run)