<<data,echo=TRUE,cache=TRUE,message=FALSE,tidy=TRUE,warning=FALSE>>=
library(RMySQL)
library(ggplot2)
library(plyr)
library(xtable)
drv <- dbDriver("MySQL")
co <- dbConnect(drv, user="2009Expo", password="R R0cks", port=3306, dbname="accidents", host="headnode.stat.iastate.edu")

dbListTables(co)
dbListFields(co, "accidents")

q = dbGetQuery(co, "select Year, DAY_WEEK, avg(DRUNK_DR>0), stddev(DRUNK_DR>0),
sum(DRUNK_DR>0) from accidents group by Year, DAY_WEEK")
names(q) = c("year", "day", "drunk_avg", "drunk_stdev", "total")
q$day = as.factor(q$day)
levels(q$day) = c("sun", "mon", "tues", "wed", "thurs", "fri", "sat", "NA")

@