关于查找drc包中LL.3三参数对数logistic函数代码的技术问询
获取drc包中三参数对数Logistic函数(LL.3)的实现代码
你可以直接在R环境中查看LL.3函数的完整实现,步骤如下:
- 先加载drc包:
library(drc)
- 直接调用
LL.3即可打印出函数的定义:
LL.3
LL.3的核心实现代码如下(来自drc包的内置定义):
LL.3 <- function(fixed = c(NA, NA, NA), names = c("b", "d", "e"), ...) { .drctemplate(fixed, names, function(x, parm) { parm[2] + (0 - parm[2])/(1 + exp(parm[1] * (log(x) - log(parm[3])))) }, ...) }
代码说明:
fixed参数用于固定模型中的某些参数(NA表示该参数需要拟合)names定义了模型的可拟合参数:b是斜率参数,d是响应的上限值,e是ED50(即响应达到中间值时的x值)- 内部的匿名函数是LL.3模型的核心计算逻辑,对应三参数对数Logistic公式:
[
f(x) = d + \frac{0 - d}{1 + e^{b(\log(x) - \log(e))}}
]
这里响应下限被固定为0,因此是三参数模型。
另外,如果你想从GitHub仓库的源码中定位它,LL.3这类模型函数都定义在仓库的R/models.R脚本文件中。
内容的提问来源于stack exchange,提问作者user20378754
相关产品推荐
相关产品推荐

