各位达人,我做了一个曲线拟合,使用的方法是nls,然后不知道怎么求拟合值的95%置信区间,哪位能告诉我方法么?谢谢
或者,类似的曲线拟合还可以使用其他什么方法吗?nls好像不能直接给出拟合后的r2值。
predict函数似乎不行。
f <-function(a,b,x1)
{a*exp(b*x1)}
result <- nls(TSI~f(a,b,S1), data=interval95,start = list(a=197.17, b=-0.063))
plot(interval95$S1,interval95$TSI)
curve(152.94886*exp(-0.04892*x),col=2,add=T)
Ypred <- predict(result, interval = "confidence", newdata= interval95, level = 0.95, type = c("terms"))