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"))

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

summarize(met, Count=n())
summarize(met, MeanCholesterol=mean(Cholesterol), Count=n())

metByStrain <- group_by(met, Strain)
summarize(metByStrain, MeanCholesterol=mean(Cholesterol), Count=n())

group_by(met, Diet) %>%
  summarize(Count=n(), MeanCholesterol=mean(Cholesterol), CholesterolSD = sd(Cholesterol))


metGrouped <- group_by(met, Strain, Diet)
metSummary <- summarize(metGrouped, 
                        Count=n(), 
                        MeanCholesterol=mean(Cholesterol), 
                        CholesterolSD = sd(Cholesterol))
write_csv(metSummary, "~/Desktop/MouseDataSummary.csv")

# All in one pipe:
read_csv("https://denvirlab.marshall.edu/BMR617-2021/data/TH-B6-metabolic.csv") %>%
  separate(MouseID, sep="-", into=c("Strain", "Diet", "Id")) %>%
  group_by(Strain, Diet) %>%
  summarize(
    Count=n(), 
    MeanCholesterol=mean(Cholesterol), 
    CholesterolSD = sd(Cholesterol)
  ) %>%
  write_csv("~/Desktop/MouseDataSummary.csv")