Re: [R] Fitting a distribution to peaks in histogram

From: Ulrik Stervbo <>
Date: Thu 20 Jul 2006 - 02:56:32 EST

On 7/19/06, hadley wickham <> wrote:
> > I would like to fit a distribution to each of the peaks in a histogram, such
> > as this:
> As a first shot, I'd try fitting a mixture of gamma distributions (say
> 3), plus a constant term for the highest bin. You could do this using
> ML. If the number of peaks is truly unknown, this will be a little
> trickier but still possible and you could use the LRT to chose between
> them.

Can you be a bit more excact? I a biologist and relatively new to R

> > Integrate the area between each two peaks, using the means and widths of the
> > distributions fitted to the two peaks. I will be using the integrate
> > function
> Why do you want to do this?

I am measureing the amount of DNA in cells, and I need to know the percentage of cells in a part of the cell cycle; that the percentage of cells in the first peak, in the second peak and so on. I want to integrate the area between to two cells, because that apparently is how its none (as far as I can tell from the literature)

> >
> > The histogram is based on approximately 15000 events, which makes Mclust and
> > pam (which both delivers the information I need) less useful.
> If you have unbinned data, it would be better (more precise/powerful)
> to use that.

It very probably is better, but mclust had no result after running for at least 2 hours (I terminated the function after two hours), and as I generate 50-100 datasets, such as the one used for the histogram, as week, I would like to find a faster solution
> Regards,
> Hadley



______________________________________________ mailing list
PLEASE do read the posting guide
and provide commented, minimal, self-contained, reproducible code.
Received on Thu Jul 20 03:05:04 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Thu 20 Jul 2006 - 06:16:20 EST.

Mailing list information is available at Please read the posting guide before posting to the list.