How can I reproduce the correlation scatterplot matrix in this link with Python? It was generated in R with ggpairs. Here is the code in R:
p <- ggpairs(iris,
