plot(som_model, type="codes", bgcol = my_palette[som_cluster])
som_cluster <- cutree(hclust(dist(som_model$codes)), 3)
add.cluster.boundaries(som_model, som_cluster)

