Perhaps kind of late but I came in this post when looking for groupedboxplot. Since other people can end here too, the answer is:
groupedboxplot(rand(["a","b"], 30), rand(30), group=rand(["C","D"], 30), legend=:bottom, bar_position = :dodge, bar_width=0.7)
