library(Rmixmod) data(birds) print(names(birds)) # les variables print(birds[1:5,]) # un aperçu des données (les cinq premiers oiseaux) train <- birds[1:50,] # l'échantillon d'apprentissage test <- birds[51:69,] # l'échantillon de test learn <- mixmodLearn(data=train[,-1],knownLabels=train[,1]) # pour l'interprétation du modèle retenu, voir ref [16] pages 6, 7. prev <- mixmodPredict(test[,-1],classificationRule=learn['bestResult']) print(prev) cat('genre estimé pour les 19 oiseaux de l\'échantillon de test','\n') print(prev[5]) cat('tableau de confusion genre réel/genre estimé','\n') print(table(test[,1],prev[5]))