OzoneQuality

library(tidyverse)

## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ dplyr     1.1.4     ✔ readr     2.1.5
## ✔ forcats   1.0.0     ✔ stringr   1.5.1
## ✔ ggplot2   3.4.4     ✔ tibble    3.2.1
## ✔ lubridate 1.9.3     ✔ tidyr     1.3.1
## ✔ purrr     1.0.2     
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag()    masks stats::lag()
## ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errors

data("airquality")

head(airquality)

##   Ozone Solar.R Wind Temp Month Day
## 1    41     190  7.4   67     5   1
## 2    36     118  8.0   72     5   2
## 3    12     149 12.6   74     5   3
## 4    18     313 11.5   62     5   4
## 5    NA      NA 14.3   56     5   5
## 6    28      NA 14.9   66     5   6

mean(airquality[,4])

## [1] 77.88235

median(airquality$Temp)

## [1] 79

sd(airquality$Wind)

## [1] 3.523001

var(airquality$Wind)

## [1] 12.41154

airquality$Month[airquality$Month == 5]<- "May"
airquality$Month[airquality$Month == 6]<- "June"
airquality$Month[airquality$Month == 7]<- "July"
airquality$Month[airquality$Month == 8]<- "August"
airquality$Month[airquality$Month == 9]<- "September"

summary(airquality$Month)

##    Length     Class      Mode 
##       153 character character

airquality$Month<-factor(airquality$Month, levels=c("May", "June","July", "August", "September"))

##Plot 5

p5 <- airquality |>
  ggplot(aes(x=Ozone, fill=Ozone)) +
  geom_histogram(position="identity", alpha=0.5, binwidth = 5, color = "white")+
  scale_fill_discrete(name = "Ozone", labels = c("May", "June","July", "August", "September")) +
  labs(x = "Average Monthly Ozone Levels from May - Sept", 
       y = "Frequency of Levels",
       title = "Histogram of Monthly Ozone Levels from May - Sept, 1973",
       caption = "New York State Department of Conservation and the National Weather Service")
p5

## Warning: Removed 37 rows containing non-finite values (`stat_bin()`).

## Warning: The following aesthetics were dropped during statistical transformation: fill
## ℹ This can happen when ggplot fails to infer the correct grouping structure in
##   the data.
## ℹ Did you forget to specify a `group` aesthetic or to convert a numerical
##   variable into a factor?

The plot I created shows a histogram of the monthly ozone levels in New York from May to September of 1973. I created this by using one of the previous plots and changed it so that it would show a different data set, while keeping the format of a histogram.