URL: http://bokeh.pydata.org/en/latest/docs/gallery/iris_splom.html
Most examples work across multiple plotting backends, this example is also available for:
import holoviews as hv
from holoviews import opts
hv.extension('matplotlib')
from bokeh.sampledata.iris import flowers
from holoviews.operation import gridmatrix
ds = hv.Dataset(flowers)
grouped_by_species = ds.groupby('species', container_type=hv.NdOverlay)
grid = gridmatrix(grouped_by_species, diagonal_type=hv.Scatter)
grid.opts(opts.Scatter(bgcolor='#efe8e2', alpha=0.4, s=16))