Most examples work across multiple plotting backends, this example is also available for:
import holoviews as hv
from holoviews import opts
hv.extension('bokeh')
from bokeh.sampledata.iris import flowers
from holoviews.operation import gridmatrix
iris_ds = hv.Dataset(flowers)
density_grid = gridmatrix(iris_ds, diagonal_type=hv.Distribution, chart_type=hv.Bivariate)
point_grid = gridmatrix(iris_ds, chart_type=hv.Points)
(density_grid * point_grid).opts(
opts.Bivariate(bandwidth=0.5, cmap='Blues'),
opts.Points(size=2, tools=['box_select']))