---
tags: real experiments
---

# Female Control 14Mar2020

Cage 0 - 18.3g - 66F

Cage 1 - 21.3g - 67F

Cage 2 - 25.2g - 62F

Cage 3 - 17g - 57F

Cage 4 - 18g - 33F

Cage 5 - 19.1g - 6F

Cage 6 - 19.7g - 31F

library(tidyverse)
library(lubridate)
library(readr)	

cageweight0 <- 18.3
cageweight1 <- 21.3
cageweight2 <- 25.2
cageweight3 <- 17
cageweight4 <- 18
cageweight5 <- 19.1
cageweight6 <- 19.7

animalID0 <- "66F"
animalID1 <- "67F"
animalID2 <- "62F"
animalID3 <- "57F"
animalID4 <- "33F"
animalID5 <- "6F"
animalID6 <- "31F"

datafile <- "~/Dropbox/Cactus_Mouse_Physiology/data/14Mar20/mar14.csv"
mar14 <- read_csv(datafile, 
    col_types = cols(Animal = col_double(), 
        StartDate = col_date(format = "%m/%d/%Y"),
        deltaCO2 = col_double(), 
        deltaH2O = col_double(),
        H2Oml = col_double(),
        VCO2 = col_double(),
        StartTime = col_time(format = "%H:%M:%S")))

'%!in%' <- function(x,y)!('%in%'(x,y))

mar14 <- mar14 %>% 
    mutate(EE = 0.06*(3.941*VO2 + 1.106*VCO2)) %>% 
    mutate(RQ = VCO2/VO2) %>%
    mutate(animal = round(Animal, digits=0)) %>%
    mutate(Animal = NULL) %>%
    unite("DateTime", StartDate:StartTime, remove = FALSE, sep =  " ") %>%
    mutate(weight = 
        ifelse(animal == 0, cageweight0, 
        ifelse(animal == 1, cageweight1,
        ifelse(animal == 2, cageweight2,
        ifelse(animal == 3, cageweight3,
        ifelse(animal == 4, cageweight4,
        ifelse(animal == 5, cageweight5,
        ifelse(animal == 6, cageweight6, NA)))))))) %>% 
    mutate(Animal_ID =
        ifelse(animal == 7, 'baseline',
        ifelse(animal == 0, animalID0, 
        ifelse(animal == 1, animalID1,
        ifelse(animal == 2, animalID2,
    ifelse(animal == 3, animalID3,
    ifelse(animal == 4, animalID4,
    ifelse(animal == 5, animalID5,
    ifelse(animal == 6, animalID6, NA))))))))) %>% 
mutate(H2Omg_edit = 
    ifelse(hour(StartTime) == 8, H2Omg - 0.001,
    ifelse(hour(StartTime) == 7, H2Omg - 0.001,
    ifelse(hour(StartTime) == 9, H2Omg - 0.002,
    ifelse(hour(StartTime) == 10, H2Omg - 0.002,
    ifelse(hour(StartTime) == 19, H2Omg + 0.001,
    ifelse(hour(StartTime) == 20, H2Omg + 0.001,
    ifelse(hour(StartTime) == 21, H2Omg + 0.001,
    ifelse(hour(StartTime) == 22, H2Omg + 0.001,
    ifelse(hour(StartTime) %!in% c(7,8,9,10,20,21,22,19), H2Omg, NA)))))))))) %>% 
mutate_at("H2Omg_edit", as.numeric) %>%
mutate(corEE = EE/weight) %>%
mutate(sex = "Female")


metric <- "corEE"

target <- c(0,1,2,4,5,6,7)
cages14mar <- mar14 %>% filter(animal %in% target)

measurement <- cages14mar %>%  select(metric)
df<-as.data.frame(measurement[[metric]])
legend_title <- "Animal ID"

p <- ggplot(data = cages14mar,aes(x=as.POSIXct(StartTime),y=measurement[[metric]]))
p <- p + geom_point(aes(group=as.factor(Animal_ID), color=as.factor(Animal_ID)), size = 3)
p <- p + theme_grey(base_size = 15)
p <- p + geom_smooth(data=df$V1, method='loess', span=.3)
p <- p + labs(x = "", y = metric)
p <- p + scale_color_brewer(legend_title, palette="Paired")
p <- p + scale_x_datetime(date_breaks = "2 hours", date_labels = "%H:%M")
#p <- p + geom_hline(yintercept = 0.8907387)
p