# NOT RUN {
  data("polyps3", package = "HSAUR")
  plot(number3m ~ age, data = polyps3, pch = as.numeric(polyps3$treatment))
  legend("topright", legend = levels(polyps3$treatment), pch = 1:2, bty = "n")
# }
Run the code above in your browser using DataLab