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

From: Charles C. Berry <cberry_at_tajo.ucsd.edu>
Date: Fri 21 Jul 2006 - 03:03:27 EST


On Wed, 19 Jul 2006, Ulrik Stervbo wrote:

[much deleted]

>
> 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)
>

I do not doubt that that is how some do it, but a better approach is takes account of errors in the values used to form the histogram.

This is reviewed in

T. Lynn Eudey
Statistical Considerations in DNA Flow Cytometry Statistical Science 1996, Vol. 11, No. 4, 320{334

A classic approach for this purpose is

John Mendelsohn; John Rice
Deconvolution of Microfluorometric Histograms with B Splines Journal of the American Statistical Association. Vol. 77, No. 380, 1982.

A FORTRAN program is mentioned there and the authors say that they will provide listings on request (but that WAS a long time ago).

I don't know if this was ever implemented in S/R, but an email to Professor John Rice (UC Berkeley Dept of Stats) might help.

HTH, Chuck

Charles C. Berry                        (858) 534-2098
                                          Dept of Family/Preventive Medicine
E mailto:cberry@tajo.ucsd.edu	         UC San Diego
http://biostat.ucsd.edu/~cberry/         La Jolla, San Diego 92093-0717

______________________________________________
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Received on Fri Jul 21 03:08:02 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 Fri 21 Jul 2006 - 04:18:53 EST.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.