onls (version 0.1-2)

plot.onls: Plotting function for 'onls' objects

Description

Plots an orthogonal nonlinear model obtained from onls.

Usage

# S3 method for onls
plot(x, fitted.nls = TRUE, fitted.onls = TRUE, 
                    segments = TRUE,...)

Value

A plot of the onls model.

Arguments

x

an object returned from onls.

fitted.nls

logical. If TRUE, the fit from the normal (vertical) nonlinear model is plotted as a blue line for comparison purposes.

fitted.onls

logical. If TRUE, the fit from the orthogonal nonlinear model is plotted as a red line.

segments

logical. If TRUE, segments connecting \((x_i, y_i)\) and \((x_{0i}, y_{0i})\) are displayed.

...

other parameters to plot.

Author

Andrej-Nikolai Spiess

Examples

## Quadratic model with 10% added noise.
## Omitting the "nls" curve,
## display orthogonal segments.
set.seed(123)
x <- 1:20
y <- 10 + 3*x^2
y <- y + rnorm(20, 0, 50)
DAT <- data.frame(x, y)
mod <- onls(y ~ a + b * x^2, data = DAT, start = list(a = 10, b = 3), extend = c(0.2, 0))
plot(mod, fitted.nls = FALSE)

## Due to different scaling, 
## orthogonality of fit is not evident.
## We need to have equal x/y-scaling for that,
## using asp = 1!
plot(mod, fitted.nls = FALSE, xlim = c(0, 10),
     ylim = c(0, 200), asp = 1)