import hvplot.pandas  # noqa

from bokeh.sampledata import sea_surface_temperature as sst

df = sst.sea_surface_temperature
df.tail()

df.hvplot.kde()

import numpy as np
import colorcet as cc

categorical_cmap = [cc.rainbow[int(i)] for i in np.linspace(0, 255, 12)]

df.hvplot.kde(by='time.month', cmap=categorical_cmap, legend='top', height=400)