/* --- CAH ----------------------------------------------------- */

/* données ASCII alimentation 1987 (food87.csv) */
clear
insheet using "food87.csv", delimiter(";")

/* CAH des pays, distance euclidienne entre les pays et distance moyenne entre les groupes */
cluster averagelinkage grcof-crisp, measure(L2) name(cahfood)

/* affichage des nouvelles variables créées, qui donnent une description complète de l'arborescence,
permettant de la construire en partant de la gauche : numéro du "nœud" (facultatif), numéro du pays
agrégé et hauteur de l'agrégation en cours */
list cahfood_id cahfood_ord cahfood_hgt

/* dessin de l'arbre, ou dendogramme */
cluster dendrogram cahfood, labels(pays)

/* partitions par troncature de l'arbre, trois classes par exemple */
cluster generate clus3=groups(3), name(cahfood) 
/* affichage des classes */
by clus3, sort : list pays clus3

cluster generate clus4=groups(4) /* quatre classes */
by clus4, sort : list pays clus4

/* ------------------------------------------------------------- */