je cherche à superposer les éboulis sur l'arbre de la CAH en indiquant exactement où les mettre (pour ne pas trop empiéter sur l'arbre).
J'ai essayé en vain :
split.screen()
layout()
par(mfrow)
sublot ne semble pas fonctionner, R me dit qu'il ne connait pas, même avec ??subplot
Dans l'idée, je fais mon arbre sur toutes l'espace de l'image ;
puis par(new=T) et là je découpe la zone graphique pour mettre le graphique des éboulis où je veux et à la taille que je souhaite...
Idéalement dans le coin en haut à droite vu que l'arbre de la CAH à l'air d'être organisé de sorte que cette zone est souvent libre.
Code reproductible :
Code : Tout sélectionner
# librairie ; je ne sais plus lesquelles il faut exactement selon la version de R utilisée...
library("ade4")
library("graphics")
# utilisation du jeu IRIS
# calcul de la matrice des distances
distance=dist(iris[,1:4])
# arbre contient la CAH
arbre=hclust(distance,meth="ward")
# longueur de l'arbe pour réduire le graphique des éboulis
length_arbre=length(arbre$height)
# plot de l'arbre
plot(arbre,cex.main=3,main="", col.main="lightskyblue3",lwd=2,cex.axis=3, yaxt="n")
rect.hclust(arbre,k=4,bor="blue")
# ajout des éboulis
par(new=T)
# découpage du nouveau graphique
layout(mat=matrix(c(4,2,1,3), 2, 2), widths = c(0.6, 0.4), heights = c(0.5,0.5) , respect = TRUE)
# graphe des eboulis
plot(arbre$height[length_arbre:(length_arbre-15)], pch=20,type="b", col=colgéné, ylab="", xlab="", main="", cex.lab=1, lwd=1.5, cex.axis=1, yaxt="n", tcl=-0.5)
abline(v=5, lty=2, lwd=2, col="gray60")
Merci ; quelle prise de tête pour faire quelque chose d'aussi simple...