n <- c(20372, 696, 23,  0,
       635, 1441, 307, 7,
       4, 204, 519, 90,
       0, 2, 54, 204)

np <- c(96.6, 3.3, 0.1, 0.0,
        26.6, 60.3, 12.8, 0.3,
        0.5, 25.0, 63.5, 11.0,
        0.0, 0.8, 20.8, 78.5)

cases <- c(218, 38, 0, 0,
           22, 96, 36, 1,
           0, 7 , 59, 24,
           0, 0, 11, 48)

casesp <- c(82.2, 14.8, 0.0, 0.0,
           14.2, 61.9, 23.2, 0.6,
           0.0, 7.8, 65.6, 26.7,
           0.0, 0.0, 18.6, 81.4)

controls <- c(19933, 642, 23, 0,
              589, 1291, 263, 6,
              4, 188, 434, 58,
              0, 2 , 38, 140)

controlsp <- c(96.8, 3.1, 0.1, 0.0,
               27.4, 60.1, 12.2, 0.3,
               0.6, 27.5, 63.4, 8.5,
               0.0, 1.1, 21.1, 77.8)

observedrisk <- c(1.3, 6.8, 0.0, NA,
                  4.4, 8.4, 14.6, 17.5,
                  0.0, 4.3, 14.3, 34.2,
                  NA, 0.0, 25.9, 29.4)

df <- data.frame(n, np, cases, casesp, controls, controlsp, observedrisk)