| esoph <- data.frame( |
| agegp = structure(ordered(c(1, 1, 1, 1, 1, 1, |
| 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, |
| 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 4, 4, |
| 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5, 5, |
| 5, 5, 5, 5, 5, 5, 5, 5, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6), levels=1:6), |
| class = c("ordered", "factor"), |
| .Label = c("25-34", "35-44", "45-54", "55-64", "65-74", "75+")), |
| alcgp = structure(ordered(c(1, 1, 1, 1, 2, 2, 2, 2, |
| 3, 3, 3, 4, 4, 4, 4, 1, 1, 1, 1, 2, 2, 2, 2, 3, 3, 3, 3, 4, 4, |
| 4, 1, 1, 1, 1, 2, 2, 2, 2, 3, 3, 3, 3, 4, 4, 4, 4, 1, 1, 1, 1, |
| 2, 2, 2, 2, 3, 3, 3, 3, 4, 4, 4, 4, 1, 1, 1, 1, 2, 2, 2, 3, 3, |
| 3, 3, 4, 4, 4, 4, 1, 1, 1, 2, 2, 2, 2, 3, 3, 4, 4), levels=1:4), |
| class = c("ordered", "factor"), |
| .Label = c("0-39g/day", "40-79", "80-119", "120+")), |
| tobgp = structure(ordered(c(1, |
| 2, 3, 4, 1, 2, 3, 4, 1, 2, 4, 1, 2, 3, 4, 1, 2, 3, 4, 1, 2, 3, |
| 4, 1, 2, 3, 4, 1, 2, 3, 1, 2, 3, 4, 1, 2, 3, 4, 1, 2, 3, 4, 1, |
| 2, 3, 4, 1, 2, 3, 4, 1, 2, 3, 4, 1, 2, 3, 4, 1, 2, 3, 4, 1, 2, |
| 3, 4, 1, 2, 3, 1, 2, 3, 4, 1, 2, 3, 4, 1, 2, 4, 1, 2, 3, 4, 1, |
| 2, 1, 2), levels=1:4), |
| class = c("ordered", "factor"), |
| .Label = c("0-9g/day", "10-19", "20-29", "30+")), |
| ncases = c(0, 0, 0, 0, 0, 0, 0, 0, |
| 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 3, 1, 0, 0, 0, 0, 0, 2, 0, |
| 2, 1, 0, 0, 0, 6, 4, 5, 5, 3, 6, 1, 2, 4, 3, 2, 4, 2, 3, 3, 4, |
| 9, 6, 4, 3, 9, 8, 3, 4, 5, 6, 2, 5, 5, 4, 2, 0, 17, 3, 5, 6, |
| 4, 2, 1, 3, 1, 1, 1, 1, 2, 1, 2, 1, 0, 1, 1, 1, 2, 1), |
| ncontrols = c(40, 10, 6, 5, 27, 7, 4, 7, 2, 1, 2, 1, 1, 1, 2, 60, 14, |
| 7, 8, 35, 23, 14, 8, 11, 6, 2, 1, 3, 3, 4, 46, 18, 10, 4, 38, 21, 15, 7, |
| 16, 14, 5, 4, 4, 4, 3, 4, 49, 22, 12, 6, 40, 21, 17, 6, 18, 15, |
| 6, 4, 10, 7, 3, 6, 48, 14, 7, 2, 34, 10, 9, 13, 12, 3, 1, 4, |
| 2, 1, 1, 18, 6, 3, 5, 3, 3, 1, 1, 1, 2, 1) |
| ) |