Loess Model for Moderate and Severe Drought Ozone Concentration

Loading and Subsetting Data

load("Ozone_Drought_Final.RData")

ozone_df <- combinedAir.final %>%
    mutate(month = as.numeric(month))

agrt_df <- ozone_df %>%
    group_by(MonitorID, USDM.categorical) %>%
    summarize(Max.Ozone = mean(Max.Ozone))

## `summarise()` has grouped output by 'MonitorID'. You can override using the
## `.groups` argument.

## No Drought
no_drt <- agrt_df %>%
    filter(USDM.categorical == "NoDrought")

## Moderate Drought
mdrt_drt <- agrt_df %>%
    filter(USDM.categorical == "ModerateDrought")

## Severe Drought
sev_drt <- agrt_df %>%
    filter(USDM.categorical == "SevereDrought")

# Severe Drought and No Drought Combined
no_sev <- merge(no_drt, sev_drt, by = "MonitorID")

# All Drought combined by MonitorID
no_mdt_sev <- merge(no_sev, mdrt_drt, by = "MonitorID")
head(no_mdt_sev, 3)

##   MonitorID USDM.categorical.x Max.Ozone.x USDM.categorical.y Max.Ozone.y
## 1 010030010          NoDrought    41.72206      SevereDrought    51.98214
## 2 010270001          NoDrought    46.29020      SevereDrought    50.11642
## 3 010331002          NoDrought    41.72557      SevereDrought    49.23077
##   USDM.categorical Max.Ozone
## 1  ModerateDrought  45.86463
## 2  ModerateDrought  48.28937
## 3  ModerateDrought  42.91655

Loess Model for Severe and Moderate Drought against No Drought

ggplot(no_mdt_sev, aes(x = Max.Ozone.x)) +
  geom_point(aes(y = Max.Ozone.y), color = "orange", alpha = 0.8, size = 1) +
  geom_point(aes(y = Max.Ozone), color = "#0d87e4", alpha = 0.8, size = 1) +
  geom_abline(intercept = 0, slope = 1, linetype = "dashed") +
  geom_smooth(aes(y = Max.Ozone), color = "#0d87e4", method = "loess", se = TRUE) +
  geom_smooth(aes(y = Max.Ozone.y), color = "orange", method = "loess", se = TRUE) +
  labs(x = "No Drought Ozone Concentration (ppb)",
      y = "Drought Ozone Concentration (ppb)",
     caption = "Orange <- Severe Drought
     Blue <- Moderate Drought", size =  16, face = "bold") +
  theme_classic() +
  theme(plot.title = element_text(size = 16, face = "bold"),
        plot.caption = element_text(hjust = 0, size = 12, face = "bold"),
        axis.title = element_text(size = 14),
        axis.text = element_text(size = 12),
        axis.line = element_line(colour = "black"),
        legend.title = element_text(size = 12, face = "bold"),
        legend.position = "bottom",
        legend.text = element_text(c("Moderate, Severe"), size = 12),
        panel.grid.major = element_line(colour = "gray90"),
        panel.grid.minor = element_blank(),
        panel.border = element_blank()) +
    xlim(15, 70) +
    ylim(15, 70)

## `geom_smooth()` using formula = 'y ~ x'
## `geom_smooth()` using formula = 'y ~ x'