1. Try the scott-factor, but the result still differ from the output of
   scipy.kde.gaussian_kde().

   I can't figure it out yet.