| ggplot(WS_all, aes(x = simpletrt, y = nodulenum_total))+ geom_boxplot(aes(fill = biochardose_pcnt, linetype = pH)) + | |
| facet_grid(.~soil, labeller = soils_labeller) + | |
| labs(y = "Nodule number per plant", | |
| x = "Treatment", | |
| fill = "Biochar dose (%)", | |
| linetype = "pH Adjusted") + | |
| scale_x_discrete(labels = c ("No amendment \ncontrol", | |
| "pH-adjusted \ncontrol", | |
| "WS biochar-\nlow dose", | |
| "WS biochar-\nhigh dose")) + | |
| theme_bw() + | |
| theme(legend.position = c(1,1), | |
| legend.box = "horizontal", | |
| legend.justification = c(1,1), | |
| legend.margin = unit(-0.4, "cm"), | |
| legend.key.size = unit(4.5, "mm"), | |
| axis.title.x = element_text(size = 11), | |
| axis.title.y = element_text(size =11, vjust = 1), | |
| axis.text.x = element_text(size = 9, angle = 45, hjust = 1), | |
| axis.text.y = element_text(size = 9), | |
| strip.text.x = element_text(size = 11), | |
| legend.title = element_text(size = 9), | |
| legend.text = element_text(size = 9), | |
| panel.grid.minor = element_blank(), panel.grid.major = element_blank()) |