{"id":129,"date":"2011-02-09T11:39:26","date_gmt":"2011-02-09T15:39:26","guid":{"rendered":"http:\/\/web.colby.edu\/mgimond\/?p=129"},"modified":"2011-10-26T12:26:24","modified_gmt":"2011-10-26T16:26:24","slug":"summarizing-meteorological-data","status":"publish","type":"post","link":"https:\/\/web.colby.edu\/mgimond\/2011\/02\/09\/summarizing-meteorological-data\/","title":{"rendered":"R: Summarizing meteorological data"},"content":{"rendered":"
\n

Dowload data for this exercise here:Waterville meteorological data for 2010<\/a><\/strong><\/p>\n

\n
\r\n\r\n# Load ggplot2 library\r\nlibrary(ggplot2)\r\n\r\n# Load data and convert date field to a date (POSIX) object.\r\n# The date\/time is in standard time.\r\nmet.raw = read.csv ( "waterville_2010.csv",\r\n                     na.strings=c("NaN","M"),header=T)\r\n# extract date\/time column. No DST.\r\nmet.date = strptime(met.raw$Date,\r\n                     format = "%m\/%d\/%Y %H:%M", tz="EST")  \r\nmet.1 = met.raw[, -1] # Remove original data\/time\r\nmet = data.frame(date = met.date,met.1) # Append date\/time formatted data \r\n# column to met.1 (overwrite met object)\r\n#\r\n# Plot monthly statistics\r\nmonth.year = function (x) {as.POSIXlt(x)$mon + 1} # get months lis\r\nmonth = month.year(met$date)\r\nMonth = factor(month, levels = 1:12, month.abb) # Convert to abbreviated months\r\nmet.month = cbind(Month, met)\r\n\r\np1 = ggplot(met.month, aes(Month, OutTemp))\r\np2 = p1 + geom_jitter(position = position_jitter(w=0.15, h=0),\r\n     alpha=I(0.1), aes(colour=(Temp scale_colour_manual(values=c("#ff9999", \r\n     "#9999ff"))+ geom_boxplot(fill="#dddddd", alpha=I(0.5), aes(line=Month),\r\n     outlier.size=1, outlier.colour="#BBBBBB", colour="#666666") +\r\n     scale_y_continuous(name="Temperature (\u00b0F)")+\r\n     opts(legend.position="none",\r\n     title="Waterville Meteorological Data Summary for 2010")\r\n<\/pre>\n
\"\"<\/a><\/p>\n
\r\n# Plot hourly temperatures\r\nhour.day = function(x){as.POSIXlt(x)$hour+1}\r\nhour = hour.day(met$date)\r\nmet.hour = cbind(hour,met)\r\n#\r\nl1 = ggplot(met.hour, aes(as.factor(hour), Temp))\r\nl2 = l1+geom_jitter(position=position_jitter(w=0.15, h=0),\r\nalpha = I(0.1), aes(colour=(Temp<32)))+\r\n        scale_colour_manual(values=c("#ff9999", "#9999ff"))+\r\n        geom_boxplot(fill="#dddddd", alpha=I(0.5), aes(line=hour),\r\n        outlier.size=1, outlier.colour="#BBBBBB",\r\n        colour="#666666")+\r\n        scale_y_continuous(name="Temperature (\u00b0F)")+\r\n        scale_x_discrete(name="Hour of day (EST)")+\r\n        opts(legend.position="none",\r\n        title="Waterville Meteorological Data Summary for 2010")\r\n<\/pre>\n
\"\"<\/a><\/p>\n
\r\n# Plot hourly temperatures by month\r\nmet.hour.month = cbind(hour,met.month)\r\n#\r\nm1 = ggplot(met.hour.month, aes(as.factor(hour), Temp))\r\nm2 = m1+geom_jitter(position=position_jitter(w=0.15, h=0), size=1,\r\nalpha = I(0.7), aes(colour=(Temp<32)))+\r\n        scale_colour_manual(values=c("#ff9999", "#9999ff"))+\r\n        facet_wrap(~Month,ncol=3)+\r\n        geom_boxplot(fill="#dddddd", alpha=I(0.5), aes(line=hour),\r\n        outlier.size=0, outlier.colour="#BBBBBB",\r\n        colour="#666666")+\r\n        scale_y_continuous(name="Temperature (\u00b0F)")+\r\n        scale_x_discrete(name="Hour of day (EST)", breaks=c(6,12,18),\r\n        labels=c(6,12,18))+\r\n        opts(legend.position="none",\r\n        title="Waterville Meteorological Data Summary for 2010")\r\n<\/pre>\n
\"\"<\/a><\/p>\n
\r\n#Save plots to png files\r\nggsave("Waterville_month_Temp.png", plot=p2, width=11, height=8)\r\nggsave("Waterville_hourly_Temp.png", plot=l2, width=11, height=8)\r\nggsave("Waterville_time_month_Temp.png", plot=m2, width=8, height=11)\r\n<\/pre>\n","protected":false},"excerpt":{"rendered":"
Dowload data for this exercise here:Waterville meteorological data for 2010<\/p>\n","protected":false},"author":1199,"featured_media":0,"comment_status":"closed","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"ngg_post_thumbnail":0,"footnotes":""},"categories":[12807],"tags":[],"_links":{"self":[{"href":"https:\/\/web.colby.edu\/mgimond\/wp-json\/wp\/v2\/posts\/129"}],"collection":[{"href":"https:\/\/web.colby.edu\/mgimond\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/web.colby.edu\/mgimond\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/web.colby.edu\/mgimond\/wp-json\/wp\/v2\/users\/1199"}],"replies":[{"embeddable":true,"href":"https:\/\/web.colby.edu\/mgimond\/wp-json\/wp\/v2\/comments?post=129"}],"version-history":[{"count":37,"href":"https:\/\/web.colby.edu\/mgimond\/wp-json\/wp\/v2\/posts\/129\/revisions"}],"predecessor-version":[{"id":167,"href":"https:\/\/web.colby.edu\/mgimond\/wp-json\/wp\/v2\/posts\/129\/revisions\/167"}],"wp:attachment":[{"href":"https:\/\/web.colby.edu\/mgimond\/wp-json\/wp\/v2\/media?parent=129"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/web.colby.edu\/mgimond\/wp-json\/wp\/v2\/categories?post=129"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/web.colby.edu\/mgimond\/wp-json\/wp\/v2\/tags?post=129"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}