R nls:将曲线拟合到数据
R nls: fitting a curve to data
我无法找到适合我的数据的正确曲线。如果比我知识渊博的人有 idea/solution 更好的拟合曲线,我将不胜感激。
数据:目的是根据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')
谢谢,F
从 y 预测 x 为 5 次多项式不是那么简洁但似乎适合:
fm <- lm(x ~ poly(y, 5), dat)
plot(x ~ y, dat)
lines(fitted(fm) ~ y, dat)
(剧情后续)
你也可以考虑drc包的UCRS.5b型号:
library(drc)
fm <- drm(x ~ y, data = dat, fct = UCRS.5b())
plot(fm)
注意: 最初,我假设你想根据 x 预测 y 并将答案写在下面。
一个立方看起来还不错:
plot(y ~ x, dat)
fm <- lm(y ~ poly(x, 3), dat)
lines(fitted(fm) ~ x, dat)
(剧情后续)
4 参数逻辑也看起来不错:
library(drc)
fm <- drm(y ~ x, data = dat, fct = LL.4())
plot(fm)
我无法找到适合我的数据的正确曲线。如果比我知识渊博的人有 idea/solution 更好的拟合曲线,我将不胜感激。
数据:目的是根据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')
谢谢,F
从 y 预测 x 为 5 次多项式不是那么简洁但似乎适合:
fm <- lm(x ~ poly(y, 5), dat)
plot(x ~ y, dat)
lines(fitted(fm) ~ y, dat)
(剧情后续)
你也可以考虑drc包的UCRS.5b型号:
library(drc)
fm <- drm(x ~ y, data = dat, fct = UCRS.5b())
plot(fm)
注意: 最初,我假设你想根据 x 预测 y 并将答案写在下面。
一个立方看起来还不错:
plot(y ~ x, dat)
fm <- lm(y ~ poly(x, 3), dat)
lines(fitted(fm) ~ x, dat)
(剧情后续)
4 参数逻辑也看起来不错:
library(drc)
fm <- drm(y ~ x, data = dat, fct = LL.4())
plot(fm)