如何解决使用 nls() 时的置信区间
我正在编写代码来模拟细菌的灭活。我能够使用 nls() 获得模型。现在,当我尝试计算和绘制置信区间和预测区间时,我无法得到适当的拟合。较低的 CI 和 PI 绘制错误,我不知道如何纠正此问题。
我的代码:
Xo <- c(0,0.461,0.978,2.087)
Y <- c(0,-2.891,-4.197,-4.489)
dat <- data.frame(Xo,Y,colClass="numeric")
f <- function(Xo,n,k) {-k*(Xo)^n}
m <- nls(Y~f(Xo,k),data=dat,start=c(n=2,k=20),trace=TRUE)
summary(m)
yfitted <- predict(m)
plot(Xo,yfitted)
lines(Xo,yfitted)
plot(dat$Xo,yfitted,xlim=c(0,5),ylim=c(-5,0))
y.conf <- predictNLS(m,interval="confidence",alpha =0.05,nsim=10000)$summary
y.pred <- predictNLS(m,interval="prediction",alpha=0.05,nsim=10000)$summary
matlines(Xo,y.conf[,c("Sim.2.5%","Sim.97.5%")],col="red",lty="dashed")
matlines(Xo,y.pred[,col="blue",lty="solid")
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。