nlsConfRegions {nlstools}R Documentation

Confidence regions

Description

Draws parameter values in the Beale's 95 percent unlinearized confidence region

Usage

nlsConfRegions (nls, length = 1000, exp = 1.5)
## S3 method for class 'nlsConfRegions':
plot (x, bounds = FALSE, ask = FALSE, ...)
## S3 method for class 'nlsConfRegions':
print (x, ...)

Arguments

nls an object of class 'nls'
length number of points to draw in the confidence region
exp expansion factor of the hypercube in which random values of parameters are drawn
x an object of class 'nlsConfRegions'
bounds logical defining whether bounds of the drawing hypercube are plotted
ask if TRUE, draw plot interactively
... further arguments passed to or from other methods

Details

A sample of points in the 95 percent confidence region is computed according to Beale's criterion (Beale, 1960). This region is also named the joint parameter likelihood region (Bates and Watts, 1988). The method used consists in a random sampling of parameters values in a hypercube centered on the least squares estimate and rejecting the parameters values whose residual sum of squares do not verify the Beale criterion. The confidence region is plotted by projection of the sampled points in each plane defined by a couple of parameters. Bounds of the hypercube in which random values of parameters are drawn may be plotted in order to check if the confidence region was totally included in the hypercube defined by default. If not the hypercube should be expanded in order to obtain the full confidence region

Value

nlsConfRegions returns a list of four objects:

cr a data frame containing the sample drawn in the Beale's confidence region
rss a vector containing the residual sums of squares corresponding to cr
rss95 the 95 percent residual sum of squares threshold according to Beale (1960)
bounds lower and upper bounds of the hypercube in which random values of parameters have been drawn

Author(s)

Florent Baty florent.baty@unibas.ch
Marie-Laure Delignette-Muller ml.delignette@vet-lyon.fr

References

Beale EML (1960) Confidence regions in non-linear estimations. Journal of the Royal Statistical Society, 22B, 41-88.

Bates DM and Watts DG (1988) Nonlinear regression analysis and its applications. Wiley, Chichester, UK.

Examples

data(growthcurve4)
nls1 <- nls(gompertzm, growthcurve4,
        list(lag = 10, mumax = 0.1, LOG10N0 = 6, LOG10Nmax = 9))
cr <- nlsConfRegions(nls1, exp = 2)
plot(cr, bounds = TRUE)

[Package nlstools version 0.0-8 Index]