# Install necessary libraries (run once)
install.packages(c("gapminder", "dplyr", "ggplot2", "ggpubr"))

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

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))