def my_kde_bandwidth(obj, fac=1./5):
    return np.power(obj.n, -1./(obj.d+4)) * fac
We use Scott's Rule, multiplied by a constant factor.
dataset/ETHPy150Open scipy/scipy/doc/source/tutorial/stats/plots/kde_plot2.py/my_kde_bandwidth