bar_plot.R

qplot(x=year, y=injuries, data=t_per_year, geom="bar", stat="identity")

parallel.R

# Needs doMC, forrach, parallel and iterators
library(doMC)
registerDoMC()

per_type.R

per_type <- ddply(storm2011, .(evtype), function(df) {
  data.frame(evtype=df$evtype[1], injuries=sum(df$INJURIES,na.rm=T), fatalities=sum(df$FATALITIES,na.rm=T), propdmg=sum(df$PROPDMG, na.rm=T), cropdmg=sum(df$CROPDMG, na.rm=T))
})

per_type_year.R

storm_type_year <- ddply(storm, .(EVTYPE), function(df) {
  ddply(df, .(year(bgn_date)), function(dg) {
    data.frame(type=dg$EVTYPE[1],
               year=year(dg$bgn_date[1]),
               injuries=sum(dg$INJURIES, na.rm=T),
               fatalities=sum(dg$FATALITIES, na.rm=T),
               propdmg=sum(dg$PROPDMG, na.rm=T),
               cropdmg=sum(dg$CROPDMG, na.rm=T))
  })
})

t_per_year.R

t_per_year <- ddply(tornado, .(year(bgn_date)), function(df) {
  data.frame(year=year(df$bgn_date[1]), injuries=sum(df$INJURIES, na.rm=T), fatalities=sum(df$FATALITIES, na.rm=T))
})