validação cruzada k-fold para o objeto de modelo de classe GLMM S4
Eu tenho um objeto GLMM apto usando oglmer
função emR
e deseja executar a validação cruzada em k-fold. Para GLMs simples eu usei oCVbinary
função doDAAG
pkg como visto abaixo.
> 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
No entanto, quando um termo aleatório para IndID é adicionado ao modelo, um erro (abaixo) resulta da classe S4 de um modelo ajustado comglmer
.
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
Tenho procurado online e não consegui encontrar uma função semelhante àCVbinary
que trabalha com objetos S4, mas queria verificar aqui antes de codificá-lo manualmente.
Em suma, (supondo que estou interpretando corretamente oR
erro) existe uma função que executa a validação cruzada de k-fold em objetos S4?