---
tags: Rscripts
---

# taking the mean for every min gam

all_noOL_M <- read.csv("~/Box Sync/UNH/metabolic chamber/Cactus_Mouse_Physiology/data/all_noOL_M.csv")

all_noOL_M$time_bucket <- as.POSIXct(round(as.numeric(as.POSIXct(paste(all_noOL_M$StartTime), format="%H:%M:%S"))/60)*60, origin='1970-01-01')

meanEE_M <- all_noOL_M[ which(all_noOL_M$experiment == "baseline"),] %>% group_by(time_bucket) %>% summarize(mean(EE))
colnames(meanEE_M)<-c("time_bucket","mean_ee")

meanH2O_M <- all_noOL_M[ which(all_noOL_M$experiment == "baseline"),] %>% group_by(time_bucket) %>% summarize(mean(H2Omg))
colnames(meanH2O_M)<-c("time_bucket","mean_h20")

meanRQ_M <- all_noOL_M[ which(all_noOL_M$experiment == "baseline"),] %>% group_by(time_bucket) %>% summarize(mean(RQ))
colnames(meanRQ_M)<-c("time_bucket","mean_rq")

meanDEGc_M <- all_noOL_M[ which(all_noOL_M$experiment == "baseline"),] %>% group_by(time_bucket) %>% summarize(mean(Deg_C))
colnames(meanDEGc_M)<-c("time_bucket","mean_degC")

gam_dataset <- inner_join(inner_join(inner_join(meanEE_M,meanH2O_M),meanRQ_M),meanDEGc_M)

gam_dataset$time <- as.POSIXct(gam_dataset$time_bucket,format="%H:%M:%S")

#gam_dataset<-separate(data=gam_dataset,col=time_bucket,into=c("date","time"),sep = " ")

ggplot(data = meanEE_M,aes(x =as.POSIXct(time_bucket), y=mean_ee))+
  geom_point()