library(tidyverse)
StudentsPerformance <- read.csv("C:/Users/hrida/Downloads/StudentsPerformance.csv")
attach(StudentsPerformance)
data<-aggregate(cbind(math.score,writing.score,reading.score)~test.preparation.course, StudentsPerformance, mean)
data3 = data %>% 
  pivot_longer(-test.preparation.course,names_to="Subject",values_to="Average_Score")
ggplot(data3,aes(fill = test.preparation.course,x=Subject,y=Average_Score))+
  geom_bar(position="dodge", stat="identity") + 
  ggtitle("Scores for preparation course students vs non-preparation students")+ 
  geom_text(aes(label=sprintf("%0.0f", round(Average_Score, digits = 0))),   position=position_dodge(width=0.9), vjust=-0.25)

Data is obtained from link which is a fictional dataset of marks secured by the students in various subjects. I wanted to present the effectiveness of the exam preparation course so compared the scores of students who completed the course vs. students who did not.