# Reproduce Fig 4.4 (left panel), p. 67
age <- mumps_uk_1986_1987$age
pos <- mumps_uk_1986_1987$pos
tot <- mumps_uk_1986_1987$tot
plot(age, pos / tot,
cex = 0.008 * tot, pch = 16, xlab = "age", ylab = "seroprevalence",
xlim = c(0, 45), ylim = c(0, 1)
)
Run the code above in your browser using DataLab