fMonth.long <- factor(air.long$Month) weather <- ggplot(air.long, aes(x=fMonth.long, y=value)) + geom_boxplot() weather <- weather + facet_wrap(~variable, nrow=2)