`arbuthnot` <-
structure(list(year = 1629:1710, boys = c(5218, 4858, 4422, 4994, 
5158, 5035, 5106, 4917, 4703, 5359, 5366, 5518, 5470, 5460, 4793, 
4107, 4047, 3768, 3796, 3363, 3079, 2890, 3231, 3220, 3196, 3441, 
3655, 3668, 3396, 3157, 3209, 3724, 4748, 5216, 5411, 6041, 5114, 
4678, 5616, 6073, 6506, 6278, 6449, 6443, 6073, 6113, 6058, 6552, 
6423, 6568, 6247, 6548, 6822, 6909, 7577, 7575, 7484, 7575, 7737, 
7487, 7604, 7909, 7662, 7602, 7676, 6985, 7263, 7632, 8062, 8426, 
7911, 7578, 8102, 8031, 7765, 6113, 8366, 7952, 8239, 8239, 7840, 
7640), girls = c(4683, 4457, 4102, 4590, 4839, 4820, 4928, 4605, 
4457, 4952, 4784, 5332, 5200, 4910, 4617, 3997, 3919, 3536, 3536, 
3181, 2746, 2722, 2840, 2908, 2959, 3179, 3349, 3382, 3289, 3013, 
2781, 3247, 4107, 4803, 4881, 5681, 4858, 4319, 5322, 5560, 5829, 
5719, 6061, 6120, 5822, 5738, 5717, 5847, 6203, 6033, 6041, 6299, 
6533, 6744, 7158, 7127, 7246, 7119, 7214, 7101, 7267, 7302, 7392, 
7316, 7483, 6647, 6713, 7229, 7767, 7626, 7452, 7061, 7514, 7656, 
7683, 5738, 7779, 7417, 7623, 7623, 7380, 7288)), .Names = c("year", 
"boys", "girls"), class = "data.frame", row.names = c(NA, -82L))