在 R 中拟合指数衰减模型

【问题标题】：Fit an exponential decay model in R在 R 中拟合指数衰减模型
【发布时间】：2020-05-11 16:53:03
【问题描述】：

我是 R 新手，非常感谢您的帮助

我有一些看起来像这样的数据。

Y 与 X 呈非线性负相关。它似乎可以通过以下形式的公式 y=1+a^x 来近似，其中 a.

如果我想在 R 中拟合这些数据以找到 a，我会使用什么函数？ NLS？

【问题讨论】：

是的，nls 值得使用。类似：nls(y ~ b+ a^x, data = mydata, start = list(b = 1, a = 0.5)) 应该可以工作。

标签： r

【解决方案1】：

下次请提供测试数据。这次我们已经为您完成了。然后我们如图所示使用nls。

set.seed(123)

# generate test data
n <- 35
x <- 1:n
a <- 0.5
y <- 1 + a^x + rnorm(n, 0, .01)

fm <- nls(y ~ 1+a^x, start = list(a = mean((y-1)^(1/x), na.rm = TRUE)))
fm

给予：

Nonlinear regression model
  model: y ~ 1 + a^x
   data: parent.frame()
     a 
0.5025 
 residual sum-of-squares: 0.003031

Number of iterations to convergence: 5 
Achieved convergence tolerance: 1.346e-06

情节

plot(y ~ x)
lines(fitted(fm) ~ x, col = "red")

【讨论】：