Set Up

knitr::opts_chunk$set(echo = TRUE)
pacman::p_load(readxl, dplyr, ggplot2)

data <- read_excel("C:/Users/DrewIvory/Downloads/HW_Data.xlsx")

## subset data to burned and control
data_burned <- subset(data, Treatment == "Burned")
data_control <- subset(data, Treatment == "Control")

Mean and Median Tree Cover Values by Treatment

mean_burned <- mean(data_burned$`Tree Cover`)
mean_burned # answer key = 25.5%
## [1] 25.9
median_burned <- median(data_burned$`Tree Cover`)
median_burned # answer key = 25%
## [1] 24
mean_control <- mean(data_control$`Tree Cover`)
mean_control # answer key = 66.5%
## [1] 66.4
median_control <- median(data_control$`Tree Cover`)
median_control # answer key = 67.5%
## [1] 68.5

Standard Deviation by Treatment

sd_burned <- sd(data_burned$`Tree Cover`)
sd_burned
## [1] 9.631545
sd_control <- sd(data_control$`Tree Cover`)
sd_control
## [1] 9.811558

Boxplot of Tree Cover by Treatment

ggplot(data, aes(x = Treatment, y = `Tree Cover`, fill = Treatment)) +
  geom_boxplot() +
  labs(title = "Boxplot of Tree Cover by Treatment",
       x = "Treatment",
       y = "Tree Cover (%)") +
  theme_minimal()