[R] relative frequency plot
Philipp Pagel
philipp.pagel.lists at t-online.de
Thu Apr 27 20:32:34 CEST 2006
On Thu, Apr 27, 2006 at 10:48:39AM -0700, nlei at sfu.ca wrote:
> Hi All,
>
> I want to use "hist" to get the relative frequency plot. But the range of
> ylab is greater than 1,which I think it should be less than 1 since it
> stands for the probability.
>
> I'm confused. Could you please help me with it?
I was pretty confused by that, too at first. The solution is that
freq=False cause hist to plot the DENSITY rather than frequency. And
density is not necesssarily the same as relative frequency. Excerpt from
?hist:
density: values f^(x[i]), as estimated density values. If
'all(diff(breaks) == 1)', they are the relative frequencies
'counts/n' and in general satisfy sum[i; f^(x[i])
(b[i+1]-b[i])] = 1, where b[i] = 'breaks[i]'.
If you want relative distance try something like this:
myhist = hist(x,breaks=52, plot=F)
myhist$counts = myhist$counts / sum(myhist$counts)
plot(myhist,main=NULL,border=TRUE,xlab="days",xlim=c(0,6),lty=2)
Not exactly clean, though -- we are messing with the myhist object...
cu
Philipp
--
Dr. Philipp Pagel Tel. +49-8161-71 2131
Dept. of Genome Oriented Bioinformatics Fax. +49-8161-71 2186
Technical University of Munich
Science Center Weihenstephan
85350 Freising, Germany
and
Institute for Bioinformatics / MIPS Tel. +49-89-3187 3675
GSF - National Research Center Fax. +49-89-3187 3585
for Environment and Health
Ingolstädter Landstrasse 1
85764 Neuherberg, Germany
http://mips.gsf.de/staff/pagel
More information about the R-help
mailing list