import numpy as np
import holoviews as hv
from holoviews import opts
hv.extension('bokeh')

xs = np.random.normal(size=100)
ys = np.random.normal(size=100) * xs
overlay = hv.Points((xs,ys)) * hv.HLine(ys.mean())
overlay.opts(
    opts.HLine(color='blue', line_width=6), 
    opts.Points(color='#D3D3D3'))