Much better base colors in R 4.0 ;-)
d <- transform(as.data.frame(replicate(2, rnorm(30))), g = gl(3, 10)) plot(V1 ~ V2, data = d, col = g, pch = 19)