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

In hoa: Higher Order Likelihood Inference

Description

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

Usage

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

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)

Value

No value is returned.

Side Effects

A plot is produced on the current graphics device.

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)

hoa documentation built on May 2, 2019, 8:56 a.m.