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