pH readings in the unit tanks and header

Subsetting and setting up data

exp$unit_ph_hach<-as.numeric(exp$unit_ph_hach)

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_ph_hach, group=treatment)) + geom_boxplot() + ggtitle("Run 1")+ylim(5,10)+ 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_ph_hach, group=treatment)) + geom_boxplot() + ggtitle("Run 2")+ylim(5,10)+ 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_ph_hach, group=treatment)) + geom_boxplot() + ggtitle("Run 3")+ylim(5,10)+ 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"))

unitph<-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(unitph,
                top = text_grob("Units", 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_pH_hach, group=treatment)) + geom_boxplot() + ggtitle("Run 1")+ylim(5,10)+ 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_pH_hach, group=treatment)) + geom_boxplot() + ggtitle("Run 2")+ylim(5,10)+ 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_pH_hach, group=treatment)) + geom_boxplot() + ggtitle("Run 3")+ylim(5,10)+ 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"))

headerph<-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(headerph,
                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))

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