From: Gabor Grothendieck <ggrothendieck_at_gmail.com>

Date: Wed, 20 Jun 2007 01:14:28 -0400

R-help_at_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 Wed 20 Jun 2007 - 05:18:26 GMT

Date: Wed, 20 Jun 2007 01:14:28 -0400

Here it is with a few improvements:

- use as.is=TRUE on read.table to get character columns rather than factor - use xaxt = FALSE rather than axes = FALSE to eliminate axis(2) - use h$mids where h is the output of hist to avoid lots of tedious calcs - eliminate all that attaching and detaching - eliminate strptime

DF <- read.table("input.dat", header = TRUE, as.is = TRUE) DF$Time <- times(paste(DF$Time, "00", sep = ":")) tod <- as.numeric(subset(DF, Trip == "m")$Time) h <- hist(tod, xaxt = "n", main = "Morning Bus Arrival Times",

xlab = "Time", col = "blue");

axis(1, h$mids, sub(":00$", "", times(h$mids)))

Although not used subsequently just in case you need it in computations that you did not show:

DF$Date <- chron(DF$Date)

DF$DateTime <- chron(DF$Date, DF$Time)

On 6/19/07, Alan Jackson <Rhelp_at_ajackson.org> wrote:

> I am wrestling with time and date data. I came up with a way to plot and

*> label a histogram with "time of day" on the x-axis, but it seemed like a
**> lot more work than should be necessary. Is there a better way to do what
**> I am trying to do?
**>
**> require(chron)
**> # read input data
**> data = read.table("input.dat", header=T)
**> # Create date-time and chron objects
**> datetime = strptime(paste(data[,3], data[,2]),format="%m/%d/%y %H:%M")
**> time = times(paste(as.vector(data[,2]), ":00",sep=""))
**> # Put it all into a data frame
**> data = data.frame(time, datetime, data$Trip)
**> names(data) = c("Time","DateTime","Trip")
**> attach(data)
**> # Create time of day array
**> times = as.numeric(chron(times = Time))
**> tod = subset(times, Trip=='m');
**> # Plot base histogram
**> hist(tod, axes=F, main="Morning Bus Arrival Times", xlab="Time", col="blue");
**> axis(2);
**> # where are the tics to be?
**> tics = seq(min(tod), max(tod), (max(tod)-min(tod))/5);
**> # build a labeled x-axis for the plot
**> axis(1, tics, labels=sub(":00$","",as.character(chron(times=tics, out.format="h:m:s"))));
**>
**> # cleanup
**> detach(data)
**>
**> --- Data ---
**> Trip Time Date
**> a 15:55 05/15/07
**> m 5:47 05/16/07
**> a 15:54 05/16/07
**> m 5:47 05/17/07
**> a 15:59 05/17/07
**> m 5:50 05/21/07
**> m 5:50 05/22/07
**> a 16:00 05/22/07
**> m 5:48 05/23/07
**> m 5:50 05/24/07
**> a 16:00 05/24/07
**> m 5:48 05/25/07
**> m 5:48 05/29/07
**> a 15:59 05/29/07
**> m 5:46 05/30/07
**> m 5:45 05/31/07
**> a 16:05 05/31/07
**> m 5:47 06/04/07
**> a 15:53 06/04/07
**> m 5:46 06/05/07
**> m 5:47 06/06/07
**> a 15:53 06/06/07
**> m 5:47 06/07/07
**> a 15:51 06/07/07
**> m 5:45 06/08/07
**> f 15:22 06/08/07
**> m 5:48 06/11/07
**> m 5:46 06/12/07
**> m 5:48 06/13/07
**> m 5:47 06/18/07
**> a 15:53 06/18/07
**> m 5:47 06/19/07
**> a 15:55 06/19/07
**>
**>
**> --
**> -----------------------------------------------------------------------
**> | Alan K. Jackson | To see a World in a Grain of Sand |
**> | alan_at_ajackson.org | And a Heaven in a Wild Flower, |
**> | www.ajackson.org | Hold Infinity in the palm of your hand |
**> | Houston, Texas | And Eternity in an hour. - Blake |
**>
**> ______________________________________________
**> R-help_at_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.
**>
*

R-help_at_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 Wed 20 Jun 2007 - 05:18:26 GMT

Archive maintained by Robert King, hosted by
the discipline of
statistics at the
University of Newcastle,
Australia.

Archive generated by hypermail 2.2.0, at Wed 20 Jun 2007 - 07:32:38 GMT.

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