# Reproduce Fig 4.4 (middle panel), p. 67
age <- rubella_uk_1986_1987$age
pos <- rubella_uk_1986_1987$pos
tot <- rubella_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