similar to: histogram of dates

Hello! Up to now I have been using hist() to display the distributions. Howevere, I noteiced strange numbers on y (vertical) axis, if I used probability = T or freq = F option. I thought it is a bug and launched the R-bug system and found some posts on that matter. Brian Ripley responded to one, that one should look at truehist() for that. Ok I can use truehist() if I want to see the ratios

Dear all, hist ( ) plots a horizontal line at y=0 when the respective bin is empty. I can deal with this by modifying the hist object before plotting it (x$density[x$density == 0] <- NA), but I'm sure I've seen a more elegant way. Perhaps this was in truehist (MASS). I have looked but can't find it. Does anyone know? Best wishes Paul -- View this message in context:

hi all : i have the data below and then below that, i call the hist function three times using the Scott method for the widths of the bins. the bin width is different for the three histograms but I would like it to always be 0.05 regfardless of the data set being histogrammed. I'm sure there must be a manual way to do this which is fine with me. i tried breaks=0.05 but it wasn't happy

Hi, This is a very basic question, but I would like to undestand hist(). I thought that the hist( , freq=FALSE) should provide the relative frequencies (probabilities), and so they should sum 1, however: set.seed(2) ah <- hist(rnorm(100), freq=F) sum(ah$intensities) [1] 2 set.seed(2) bh <- hist(rlnorm(100), freq=F) sum(bh$intensities) [1] 0.4999996 I'm getting similar figures with

Hi there, # Consider the following example: A = 19 B = 20 A< B A==B hist(c(1:15,B,50),breaks=c(0,15,A,50),plot=F)$count hist(c(1:15,A,50),breaks=c(0,15,B,50),plot=F)$count # I was expecting the same results with the following values of A and B: A = 19.6019203953960 B = 19.6019204365543 A< B A==B hist(c(1:15,B,50),breaks=c(0,15,A,50),plot=F)$count# wrong

Hi, I am facing the problem that I want to plot a histogram chart set freq to true and overlay with normal or weibull or exponential distribution curve. The sample code is shown as below: >samp<-c(-8.2262,-8.2262,-8.2262,-8.20209,-8.09294,-8.07321,-8.07321, -8.07321,-8.07175,-8.04948,-8.04948,-8.04948,-8.03848,-8.03848,

rdevel from yesterday, also in 1.9.1 (Debian, -3). The man page of 'hist' (from the graphics package, obtained by "?hist") suggests 'truehist' under "see also", but 'truehist' isn't loaded by default (a sensible "of course"), so that "?truehist" fails. One can do help.search(), of course, which points what to do, but leads to my

Dear People, I have a data set of data x from a probability distribution, and I have a function, mydensity, of the pdf of that distribution. I'm asking for help in superimposing the histogram of x and the plot of mydensity. In the function below, I call truehist and curve, but these are plotted in different figures. I'd like them to be plotted on the same figure, and to use common

>selectiris<-sample(1:5, 150, replace=T) > hist(selectiris, plot=TRUE) #freq=TRUE, breaks=5, etc all do the same thing >truehist(selectiris) Is the behavior of "hist" in this sample of code correct. On my system (Windows 2K, SP2) the first and second bars are pushed together, while the 3rd, 4th, and 5th bars are correctly positioned. The function "truehist"

Hi, This might be a simple problem but I don't know how to calculate a random variable density the way panel.histogram does it before it creates the actual density rectangles. The documentation says that it uses the density function but the actual code suggests that the hist.constructor function (which does not seem to be easily accessible). Any suggestion for computing the density values of

Dear All, According to the Sturges rule, the number of classes of a histogram is the closest integer to 1 + logb(n,base=2) where n is the number of observations. The function hist(), by default, uses the Sturges rule. However, the code x <- 1:200 hist(x) produces a histogram with 10 classes and not 9 classes as determined by the Sturges rule. What am I missing? Thanks in advance, Paul

Hi, Is this a bug in truehist()? > library(MASS) > x <- rep(1, 10) > truehist(x) Error in pretty(data, nbins) : invalid 'n' value Thanks, Gad > R.version platform i486-pc-linux-gnu arch i486 os linux-gnu system i486, linux-gnu status major 2 minor 4.1 year 2006 month 12 day 18 svn

Hi, Is this a bug in truehist()? > library(MASS) > x <- rep(1, 10) > truehist(x) Error in pretty(data, nbins) : invalid 'n' value Thanks, Gad > R.version platform i486-pc-linux-gnu arch i486 os linux-gnu system i486, linux-gnu status major 2 minor 4.1 year 2006 month 12 day 18 svn

Hello again, I'm interested in manipulating some date objects, and have been playing with the difftime function. I'm not sure that this is the desired behavior: > pa$days.before.apr30 <- difftime(may01,as.POSIXct(pa$training.max,origin=as.POSIXct("1969-12-31 16:00:00",tz="PDT")),units="days") > units(pa$days.before.apr30) [1] "days" >

To whom it may concern, I have a problem concerning the design of a histogram. How do I change the border widths of the bars of a histogram. The initial command is: hist (punkte,breaks=30, xlab="Punkte", ylab="H?ufigkeit", main="Histogramm", col= heat.colors(30), border= "red") I suspect that it has to do with the "lwd" command but can't

I have problems getting a histogram with relative frequencies on the y-axis. Here is an example data set: > a <- c(4.626, 4.627, 4.627, 4.628, 4.629, 4.629, 4.630, 4.631, 4.632, 4.632) > d = hist(a,freq=F) > d$density [1] 299.9999 100.0000 200.0000 100.0000 100.0000 200.0000 The obtained densities are given by counts/(total n * bin width), with bin width being

Hello, how do you create a histogram with a data frame? year snow.cover 1970 6.5 1971 12.0 1972 14.9 1973 10.0 1974 10.7 1975 7.9 ... mydata=data.frame(year=c(1970,...),snow.cover=c(6.5,...)) hist(mydata) does not work. Many thanks. PR

Is there any neat way to add a curve (frequency function or the like) to a histogram or other plots? I haven't found one yet... Robert

Dear People, I am sure I am missing something obvious as usual, but in the following graph I can't plot ylab. Ignoring unimportant details, I am plotting one instance of truehist() and one instance of curve() on the same graph. Truehist() won't let me pass the ylab argument. It gives me the error Error in plot.default(xlim, c(0, ymax), type = "n", xlab = xlab, ylab =

Hello - This is not a real problem, just an annoyance. Sometimes, but not always, I get a set of strange warnings from hist(). Example follows. ############################ #Produce a histogram of start dates for a set of field measurements. # I didn't reproduce all the dates, because not sure it's relevant, but here's a sample. # Note the invalid date in the fifth element. >