Calculate and graph percent reproductive dry weight= ((reproductive dry weight)/(total dry weight))*100

Load packages and read in data

Convert date column to date

exp$date<-as.Date(exp$date, format = c("%m/%d/%Y"))

Calculate percent reproductive dry weight

exp$dwt<-exp$dwr+exp$dwv

exp$perdwr<-(exp$dwr/exp$dwt)*100

Subset data by site

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

Graph

a<-ggplot(data = run1, mapping = aes(x = treatment, y = perdwr, group=treatment)) + geom_boxplot() + ggtitle("Run 1")+ylim(40,80)+ 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")) 


b<-ggplot(data = run2, mapping = aes(x = treatment, y = perdwr, group=treatment)) + geom_boxplot() + ggtitle("Run 2")+ylim(40,80)+ 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")) 


c<-ggplot(data = run3, mapping = aes(x = treatment, y = perdwr, group=treatment)) + geom_boxplot() + ggtitle("Run 3")+ylim(40,80)+ 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")) 



figure<-ggarrange(a, b, c, ncol = 3, nrow = 1)
## Warning: Removed 41 rows containing non-finite values (stat_boxplot).
## Warning: Removed 42 rows containing non-finite values (stat_boxplot).
## Warning: Removed 40 rows containing non-finite values (stat_boxplot).
annotate_figure(figure,
                top = text_grob("Average percent reproductive dry weight per thalli", color = "black", face = "plain", size = 14), left = text_grob("Average percent reproductive dry weight", color = "black", rot = 90, size=12))