[R] Visualization of histograms
(Ted Harding)
ted.harding at wlandres.net
Mon Nov 8 14:34:13 CET 2010
On 08-Nov-10 12:56:26, Etienne Stockhausen wrote:
> Dear r-user,
> I've wrote a small function to visualize the creation of a histogram.
> Therefore I create a vector with random values and let them plot in
> histograms. I better give a small example in R-Code:
>
> i=0
> funVec=rnorm(1000)
> temp=hist(funVec,plot=FALSE)
>
> plot(1,xlim=c(max(temp$breaks),min(temp$breaks)),
> ylim=c(0,max(temp$counts)+0.1*max(temp$counts)),col="white")
> for(i in 1:1000)
> {
> hist(funVec[1:i],add=TRUE,col="white",)
> Sys.sleep(1/25)
> }
>
> I don't find any mistake in the example above and nevertheless there
> something strange about the result. R should plot 1000 histograms
> with time breaks of 1/25 seconds between the single plots. But the
> program visualize not even 100 steps. I don't understand that
> circumstance. If I don't use the add parameter, R visualize all 1000
> steps.
>
> I hope somebody could help me or explain me what ist happening there.
> I m looking forward to any hints or clues.
>
> Greetings and thanks in advance
> Etienne
I tried your code, and observed what you expected to see (not what
you say you saw). It may be a problem with your system. As a check,
I added a couple of lines to your code, as follows:
i=0
funVec=rnorm(1000)
temp=hist(funVec,plot=FALSE)
plot(1,xlim=c(max(temp$breaks),min(temp$breaks)),
ylim=c(0,max(temp$counts)+0.1*max(temp$counts)),col="white")
maxx<-max(temp$breaks) ; maxy<-max(temp$counts) ## upper x & y lims
for(i in 1:1000)
{
hist(funVec[1:i],add=TRUE,col="white",)
points(i*maxx/1000,i*maxy/1000,pch="+") ## plot a "tracker" point
Sys.sleep(1/25)
}
and saw the "tracker" move steadily from start to finish, with
clearly far more than a mere 100 steps, and at each new tracker
point there is also a change in the histogram. So in my case the
results were plotted as you wanted.
My system:
sessionInfo()
R version 2.11.0 (2010-04-22)
i486-pc-linux-gnu
locale:
[1] LC_CTYPE=en_GB.UTF-8 LC_NUMERIC=C
[3] LC_TIME=en_GB.UTF-8 LC_COLLATE=en_GB.UTF-8
[5] LC_MONETARY=C LC_MESSAGES=en_GB.UTF-8
[7] LC_PAPER=en_GB.UTF-8 LC_NAME=C
[9] LC_ADDRESS=C LC_TELEPHONE=C
[11] LC_MEASUREMENT=en_GB.UTF-8 LC_IDENTIFICATION=C
attached base packages:
[1] stats graphics grDevices utils datasets methods base
loaded via a namespace (and not attached):
[1] tools_2.11.0
--------------------------------------------------------------------
E-Mail: (Ted Harding) <ted.harding at wlandres.net>
Fax-to-email: +44 (0)870 094 0861
Date: 08-Nov-10 Time: 13:34:10
------------------------------ XFMail ------------------------------
More information about the R-help
mailing list