Specificaly I have the follwing function where the variable "edu" has the values 1-5:
library(mosaic) gf_boxplot (di01a~edu, data=ALLBUSNEU)