# Reproduce Fig 4.5 (left upper panel), p. 68
library(dplyr)
df <- parvob19_be_2001_2003 %>%
group_by(age) %>%
summarise(pos = sum(seropositive), tot = n())
plot(df$age, df$pos / df$tot,
cex = 0.02 * df$tot, pch = 16, xlab = "age", ylab = "seroprevalence",
xlim = c(0, 82), ylim = c(0, 1)
)
Run the code above in your browser using DataLab