# NOT RUN {
data("BCG", package = "HSAUR")
boxplot(BCG$BCGTB/BCG$BCGVacc, BCG$NoVaccTB/BCG$NoVacc,
names = c("BCG Vaccination", "No Vaccination"),
ylab = "Percent BCG cases")
# }
Run the code above in your browser using DataLab