# 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