来自 medrm包 medrc的 str2lang(x) 错误与混合效应 Brain-Cousens hormesis 模型

问题描述

我有一个实验，在 5 次不同的消化持续时间（0 次和 4 次非 0 次）后检查种子发芽情况。有四个重复，四个块中的一个。我想使用 medrm 包中的 medrc 将结果建模为非线性混合效应模型，并比较几个模型的拟合，包括 Brain-Cousens (1989) hormesis 模型，以便选择最合适。这是一个示例数据集：

time <- c(0,12,24,48,96,96)
block <- c(1,2,3,4,1,4)
block <- as.factor(block)
germination <- c(.78,.80,.82,.75,.90,.95,.89,.91,.68,.71,.65,.69,.08,.05,.06,.02,.01,0)
data <- data.frame(block,time,germination)

我的模型

m <- medrm(germination ~ time,data = data,fct = BC.4(),random = b + e ~ 1|block)

该代码非常适合我想比较的其他模型类型：log-logistic (fct = LL.3())、Weibull I (fct = W1.3()) 和 Weibull II (fct = W2.3())，但是，当我使用 Brain-Cousens 模型 (fct = BC.4()) 时，出现以下错误：

Error in str2lang(x) : <text>:1:18: unexpected ','
1: germination~(time,^

这个错误是什么意思，我该怎么办？我查看了一些源代码（可在此处获得：https://rdrr.io/github/DoseResponse/medrc/src/R/medrm.R）（并不是我了解其中 99% 的内容......）并注意到 BC.4 似乎没有被包含在内，这让我想知道 medrc 是否完全支持 Brain-Cousens 模型？如果 medrc 不支持它，那么我使用的最好/最简单的功能是什么？ （仅供参考，我的编程经验有限）

在大多数 fct 在线教程中提供的 BC.4() 示例中，通过将 vinclozolin 参数设置为 medrc 也可以生成相同的错误。

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）