Something I really like about R is that the plot() command lets you quickly see a lot of what is going on in a data frame
e.g.
library(datasets)
iris_df<-data.frame(iris)
plot(iris_df)
Whilst this set of graphs is nice it doesn't deal well with categorical variables. Column 5 of our data is species, and whilst the last row of the graphs tells us how species relates to the other variables it would be nice to see how the clustering in the other graphs relates to species.
You can do this via
plot(iris_df[is.element(iris_df$Species, "versicolor"),])
plot(iris_df[is.element(iris_df$Species, "setosa"),])
plot(iris_df[is.element(iris_df$Species, "virginica"),])
but this gives you three separate plots. I'd like to combine them and visualise them with different colours so you can see how the clustering works in one image.
(I imagine this might be equivalent to using hold on in matlab )

