{"id":32,"date":"2015-01-19T11:52:43","date_gmt":"2015-01-19T18:52:43","guid":{"rendered":"https:\/\/blogs.ubc.ca\/wtamagi\/?p=32"},"modified":"2015-01-19T11:54:08","modified_gmt":"2015-01-19T18:54:08","slug":"r-stat-frequency-plot-adding-titles-to-graph-plot","status":"publish","type":"post","link":"https:\/\/blogs.ubc.ca\/wtamagi\/2015\/01\/19\/r-stat-frequency-plot-adding-titles-to-graph-plot\/","title":{"rendered":"R – Stat Frequency Plot – adding titles to graph plot"},"content":{"rendered":"

> AirQualityIndex <- read.csv(“h:\\\\AirQualityIndex.txt”, header=TRUE, stringsAsFactors=FALSE)<\/p>\n

Read in data into table AirQualityIndex<\/p>\n

> AQI=AirQualityIndex$AQI.O3<\/p>\n

Set AQI array to data in column AQI.03<\/p>\n

> range(AQI)
\n[1]\u00a0 8 49<\/p>\n

Range of values for AQI<\/p>\n

> breaks=seq(8,50,by=6)<\/p>\n

Create an array of values broken into intervals of 6 beginning at 8, then 14,…<\/p>\n

> AQI.cut=cut(AQI,breaks,right=F)<\/p>\n

Data is classified into intervals<\/p>\n

> AQI.freq=table(AQI.cut)<\/p>\n

Summarizes data, counting the frequency of data per interval<\/p>\n

> AQI.freq
\nAQI.cut
\n[8,14) [14,20) [20,26) [26,32) [32,38) [38,44) [44,50)
\n5\u00a0\u00a0\u00a0\u00a0\u00a0 12\u00a0\u00a0\u00a0\u00a0\u00a0 10\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 4\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 3\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 1\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 1<\/p>\n

Outputs the the freqency of values per interval<\/p>\n

> nrow(AirQualityIndex)
\n[1] 36<\/p>\n

Row count of values.<\/p>\n

> AQI.relfreq=AQI.freq\/(nrow(AirQualityIndex))<\/p>\n

Calculate the percentage of values per interval divide the interval count by the total number of rows<\/p>\n

> AQI.relfreq
\nAQI.cut
\n[8,14)\u00a0\u00a0\u00a0 [14,20)\u00a0\u00a0\u00a0 [20,26)\u00a0\u00a0\u00a0 [26,32)\u00a0\u00a0\u00a0 [32,38)\u00a0\u00a0\u00a0 [38,44)\u00a0\u00a0\u00a0 [44,50)
\n0.13888889 0.33333333 0.27777778 0.11111111 0.08333333 0.02777778 0.02777778
\n> nrow(AirQualityIndex)
\n[1] 36
\n> midpoints.breaks=seq(11,47,by=6)
\n> midpoints.breaks
\n[1] 11 17 23 29 35 41 47
\n> plot(midpoints.breaks,AQI.relfreq,\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 #plot the data x is midpoints.breaks, y is AQI.relfreq
\n+\u00a0 main=”Air Quality Index, Sudbury Ontario”,\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 #+ main title
\n+\u00a0 xlab=”Daily Ozone Readings”,\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 #+ label – horizontal axis
\n+\u00a0 ylab=”Relative Frequency”)
\n> lines(midpoints.breaks,AQI.relfreq)\u00a0 # connects the dots with lines
\n> axis(side=2) # display y-axis values<\/p>\n

\"frequency<\/a><\/p>\n","protected":false},"excerpt":{"rendered":"

> AirQualityIndex <- read.csv(“h:\\\\AirQualityIndex.txt”, header=TRUE, stringsAsFactors=FALSE) Read in data into table AirQualityIndex > AQI=AirQualityIndex$AQI.O3 Set AQI array to data in column AQI.03 > range(AQI) [1]\u00a0 8 49 Range of values for AQI > breaks=seq(8,50,by=6) Create an array of values broken into intervals of 6 beginning at 8, then 14,… > AQI.cut=cut(AQI,breaks,right=F) Data is classified into […]<\/p>\n","protected":false},"author":25150,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[1],"tags":[],"class_list":["post-32","post","type-post","status-publish","format-standard","hentry","category-uncategorized"],"_links":{"self":[{"href":"https:\/\/blogs.ubc.ca\/wtamagi\/wp-json\/wp\/v2\/posts\/32","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/blogs.ubc.ca\/wtamagi\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/blogs.ubc.ca\/wtamagi\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/blogs.ubc.ca\/wtamagi\/wp-json\/wp\/v2\/users\/25150"}],"replies":[{"embeddable":true,"href":"https:\/\/blogs.ubc.ca\/wtamagi\/wp-json\/wp\/v2\/comments?post=32"}],"version-history":[{"count":2,"href":"https:\/\/blogs.ubc.ca\/wtamagi\/wp-json\/wp\/v2\/posts\/32\/revisions"}],"predecessor-version":[{"id":35,"href":"https:\/\/blogs.ubc.ca\/wtamagi\/wp-json\/wp\/v2\/posts\/32\/revisions\/35"}],"wp:attachment":[{"href":"https:\/\/blogs.ubc.ca\/wtamagi\/wp-json\/wp\/v2\/media?parent=32"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/blogs.ubc.ca\/wtamagi\/wp-json\/wp\/v2\/categories?post=32"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/blogs.ubc.ca\/wtamagi\/wp-json\/wp\/v2\/tags?post=32"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}