def my_kde_bandwidth(obj, fac=(1.0 / 5)):
 return (np.power(obj.n, ((-1.0) / (obj.d + 4))) * fac)

def my_kde_bandwidth ( obj , fac = ( 1.0 / 5 ) ) : return ( np . power ( obj . n , ( ( - 1.0 ) / ( obj . d + 4 ) ) ) * fac )
we use scotts rule .