---
tags: Rscripts
---

# heatmap for t tests
note not publication ready and still need work to look nice. Could not get loop or function to work so I copy and pasted

VO2

vo2 <- read.csv("/Users/danielleblumstein/Documents/UNH/metabolic chamber/stats/temp_experiment/VO2_t_tests.csv")

#convert to long table
VO2_long <- reshape(data=vo2, idvar="VO2",
                    varying = c("BL_day_noOL_F","BL_night_noOL_F","cold_day_noOL_F","cold_night_noOL_F","hot_day_noOL_F","hot_night_noOL_F"),
                    v.name=c("p_val"),
                    times=c("BL_day_noOL_F","BL_night_noOL_F","cold_day_noOL_F","cold_night_noOL_F","hot_day_noOL_F","hot_night_noOL_F"),
                    new.row.names = 1:1000,
                    direction="long")

#rename columns
colnames(VO2_long)<-c("VO2_1","VO2_2","p_val")
#set pval as numeric for color coding
VO2_long$p_val<-as.numeric(VO2_long$p_val)

#plot using geom_raster
vo2_plot <- ggplot(data=VO2_long)+
  geom_raster(aes(x=VO2_1,y=VO2_2,fill=p_val))+
  #scale_fill_viridis(option="virdis", alpha = 0.75) +
  geom_text(aes(x=VO2_1,y=VO2_2,label=p_val),size=3,color="black")+
  theme(axis.text.x.top = element_text(angle = 45, hjust = 0, vjust=0))+
  scale_x_discrete(position="top")+ylab("VO2_2")


VCO2

vco2 <- read.csv("/Users/danielleblumstein/Documents/UNH/metabolic chamber/stats/temp_experiment/VCO2_t_tests.csv")

#convert to long table
VCO2_long <- reshape(data=vco2, idvar="VCO2",
                    varying = c("BL_day_noOL_F","BL_night_noOL_F","cold_day_noOL_F","cold_night_noOL_F","hot_day_noOL_F","hot_night_noOL_F"),
                    v.name=c("p_val"),
                    times=c("BL_day_noOL_F","BL_night_noOL_F","cold_day_noOL_F","cold_night_noOL_F","hot_day_noOL_F","hot_night_noOL_F"),
                    new.row.names = 1:1000,
                    direction="long")

#rename columns
colnames(VCO2_long)<-c("VCO2_1","VCO2_2","p_val")
#set pval as numeric for color coding
VCO2_long$p_val<-as.numeric(VCO2_long$p_val)

#plot using geom_raster
vco2_plot <- ggplot(data=VCO2_long)+
  geom_raster(aes(x=VCO2_1,y=VCO2_2,fill=p_val))+
  #scale_fill_viridis(option="virdis", alpha = 0.75) +
  geom_text(aes(x=VCO2_1,y=VCO2_2,label=p_val),size=3,color="black")+
  theme(axis.text.x.top = element_text(angle = 45, hjust = 0, vjust=0))+
  scale_x_discrete(position="top")+ylab("VCO2_2")


H2Omg

h2o <- read.csv("/Users/danielleblumstein/Documents/UNH/metabolic chamber/stats/temp_experiment/H2O_t_tests.csv")

#convert to long table
H2O_long <- reshape(data=h2o, idvar="H2OMG",
                    varying = c("BL_day_noOL_F","BL_night_noOL_F","cold_day_noOL_F","cold_night_noOL_F","hot_day_noOL_F","hot_night_noOL_F"),
                    v.name=c("p_val"),
                    times=c("BL_day_noOL_F","BL_night_noOL_F","cold_day_noOL_F","cold_night_noOL_F","hot_day_noOL_F","hot_night_noOL_F"),
                    new.row.names = 1:1000,
                    direction="long")

#rename columns
colnames(H2O_long)<-c("H2O_1","H2O_2","p_val")
#set pval as numeric for color coding
H2O_long$p_val<-as.numeric(H2O_long$p_val)

#plot using geom_raster
h2o_plot <- ggplot(data=H2O_long)+
  geom_raster(aes(x=H2O_1,y=H2O_2,fill=p_val))+
  #scale_fill_viridis(option="virdis", alpha = 0.75) +
  geom_text(aes(x=H2O_1,y=H2O_2,label=p_val),size=3,color="black")+
  theme(axis.text.x.top = element_text(angle = 45, hjust = 0, vjust=0))+
  scale_x_discrete(position="top")+ylab("H2O_2")


RQ

rq <- read.csv("/Users/danielleblumstein/Documents/UNH/metabolic chamber/stats/temp_experiment/RQ_t_tests.csv")

#convert to long table
RQ_long <- reshape(data=rq, idvar="H2OMG",
                    varying = c("BL_day_noOL_F","BL_night_noOL_F","cold_day_noOL_F","cold_night_noOL_F","hot_day_noOL_F","hot_night_noOL_F"),
                    v.name=c("p_val"),
                    times=c("BL_day_noOL_F","BL_night_noOL_F","cold_day_noOL_F","cold_night_noOL_F","hot_day_noOL_F","hot_night_noOL_F"),
                    new.row.names = 1:1000,
                    direction="long")

#rename columns
colnames(RQ_long)<-c("RQ_1","RQ_2","p_val")
#set pval as numeric for color coding
RQ_long$p_val<-as.numeric(RQ_long$p_val)

#plot using geom_raster
rq_plot <- ggplot(data=RQ_long)+
  geom_raster(aes(x=RQ_1,y=RQ_2,fill=p_val))+
  #scale_fill_viridis(option="virdis", alpha = 0.75) +
  geom_text(aes(x=RQ_1,y=RQ_2,label=p_val),size=3,color="black")+
  theme(axis.text.x.top = element_text(angle = 45, hjust = 0, vjust=0))+
  scale_x_discrete(position="top")+ylab("RQ_2")


EE

ee <- read.csv("/Users/danielleblumstein/Documents/UNH/metabolic chamber/stats/temp_experiment/EE_t_tests.csv")

#convert to long table
EE_long <- reshape(data=ee, idvar="EE",
                    varying = c("BL_day_noOL_F","BL_night_noOL_F","cold_day_noOL_F","cold_night_noOL_F","hot_day_noOL_F","hot_night_noOL_F"),
                    v.name=c("p_val"),
                    times=c("BL_day_noOL_F","BL_night_noOL_F","cold_day_noOL_F","cold_night_noOL_F","hot_day_noOL_F","hot_night_noOL_F"),
                    new.row.names = 1:1000,
                    direction="long")

#rename columns
colnames(EE_long)<-c("EE_1","EE_2","p_val")
#set pval as numeric for color coding
EE_long$p_val<-as.numeric(EE_long$p_val)

#plot using geom_raster
ee_plot <- ggplot(data=EE_long)+
  geom_raster(aes(x=EE_1,y=EE_2,fill=p_val))+
  #scale_fill_viridis(option="virdis", alpha = 0.75) +
  geom_text(aes(x=EE_1,y=EE_2,label=p_val),size=3,color="black")+
  theme(axis.text.x.top = element_text(angle = 45, hjust = 0, vjust=0))+
  scale_x_discrete(position="top")+ylab("EE_2")