PlotPairs(iris[1:4], g=iris$Species, main = "Anderson's Iris Data -- 3 species",
col=c(DescTools::hred, DescTools::hblue, DescTools::hgreen),
col.smooth=c("black", DescTools::hred, DescTools::hblue, DescTools::hgreen))
Run the code above in your browser using DataLab