View source for DScience:Histograms

{{sidebar box|[[DScience:Start|Table of contents]]}}


{{Progression|70}}



__TOC__

A [[Histogram|histogram]] is a summary graph that shows counts of data points falling into various ranges, thus it gives an approximation of the frequency distribution of data. It is an elegant tool to project multidimensional data to lower dimensions and display such projections for visual inspection. The histogram shows data in the form of a bar graph in which the bar heights display event frequencies. Events are measured on the horizontal axis, "X", which has to be binned. The larger number of bins, the higher chances that a fine structure of data can be resolved. Obviously, the binning destroys the fine-grain information of original data.

In this respect, the histogram representation is useful if one needs to create a statistical snapshot of a large data sample in a compact form. Let us illustrate this: Assume we have <math>N</math> numbers, each representing a single measurement. One can store such data in the form of Java or Python lists. Thus, one needs to store <math>8\times N</math> bytes (assuming 8 bytes to keep one number). In the case of a large number measurements, we need to be ready to store a very big output file, as the size of this file is proportional to the number of events. Instead, one can keep for future use only the most important statistical summary of data, such as the shape of the frequency distribution and the total numbers of events. The information which needs to be stored is proportional to the number of bins, thus the file storage has nothing to do with the size of the original data.
An example of a histogram is shown in

<jput>
image 84227025
</jput>

Each bin of the histogram has associated statistical uncertainties showed with vertical lines. 
The size of the error bars is equal to the square root of the number of  entries.


= 1D histogram =

To create a histogram in one dimension (1D), one needs to define the number of bins, <span>Nbins</span>, and the minimum (<span>Min</span>) and the maximum (<span>Max</span>) value for a certain variable. The bin width is given by <span>(Max-Min)/Nbins</span>. If the bins are too wide (<span>Nbins</span> is small), important information might get lost. On the other hand, if the bins are too narrow (<span>Nbins</span> is large), what may appear to be meaningful information really may be due to statistical variations of data entries in bins. To determine whether the bin width is set to an appropriate size, different bin sizes should be tried.

DMelt histograms are designed on the bases of the JAIDA FreeHEP library . For one-dimensional histograms, use the class <javadoc>jhplot.H1D</javadoc>. To initialize an empty histogram, the following constructor can be used:

<pycode>
jhplot import *
h1=H1D("data", 100, 0, 20)
</pycode>

The code above  creates a 1D histogram with the title <span>“data”</span>, the number of bins <span>Nbins=100</span>, and the range of axis <math>X</math> to be binned, which is defined by the minimum and the maximum values, <span>Min=0</span> and <span>Max=20</span>. Thus, the bin width of this histogram is fixed to 0.2. The bin size and the number of bins are given by the following methods:

<pre>
d=h1.getBinSize() 
i=h1.getBins() </pre>

We should note that a fixed-size binning is used. In the following sections, we will consider a more general case when the histogram bin size is not fixed to a single value.

The method <span>fill(d)</span> fills a histogram with a single value, where <span>“d”</span> is a double number. The histograms can be displayed on the <javadoc>jhplot.HPlot</javadoc> canvas using the standard <span>draw(h1)</span> method.

Let us give a complete example of how to fill a histogram with Gaussian random numbers:
<pycode>
from java.util import Random
r=Random()
for i in range(Max):
   h1.fill(m+sd*r.nextGaussian())
</pycode>

that fills <span>Max</span> Gaussian random numbers with a mean <span>m</span> and a standard deviation <span>sd</span>. To do this, we use the package <span>java.util.Random</span>. 

In the simplest case, that code that makes the figure shown above looks like this:

<jcode>
dmelt  84227025.py
</jcode>

We have used two new methods of the <javadoc>jhplot.H1D</javadoc> class in the script shown above. Due to their importance, we will discuss them here:

; <span>setFill(b)</span>
: fill a histogram area when <math>b=1</math> (boolean “True”). When <math>b=0</math> (boolean “False”), the area is not filled;
; <span>setFillColor(c)</span>
: color (Java AWT class) for filling the histogram area.

The histograms have the following features:

* The height of bins for each histogram depends on the bin size. Even when the number of entries is the same, histograms are difficult to compare in shape when the histogram bins are different, see Fig. [h1d] (left).
* Relative size of errors decreases with increasing the number of entries.

We will show a few basic manipulations useful for examining the shapes of histograms, assuming that the underlying mechanism for occurrence of events reveals itself in shapes of event distributions, rather than in the overall statistics or chosen bin size. The shape of the distributions is very important as it conveys information about the probability distribution of event samples.

First of all, let us get rid of the bin dependence. To do this, we will divide each bin height by the bin size. Assuming that <span>“h1”</span> represents a histogram, this can be done as:

<pre>
width=h1.getBinSize()
h1.scale(1/width) </pre>

After this operation, all histogram entries (including statistical uncertainties) will be divided by the bin width. You may still want to keep a copy of the original histogram using the method <span>h2=h1.copy()</span>, so one can come back to the original histogram if needed.

Different histograms can have different normalization, thus a visual comparison of histograms might be difficult. In the case if we are interested in histogram shapes, one can divide each bin height by the total number of histogram entries. This involves another scaling:

<pre>
entries=h1.allEntries()
h1.scale(1/entries)</pre>
Obviously, both operations above can be done using one line:

<pre>
h1.scale(1/(h1.getBinSize()*h1.allEntries())) </pre>
The second step in comparing our histograms would be to shift the bins of the third histogram. Normally, we do not know the exact shift (what should be don in this case will be considered later). At this moment, for the sake of simplicity, we assume that the shift is known and equals to -2. There is a special histogram operation which does such shift:

<pre>h2.shift(-2) # shift all bins by -2</pre>
Now we are ready to modify all the histograms and to compare them. Look at the example below:

After execution of this script, you will find three overlaid histograms. The shapes of all histograms will be totally consistent with each other, i.e. all bin heights will agree within their statistical uncertainties.

One can find many situations in which you may be interested in how well histogram shapes agree to each other. For example, let us assume that each histogram represents the number of days with rainfall measured during one year for one state. If the distributions are shown as histograms, it is obvious that bigger states have a larger number of days with rainfalls compare to small states. This means that all histograms are shifted (roughly by a factor proportional to the area of states, ignoring other geographical differences). The measurements could be done by different weather stations and the bin widths could be rather different, assuming that there is no agreement between the weather stations about how the histograms should be defined. Moreover, the measurements could be done during different time intervals, therefore, the histograms could have rather different numbers of entries. How one can compare the results from different weather stations, if we are only interested in some regularities in the rainfall distributions? The answer to this question is in the example above: all histograms have to be: 1) normalized; 2) shifted; 3) a bin dependence should be removed.

The only unclear question is how to find the horizontal shifts, since the normalization issue is rather obvious and can be done with the method discussed above. This problem will be addressed in the following chapters when we will discuss a statistical test that evaluates the “fit” of a hypothesis to a sample.

= Probability distribution and probability density =

The examples above tell that there are several quantities which can be derived from a histogram. One can extract a <span>''probability distribution''</span> by dividing histogram entries by the total number of entries. The second important quantity is a <span>''probability density''</span>, when the probability distribution is divided by the bin width, so that the total area of the rectangles sums to one (which is, in fact, the definition of the probability density).

Both the probability distribution and the density distribution can be obtained after dividing histogram entries as discussed above. However, these two characteristics can be obtained easier by calling the following methods:

<pre>h2=h1.getProbability()
h2=h1.getDensity() </pre>
which return two new <span>H1D</span> objects: the first represents the probability distribution and the second returns the probability density. In addition to the obvious simplicity, such methods are very useful for variable-bin-size histograms, since this case is taken into account automatically during the division by bin widths. You can check the density distribution using this statement:

<pre>
print h1.integral()</pre>
which prints “1” if the histogram is properly normalized.

Note the following: one can save computation time in the case of calculation of the probability distributions if you know the total number of events (or entries) <math>N_{tot}</math> beforehand. In this case, one can obtain the probability distribution using the weight <math>w1 = 1.0/N_{tot}</math> in the method <span>fill(x,w1)</span>, without subsequent call to the method <span>getProbability()</span>. After the end of the fill, the histogram will represent the probability distribution normalized to unity by definition. In addition, one can remove the bin dependence by specifying another weight as <math>w2 = 1.0/bsize</math>, where <math>bsize</math> is the size of the bin. Finally, the density distribution can be obtained using the weight <math>w3=w1*w2</math>.

= Histogram characteristics =

This section continues our discussion of important characteristics of the <javadoc>jhplot.H1D</javadoc> histogram class. The most popular characteristics of a histogram are the median and the standard deviations (RMS). Assuming that <span>h1</span> represents a <span>H1D</span> histogram, both (double) values can be obtained as:

<pre>
d=h1.mean()
d=h1.rms() </pre>
We already know that one can obtain the number of entries with the method <span>allEntries()</span>. However, some values could fall outside of the selected range during the <span>fill()</span> method. Luckily, the histogram class has the following list of methods for to access the number of entries:

<pre>
gt; i=h1.allEntries()   # all entries 
i=h1.entries()      # number entries in the range
i=h1.extraEntries() # under and overflow entries 
i=h1.getUnderflow() # underflow entries 
i=h1.getOverflow()  # overflow entries </pre>
All methods above return integer numbers.

Another useful characteristics is the histogram entropy. It is defined to be the negation of the sum of the products of the probability associated with each bin with the base-2 log of the probability. One can get the value of the entropy with the method:

<pre>&gt;&gt;&gt; print &quot;Entropy=&quot;,h1.getEntropy() </pre>

= Initialization and filling methods =

Previously, it has been shown how to initialize a histogram with fixed bin sizes. One can also create a histogram using a simpler constructor, and then using a sequence of methods to set histogram characteristics:

<pre>
h1=H1D(&quot;Title&quot;)
h1.setMin(min); h1.setMax(max); h1.setBins(bins) </pre>
which are used to set the minimum, maximum and the number of bins. These methods can also be useful to redefine these histogram characteristics after the histogram was created using the usual approach.

One can also build a variable bin-size histogram by passing a list with the bin edges as shown in this example:

<pre>
bins=[0,10,100,1000]
h1=H1D(&quot;Title&quot;,bins) </pre>
This creates a histogram with three bins with the bin edges given by the input list. This constructor is handy when a data is represented by a bell-shaped or falling distribution; in this case it is important to increase the bin size in certain regions (tails) in order to reduce statistical fluctuations.

As we already know, to fill a histogram with numbers, use the method <span>fill(d)</span>. More generally, one can assign a weight <span>“w”</span> to each value as

<pre>h1.fill(d, w) </pre>
where <span>“w”</span> is any arbitrary number representing a weight for a value <span>“d”</span>. The original method <span>fill(d)</span> assumes that all weights are 1.

But why do we need these weights? We have already discussed in Sect. [h1d:prob] that the weights are useful to reduce the computational time when the expected final answer should be either a probability distribution or density distribution. There are other cases when the weights are useful: We should note again that a histogram object stores the sum of all weights in each bin. This sum runs over the number of entries in a bin only when the weights are set to 1. Events may have smaller weights if they are relatively unimportant compared to other events. It is up to you to make this decision since this depends on a concrete situation.

The method <span>fill(d)</span> is slow in Jython when it used inside loops, therefore, it is more efficient to fill a histogram at once using the method <span>fill(list)</span>, where <span>list</span> is an array with numbers passed from another program or file. As before, <span>fill(list, wlist)</span> can be used to fill a histogram from two lists. Each number in <span>list</span> has an appropriate weight given by the second argument.

Instead of Jython (or Java) lists, one can pass a <javadoc>jhplot.P0D</javadoc> Java array:

<pre>h1.fill(p0d)</pre>
where <span>p0d</span> represents the <span>P0D</span> class.

Analogously, one can fill a histogram by passing a <javadoc>jhplot.PND</javadoc> multidimensional array. This can be done again with the method <span>fill(pnd)</span>, where <span>pnd</span> is an array with any size or dimension. One can specify also weights in the form of an additional <span>PND</span> object passed as a second argument to the method <span>fill(pnd, w)</span>.

Histograms can be filled with weights which are inversely proportional to the bin size - as it was shown in the previous section, removing the bin size dependence is one the most common operations:

<pre>h1.fillInvBinSizeWeight(d)</pre>
It should be noted that this method works even when histograms have irregular binning.

Finally, one can set the bin contents (bin heights and their errors) from an external source as shown below:

<pre>
h1.setContents(values, errors)
h1.setMeanAndRms(mean,rms)</pre>
where <span>values</span> and <span>errors</span> are input arrays. Together with the settings for the bin content, the second line of the above example shows how to set the global histogram characteristics, such as the mean and the standard deviation. There are more methods dealing with external arrays; advanced users can find appropriate methods in the API documentation of the class <span>H1D</span> or using the code assist.

One can create a histogram object from a function. Here is a simple example that shows how to do this:

= Accessing histogram values =

One-dimensional histograms based on the <javadoc>jhplot.H1D</javadoc> class can easily be viewed using the following convenient methods designed for visual inspection:

; <span>toString()</span>
: - convert a <span>H1D</span> histogram into a string
; <span>print()</span>
: - print a histogram
; <span>toTable()</span>
: - show a histogram as a table

= Integration =

Histogram integration is similar to the <span>F1D</span> functions considered in the previous chapters: We simply sum up all bin heights. This can be done using the method <span>integral()</span>. More often, however, it is necessary to sum up heights in a certain bin region, say between a bin <span>“i1”</span> and <span>“i2”</span>. Then use this method:

<pre>sum=h1.integral(i1,i2) </pre>
We should note that the integral is not just the number of events between these two bins: the summation is performed using the bin heights. However, if the weights for the method <span>fill()</span> are all set to one, then the integral is equivalent to the summation of numbers of events.

The integration shown above does include multiplication by a bin width. If one needs to calculate an integral over all bins by multiplying the bin content by the bin width (which can be either fixed or variable), use the method:

<pre>
sum=h1.integral(i1,i2,1)</pre>
where the last parameter should be set to 1 (or to “true” in case of Java codding).

The next question is how to integrate a region in <math>X</math> by translating <math>X</math>-coordinates into the bin indexes. This can be done by calling the method <span>findBin(x)</span>, which returns an index of the bin corresponding to a coordinate <math>X</math>. One can call this method every time when you need to identify the bin index before calling the method <span>integrate()</span>. Alternatively, this can be done in one line as:

<pre>
sum=integralRegion(xmin,xmax,b)</pre>
The method returns a value of the integral between two coordinates, <span>xmin</span> and <span>xmax</span>. The bin content will be multiplied by the bin width if the boolean value <span>b</span> is set to 1 (boolean “true” in Java).

= Histogram operations =

Histograms can be added, subtracted, multiplied and divided. Assuming that we have filled two histograms, <span>h1</span> and <span>h2</span>, all operations can be done using the following generic method:

<pre>h1.oper(h2,&quot;NewTitle&quot;,&quot;operation&quot;)</pre>
where <span>“operation”</span> is a string which takes the following values: “+” (add), “-” (subtract), “*” (multiply) and “/” (divide). The operations are applied to the histogram <span>h1</span> using the histogram <span>h2</span> as an input. One can skip the string with a new title if one has to keep the same title as for the original histogram. In this case, the additive operation will look as <span>h1.oper(h1,“+”)</span>

To create an exact copy of a histogram, use the method <span>copy()</span>. Previously, we have already discussed the <span>scale(d)</span> and <span>shift(d)</span> operations.

A histogram can be smoothed. This topic will be described in Sect. [sec_smoothing], since smoothing and interpolation are widely used data-analysis techniques. Here we just mention that a smoothing of histogram can be done using this method:

<pre>h1=h1.operSmooth(b,k)</pre>
This is done by averaging over a moving window. If <span>“b=1”</span> then the bins will be weighted using a triangular weighting scheme favoring bins near the central bin (<span>“b=0”</span> for the ordinary smoothing) One should specify the integer parameter <span>“k”</span> which defines the window as “2*k + 1”. The smoothing may be weighted to favor the central value using a “triangular” weighting. For instance, for “k=2”, the central bin would have weight 1/3, the adjacent bins 2/9, and the next adjacent bins 1/9. For all these operations, errors are kept to be the same as for the original (non-smoothed) histogram.

One can also create a Gaussian smoothed version of a <span>H1D</span> histogram. Each band of the histogram is smoothed by a discrete convolution with a kernel approximating a Gaussian impulse response with the specified standard deviation.

<pre>h2=h1.operSmoothGauss(rms)</pre>
where <span>rms</span> is a double value representing a standard deviation of the Gaussian smoothing kernel (must be non-negative).

One useful technique is histogram re-binning, i.e. when groups of bins are joined together. This approach could be used if statistics in bins is low; in this case, it makes sense to make bins larger in order to reduce relative statistical uncertainty for entries inside bins (we remind that in case of counting experiments, such uncertainty is <math>\sqrt{N}</math>, where <math>N</math> is a number of entries). The method which implements this operation is called <span>rebin(group)</span>, where <span>group</span> defines how many bins should merged together. This method returns a new histogram with a smaller number of bins. However, there is one restriction: the method <span>rebin</span> cannot be used for histograms with non-constant bin sizes.

= Accessing low-level Jaida classes =

The <span>H1D</span> class is based on the two classes, <span>IAxis</span> and <span>Histogram1D</span> of the Jaida FreeHep library. Assuming <span>h1</span> represents a <span>H1D</span> object, these two Jaida classes can be obtained as:

<pre>
a=h1.getAxis() # get IAxis object
h=h1.get()     # get Histogram1D class </pre>
Both objects are rather useful. Although they do not contain graphical attributes, they have many methods for histogram manipulations, which are not present for the higher-level <span>H1D</span> class. The description of these Jaida classes is beyond the scope of this book. Please look at the Java documentation of these classes or use the code assist.

= Graphical attributes =

Histograms can be shown either by lines (default) or by using symbols. For the default option (lines), one can consider either to fill histogram area or keep this area empty. The following methods below can be useful:

<pre>h1.setFill(b) 
h1.setFillColor(c)</pre>
For the first method, Jython boolean <span>“b=1”</span> means to fill the histogram, while <span>“b=0”</span> (false) keeps the histogram empty. If the histogram area has to be filled, you may consider to select an appropriate color by specifying Java AWT Color object <span>“c”</span>. How to find appropriate color has been discussed in Sect. [hplot_func].

Histograms can be shown using symbols as:

<pre>h1.setStyle(&quot;p&quot;) </pre>

The style can be set back to the default style (histogram bars). This can be done by passing the string <span>“h”</span> instead of <span>“p”</span>. One can also use symbols connected by lines; in this case, use the character <span>“l”</span> (draw lines) or the string <span>“lp”</span> (draw lines and symbols).