I have a data set with N values residing in the range [-1, 1]. Assume that N is quite large (1,000 to 10,000). I would like to compute a histogram for this data set with the usual binning requirements of not too few bins (at risk of hiding points of interest) and not too many bins (at risk of introducing a lot of noise). Currently, I'd like to use the Freedman-Diaconis rule for determine bin widths: $h=2\times\text{IQR}\times N^{-1/3}$ (where $h$ is the bin width).
When I compute this value, I get quite a large bin size (~600-900 bins). This number of bins is difficult to fit within a histogram (my histogram must be rendered on a webpage).
I think the problem is my domain data is unscaled. My reasoning for this is as such:
Imagine you have two datasets with gaussian distributions. However, assume the standard deviation for the first is 1, while the standard deviation for the second is 10. Still, they would have an identical shape. Presumably you'd want the histograms for each of these to look the same, with the only difference being the x-axis interval labels. The problem with this is that each dataset would have a different IQR, and thus a different recommendation for the bin width according to Freedman-Diaconis.
It seems like I need to scale my dataset to get a reasonable bin width (that captures the distribution well but isn't affected by its scale). What strategy would I use to scale the IQR?