Freedman–Diaconis rule

(Redirected from Freedman-Diaconis rule)

In statistics, the Freedman–Diaconis rule can be used to select the width of the bins to be used in a histogram.[1] It is named after David A. Freedman and Persi Diaconis.

For a set of empirical measurements sampled from some probability distribution, the Freedman–Diaconis rule is designed approximately minimize the integral of the squared difference between the histogram (i.e., relative frequency density) and the density of the theoretical probability distribution.

In detail, the Integrated Mean Squared Error (IMSE) is

where is the histogram approximation of on the interval computed with data points sampled from the distribution . denotes the expectation across many independent draws of data points. Under mild conditions, namely that and its first two derivatives are , Freedman and Diaconis show that the integral is minimised by choosing the bin width

A formula which was derived earlier by Scott.[2] Swapping the order of the integration and expectation is justified by Fubini's Theorem. The Freedman–Diaconis rule is derived by assuming that is a Normal distribution, making it an example of a normal reference rule. In this case .[3]

Freedman and Diaconis use the interquartile range to estimate the standard deviation: [4] where is the cumulative distribution function for a normal density. This gives the rule

where is the interquartile range of the data and is the number of observations in the sample . In fact if the normal density is used the factor 2 in front comes out to be ,[4] but 2 is the factor recommended by Freedman and Diaconis.

Other approaches

edit
 
10000 samples from a normal distribution data binned using different rules. The Freedman-Diaconis rule results in 61 bins, the Scott rule 48 and Sturges' rule 15.

With the factor 2 replaced by approximately 2.59, the Freedman–Diaconis rule asymptotically matches Scott's Rule for data sampled from a normal distribution.

Another approach is to use Sturges's rule: use a bin width so that there are about   non-empty bins, however this approach is not recommended when the number of data points is large.[4] For a discussion of the many alternative approaches to bin selection, see Birgé and Rozenholc.[5]

References

edit
  1. ^ Freedman, David; Diaconis, Persi (December 1981). "On the histogram as a density estimator: L2 theory". Probability Theory and Related Fields. 57 (4): 453–476. CiteSeerX 10.1.1.650.2473. doi:10.1007/BF01025868. ISSN 0178-8051. S2CID 14437088.
  2. ^ D.W. Scott (1979). "On optimal and data-based histograms". Biometrika. 66 (3): 605–610. doi:10.1093/biomet/66.3.605. JSTOR 2335182.
  3. ^ Scott, D.W. (2009). "Sturges' rule". WIREs Computational Statistics. 1 (3): 303–306. doi:10.1002/wics.35. S2CID 197483064.
  4. ^ a b c D.W. Scott (2010). "Scott's Rule". Wiley Interdisciplinary Reviews: Computational Statistics. 2 (4). Wiley: 497–502. doi:10.1002/wics.103.
  5. ^ Birgé, L.; Rozenholc, Y. (2006). "How many bins should be put in a regular histogram". ESAIM: Probability and Statistics. 10: 24–45. CiteSeerX 10.1.1.3.220. doi:10.1051/ps:2006001.