x <- read.table("samplefile.txt", header=TRUE) # see the "samplefile.txt" file above
with(x, stripchart(age ~ race, vertical=TRUE, method="jitter", pch=19, jitter=0.3, col="gray"))
with(x, boxplot(age ~ race, add=TRUE, varwidth=TRUE, notch=TRUE))
# use the "add = TRUE" argument of the boxplot() function to force the
# boxplot() function (usually a high-level plotting command) to act as a
# low-level plotting command - i.e., to add the boxplot on top of the stripchart
# You can add more to the plot using additional arguments - such as the labels.