拟合 1 - exp(x) 赋予第一个值更高的权重

Question

我想将 1 - exp(x) 函数拟合到数据集，但给予第一个值更高的权重。但是，以下代码无法以这种方式工作：

x <-sqrt((0.05)^2+(0.05)^2)*seq(from = 1, to = 20, by = 1)
y <- c(11,20,27,32,35,36,36.5,25,16,9,4,1,7.87e-16,2.07e-15,-9.36e-16,1.61e-15,-3.81e-16,3.92e-16,7.65e-16,-8.26e-16)
temp <- data.frame(cbind(x,y))
we <- 1/(log1p(seq_along(x)))
# fit non-linear model
mod <- nls(y ~ (1 - exp(a + b * x)), data = temp, start = list(a = 0, b = 0), weights = we)

#add fitted curve
lines(temp$x, predict(mod, list(x = temp$x)))

这是输出：

Answer 1

您指定的重量是正确的。您获得的不合适是由于您错误的模型假设。你假设：

y ~ 1 - exp(a + b * x)

请注意，exp() 给出严格的正值，因此 y 不会大于 1。但是，y 数据范围内的值最多为 35。

我的想法并不完美，但它可能会给你一个更好的起点。考虑：

y ~ a * x * exp(b * x * x + c * x)

使用您的数据：

x <- c(0, sqrt((0.05)^2+(0.05)^2)*seq(from = 1, to = 20, by = 1))
y <- c(0, 11,20,27,32,35,36,36.5,25,16,9,4,1,7.87e-16,2.07e-15,-9.36e-16,1.61e-15,-3.81e-16,3.92e-16,7.65e-16,-8.26e-16)
fit <- nls(y ~ a * x * exp(b * x * x + c * x), start = list(a = 30, b= -1, c = -1))
plot(x, y)
lines(x, predict(fit, list(x)))

拟合 1 - exp(x) 赋予第一个值更高的权重

Fiting 1 - exp(x) giving higher weight to the first values

r

nls