Bonjour,

Je souhaiterais faire une validation croisée avec la méthode leave-one-out mais plutot que la somme des erreurs au carré, je voudrais afficher le coefficient de détermination r2. J'ai un prédictant et des prédicteurs, je construit une fonction qui calcule le r2, je calcule mon modèle avec glm et je fais la cross validation avec cv.glm (les deux dernières fonctions sont du package boot):

Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
cost<-function(x,y) (cor(x,y))^2
yields.reg<-glm(Predictant~., data = datayields, family=gaussian)
yields.loo<-cv.glm(datayields,yields.reg,cost)
print(yields.loo)
Tout fonctionne à part que l'estimation delta en sortie de cv.glm est égale à NA alors que cela marche si je prends la fonction cost par défaut.

Merci d'avance si quelqu'un a déjà utilisé ces fonctions et peut m'aider.