k-fache Kreuzvalidierung für GLMM S4-Klassenmodellobjekt
Ich habe ein GLMM-Objekt, das mit dem Befehl "glmer
Funktion inR
und wollen eine k-fache Kreuzvalidierung durchführen. Für einfache GLMs habe ich das verwendetCVbinary
Funktion von derDAAG
pkg wie unten gezeigt.
> 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
Wenn dem Modell jedoch ein zufälliger Ausdruck für IndID hinzugefügt wird, ergibt sich ein Fehler (siehe unten) aus der S4-Klasse eines Modells, zu dem es passtglmer
.
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
Ich habe online gesucht und konnte keine ähnliche Funktion findenCVbinary
das funktioniert mit s4 objekten, wollte aber hier nochmal nachschauen bevor ich es manuell codiere.
Kurz gesagt (vorausgesetzt, ich interpretiere das richtig)R
Fehler) Gibt es eine Funktion, die eine k-fache Kreuzvalidierung für S4-Objekte durchführt?