Lukasz Stasielowicz wrote me a message to indicate that there is a way to plot different colours according to sex in the "traditional" framework in R: plot(df$population,df$suicides_no,col=df$sex) plot(df$population,df$suicides_no,col=(1+df$sex)) plot(df$population,df$suicides_no,col=c("red","green")[df$sex]) I wasn't aware of these, but my general point still stands since I think in other languages this is trickier.