我使用nlmeR包和包含的函数构建了几个广义非线性最小二乘模型（指数衰减）。我不简单地使用基函数构建非线性最小二乘模型的原因是因为我希望能够对异方差进行建模以避免变换。我的模型看起来像这样：gnls()nls()
model <- gnls(Response ~ C * exp(k * Explanatory1) + A,
start = list(C = c(C1,C1), k = c(k1,k1), A = c(A1,A1)),
params = list(C ~ Explanatory2, k ~ Explanatory2,
A ~ Explanatory2),
weights = varPower(),
data = Data)
与简单模型的主要区别nls()在于weights参数，它可以通过解释变量对异方差性进行建模。的线性等效项gnls()是广义最小二乘法，使用nlmegls()函数运行运行。
现在我想计算置信区间R并将它们与我的模型拟合ggplot()（ggplot2包）一起绘制。我对一个对象执行此操作的方法gls()是这样的：
NewData <- data.frame(Explanatory1 = c(...), Explanatory2 = c(...))
NewData$fit <- predict(model, newdata = NewData)
到目前为止，一切正常，我的模型适合我。
modmat …