Temperature hach readings in the unit tanks and header

Subsetting and setting up data

run1<-subset(exp, run == "1")
run2<-subset(exp, run == "2")
run3<-subset(exp, run == "3")

Graph units

r1<-ggplot(data = run1, mapping = aes(x = treatment, y = unit_temp_hach, group=treatment)) + geom_boxplot() + ggtitle("Run 1")+ylim(14,18)+ theme(plot.title = element_text(size = 12, face = "plain")) +theme(axis.text=element_text(size=12),axis.title=element_text(size=12,face="plain")) + ylab("") + xlab("") + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank(), panel.background = element_blank(), axis.line = element_line(colour = "black")) 

r2<-ggplot(data = run2, mapping = aes(x = treatment, y = unit_temp_hach, group=treatment)) + geom_boxplot() + ggtitle("Run 2")+ylim(14,18)+ theme(plot.title = element_text(size = 12, face = "plain")) +theme(axis.text=element_text(size=12),axis.title=element_text(size=12,face="plain")) + ylab("") + xlab("") + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank(), panel.background = element_blank(), axis.line = element_line(colour = "black"))

r3<-ggplot(data = run3, mapping = aes(x = treatment, y = unit_temp_hach, group=treatment)) + geom_boxplot() + ggtitle("Run 3")+ylim(14,18)+ theme(plot.title = element_text(size = 12, face = "plain")) +theme(axis.text=element_text(size=12),axis.title=element_text(size=12,face="plain")) + ylab("") + xlab("") + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank(), panel.background = element_blank(), axis.line = element_line(colour = "black"))

unittemp<-ggarrange(r1, r2, r3, ncol = 3, nrow = 1)
## Warning: Removed 20 rows containing non-finite values (stat_boxplot).

## Warning: Removed 20 rows containing non-finite values (stat_boxplot).
unit<-annotate_figure(unittemp,
                top = text_grob("Unit", color = "black", face = "plain", size = 14), bottom=text_grob("", color = "black", rot = 0, size=12), left = text_grob("", color = "black", rot = 90, size=12))

Graph headers

r1h<-ggplot(data = run1, mapping = aes(x =treatment, y = header_temp_hach, group=treatment)) + geom_boxplot() + ggtitle("Run 1")+ylim(14,18)+ theme(plot.title = element_text(size = 12, face = "plain")) +theme(axis.text=element_text(size=12),axis.title=element_text(size=12,face="plain")) + ylab("") + xlab("") + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank(), panel.background = element_blank(), axis.line = element_line(colour = "black")) 

r2h<-ggplot(data = run2, mapping = aes(x = treatment, y = header_temp_hach, group=treatment)) + geom_boxplot() + ggtitle("Run 2")+ylim(14,18)+ theme(plot.title = element_text(size = 12, face = "plain")) +theme(axis.text=element_text(size=12),axis.title=element_text(size=12,face="plain")) + ylab("") + xlab("") + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank(), panel.background = element_blank(), axis.line = element_line(colour = "black"))

r3h<-ggplot(data = run3, mapping = aes(x = treatment, y = header_temp_hach, group=treatment)) + geom_boxplot() + ggtitle("Run 3")+ylim(14,18)+ theme(plot.title = element_text(size = 12, face = "plain")) +theme(axis.text=element_text(size=12),axis.title=element_text(size=12,face="plain")) + ylab("") + xlab("") + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank(), panel.background = element_blank(), axis.line = element_line(colour = "black"))

headertemp<-ggarrange(r1h, r2h, r3h, ncol = 3, nrow = 1)
## Warning: Removed 52 rows containing non-finite values (stat_boxplot).

## Warning: Removed 52 rows containing non-finite values (stat_boxplot).
## Warning: Removed 48 rows containing non-finite values (stat_boxplot).
header<-annotate_figure(headertemp,
                top = text_grob("Header Tanks", color = "black", face = "plain", size = 14), bottom=text_grob("", color = "black", rot = 0, size=12), left = text_grob("", color = "black", rot = 90, size=12))

alltemp<-ggarrange(unit, header, ncol = 1, nrow = 2)
annotate_figure(alltemp,
                top = text_grob("Temperature (C)(Hach)", color = "black", face = "plain", size = 16), bottom=text_grob("Treatment", color = "black", rot = 0, size=15), left = text_grob("Temperature (C)", color = "black", rot = 90, size=15))