# Reproduce Fig 4.1 (lower panel), p. 63
age <- hav_bg_1964$age
pos <- hav_bg_1964$pos
tot <- hav_bg_1964$tot
plot(
age, pos / tot,
pty = "s", cex = 0.08 * tot, pch = 16, xlab = "age",
ylab = "seroprevalence", xlim = c(0, 86), ylim = c(0, 1)
)
Run the code above in your browser using DataLab