Mam dopasowanie obiektu GLMM za pomocąglmer funkcja wR i chcesz wykonać k-krotną walidację krzyżową. Dla prostych GLM użyłemCVbinary funkcja zDAAG pkg, jak widać poniżej.

> 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

Jednak gdy do modelu dodany zostanie losowy termin dla IndID, błąd (poniżej) wynika z klasy S4 dopasowania modelu zglmer.

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

Szukałem online i nie byłem w stanie znaleźć podobnej funkcjiCVbinary to działa z obiektami S4, ale chciałem podwójnie sprawdzić tutaj, zanim zakoduję je ręcznie.

W skrócie (zakładając, że poprawnie interpretujęR błąd) Czy istnieje funkcja, która wykonuje k-krotną walidację krzyżową na obiektach S4?