library(gapminder)
library(dplyr)
library(ggplot2)
data("gapminder")
head(gapminder)

anova_model <- aov(lifeExp ~ continent, data = gapminder)
summary(anova_model)

ggplot(gapminder, aes(x= continent, y = lifeExp)) +
  geom_boxplot() +
  theme_minimal()

plot(TukeyHSD(anova_model))