r - R drm 功率曲线 (DRC.powerCurve) 二阶导数

Question

我正在尝试使用基于 https://www.statforbiology.com/nonlinearregression/usefulequations的 drm 计算 R 中的功率曲线（Freundlich 方程） 。这是代码：

  model <- drm(numSpecies ~ Area, fct = DRC.powerCurve(),
                 data = speciesArea)

DRC.powerCurve() 提供一阶导数（即斜率）的估计值，但是，我想知道是否有可能获得二阶导数？根据 R 文档（http://ftp.uni-bayreuth.de/math/statlib/R/CRAN/doc/packages/drc.pdf），您可以在“固定”语句中指定要估计的参数

LL.4(fixed = c(NA, NA, NA, NA), names = c("b", "c", "d", "e"), ...)

但是，当我尝试在 R 中使用 DRC.PowerCurve 将参数估计值设置为零时，我收到错误“不正确的'固定'参数”。我也不确定哪个初始对应于二阶导数。

我几乎找不到关于 DRC.PowerCurve 的信息，更不用说获得二阶导数了。事实上，DRC.PowerCurve 甚至没有在 R 文档或 getMeanFunction() 中列为潜在函数。