R nls:拟合数据曲线

时间:2015-11-02 13:34:42

标签: r nls

我无法找到适合我数据的正确曲线。如果比我更了解情况的人有更好的拟合曲线的想法/解决方案,我将非常感激。

数据:目的是从y

预测x
dat <- data.frame(x = c(15,25,50,100,150,200,300,400,500,700,850,1000,1500),
                  y = c(43,45.16,47.41,53.74,59.66,65.19,76.4,86.12,92.97,
                        103.15,106.34,108.21,113) ) 

这就是我走了多远:

model <- nls(x ~ a * exp( (log(2) / b ) * y),
             data = dat, start = list(a = 1, b = 15 ), trace = T)

哪个不合适:

dat$pred <- predict(model, list(y = dat$y))
plot( dat$y, dat$x, type = 'o', lty = 2)
points( dat$y, dat$pred, type = 'o', col = 'red')

fit plot

谢谢,F

1 个答案:

答案 0 :(得分:5)

从5度多项式预测x不是那么简约,但似乎很合适:

fm <- lm(x ~ poly(y, 5), dat)
plot(x ~ y, dat)
lines(fitted(fm) ~ y, dat)

(情节后继续)

screenshot

您还可以考虑drc包的UCRS.5b模型:

library(drc)
fm <- drm(x ~ y, data = dat, fct = UCRS.5b())
plot(fm)

screenshot

注意:最初,我认为你想从x预测y,并在下面写了答案。

立方体看起来很不错:

plot(y ~ x, dat)
fm <- lm(y ~ poly(x, 3), dat)
lines(fitted(fm) ~ x, dat)

(情节后继续)

screenshot

4参数逻辑看起来也很好:

library(drc)
fm <- drm(y ~ x, data = dat, fct = LL.4())
plot(fm)

screenshot