Comparison

This document is a comparison of the BSCS Program between the academic years 2024-2025 and 2025-2026. It displays the ABET scores for the two years to facilitate a discussion of progress in the program.

Setup

library(tidyverse)
library(ggbeeswarm)
For_Graphics2526 <- read_delim("For Graphics2526.tsv", 
    delim = "\t", escape_double = FALSE, 
    trim_ws = TRUE)
For_Graphics2425 <- read_delim("~/Library/CloudStorage/Dropbox/Documents/SMU/ABET 24-25/For Graphics2425.tsv", 
    delim = "\t", escape_double = FALSE, 
    trim_ws = TRUE)
View(For_Graphics2425)

Mark and Combine

Combine the data for the two academic years after marking the rows in each of the dataframes with a label for the academic year.

For_Graphics2425 = For_Graphics2425 |> 
  select(program,outcome,score) |> 
  mutate(AY = "24-25") |> 
  filter(program == "CS")
For_Graphics2526 = For_Graphics2526 |> 
  select(program,outcome,score) |>
  mutate(AY = "25-26") |> 
  filter(program == "CS")
combined = rbind(For_Graphics2425,For_Graphics2526)

Overall Scores

Get a statistical summary of the scores by academic year.

tapply(combined$score,combined$AY, summary)

## $`24-25`
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   3.830   4.400   4.550   4.526   4.750   5.000 
## 
## $`25-26`
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   2.000   3.950   4.450   4.282   4.650   5.000

The statistics suggest that the average score for has declined between 2024-25 and 2025-26.

A formal test of this hypothesis follows.

t.test(combined$score~combined$AY)

## 
##  Welch Two Sample t-test
## 
## data:  combined$score by combined$AY
## t = 4.403, df = 180.04, p-value = 1.827e-05
## alternative hypothesis: true difference in means between group 24-25 and group 25-26 is not equal to 0
## 95 percent confidence interval:
##  0.1347907 0.3537173
## sample estimates:
## mean in group 24-25 mean in group 25-26 
##            4.526014            4.281760

The conclusion, based on an extremely low p-value is that the average score has actually declined.

Overall Graphic

The following graphic shows how the average score for each outcome measure has changed.

by_outcome_AY = combined |> 
  group_by(outcome,AY) |> 
  summarize(mean_score = mean(score)) |> 
  ungroup()
#View(by_outcome_AY)

by_outcome_AY |> 
  ggplot(aes(mean_score,outcome,color = AY)) +
  geom_point()

Detail by Outcome and Academic Year

The following graphic displays the individual scores for each outcome measure and academic year.

combined %>% 
  ggplot(aes(factor(AY),score,color = AY)) +
  geom_beeswarm(cex=10) +
  facet_wrap(~outcome)+
  labs(x = "Academic Year",
       y = "Score",
       title = "Detail by Academic Year")