Я пытаюсь обработать пример кода из опубликованной статьи и загруженного набора данных.Данные и код взяты из пакета 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» содержат кумулятивные пропорции проросших семян, и это не является обязательным для моделей времени до события.Модель HTE подходит с помощью функции HTE1 ().Никаких начальных значений не требуется, поскольку в определение модели встроена самозапускающаяся процедура.
Я что-то упустил для запуска кода drm?Я запутался, почему head (rape) или colnames (rape) будут показывать мне имена столбцов в данных, но код их не распознает?
colnames (rape) [1]"Psi" "Блюдо" "timeBef" "timeAf" "nSeeds" "nCum" "propCum"
Спасибо!