import numpy as np
import holoviews as hv

hv.extension('matplotlib')

normal = np.random.randn(1000, 2)
hv.Bivariate(normal)

hv.Bivariate(normal).opts(colorbar=True, cmap='Blues', filled=True)

hv.NdLayout({bw: hv.Bivariate(normal).opts(axiswise=True, bandwidth=bw)
             for bw in [0.05, 0.1, 0.5, 1]}, 'Bandwidth')

from holoviews.operation.stats import bivariate_kde
dist = hv.Bivariate(normal)
kde = bivariate_kde(dist, x_range=(-4, 4), y_range=(-4, 4), bw_method='silverman', n_samples=20)
kde