library(tidyverse)

met <- read_csv("https://denvirlab.marshall.edu/BMR617-2021/data/TH-B6-metabolic.csv") %>%
  separate(MouseID, sep="-", into=c("Strain","Diet","Id"))
met %>% group_by(Strain, Diet) %>% summarise(Chol=mean(Cholesterol))

met_summary <- met %>% group_by(Strain, Diet) %>% 
  summarise(Chol=mean(Cholesterol), CholSD = sd(Cholesterol), n=n(), sem=CholSD/sqrt(n))
met_summary


th_summary <- met_summary %>% filter(Strain=="TH")

ggplot(th_summary, aes(x=Diet, y=Chol)) +
  geom_bar(stat="identity", fill="#00b140") +
  geom_errorbar(aes(ymin=Chol - sem, ymax=Chol + sem), width=0.2)

ggplot(met_summary, aes(x=Strain, y=Chol, fill=Diet)) +
  geom_bar(stat="identity", position=position_dodge()) +
  geom_errorbar(aes(ymin=Chol - sem, ymax=Chol + sem), width=0.2, position=position_dodge(0.9))