library(tidyverse)
met <- read_csv("https://denvirlab.marshall.edu/BMR617-2022/data/TH-B6-metabolic.csv") %>%
  separate(MouseID, sep="-", into=c("Strain","Diet","ID"))

b6_bw <- pull((filter(met, Strain=="B6")), BodyWeight)
th_bw <- pull((filter(met, Strain=="TH")), BodyWeight)

shapiro.test(b6_bw)
shapiro.test(th_bw)

var.test(BodyWeight ~ Strain, data = met)

t.test(BodyWeight ~ Strain, data = met, var.equal = TRUE)