R plot.profile.nls 绘制 profile.nls 对象

说明

显示非线性回归模型中参数的轮廓 t 函数和插值置信区间的一系列图，该模型已使用 nls 拟合并使用 profile.nls 进行分析。

用法

## S3 method for class 'profile.nls'
plot(x, levels, conf = c(99, 95, 90, 80, 50)/100,
     absVal = TRUE, ylab = NULL, lty = 2, ...)

参数

x	类 "profile.nls" 的对象
levels	级别，以 t 统计量的绝对值为尺度，在该级别上插值间隔。通常使用 conf 而不是显式给出 levels。
conf	参数的基于配置文件的置信区间的置信水平的数值向量。默认为c(0.99, 0.95, 0.90, 0.80, 0.50).
absVal	一个逻辑值，指示绘图是否应采用轮廓 t 的绝对值的比例。默认为 TRUE 。
lty	用于轴和删除线的线类型。
ylab, ...	plot.default 函数的其他参数可以在此处传递(但不能传递 xlab 、 xlim 、 ylim 或 type )。

细节

这些图以一组 hard-coded 颜色生成，但由于这些颜色是按数字编码的，因此可以通过设置 palette 来更改其效果。颜色 1 用于轴，颜色 4 用于轮廓本身。颜色 3 和 6 用于零处的轴线以及落在轴上的水平/垂直线。

例子

require(graphics)

# obtain the fitted object
fm1 <- nls(demand ~ SSasympOrig(Time, A, lrc), data = BOD)
# get the profile for the fitted model
pr1 <- profile(fm1, alphamax = 0.05)
opar <- par(mfrow = c(2,2), oma = c(1.1, 0, 1.1, 0), las = 1)
plot(pr1, conf = c(95, 90, 80, 50)/100)
plot(pr1, conf = c(95, 90, 80, 50)/100, absVal = FALSE)
mtext("Confidence intervals based on the profile sum of squares",
      side = 3, outer = TRUE)
mtext("BOD data - confidence levels of 50%, 80%, 90% and 95%",
      side = 1, outer = TRUE)
par(opar)

作者

Douglas M. Bates and Saikat DebRoy

参考

Bates, D.M. and Watts, D.G. (1988), Nonlinear Regression Analysis and Its Applications, Wiley (chapter 6)

也可以看看

nls , profile , profile.nls