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"))
th_chow <- filter(met, Strain=="TH" & Diet=="Chow")
th_chow_chol <- pull(th_chow, Cholesterol)
mean(th_chow_chol)
median(th_chow_chol)

th_hf <- filter(met, Strain=="TH" & Diet=="HF")
th_hf_chol <- pull(th_hf, Cholesterol)
mean(th_hf_chol)
median(th_hf_chol)

sd(th_chow_chol)
IQR(th_chow_chol)