错误消息“已选择未定义的列”，但标题（数据）显示名称？

Question

我正在尝试研究已发表论文和上载数据集的一些示例代码。数据和代码来自软件包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）会向我显示数据中的列名，但是代码无法识别它们？

  

姓氏（强奸）   [1]“ Psi”“ Dish”“ timeBef”“ timeAf”“ nSeeds”“ nCum”“ propCum”

谢谢！