r - 如何在R中平滑非线性回归曲线

Question

所以我被要求获得 MASS 包中变量 Length 的估计值。我使用的代码如下所示，以及生成的曲线。不知何故，我最终没有得到一条平滑的曲线，而是一条非常“块状”的曲线，以及曲线上各点之间的一些线。谁能帮我得到一个平滑的曲线？

utils::data(muscle,package = "MASS")
Length.fit<-nls(Length~t1+t2*exp(-Conc/t3),muscle,
                start=list(t1=3,t2=-3,t3=1))
plot(Length~Conc,data=muscle)
lines(muscle$Conc, predict(Length.fit))

情节图像： 。

编辑：作为后续问题：如果我想更准确地预测曲线，我使用非线性回归来预测 21 个物种中的每一个的曲线。这给了我一个向量

theta=(T11,T12,...,T21,T22,...,T3).

我可以创建一个绘制所有图形的 for 循环，但像以前一样，我最终得到了块状曲线。但是，看到我必须将这些曲线绘制如下：

for(i in 1:21) {
  lines(muscle$Conc,theta[i]+theta[i+21]*
        exp(-muscle$Conc/theta[43]), col=color[i])
  i = i+1
}

我不知道如何使用相同的技巧来平滑这些曲线，因为肌肉$Conc 仍然只有 4 个值。

编辑2：我想通了，并将其更改为以下内容：

lines(seq(0,4,0.1),theta[i]+theta[i+21]*exp(-seq(0,4,0.1)/theta[43]), col=color[i])

Answer 1

如果您查看 的输出cbind(muscle$Conc, predict(Length.fit))，您会看到许多点是重复的，并且它们没有按 的顺序排序Conc。lines只需按顺序绘制点并连接点，为您提供多条来回线。下面的代码predict在一组唯一的有序值上运行Conc。

plot(Length ~ Conc,data=muscle)
lines(seq(0,4,0.1),
      predict(Length.fit, newdata=data.frame(Conc=seq(0,4,0.1))))