Discussion :

[iwgaqm]

Bonjour à tous,

Je travaille actuellement sur les arbres décisionnels et j'essaie de représenter graphiquement mon arbre C5.0.

Mon arbre semble plutôt bien fonctionner :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
cnt_c50 <- C5.0Control(CF=0.25,minCases=30,sample=0.7)
myTree <- C5.0(Y ~ X1+X2+X3+....+X30, data=data, control= cnt_c50,trials=100)
summary(myTree)

ne renvoie pas d'erreur.
Par contre, quand j'essaie de le visualiser :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

plot(myTree)

J'ai ce message d'erreur :

Error in partysplit(varid = as.integer(i), index = index, info = k, prob = NULL) : minimum of ‘index’ is not equal to 1
In addition: Warning message:
In min(index, na.rm = TRUE) :
no non-missing arguments to min; returning Inf

J'ai aussi essayé avec :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
library(rattle)
fancyRpartplot(myTree)

qui renvoit cette erreur :

Error in if (model$method == "class") { : argument is of length zero

J'ai lu que cela pouvait être lié à des catégories comprenant trop peu d'individus donc j'ai essayé en supprimant toutes les variables qui pourraient être concernées et j'ai les mêmes erreurs.

J'ai aussi essayé cette fonction : http://r-project-thanos.blogspot.fr/...rees-in-r.html mais je suis vraiment une bille en programmation donc j'ai des erreurs mais je ne saurai pas dire d'où elles viennent.

Quelqu'un a-t-il déjà eu ce problème ?

Par avance merci beaucoup !

[olivier.decourt]

Bonjour.
Pour l'utilisation de la fonction rattle::fancyRpartplot, le message est dû au fait qu'elle attend un objet produit par la fonction rpart (arbres CART) et pas C5.0.
Concernant ton C50::plot, sans certitude, est-ce que l'erreur ne viendrait pas du boosting que tu demandes avec trials=100 ? Si tu produis 100 arbres, lequel veux-tu afficher ?

[iwgaqm]

Bonjour Olivier,

Merci pour ta réponse. J'avais lu qu'on pouvait utiliser fancyRpartplot sur Stackoverflow en réponse à quelqu'un qui avait le même problème que moi mais en effet il est bien plus logique d'utiliser cette fonction avec un arbre CART !

Pour le boosting, en effet si je note trials=1 et non trials=100, je n'obtiens pas le même message d'erreur ! Mais toujours un message d'erreur

Error in if (!n.cat[i]) { : missing value where TRUE/FALSE needed

Et si je mets trials=2 :

Error in xval[j, 1] : subscript out of bounds
In addition: Warning message:
In 1:which(out == "Decision tree:") :
numerical expression has 2 elements: only the first used

C'est à n'y rien comprendre

Je pense d'ailleurs qu'il doit être possible de choisir quel arbre des X réalisés on veut afficher, mais je ne connais pas la manip...

Edit : on peut le faire en ajoutant l'argument trial

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

plot(myTree,trial=9)

Et le message d'erreur change encore

Error in !all.equal(diff(sort(unique(index))), rep(1, max(index, na.rm = TRUE) - :
invalid argument type

Je me remets donc à mes recherches...