文章/答案/技术大牛

发布

社区首页 >问答首页 >在R中绘制非线性回归conf区间

问在R中绘制非线性回归conf区间
EN

Stack Overflow用户

提问于 2015-03-05 05:14:59

回答 1查看 3.6K关注 0票数 1

请帮帮忙。我已经成功地将一个相当不错的非线性回归曲线拟合到我的数据中，使用：

model2 <- nls(urine~a*(1-exp(-c*water)), data = eagle, start = list(a=550, c=0.00385))

然后我使用以下命令生成我的图：

av <- seq(0,1200, 0.5)
bv <- predict (model2, data=eagle, list(water=av))
plot (urine~water, data = eagle)
lines (av,bv)

我想画一些虚线表示上下限95%的置信区间。我已经设法为我的线性回归分析编写了代码，但我正在与当前的非线性回归作斗争。以下代码：

a <- predict (model2, newdata=data.frame(water=av), interval="confidence")

对于线性回归，这给出了一个包含fit、lwr和upr三列的数据框架。对于我的非线性回归，它只给了我一个我不能做任何事情的数字列表。有什么建议吗？有没有一种简单的方法来获得非线性回归的上线和置信线？

发布于 2015-03-06 01:24:25

se.fit显然没有运行。找到了很多论坛，讨论了R中非线性回归拟合置信区间的困难。经过一整天的搜索和在互联网上的猛烈抨击，我终于找到了一个非常有用的公式。我不是100%确定它是如何工作的，但它是工作的，所以这里是为任何其他需要这样做的人…

http://www.r-bloggers.com/predictnls-part-1-monte-carlo-simulation-confidence-intervals-for-nls-models/

根据一系列x值生成预测值的统计表后，可以访问输出的第6列和第7列中的lwr和upr CI值数据，并以这种方式绘制。

票数 1

页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持

原文链接：

https://stackoverflow.com/questions/28865097

复制

相似问题

问在R中绘制非线性回归conf区间EN