now using Freedman-Diaconis method to calculate histogram bin width