[R] Comparing histograms?

(Ted Harding) Ted.Harding at manchester.ac.uk
Thu May 13 21:14:54 CEST 2010

On 13-May-10 16:00:51, Jonathan Greenberg wrote:
> Rhelpers:
> I'm curious what the appropriate analysis to use for testing the
> hypothesis that two histograms are statistically different from one
> another?  Thanks!
> --j

That's potentially several questions in one, and each question
might have more than one answer!

If you have two histograms with the same number (k) of bins,
then you could think of it as a two-way (2xk) table, and use
a chi-squared or likelihood-ratio test. But then you may well
have to take account of what reasoning went into the choice
of breaks to define the bins.

If the two histograms have different bumbers of bins (and
possibly different break-points), then that approach would
not work. You could consider fitting a distribution (of the
same type, e.g. Normal) to each, and to them jointly, and
then doing a likelihood ratio test. Basically on the lines of

Hist1: m1.1 out of M in (a0,a1), m1.2 out of M in (a1,a2),
   ... m1.k1 out of M in (a[k1-1],ak1)

Hist2: n11 out of N in (b0,b1), n12 out of N in (b1,b2),
   ... n1.k2 out of N in (b[k2-1],ak2)

A: Fit a normal distribution by maximum likelihood (grouped
data) to the joint Hist1 and Hist2 data (same mean mu and
variance V for each); Watch out that you need to condition
on M = N.

B: Fit a Normal distribution to Hist 1 (mean mu1, variance V1)
to the Hist1 data. Fit one (mean mu2m variance V2) to Hist2.

In fit B you have fitted 4 parameters, in fit A you have fitted 2.
Hence 2*(log likelihood from B - log likelihood from A) will
have a chi-squared distribution with 4-2 = 2 degrees of freedom.

You might want to impose that V1 = V2 (or not).

Similar for distributions of other types that you might consider
fitting to the data.

And so on!

E-Mail: (Ted Harding) <Ted.Harding at manchester.ac.uk>
Fax-to-email: +44 (0)870 094 0861
Date: 13-May-10                                       Time: 20:14:51
------------------------------ XFMail ------------------------------

More information about the R-help mailing list