来自`nls`的R表达式适合吗？

Question

我已经使用nls拟合了参数函数，现在我想打印出该函数的表达式，并将学习到的参数替换回去。例如：

x <- runif(100, 0, 100)
m <- 13 * exp(-0.05 * x^2) + 0.1 + runif(100,0,0.1)
mod <- nls(m ~ a*exp(-b*x^2)+c, start=list(a=10,b=0.1,c=0.1))

我可以像这样提取公式和系数：

formula(mod)
# m ~ a * exp(-b * x^2) + c
coef(mod)
#           a           b           c 
# 13.00029360  0.04975388  0.14457936 

但是我看不到直接替代它们的方法。 我似乎唯一能做的就是再次写出公式：

substitute(m ~ a * exp(-b * x^2) + c, as.list(round(coef(mod), 4)))
# m ~ 13.0003 * exp(-0.0498 * x^2) + 0.1446

我的最终目标是从磁盘上的RDS文件读取一个合适的nls对象，并在组织模式文档中显示其功能表达式。

Answer 1

这是您要找的东西吗？

do.call(substitute, args=list(formula(mod), as.list(round(coef(mod),4))))
# m ~ 13.0097 * exp(-0.0501 * x^2) + 0.1536

之所以起作用，是因为do.call首先对args两个参数求值，然后才使用substitute()将系数代入公式表达式中。 即， do.call()最终求值的表达式如下所示：

as.call(list(substitute, formula(mod), as.list(round(coef(mod),4))))
# .Primitive("substitute")(m ~ a * exp(-b * x^2) + c, list(a = 13.0097, 
# b = 0.0501, c = 0.1536))