som_cluster <- cutree(hclust(dist(som_model$codes)), 3)
my_palette <- c("#ffffff", '#dddddd', '#ccccff', '#aaaaaa')
plot(som_model, type="codes", bgcol = my_palette [som_cluster])
add.cluster.boundaries(som_model, som_cluster)

