library(tidyverse)
met <- read_csv("https://denvirlab.marshall.edu/BMR617-2021/data/TH-B6-metabolic.csv")
met <- separate(met, MouseID, sep="-", into=c("Strain", "Diet", "ID"))
chow_only <- filter(met, Diet=="Chow")
ggplot(chow_only, aes(x=Strain, y=BodyWeight)) + geom_boxplot() +  geom_point()

t.test(BodyWeight ~ Strain, data=chow_only)

