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