Hi! Sorry for bothering you again but I can't seem to get it right.
When i multiply the density with the number of observations it seems to be way to high, The reference curve is drawn at maybe 20 times higher frequency count than it should be. I use the following code where "weights$Weight" is my weights data and "reference" is my reference dataset. # calucate the right breakpoints breakpoints <- seq(min(weights), max(weights), by=binwidth) #scale density dens <- density(reference) dens$y <- dens$y * (length(weights$Weight)) #graph it hist(weights$Weight, freq=TRUE, breaks=breakpoints, main=wfiles[i]) lines(dens) Any ideas are greatly appreciated. /Thomas On Fri, Jun 27, 2008 at 10:54 PM, Daniel Folkinshteyn <[EMAIL PROTECTED]> wrote: > if you want the "frequency" scale rather than density scale, then leave hist > as is (by default it uses the frequency scale), and rescale the density by > multiplying it by the appropriate NOBS. > > on 06/27/2008 01:16 PM Thomas Frööjd said the following: >> >> Hi >> >> Thank you very much for taking time to answer. >> >> The solution of using hist(data) for the main dataset and adding >> lines(density(refdata)) for the reference data seem to work great. I >> forgot to mention one thing however, I need to have frequency on the y >> azis instead of density as now. >> >> I know this is not a "real" histogram but since the audience is not >> very statistically experienced I would prefer to do it this way. >> Anyone have an idea? >> >> Thanks again for your help. >> >> Thomas Fröjd >> >> On Wed, Jun 25, 2008 at 6:16 PM, Daniel Folkinshteyn <[EMAIL PROTECTED]> >> wrote: >>>> >>>> I don't understand this. Why not just get hist() to plot on the >>>> density scale, >>>> thereby making its output commensurate with the output of density()? >>>> The hist() function will plot on the density scale if you ask it to. >>>> Set freq=FALSE >>>> (or prob=TRUE) in the call to hist. >>> >>> ehrm... because i didn't realize that option existed :) that certainly is >>> easier than manually scaling hist output by NOBS - thanks for the tip! >>> >> > > ______________________________________________ R-help@r-project.org 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.
