K-кратная перекрестная проверка для объекта модели класса GLMM S4
У меня есть объект GLMM, подходящий с помощьюglmer
функция вR
и хотите выполнить к-кратную перекрестную проверку. Для простых GLM я использовалCVbinary
функция отDAAG
ПКГ, как показано ниже.
> SimpleGLM <- glm(Res ~ Var1 + Var2, data = Data, family=binomial)
> CVbinary(SimpleGLM, nfolds=10, print.details=TRUE)
Fold: 3 2 4 1 7 10 6 9 5 8
Internal estimate of accuracy = 0.828
Cross-validation estimate of accuracy = 0.827
Однако, когда в модель добавляется случайный член для IndID, ошибка (ниже) получается из класса S4 модели, соответствующейglmer
.
GLMMod <- glmer(Res ~ Var1 + Var2 + (1|IndID), data = Data, family=binomial)
> CVbinary(GLMMod , nfolds=10, print.details=TRUE)
Error in obj$data : $ operator not defined for this S4 class
Я искал в Интернете и не смог найти функцию, аналогичнуюCVbinary
это работает с объектами S4, но я хотел бы проверить здесь дважды, прежде чем кодировать его вручную.
Короче говоря, (при условии, что я правильно интерпретируюR
ошибка) есть ли функция, которая выполняет k-кратную перекрестную проверку на объектах S4?