# 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