0

我正在尝试从已发表的论文和上传的数据集中处理一些示例代码。数据和代码来自包 drcSeedGerm

我正在尝试运行此示例代码：

library(devtools)
install_github("OnofriAndreaPG/drcSeedGerm")

library(drc)
library(drcSeedGerm)
library(lmtest)
library(sandwich)
data(rape)
head(rape)

运行没有问题。下一个例子是这样的：

modHTE<-drm(nSeeds~timeBef+timeAf+Psi,data=rape,fct=HTE1(),type="event")
colnames(rape)

我收到错误消息：

[.data.frame`(temp, , 3) 中的错误：选择了未定义的列

在论文中，他们说： 我们可以看到数据需要按评估间隔分组：'timeBef'和'timeAf'分别是评分间隔的开始和结束（以天为单位），'nSeeds'是发芽的种子。'propCum' 列包含发芽种子的累积比例，对于时间事件模型来说不是必需的。使用函数 HTE1() 拟合 HTE 模型。不需要起始值，因为模型定义中内置了自启动例程。

我错过了运行drm代码所需的东西吗？我很困惑为什么 head(rape) 或 colnames(rape) 会显示数据中列的名称，但代码无法识别它们？

colnames(rape) [1] "Psi" "Dish" "timeBef" "timeAf" "nSeeds" "nCum" "propCum"

谢谢！

4

0 回答 0