I need to make boxplot of two way ANOVA using ggplot. I tried using boxplot, however the visualization of plot it not appropriate especially x-axis is not according bars and so distance between bars. I am wondering if I can use ggplot, please have a look into below data and current plot. I need to arrange all factors (whc) on x-axis according different fields in separate color.
boxplot(numbers ~ whc * field, data = fac, cex.main=1.2, xlab="sites", ylab="mean", outline=FALSE, font.lab=3, ylim = c(0, 2000))
File 1: Example
field whc numbers CL W1.100 43236 CL W1.100.R 15961 CL W1.50 14284 CL W1.50.R 1465 SiC W4.100 344 SiC W4.100.R 432 SiC W4.50 56 SiC W4.50.R 456 L W5.100 4278 L W5.100.R 3743 L W5.50 1844 L W5.50.R 1471