Tengo un objeto GLMM ajustado usando elglmer funcionar enR y quieren realizar la validación cruzada de k-fold. Para GLMs simples he usado elCVbinary función de laDAAG pkg como se ve abajo.

> 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

Sin embargo, cuando se agrega un término aleatorio para IndID al modelo, aparece un error (a continuación) de la clase S4 de un modelo que encaja conglmer.

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

He estado buscando en línea y no he podido encontrar una función similar aCVbinary eso funciona con los objetos S4, pero quería verificar dos veces aquí antes de codificarlo manualmente.

En resumen, (asumiendo que estoy interpretando correctamente elR error) ¿existe una función que realice la validación cruzada de k-fold en objetos S4?