# Install packages (run once)
install.packages("gapminder")
install.packages("dplyr")
install.packages("ggplot2")

# Load libraries
library(gapminder)
library(dplyr)
library(ggplot2)

# Load and inspect data
data("gapminder")
head(gapminder)

# One-way ANOVA
anova_model <- aov(lifeExp ~ continent, data = gapminder)
summary(anova_model)

# Boxplot of life expectancy by continent
ggplot(gapminder, aes(x = continent, y = lifeExp)) +
  geom_boxplot() +
  theme_minimal()

# Tukey post-hoc test
TukeyHSD(anova_model)

# Plot Tukey results
plot(TukeyHSD(anova_model))

