diff --git a/DESCRIPTION b/DESCRIPTION index 1e8126a..734e502 100644 --- a/DESCRIPTION +++ b/DESCRIPTION @@ -8,7 +8,7 @@ Authors@R: c( person("Brandon", "Bertelsen", , "brandon@bertelsen.ca", role = c("aut", "cre")), person("Santiago", "Olivella", , "olivella@unc.edu", role = "aut"), person("Evan", "Rosenman", , "etrrosenman@gmail.com", role = "aut"), - person("Alex", "Rossell Hayes", , "alexander@rossellhayes.com", role = "aut"), + person("Alexander", "Rossell Hayes", , "alexander@rossellhayes.com", role = "aut"), person("Kosuke", "Imai", , "imai@harvard.edu", role = "aut") ) Description: Predicts individual race/ethnicity using surname, first name, diff --git a/R/census_geo_api_names.R b/R/census_geo_api_names.R index 3e9514a..0abc53e 100644 --- a/R/census_geo_api_names.R +++ b/R/census_geo_api_names.R @@ -33,36 +33,78 @@ census_geo_api_names <- function( prefix <- "P12" separator <- "_" suffix <- "N" + sex_codes <- c("_mal" = 2, "_fem" = 26) + age_codes <- 1:23 + names(age_codes) <- paste0("_", age_codes) } else if (year %in% c("2010", "2000")) { - prefix <- "P012" + prefix <- "PCT012" separator <- "" suffix <- "" + sex_codes <- c("_mal" = 2, "_fem" = 106) + age_codes <- list( + "_1" = 1:5, + "_2" = 6:10, + "_3" = 11:15, + "_4" = 16:18, + "_5" = 19:20, + "_6" = 21, + "_7" = 22, + "_8" = 23:25, + "_9" = 26:30, + "_10" = 31:35, + "_11" = 36:40, + "_12" = 41:45, + "_13" = 46:50, + "_14" = 51:55, + "_15" = 56:60, + "_16" = 61:62, + "_17" = 63:65, + "_18" = 66:67, + "_19" = 68:70, + "_20" = 71:75, + "_21" = 76:80, + "_22" = 81:85, + "_23" = 86:103 + ) } race_codes <- list( "_whi" = "I", - "_bla" = "B", + "_bla" = "J", "_his" = "H", - "_asi" = c("D", "E"), - "_oth" = c("C", "F", "G") + "_asi" = c("L", "M"), + "_oth" = c("K", "N", "O") ) - sex_codes <- c("_mal" = 2, "_fem" = 26) - - age_codes <- 1:23 - names(age_codes) <- paste0("_", age_codes) - numeric_codes <- if (age && sex) { - age_sex_codes <- unlist( - purrr::map(sex_codes, function(x) x + age_codes) + age_sex_codes <- purrr::imap( + sex_codes, + function(sex_code, name) { + codes <- purrr::map( + age_codes, + function(age_code) { + str_pad(age_code + sex_code, 3, "left", pad = "0") + } + ) + names(codes) <- paste0(name, names(codes)) + codes + } ) - names(age_sex_codes) <- sub(".", "", names(age_sex_codes), fixed = TRUE) - age_sex_codes[] <- str_pad(age_sex_codes, 3, "left", pad = "0") - as.list(age_sex_codes) + + do.call(c, unname(age_sex_codes)) } else if (age) { purrr::map( age_codes, - function(x) str_pad(x + sex_codes, 3, "left", pad = "0") + function(age_code) { + unlist( + purrr::map( + sex_codes, + function(sex_code) { + str_pad(age_code + sex_code, 3, "left", pad = "0") + } + ) + ) + } ) } else if (sex) { sex_codes[] <- str_pad(sex_codes, 3, "left", pad = "0") diff --git a/inst/CITATION b/inst/CITATION index 7967551..159fa69 100644 --- a/inst/CITATION +++ b/inst/CITATION @@ -6,8 +6,9 @@ bibentry( person("Brandon", "Bertelsen", , "brandon@bertelsen.ca", role = c("aut", "cre")), person("Santiago", "Olivella", , "olivella@unc.edu", role = "aut"), person("Evan", "Rosenman", , "etrrosenman@gmail.com", role = "aut"), + person("Alexander", "Rossell Hayes", , "alexander@rossellhayes.com", role = "aut"), person("Kosuke", "Imai", , "imai@harvard.edu", role = "aut") ), - year = 2023, + year = 2024, url = "https://CRAN.R-project.org/package=wru" ) diff --git a/tests/testthat/_snaps/census_geo_api.md b/tests/testthat/_snaps/census_geo_api.md index faf5502..08d8ba2 100644 --- a/tests/testthat/_snaps/census_geo_api.md +++ b/tests/testthat/_snaps/census_geo_api.md @@ -1,359 +1,918 @@ # snapshot structure(list(state = c("DE", "DE", "DE"), county = c("003", - "005", "001"), P12I_001N = c(303265, 171741, 104845), P12B_001N = c(146544, - 25357, 46998), P12H_001N = c(63516, 26793, 13981), P12D_001N = c(35200, - 3070, 4429), P12E_001N = c(175, 112, 125), P12C_001N = c(2156, - 1843, 1149), P12F_001N = c(28714, 14137, 5658), P12G_001N = c(43699, - 17012, 15807), r_whi = c(0.523005047848499, 0.296181260358264, - 0.180813691793237), r_bla = c(0.669459431061814, 0.115838811506677, - 0.21470175743151), r_his = c(0.609032505513472, 0.256908620193691, - 0.134058874292837), r_asi = c(0.820556238547007, 0.0738094685811046, - 0.105634292871889), r_oth = c(0.572836566160937, 0.253443441521029, - 0.173719992318033)), class = c("tbl_df", "tbl", "data.frame"), row.names = c(NA, + "005", "001"), P12I_001N = c(303265, 171741, 104845), P12J_001N = c(142388, + 24835, 45737), P12H_001N = c(63516, 26793, 13981), P12L_001N = c(35022, + 3014, 4362), P12M_001N = c(127, 70, 107), P12K_001N = c(893, + 764, 864), P12N_001N = c(2555, 1018, 1028), P12O_001N = c(22953, + 9143, 10927), r_whi = c(0.523005047848499, 0.296181260358264, + 0.180813691793237), r_bla = c(0.668613824192337, 0.116618144252442, + 0.214768031555222), r_his = c(0.609032505513472, 0.256908620193691, + 0.134058874292837), r_asi = c(0.823123038733549, 0.0722214416186596, + 0.104655519647792), r_oth = c(0.526493169807558, 0.217868182271413, + 0.255638647921029)), class = c("tbl_df", "tbl", "data.frame"), row.names = c(NA, -3L)) --- structure(list(state = c("DE", "DE", "DE"), county = c("001", - "005", "003"), P012I001 = c(105891, 149025, 331836), P012B001 = c(38913, - 25115, 127786), P012H001 = c(9346, 16954, 46921), P012D001 = c(3306, - 1943, 23300), P012E001 = c(91, 150, 159), P012C001 = c(1043, - 1564, 1574), P012F001 = c(3235, 8169, 19115), P012G001 = c(5723, - 4541, 13590), r_whi = c(0.180469772578534, 0.253982943390052, - 0.565547284031414), r_bla = c(0.202868403766148, 0.130934134109085, - 0.666197462124767), r_his = c(0.127640977315251, 0.23154559484301, - 0.64081342784174), r_asi = c(0.117344295139728, 0.0722995612974541, - 0.810356143562817), r_oth = c(0.170799603784541, 0.243774976944359, - 0.5854254192711)), class = c("tbl_df", "tbl", "data.frame"), row.names = c(NA, + "005", "003"), PCT012I001 = c(105891, 149025, 331836), PCT012J001 = c(37812, + 24544, 124426), PCT012H001 = c(9346, 16954, 46921), PCT012L001 = c(3266, + 1910, 23132), PCT012M001 = c(74, 62, 102), PCT012K001 = c(916, + 924, 984), PCT012N001 = c(347, 305, 873), PCT012O001 = c(4658, + 3421, 10205), r_whi = c(0.180469772578534, 0.253982943390052, + 0.565547284031414), r_bla = c(0.20243920720412, 0.131404525061301, + 0.666156267734578), r_his = c(0.127640977315251, 0.23154559484301, + 0.64081342784174), r_asi = c(0.117004133678974, 0.0690814825194423, + 0.813914383801583), r_oth = c(0.261609154774003, 0.205452215791101, + 0.532938629434896)), class = c("tbl_df", "tbl", "data.frame"), row.names = c(NA, -3L)) --- structure(list(state = c("DE", "DE", "DE"), county = c("001", - "005", "003"), P012I001 = c(105891, 149025, 331836), P012B001 = c(38913, - 25115, 127786), P012H001 = c(9346, 16954, 46921), P012D001 = c(3306, - 1943, 23300), P012E001 = c(91, 150, 159), P012C001 = c(1043, - 1564, 1574), P012F001 = c(3235, 8169, 19115), P012G001 = c(5723, - 4541, 13590), P012I002 = c(51574, 72032, 160694), P012B002 = c(18194, - 12160, 60237), P012H002 = c(4579, 9154, 24193), P012D002 = c(1355, - 889, 11467), P012E002 = c(44, 83, 71), P012C002 = c(509, 828, - 767), P012F002 = c(1660, 4492, 10111), P012G002 = c(2705, 2207, - 6433), P012I026 = c(54317, 76993, 171142), P012B026 = c(20719, - 12955, 67549), P012H026 = c(4767, 7800, 22728), P012D026 = c(1951, - 1054, 11833), P012E026 = c(47, 67, 88), P012C026 = c(534, 736, - 807), P012F026 = c(1575, 3677, 9004), P012G026 = c(3018, 2334, - 7157), r_whi = c(0.180469772578534, 0.253982943390052, 0.565547284031414 - ), r_bla = c(0.202868403766148, 0.130934134109085, 0.666197462124767 - ), r_his = c(0.127640977315251, 0.23154559484301, 0.64081342784174 - ), r_asi = c(0.117344295139728, 0.0722995612974541, 0.810356143562817 - ), r_oth = c(0.170799603784541, 0.243774976944359, 0.5854254192711 - ), r_mal_whi = c(0.0878974421902269, 0.122763961605585, 0.273870391579407 - ), r_mal_bla = c(0.0948523048369775, 0.0633947469944842, 0.314038599893647 - ), r_mal_his = c(0.0625367039510523, 0.125018778765655, 0.330410674533262 - ), r_mal_asi = c(0.0483263670593112, 0.0335762893364192, 0.398562990086013 - ), r_mal_oth = c(0.0832394029442907, 0.128548006967927, 0.295641629948424 - ), r_fem_whi = c(0.0925723303883072, 0.131218981784468, 0.291676892452007 - ), r_fem_bla = c(0.108016098929171, 0.0675393871146006, 0.35215886223112 - ), r_fem_his = c(0.0651042733641988, 0.106526816077355, 0.310402753308477 - ), r_fem_asi = c(0.0690179280804173, 0.0387232719610349, 0.411793153476804 - ), r_fem_oth = c(0.08756020084025, 0.115226969976432, 0.289783789322677 - )), class = c("tbl_df", "tbl", "data.frame"), row.names = c(NA, + "005", "003"), PCT012I001 = c(105891, 149025, 331836), PCT012J001 = c(37812, + 24544, 124426), PCT012H001 = c(9346, 16954, 46921), PCT012L001 = c(3266, + 1910, 23132), PCT012M001 = c(74, 62, 102), PCT012K001 = c(916, + 924, 984), PCT012N001 = c(347, 305, 873), PCT012O001 = c(4658, + 3421, 10205), PCT012I002 = c(51574, 72032, 160694), PCT012J002 = c(17702, + 11860, 58678), PCT012H002 = c(4579, 9154, 24193), PCT012L002 = c(1339, + 875, 11390), PCT012M002 = c(33, 29, 42), PCT012K002 = c(448, + 481, 473), PCT012N002 = c(163, 152, 431), PCT012O002 = c(2197, + 1650, 4770), PCT012I106 = c(54317, 76993, 171142), PCT012J106 = c(20110, + 12684, 65748), PCT012H106 = c(4767, 7800, 22728), PCT012L106 = c(1927, + 1035, 11742), PCT012M106 = c(41, 33, 60), PCT012K106 = c(468, + 443, 511), PCT012N106 = c(184, 153, 442), PCT012O106 = c(2461, + 1771, 5435), r_whi = c(0.180469772578534, 0.253982943390052, + 0.565547284031414), r_bla = c(0.20243920720412, 0.131404525061301, + 0.666156267734578), r_his = c(0.127640977315251, 0.23154559484301, + 0.64081342784174), r_asi = c(0.117004133678974, 0.0690814825194423, + 0.813914383801583), r_oth = c(0.261609154774003, 0.205452215791101, + 0.532938629434896), r_mal_whi = c(0.0878974421902269, 0.122763961605585, + 0.273870391579407), r_mal_bla = c(0.0947735863198809, 0.0634964825304365, + 0.314152327312054), r_mal_his = c(0.0625367039510523, 0.125018778765655, + 0.330410674533262), r_mal_asi = c(0.0480627758705248, 0.0316681846843691, + 0.400476424017375), r_mal_oth = c(0.124066628374497, 0.100870410462599, + 0.250695886537357), r_fem_whi = c(0.0925723303883072, 0.131218981784468, + 0.291676892452007), r_fem_bla = c(0.107665620884239, 0.0679080425308649, + 0.352003940422525), r_fem_his = c(0.0651042733641988, 0.106526816077355, + 0.310402753308477), r_fem_asi = c(0.0689413578084495, 0.0374132978350732, + 0.413437959784208), r_fem_oth = c(0.137542526399505, 0.104581805328503, + 0.282242742897539)), class = c("tbl_df", "tbl", "data.frame"), row.names = c(NA, -3L)) --- structure(list(state = c("DE", "DE", "DE"), county = c("001", - "005", "003"), P012I001 = c(105891, 149025, 331836), P012B001 = c(38913, - 25115, 127786), P012H001 = c(9346, 16954, 46921), P012D001 = c(3306, - 1943, 23300), P012E001 = c(91, 150, 159), P012C001 = c(1043, - 1564, 1574), P012F001 = c(3235, 8169, 19115), P012G001 = c(5723, - 4541, 13590), P012I003 = c(3159, 3151, 7793), P012I027 = c(2913, - 3035, 7497), P012B003 = c(1488, 1030, 4669), P012B027 = c(1422, - 991, 4637), P012H003 = c(600, 1253, 2735), P012H027 = c(616, - 1198, 2630), P012D003 = c(86, 51, 917), P012E003 = c(1, 11, 6 - ), P012D027 = c(104, 54, 872), P012E027 = c(3, 4, 4), P012C003 = c(26, - 78, 64), P012F003 = c(178, 587, 1029), P012G003 = c(538, 443, - 1203), P012C027 = c(22, 73, 62), P012F027 = c(195, 570, 1008), - P012G027 = c(541, 493, 1252), P012I004 = c(3265, 3384, 8494 - ), P012I028 = c(3030, 3125, 8084), P012B004 = c(1484, 952, - 4917), P012B028 = c(1548, 966, 4743), P012H004 = c(525, 986, - 2529), P012H028 = c(504, 911, 2426), P012D004 = c(111, 74, - 796), P012E004 = c(6, 10, 5), P012D028 = c(118, 74, 863), - P012E028 = c(3, 6, 10), P012C004 = c(24, 68, 60), P012F004 = c(163, - 475, 1012), P012G004 = c(458, 355, 1030), P012C028 = c(25, - 48, 64), P012F028 = c(156, 450, 947), P012G028 = c(488, 374, - 1031), P012I005 = c(3292, 3462, 9316), P012I029 = c(3044, - 3367, 8870), P012B005 = c(1584, 989, 5262), P012B029 = c(1522, - 920, 5015), P012H005 = c(381, 666, 2147), P012H029 = c(440, - 695, 2124), P012D005 = c(106, 54, 672), P012E005 = c(5, 1, - 10), P012D029 = c(93, 85, 645), P012E029 = c(3, 3, 9), P012C005 = c(34, - 42, 56), P012F005 = c(131, 316, 845), P012G005 = c(393, 283, - 769), P012C029 = c(31, 38, 70), P012F029 = c(124, 331, 850 - ), P012G029 = c(413, 283, 853), P012I006 = c(2040, 2352, - 6242), P012I030 = c(1942, 2303, 5956), P012B006 = c(965, - 652, 3416), P012B030 = c(974, 653, 3273), P012H006 = c(271, - 411, 1290), P012H030 = c(243, 332, 1172), P012D006 = c(59, - 41, 399), P012E006 = c(2, 1, 2), P012D030 = c(52, 49, 362 - ), P012E030 = c(3, 3, 8), P012C006 = c(28, 34, 50), P012F006 = c(79, - 185, 544), P012G006 = c(185, 154, 453), P012C030 = c(20, - 30, 49), P012F030 = c(86, 163, 495), P012G030 = c(183, 161, - 404), P012I007 = c(1389, 1421, 5291), P012I031 = c(1311, - 1311, 5727), P012B007 = c(1018, 447, 2272), P012B031 = c(1174, - 382, 2093), P012H007 = c(174, 276, 945), P012H031 = c(207, - 194, 944), P012D007 = c(49, 18, 362), P012E007 = c(0, 1, - 4), P012D031 = c(29, 26, 354), P012E031 = c(0, 4, 0), P012C007 = c(14, - 20, 28), P012F007 = c(71, 147, 362), P012G007 = c(120, 81, - 279), P012C031 = c(16, 17, 26), P012F031 = c(66, 92, 344), - P012G031 = c(150, 79, 289), P012I008 = c(678, 708, 2730), - P012I032 = c(643, 641, 3067), P012B008 = c(463, 194, 1040 - ), P012B032 = c(575, 170, 1003), P012H008 = c(101, 163, 508 - ), P012H032 = c(77, 99, 457), P012D008 = c(13, 12, 217), - P012E008 = c(3, 1, 2), P012D032 = c(24, 9, 173), P012E032 = c(3, - 1, 2), P012C008 = c(8, 13, 14), P012F008 = c(37, 88, 215), - P012G008 = c(54, 19, 126), P012C032 = c(4, 10, 13), P012F032 = c(25, - 40, 181), P012G032 = c(42, 23, 140), P012I009 = c(738, 577, - 2685), P012I033 = c(735, 589, 2939), P012B009 = c(419, 170, - 934), P012B033 = c(463, 170, 906), P012H009 = c(91, 185, - 485), P012H033 = c(97, 124, 429), P012D009 = c(15, 11, 139 - ), P012E009 = c(1, 2, 3), P012D033 = c(16, 6, 185), P012E033 = c(1, - 1, 3), P012C009 = c(5, 9, 10), P012F009 = c(36, 102, 188), - P012G009 = c(41, 34, 102), P012C033 = c(13, 13, 16), P012F033 = c(35, - 64, 172), P012G033 = c(59, 30, 131), P012I010 = c(2122, 1951, - 6605), P012I034 = c(2010, 1970, 6403), P012B010 = c(945, - 500, 2657), P012B034 = c(955, 559, 2756), P012H010 = c(302, - 623, 1397), P012H034 = c(260, 482, 1156), P012D010 = c(74, - 35, 482), P012E010 = c(4, 5, 5), P012D034 = c(67, 33, 510 - ), P012E034 = c(3, 2, 6), P012C010 = c(25, 38, 27), P012F010 = c(114, - 329, 627), P012G010 = c(115, 72, 236), P012C034 = c(19, 34, - 27), P012F034 = c(99, 252, 466), P012G034 = c(136, 79, 311 - ), P012I011 = c(3277, 3402, 10187), P012I035 = c(3380, 3429, - 10014), P012B011 = c(1158, 837, 4120), P012B035 = c(1340, - 842, 4478), P012H011 = c(431, 1169, 2415), P012H035 = c(426, - 835, 2092), P012D011 = c(87, 49, 1151), P012E011 = c(5, 15, - 6), P012D035 = c(125, 64, 1355), P012E035 = c(2, 11, 6), - P012C011 = c(28, 64, 50), P012F011 = c(186, 638, 1181), P012G011 = c(167, - 129, 370), P012C035 = c(34, 57, 55), P012F035 = c(146, 413, - 910), P012G035 = c(187, 147, 465), P012I012 = c(2963, 3214, - 9253), P012I036 = c(3226, 3168, 9400), P012B012 = c(922, - 699, 4038), P012B036 = c(1228, 767, 4581), P012H012 = c(366, - 930, 2184), P012H036 = c(345, 777, 2003), P012D012 = c(91, - 71, 1283), P012E012 = c(1, 4, 8), P012D036 = c(141, 68, 1221 - ), P012E036 = c(5, 5, 4), P012C012 = c(37, 59, 49), P012F012 = c(161, - 461, 1012), P012G012 = c(121, 119, 386), P012C036 = c(28, - 48, 61), P012F036 = c(147, 385, 892), P012G036 = c(150, 113, - 420), P012I013 = c(3048, 3523, 9894), P012I037 = c(3327, - 3718, 9876), P012B013 = c(1045, 707, 4174), P012B037 = c(1282, - 712, 4821), P012H013 = c(297, 767, 1866), P012H037 = c(332, - 657, 1784), P012D013 = c(123, 70, 1173), P012E013 = c(2, - 8, 3), P012D037 = c(163, 93, 1091), P012E037 = c(2, 9, 10 - ), P012C013 = c(32, 48, 66), P012F013 = c(129, 431, 837), - P012G013 = c(110, 82, 320), P012C037 = c(37, 64, 69), P012F037 = c(127, - 337, 764), P012G037 = c(132, 90, 380), P012I014 = c(3603, - 4282, 11235), P012I038 = c(3657, 4435, 11648), P012B014 = c(1142, - 842, 4436), P012B038 = c(1347, 937, 5499), P012H014 = c(240, - 565, 1637), P012H038 = c(265, 448, 1581), P012D014 = c(126, - 83, 864), P012E014 = c(3, 4, 4), P012D038 = c(171, 86, 896 - ), P012E038 = c(2, 2, 5), P012C014 = c(28, 50, 52), P012F014 = c(83, - 288, 722), P012G014 = c(85, 90, 297), P012C038 = c(48, 46, - 53), P012F038 = c(95, 202, 658), P012G038 = c(119, 72, 350 - ), P012I015 = c(4108, 5377, 13313), P012I039 = c(4280, 5636, - 13731), P012B015 = c(1261, 918, 4578), P012B039 = c(1413, - 973, 5489), P012H015 = c(229, 353, 1290), P012H039 = c(245, - 311, 1190), P012D015 = c(106, 68, 752), P012E015 = c(1, 4, - 6), P012D039 = c(178, 100, 806), P012E039 = c(4, 3, 8), P012C015 = c(42, - 76, 64), P012F015 = c(105, 161, 548), P012G015 = c(70, 92, - 252), P012C039 = c(40, 48, 62), P012F039 = c(94, 137, 470 - ), P012G039 = c(105, 101, 313), P012I016 = c(3831, 5618, - 13725), P012I040 = c(4119, 6144, 14291), P012B016 = c(1073, - 875, 4074), P012B040 = c(1320, 982, 4946), P012H016 = c(172, - 259, 957), P012H040 = c(193, 212, 816), P012D016 = c(81, - 65, 622), P012E016 = c(6, 4, 2), P012D040 = c(137, 80, 646 - ), P012E040 = c(9, 4, 2), P012C016 = c(47, 61, 57), P012F016 = c(64, - 98, 379), P012G016 = c(79, 68, 180), P012C040 = c(46, 36, - 60), P012F040 = c(65, 85, 285), P012G040 = c(90, 64, 252), - P012I017 = c(3198, 5713, 11775), P012I041 = c(3583, 6570, - 12620), P012B017 = c(913, 703, 3013), P012B041 = c(1161, - 801, 3955), P012H017 = c(129, 180, 679), P012H041 = c(113, - 161, 626), P012D017 = c(51, 44, 490), P012E017 = c(1, 3, - 3), P012D041 = c(133, 60, 552), P012E041 = c(2, 1, 3), P012C017 = c(37, - 49, 50), P012F017 = c(45, 78, 247), P012G017 = c(45, 40, - 135), P012C041 = c(45, 41, 39), P012F041 = c(27, 66, 201), - P012G041 = c(56, 52, 200), P012I018 = c(1273, 2551, 4128), - P012I042 = c(1382, 2830, 4509), P012B018 = c(320, 227, 1110 - ), P012B042 = c(422, 275, 1377), P012H018 = c(42, 59, 202 - ), P012H042 = c(61, 52, 200), P012D018 = c(28, 20, 165), - P012E018 = c(0, 1, 0), P012D042 = c(49, 27, 182), P012E042 = c(0, - 3, 2), P012C018 = c(13, 11, 7), P012F018 = c(11, 23, 65), - P012G018 = c(18, 20, 46), P012C042 = c(18, 14, 9), P012F042 = c(13, - 22, 57), P012G042 = c(21, 21, 55), P012I019 = c(1925, 4078, - 5826), P012I043 = c(2108, 4596, 6378), P012B019 = c(450, - 319, 1340), P012B043 = c(550, 350, 1776), P012H019 = c(51, - 71, 233), P012H043 = c(79, 68, 272), P012D019 = c(29, 29, - 217), P012E019 = c(0, 1, 2), P012D043 = c(75, 29, 289), P012E043 = c(0, - 2, 2), P012C019 = c(19, 24, 15), P012F019 = c(21, 27, 76), - P012G019 = c(30, 22, 63), P012C043 = c(20, 21, 15), P012F043 = c(15, - 18, 83), P012G043 = c(28, 34, 79), P012I020 = c(1123, 2492, - 3106), P012I044 = c(1247, 2728, 3446), P012B020 = c(230, - 160, 721), P012B044 = c(296, 208, 933), P012H020 = c(17, - 42, 111), P012H044 = c(31, 33, 140), P012D020 = c(24, 7, - 116), P012E020 = c(0, 2, 0), P012D044 = c(50, 28, 144), P012E044 = c(0, - 0, 0), P012C020 = c(9, 13, 7), P012F020 = c(2, 14, 39), P012G020 = c(21, - 16, 37), P012C044 = c(6, 4, 7), P012F044 = c(7, 6, 46), P012G044 = c(19, - 13, 26), P012I021 = c(1507, 3532, 4076), P012I045 = c(1706, - 3896, 4876), P012B021 = c(330, 256, 923), P012B045 = c(397, - 270, 1213), P012H021 = c(38, 50, 150), P012H045 = c(55, 49, - 169), P012D021 = c(25, 26, 200), P012E021 = c(1, 0, 0), P012D045 = c(65, - 20, 180), P012E045 = c(0, 1, 0), P012C021 = c(13, 21, 8), - P012F021 = c(7, 12, 47), P012G021 = c(13, 20, 40), P012C045 = c(10, - 17, 14), P012F045 = c(11, 12, 47), P012G045 = c(36, 27, 52 - ), P012I022 = c(1923, 4552, 5078), P012I046 = c(2285, 4883, - 6167), P012B022 = c(448, 262, 1128), P012B046 = c(497, 331, - 1537), P012H022 = c(54, 71, 184), P012H046 = c(65, 66, 204 - ), P012D022 = c(40, 34, 234), P012E022 = c(0, 1, 0), P012D046 = c(61, - 28, 245), P012E046 = c(0, 0, 0), P012C022 = c(16, 27, 15), - P012F022 = c(20, 17, 68), P012G022 = c(24, 35, 43), P012C046 = c(20, - 31, 12), P012F046 = c(15, 10, 59), P012G046 = c(26, 27, 42 - ), P012I023 = c(1563, 3347, 4147), P012I047 = c(1778, 3561, - 5592), P012B023 = c(286, 186, 724), P012B047 = c(397, 280, - 1074), P012H023 = c(27, 39, 114), P012H047 = c(44, 40, 142 - ), P012D023 = c(14, 17, 122), P012E023 = c(0, 1, 0), P012D047 = c(57, - 21, 129), P012E047 = c(1, 1, 3), P012C023 = c(13, 11, 10), - P012F023 = c(5, 9, 28), P012G023 = c(6, 17, 22), P012C047 = c(12, - 23, 9), P012F047 = c(9, 12, 33), P012G047 = c(24, 23, 45), - P012I024 = c(929, 2123, 3204), P012I048 = c(1259, 2540, 4733 - ), P012B024 = c(141, 120, 395), P012B048 = c(200, 211, 747 - ), P012H024 = c(22, 22, 84), P012H048 = c(39, 29, 88), P012D024 = c(11, - 8, 53), P012E024 = c(1, 3, 0), P012D048 = c(26, 10, 72), - P012E048 = c(1, 0, 0), P012C024 = c(8, 8, 5), P012F024 = c(7, - 4, 26), P012G024 = c(6, 9, 28), P012C048 = c(15, 9, 5), P012F048 = c(12, - 7, 19), P012G048 = c(9, 12, 26), P012I025 = c(620, 1222, - 2591), P012I049 = c(1352, 2518, 5318), P012B025 = c(109, - 115, 296), P012B049 = c(236, 205, 697), P012H025 = c(19, - 14, 51), P012H049 = c(30, 27, 83), P012D025 = c(6, 2, 41), - P012E025 = c(1, 0, 0), P012D049 = c(17, 4, 61), P012E049 = c(0, - 1, 1), P012C025 = c(3, 4, 3), P012F025 = c(5, 2, 14), P012G025 = c(6, - 7, 16), P012C049 = c(5, 14, 10), P012F049 = c(6, 3, 17), - P012G049 = c(4, 16, 41), r_whi = c(0.180469772578534, 0.253982943390052, - 0.565547284031414), r_bla = c(0.202868403766148, 0.130934134109085, - 0.666197462124767), r_his = c(0.127640977315251, 0.23154559484301, - 0.64081342784174), r_asi = c(0.117344295139728, 0.0722995612974541, - 0.810356143562817), r_oth = c(0.170799603784541, 0.243774976944359, - 0.5854254192711), r_1_whi = c(0.0103484947643979, 0.0105427846858639, - 0.0260587096422339), r_1_bla = c(0.0151709468547656, 0.0105362486575537, - 0.0485157496324564), r_1_his = c(0.0166072574807774, 0.0334740033596919, - 0.0732713292634627), r_1_asi = c(0.00670144046426474, 0.00414522090573077, - 0.0621437700784138), r_1_oth = c(0.0256173788297981, 0.038323598729378, - 0.0788673702906719), r_2_whi = c(0.0107285531195462, 0.0110932727966841, - 0.028253844895288), r_2_bla = c(0.01580697967823, 0.00999927012626815, - 0.050361287497263), r_2_his = c(0.0140533453517434, 0.0259078679613772, - 0.0676718427773453), r_2_asi = c(0.00822135479636602, 0.00566513523783205, - 0.0578258316349442), r_2_oth = c(0.0224408238549032, 0.0302285070191618, - 0.0707722785804556), r_3_whi = c(0.0107984293193717, 0.0116386480148342, - 0.0309943553664921), r_3_bla = c(0.0161927700793477, 0.00995234967207816, - 0.0535779453011772), r_3_his = c(0.0112126302563472, 0.0185875636770872, - 0.058330260444408), r_3_asi = c(0.00715050606238558, 0.00493972157932916, - 0.0461501260838025), r_3_oth = c(0.0192301123749018, 0.022082180551286, - 0.0588004235406633), r_4_whi = c(0.00678651287085515, 0.00793350512652705, - 0.0207890215968586), r_4_bla = c(0.0101087511860448, 0.0068034658575495, - 0.0348723242307652), r_4_his = c(0.00701984403381544, 0.0101473620955737, - 0.0336242334849292), r_4_asi = c(0.00400704687553974, 0.0032470897094891, - 0.0266330443193202), r_4_oth = c(0.00992246473340848, 0.0124158896061755, - 0.0340711138436315), r_5_whi = c(0.00460160340314136, 0.00465614092495637, - 0.0187779504799302), r_5_bla = c(0.0114277372871636, 0.00432189516927857, - 0.0227564202821483), r_5_his = c(0.00520342524685541, 0.00641892353286625, - 0.0257986096884773), r_5_asi = c(0.002694393588725, 0.00169263186984006, - 0.0248713254343846), r_5_oth = c(0.00746319636574786, 0.00744611811319466, - 0.0226799193906479), r_6_whi = c(0.00225137707242583, 0.00229909740401396, - 0.00987981293630017), r_6_bla = c(0.00541149238324627, 0.00189767170279542, - 0.0106509431011292), r_6_his = c(0.00243099657202169, 0.00357820843747012, - 0.0131792791685445), r_6_asi = c(0.00148537082455353, 0.000794500673598397, - 0.013610141973816), r_6_oth = c(0.00290330293404379, 0.00329610274276736, - 0.0117669160091539), r_7_whi = c(0.00251043030104712, 0.00198721095113438, - 0.00958496945898778), r_7_bla = c(0.00459820451061966, 0.00177255049162209, - 0.00959262618995485), r_7_his = c(0.0025675694131465, 0.00422010079075675, - 0.012482757678808), r_7_asi = c(0.00113993574907596, 0.000690870150955128, - 0.0113993574907596), r_7_oth = c(0.00322778973255457, 0.00430371964340609, - 0.01057143833043), r_8_whi = c(0.007042157504363, 0.00668255071989529, - 0.022169502617801), r_8_bla = c(0.00990542921788816, 0.00552097344302293, - 0.0282200465033835), r_8_his = c(0.00767539367121454, 0.0150912989442919, - 0.034867046339165), r_8_asi = c(0.00511243911706795, 0.00259076306608173, - 0.0346471380703997), r_8_oth = c(0.00867575229702497, 0.0137309150527718, - 0.0289305598251187), r_9_whi = c(0.0113455088350785, 0.0116420566099476, - 0.034428514943281), r_9_bla = c(0.0130230327296235, 0.0087532713983338, - 0.0448246739028434), r_9_his = c(0.0117042924843966, 0.0273691973614127, - 0.0615533794949536), r_9_asi = c(0.00756502815295865, 0.00480154754913814, - 0.0869805520052506), r_9_oth = c(0.0127745329097927, 0.0247293096970318, - 0.0517641834887454), r_10_whi = c(0.010547897578534, 0.0108768270069808, - 0.0317902623254799), r_10_bla = c(0.0112087751676103, 0.00764282064917055, - 0.04493415496262), r_10_his = c(0.00971032900397427, 0.0233129839800057, - 0.0571830485789596), r_10_asi = c(0.00822135479636602, 0.00511243911706795, - 0.0869114649901551), r_10_oth = c(0.01099839464426, 0.0202377292755405, - 0.0481606722000205), r_11_whi = c(0.0108648969240838, 0.0123408186082024, - 0.0336939626963351), r_11_bla = c(0.0121315441000136, 0.00739779161062279, - 0.0468943872710021), r_11_his = c(0.0085904317067508, 0.0194479725761735, - 0.0498490870105571), r_11_asi = c(0.0100176171888494, 0.00621783135859615, - 0.0786555666862413), r_11_oth = c(0.0096833691976637, 0.0179663216859651, - 0.0416026232195922), r_12_whi = c(0.0123732002617801, 0.0148563618019197, - 0.0389994409904014), r_12_bla = c(0.0129761122754335, 0.00927460977822265, - 0.0517949680419573), r_12_his = c(0.0068969284768031, 0.0138348288059437, - 0.0439491402739651), r_12_asi = c(0.0104321392794224, 0.00604511382085737, - 0.0611074648519811), r_12_oth = c(0.00782183966936503, 0.0127745329097927, - 0.0364108344434197), r_13_whi = c(0.0142956479057592, 0.0187694289921466, - 0.0460910231239093), r_13_bla = c(0.0139405882782279, 0.00985850876369817, - 0.0524831347034106), r_13_his = c(0.00647355266931618, 0.00906843665068764, - 0.0338700645989539), r_13_asi = c(0.0099830736813016, 0.00604511382085737, - 0.0543023938650731), r_13_oth = c(0.00778768316425863, 0.0105031253202172, - 0.0291867336134167), r_14_whi = c(0.0135491655759162, 0.0200459478621291, - 0.0477476003490401), r_14_bla = c(0.0124756274307402, 0.00968125371453596, - 0.0470247218659743), r_14_his = c(0.00498490870105571, 0.00643258081697874, - 0.0242143647314295), r_14_asi = c(0.00804863725862724, 0.00528515665480673, - 0.0439393416007461), r_14_oth = c(0.00667759674830071, 0.00703624005191789, - 0.0207159203470301), r_15_whi = c(0.0115568417321117, 0.0209338868891798, - 0.0415763388961606), r_15_bla = c(0.0108125579988948, 0.00784092923352831, - 0.0363268583106551), r_15_his = c(0.0033050627552205, 0.00465713388235616, - 0.0178227557667882), r_15_asi = c(0.00645963591143045, 0.00373069881515769, - 0.0362015959100487), r_15_oth = c(0.00435495440106568, 0.0055675103323428, - 0.0148922362263893), r_16_whi = c(0.00452491001308901, 0.00917082515270506, - 0.0147200179973822), r_16_bla = c(0.00386833077877527, 0.00261711866704203, - 0.0129656855078357), r_16_his = c(0.00140670026358558, 0.00151595853648543, - 0.00549022821321752), r_16_asi = c(0.00265985008117724, 0.00176171888493558, - 0.012055684134167), r_16_oth = c(0.00160535574000068, 0.00189568603340506, - 0.0040817023602145), r_17_whi = c(0.00687343204624782, 0.0147830770069808, - 0.020799247382199), r_17_bla = c(0.00521338379888851, 0.00348775376145641, - 0.0162449039173366), r_17_his = c(0.00177544693462258, 0.00189836249163491, - 0.0068969284768031), r_17_asi = c(0.00359252478496667, 0.00210715396041314, - 0.0176171888493558), r_17_oth = c(0.00227140758957543, 0.00249342487276702, - 0.00565290159510879), r_18_whi = c(0.00403918520942408, 0.0088964332460733, - 0.011166557591623), r_18_bla = c(0.00274223987821535, 0.00191852523799097, - 0.00862293680336159), r_18_his = c(0.000655549637399107, - 0.0010242963084361, 0.00342797831223283), r_18_asi = c(0.00255621955853397, - 0.00127810977926699, 0.00898131196241666), r_18_oth = c(0.00109300816340472, - 0.00112716466851112, 0.0027666769136182), r_19_whi = c(0.00547590804973822, - 0.0126595222513089, 0.0152568717277487), r_19_bla = c(0.00379013002179194, - 0.00274223987821535, 0.0111357877944258), r_19_his = c(0.00127012742246077, - 0.00135207112713566, 0.00435667363188156), r_19_asi = c(0.00314345918684583, - 0.00162354485474455, 0.0131265328681474), r_19_oth = c(0.00153704272978789, - 0.00186152952829866, 0.00355227653106534), r_20_whi = c(0.00717168411867365, - 0.016080047447644, 0.0191648260253054), r_20_bla = c(0.00492664768994964, - 0.00309153659274088, 0.0138936678240379), r_20_his = c(0.00162521680938529, - 0.00187104792340995, 0.00529902623564278), r_20_asi = c(0.0034888942623234, - 0.00217624097550865, 0.0165463401153753), r_20_oth = c(0.00206646855893705, - 0.00251050312532022, 0.0040817023602145), r_21_whi = c(0.00569405813699825, - 0.011773287521815, 0.0165981539048866), r_21_bla = c(0.00356074113464085, - 0.00242943685028204, 0.00937366407040154), r_21_his = c(0.000969667171986179, - 0.00107892544488603, 0.00349626473279524), r_21_asi = c(0.00248713254343846, - 0.00138174030191026, 0.00877405091713013), r_21_oth = c(0.00117839942617071, - 0.00162243399255388, 0.00251050312532022), r_22_whi = c(0.00372900305410122, - 0.0079471395069808, 0.0135270097076789), r_22_bla = c(0.00177776387542098, - 0.0017256300374321, 0.00595368429833067), r_22_his = c(0.000833094330861365, - 0.000696521489736551, 0.0023490528673468), r_22_asi = c(0.0013471967943625, - 0.000725413658502884, 0.00431793844346955), r_22_oth = c(0.000973460395532329, - 0.000836834375106739, 0.00186152952829866), r_23_whi = c(0.00336087478184991, - 0.00637407286212914, 0.0134792893760908), r_23_bla = c(0.00179861741061653, - 0.00166828281564432, 0.00517689011229629), r_23_his = c(0.000669206921511588, - 0.000559948648611737, 0.00183007607107251), r_23_asi = c(0.000829044181146154, - 0.000241804552834295, 0.00355798127741891), r_23_oth = c(0.000495269324042764, - 0.000785599617447143, 0.00172490350787307)), class = c("tbl_df", + "005", "003"), PCT012I001 = c(105891, 149025, 331836), PCT012J001 = c(37812, + 24544, 124426), PCT012H001 = c(9346, 16954, 46921), PCT012L001 = c(3266, + 1910, 23132), PCT012M001 = c(74, 62, 102), PCT012K001 = c(916, + 924, 984), PCT012N001 = c(347, 305, 873), PCT012O001 = c(4658, + 3421, 10205), PCT012I003 = c(617, 593, 1507), PCT012I004 = c(631, + 603, 1546), PCT012I005 = c(613, 647, 1540), PCT012I006 = c(636, + 683, 1573), PCT012I007 = c(662, 625, 1627), PCT012I107 = c(563, + 605, 1377), PCT012I108 = c(607, 590, 1424), PCT012I109 = c(601, + 575, 1553), PCT012I110 = c(551, 618, 1548), PCT012I111 = c(591, + 647, 1595), PCT012J003 = c(280, 216, 872), PCT012J004 = c(272, + 212, 897), PCT012J005 = c(260, 182, 945), PCT012J006 = c(309, + 182, 881), PCT012J007 = c(277, 204, 857), PCT012J107 = c(260, + 197, 865), PCT012J108 = c(256, 177, 862), PCT012J109 = c(273, + 207, 928), PCT012J110 = c(273, 188, 875), PCT012J111 = c(275, + 183, 868), PCT012H003 = c(116, 254, 523), PCT012H004 = c(121, + 259, 530), PCT012H005 = c(138, 252, 551), PCT012H006 = c(113, + 223, 558), PCT012H007 = c(112, 265, 573), PCT012H107 = c(126, + 220, 520), PCT012H108 = c(95, 243, 478), PCT012H109 = c(137, + 252, 525), PCT012H110 = c(124, 242, 584), PCT012H111 = c(134, + 241, 523), PCT012L003 = c(18, 7, 175), PCT012M003 = c(0, 0, 0 + ), PCT012L004 = c(9, 11, 167), PCT012M004 = c(0, 0, 2), PCT012L005 = c(17, + 14, 184), PCT012M005 = c(0, 3, 0), PCT012L006 = c(25, 8, 181), + PCT012M006 = c(0, 0, 0), PCT012L007 = c(15, 9, 200), PCT012M007 = c(0, + 1, 1), PCT012L107 = c(17, 9, 171), PCT012M107 = c(1, 1, 0 + ), PCT012L108 = c(18, 9, 175), PCT012M108 = c(0, 0, 1), PCT012L109 = c(23, + 11, 178), PCT012M109 = c(0, 0, 0), PCT012L110 = c(21, 9, + 190), PCT012M110 = c(1, 0, 0), PCT012L111 = c(24, 15, 148 + ), PCT012M111 = c(0, 1, 0), PCT012K003 = c(7, 2, 10), PCT012N003 = c(7, + 1, 13), PCT012O003 = c(91, 66, 176), PCT012K004 = c(6, 4, + 7), PCT012N004 = c(2, 11, 7), PCT012O004 = c(84, 73, 191), + PCT012K005 = c(3, 4, 2), PCT012N005 = c(3, 7, 11), PCT012O005 = c(67, + 69, 180), PCT012K006 = c(2, 6, 4), PCT012N006 = c(9, 5, 18 + ), PCT012O006 = c(88, 67, 168), PCT012K007 = c(4, 2, 3), + PCT012N007 = c(4, 6, 14), PCT012O007 = c(88, 67, 172), PCT012K107 = c(2, + 2, 6), PCT012N107 = c(6, 6, 13), PCT012O107 = c(87, 69, 184 + ), PCT012K108 = c(2, 3, 8), PCT012N108 = c(8, 4, 10), PCT012O108 = c(84, + 78, 190), PCT012K109 = c(4, 1, 6), PCT012N109 = c(6, 1, 9 + ), PCT012O109 = c(80, 79, 188), PCT012K110 = c(2, 7, 9), + PCT012N110 = c(11, 3, 11), PCT012O110 = c(87, 68, 168), PCT012K111 = c(4, + 4, 5), PCT012N111 = c(6, 3, 17), PCT012O111 = c(84, 76, 171 + ), PCT012I008 = c(645, 646, 1592), PCT012I009 = c(669, 672, + 1751), PCT012I010 = c(646, 643, 1717), PCT012I011 = c(624, + 692, 1646), PCT012I012 = c(681, 731, 1788), PCT012I112 = c(591, + 609, 1630), PCT012I113 = c(612, 595, 1604), PCT012I114 = c(609, + 644, 1588), PCT012I115 = c(607, 617, 1564), PCT012I116 = c(611, + 660, 1698), PCT012J008 = c(265, 200, 913), PCT012J009 = c(286, + 199, 901), PCT012J010 = c(279, 167, 951), PCT012J011 = c(265, + 183, 985), PCT012J012 = c(317, 161, 953), PCT012J112 = c(298, + 192, 851), PCT012J113 = c(325, 170, 904), PCT012J114 = c(267, + 181, 874), PCT012J115 = c(287, 187, 899), PCT012J116 = c(288, + 206, 988), PCT012H008 = c(104, 235, 523), PCT012H009 = c(112, + 201, 512), PCT012H010 = c(111, 204, 521), PCT012H011 = c(113, + 189, 478), PCT012H012 = c(85, 157, 495), PCT012H112 = c(111, + 234, 506), PCT012H113 = c(96, 169, 530), PCT012H114 = c(123, + 200, 496), PCT012H115 = c(86, 165, 448), PCT012H116 = c(88, + 143, 446), PCT012L008 = c(23, 12, 158), PCT012M008 = c(1, + 0, 0), PCT012L009 = c(21, 9, 163), PCT012M009 = c(1, 0, 1 + ), PCT012L010 = c(18, 13, 166), PCT012M010 = c(0, 1, 2), + PCT012L011 = c(22, 26, 148), PCT012M011 = c(0, 1, 0), PCT012L012 = c(25, + 12, 150), PCT012M012 = c(0, 0, 0), PCT012L112 = c(17, 17, + 209), PCT012M112 = c(0, 0, 1), PCT012L113 = c(19, 20, 196 + ), PCT012M113 = c(1, 2, 2), PCT012L114 = c(25, 17, 166), + PCT012M114 = c(1, 0, 1), PCT012L115 = c(23, 11, 157), PCT012M115 = c(0, + 0, 1), PCT012L116 = c(31, 8, 120), PCT012M116 = c(1, 0, 2 + ), PCT012K008 = c(5, 7, 10), PCT012N008 = c(3, 4, 9), PCT012O008 = c(73, + 72, 170), PCT012K009 = c(4, 3, 6), PCT012N009 = c(7, 0, 9 + ), PCT012O009 = c(68, 57, 161), PCT012K010 = c(0, 9, 7), + PCT012N010 = c(2, 5, 10), PCT012O010 = c(65, 51, 160), PCT012K011 = c(4, + 8, 5), PCT012N011 = c(7, 5, 5), PCT012O011 = c(80, 58, 144 + ), PCT012K012 = c(2, 5, 8), PCT012N012 = c(13, 4, 12), PCT012O012 = c(78, + 51, 130), PCT012K112 = c(3, 6, 7), PCT012N112 = c(1, 3, 8 + ), PCT012O112 = c(86, 50, 160), PCT012K113 = c(1, 7, 3), + PCT012N113 = c(7, 7, 14), PCT012O113 = c(82, 56, 173), PCT012K114 = c(8, + 5, 7), PCT012N114 = c(6, 3, 13), PCT012O114 = c(64, 57, 139 + ), PCT012K115 = c(3, 2, 8), PCT012N115 = c(9, 6, 8), PCT012O115 = c(77, + 60, 154), PCT012K116 = c(5, 2, 9), PCT012N116 = c(3, 8, 7 + ), PCT012O116 = c(90, 65, 172), PCT012I013 = c(668, 699, + 1828), PCT012I014 = c(670, 682, 1832), PCT012I015 = c(663, + 684, 1872), PCT012I016 = c(639, 678, 1843), PCT012I017 = c(652, + 719, 1941), PCT012I117 = c(631, 646, 1736), PCT012I118 = c(625, + 678, 1701), PCT012I119 = c(617, 651, 1785), PCT012I120 = c(573, + 683, 1785), PCT012I121 = c(598, 709, 1863), PCT012J013 = c(320, + 196, 1000), PCT012J014 = c(293, 211, 1008), PCT012J015 = c(336, + 187, 997), PCT012J016 = c(282, 179, 1013), PCT012J017 = c(310, + 178, 1035), PCT012J117 = c(265, 168, 1015), PCT012J118 = c(308, + 164, 965), PCT012J119 = c(305, 181, 944), PCT012J120 = c(299, + 195, 929), PCT012J121 = c(282, 169, 972), PCT012H013 = c(78, + 165, 482), PCT012H014 = c(79, 123, 410), PCT012H015 = c(78, + 138, 444), PCT012H016 = c(72, 114, 413), PCT012H017 = c(74, + 126, 398), PCT012H117 = c(109, 169, 464), PCT012H118 = c(94, + 139, 458), PCT012H119 = c(84, 124, 419), PCT012H120 = c(86, + 130, 393), PCT012H121 = c(67, 133, 390), PCT012L013 = c(22, + 6, 151), PCT012M013 = c(2, 0, 2), PCT012L014 = c(17, 22, + 119), PCT012M014 = c(0, 0, 0), PCT012L015 = c(27, 8, 113), + PCT012M015 = c(0, 0, 0), PCT012L016 = c(15, 12, 122), PCT012M016 = c(1, + 0, 2), PCT012L017 = c(22, 5, 156), PCT012M017 = c(1, 0, 0 + ), PCT012L117 = c(17, 19, 145), PCT012M117 = c(0, 0, 3), + PCT012L118 = c(22, 14, 133), PCT012M118 = c(1, 1, 1), PCT012L119 = c(17, + 14, 130), PCT012M119 = c(1, 0, 1), PCT012L120 = c(17, 20, + 123), PCT012M120 = c(0, 0, 0), PCT012L121 = c(16, 17, 109 + ), PCT012M121 = c(0, 0, 1), PCT012K013 = c(5, 6, 8), PCT012N013 = c(5, + 7, 4), PCT012O013 = c(70, 61, 134), PCT012K014 = c(8, 3, + 6), PCT012N014 = c(5, 4, 11), PCT012O014 = c(59, 47, 100), + PCT012K015 = c(4, 5, 6), PCT012N015 = c(3, 5, 14), PCT012O015 = c(72, + 40, 151), PCT012K016 = c(6, 4, 5), PCT012N016 = c(4, 2, 5 + ), PCT012O016 = c(68, 30, 118), PCT012K017 = c(5, 5, 7), + PCT012N017 = c(6, 3, 12), PCT012O017 = c(59, 45, 92), PCT012K117 = c(4, + 1, 12), PCT012N117 = c(4, 4, 8), PCT012O117 = c(64, 67, 143 + ), PCT012K118 = c(5, 3, 5), PCT012N118 = c(7, 1, 4), PCT012O118 = c(72, + 46, 139), PCT012K119 = c(10, 3, 3), PCT012N119 = c(6, 9, + 11), PCT012O119 = c(62, 46, 109), PCT012K120 = c(5, 2, 8), + PCT012N120 = c(5, 5, 7), PCT012O120 = c(67, 37, 124), PCT012K121 = c(2, + 6, 3), PCT012N121 = c(1, 3, 9), PCT012O121 = c(71, 33, 136 + ), PCT012I018 = c(701, 810, 1931), PCT012I019 = c(673, 754, + 2123), PCT012I020 = c(666, 788, 2188), PCT012I122 = c(640, + 719, 1958), PCT012I123 = c(644, 799, 1911), PCT012I124 = c(658, + 785, 2087), PCT012J018 = c(266, 204, 1087), PCT012J019 = c(312, + 230, 1084), PCT012J020 = c(336, 196, 1157), PCT012J122 = c(269, + 204, 1031), PCT012J123 = c(332, 222, 1054), PCT012J124 = c(335, + 213, 1085), PCT012H018 = c(86, 145, 420), PCT012H019 = c(99, + 123, 414), PCT012H020 = c(86, 143, 456), PCT012H122 = c(88, + 102, 373), PCT012H123 = c(75, 122, 416), PCT012H124 = c(80, + 108, 383), PCT012L018 = c(19, 11, 148), PCT012M018 = c(2, + 0, 1), PCT012L019 = c(25, 11, 128), PCT012M019 = c(0, 0, + 1), PCT012L020 = c(15, 15, 122), PCT012M020 = c(0, 1, 0), + PCT012L122 = c(19, 13, 118), PCT012M122 = c(0, 1, 2), PCT012L123 = c(19, + 17, 111), PCT012M123 = c(1, 0, 1), PCT012L124 = c(12, 18, + 130), PCT012M124 = c(2, 0, 1), PCT012K018 = c(10, 8, 10), + PCT012N018 = c(2, 4, 5), PCT012O018 = c(58, 34, 136), PCT012K019 = c(10, + 8, 10), PCT012N019 = c(4, 6, 10), PCT012O019 = c(49, 36, + 95), PCT012K020 = c(5, 4, 7), PCT012N020 = c(4, 6, 8), PCT012O020 = c(38, + 43, 103), PCT012K122 = c(3, 4, 2), PCT012N122 = c(4, 4, 8 + ), PCT012O122 = c(49, 43, 109), PCT012K123 = c(10, 8, 15), + PCT012N123 = c(4, 4, 7), PCT012O123 = c(59, 46, 114), PCT012K124 = c(4, + 3, 10), PCT012N124 = c(5, 3, 4), PCT012O124 = c(46, 36, 93 + ), PCT012I021 = c(663, 764, 2448), PCT012I022 = c(726, 657, + 2843), PCT012I125 = c(636, 721, 2703), PCT012I126 = c(675, + 590, 3024), PCT012J021 = c(518, 214, 1123), PCT012J022 = c(475, + 222, 1064), PCT012J125 = c(538, 194, 1060), PCT012J126 = c(587, + 180, 954), PCT012H021 = c(69, 129, 466), PCT012H022 = c(105, + 147, 479), PCT012H125 = c(104, 100, 462), PCT012H126 = c(103, + 94, 482), PCT012L021 = c(29, 10, 162), PCT012M021 = c(0, + 0, 1), PCT012L022 = c(18, 8, 197), PCT012M022 = c(0, 1, 1 + ), PCT012L125 = c(18, 12, 167), PCT012M125 = c(0, 1, 0), + PCT012L126 = c(11, 13, 182), PCT012M126 = c(0, 0, 0), PCT012K021 = c(7, + 3, 9), PCT012N021 = c(6, 1, 5), PCT012O021 = c(46, 27, 93 + ), PCT012K022 = c(6, 3, 8), PCT012N022 = c(7, 2, 4), PCT012O022 = c(57, + 30, 109), PCT012K125 = c(7, 4, 6), PCT012N125 = c(1, 2, 13 + ), PCT012O125 = c(61, 23, 101), PCT012K126 = c(6, 6, 6), + PCT012N126 = c(2, 0, 9), PCT012O126 = c(65, 35, 118), PCT012I023 = c(678, + 708, 2730), PCT012I127 = c(643, 641, 3067), PCT012J023 = c(445, + 188, 1016), PCT012J127 = c(559, 166, 961), PCT012H023 = c(101, + 163, 508), PCT012H127 = c(77, 99, 457), PCT012L023 = c(12, + 12, 216), PCT012M023 = c(3, 0, 1), PCT012L127 = c(24, 8, + 173), PCT012M127 = c(3, 1, 1), PCT012K023 = c(5, 6, 6), PCT012N023 = c(4, + 3, 9), PCT012O023 = c(49, 13, 96), PCT012K127 = c(4, 6, 8 + ), PCT012N127 = c(3, 1, 8), PCT012O127 = c(32, 18, 109), + PCT012I024 = c(738, 577, 2685), PCT012I128 = c(735, 589, + 2939), PCT012J024 = c(408, 165, 904), PCT012J128 = c(443, + 166, 864), PCT012H024 = c(91, 185, 485), PCT012H128 = c(97, + 124, 429), PCT012L024 = c(15, 11, 138), PCT012M024 = c(1, + 0, 3), PCT012L128 = c(15, 5, 183), PCT012M128 = c(1, 0, 2 + ), PCT012K024 = c(4, 4, 4), PCT012N024 = c(4, 2, 7), PCT012O024 = c(35, + 26, 67), PCT012K128 = c(12, 6, 11), PCT012N128 = c(4, 4, + 12), PCT012O128 = c(49, 20, 103), PCT012I025 = c(716, 597, + 2322), PCT012I026 = c(712, 649, 2159), PCT012I027 = c(694, + 705, 2124), PCT012I129 = c(650, 609, 2423), PCT012I130 = c(655, + 659, 1934), PCT012I131 = c(705, 702, 2046), PCT012J025 = c(385, + 151, 868), PCT012J026 = c(310, 170, 885), PCT012J027 = c(230, + 161, 819), PCT012J129 = c(379, 179, 886), PCT012J130 = c(302, + 172, 835), PCT012J131 = c(244, 193, 929), PCT012H025 = c(94, + 201, 446), PCT012H026 = c(106, 192, 496), PCT012H027 = c(102, + 230, 455), PCT012H129 = c(87, 162, 394), PCT012H130 = c(94, + 155, 358), PCT012H131 = c(79, 165, 404), PCT012L025 = c(16, + 12, 151), PCT012M025 = c(2, 0, 1), PCT012L026 = c(32, 13, + 149), PCT012M026 = c(1, 1, 0), PCT012L027 = c(24, 9, 180), + PCT012M027 = c(0, 1, 3), PCT012L129 = c(19, 7, 151), PCT012M129 = c(0, + 1, 1), PCT012L130 = c(29, 13, 148), PCT012M130 = c(2, 0, + 2), PCT012L131 = c(19, 11, 206), PCT012M131 = c(1, 0, 1), + PCT012K025 = c(10, 2, 4), PCT012N025 = c(2, 5, 7), PCT012O025 = c(38, + 16, 64), PCT012K026 = c(7, 1, 3), PCT012N026 = c(1, 2, 3), + PCT012O026 = c(25, 11, 57), PCT012K027 = c(7, 6, 5), PCT012N027 = c(4, + 1, 10), PCT012O027 = c(27, 14, 56), PCT012K129 = c(5, 7, + 4), PCT012N129 = c(1, 2, 8), PCT012O129 = c(36, 20, 91), + PCT012K130 = c(5, 1, 4), PCT012N130 = c(2, 4, 6), PCT012O130 = c(40, + 22, 58), PCT012K131 = c(5, 7, 8), PCT012N131 = c(2, 0, 8), + PCT012O131 = c(37, 15, 83), PCT012I028 = c(656, 657, 2154 + ), PCT012I029 = c(702, 691, 1974), PCT012I030 = c(646, 684, + 2106), PCT012I031 = c(626, 680, 1989), PCT012I032 = c(647, + 690, 1964), PCT012I132 = c(678, 668, 1981), PCT012I133 = c(682, + 688, 1935), PCT012I134 = c(666, 666, 2019), PCT012I135 = c(671, + 689, 1973), PCT012I136 = c(683, 718, 2106), PCT012J028 = c(229, + 188, 811), PCT012J029 = c(225, 155, 802), PCT012J030 = c(233, + 160, 842), PCT012J031 = c(240, 141, 767), PCT012J032 = c(192, + 162, 797), PCT012J132 = c(266, 164, 867), PCT012J133 = c(241, + 154, 848), PCT012J134 = c(277, 154, 826), PCT012J135 = c(250, + 177, 920), PCT012J136 = c(255, 169, 881), PCT012H028 = c(79, + 260, 483), PCT012H029 = c(82, 246, 470), PCT012H030 = c(102, + 227, 496), PCT012H031 = c(83, 219, 481), PCT012H032 = c(85, + 217, 485), PCT012H132 = c(70, 161, 394), PCT012H133 = c(77, + 163, 424), PCT012H134 = c(96, 174, 440), PCT012H135 = c(83, + 187, 433), PCT012H136 = c(100, 150, 401), PCT012L028 = c(13, + 7, 180), PCT012M028 = c(0, 1, 1), PCT012L029 = c(17, 9, 196 + ), PCT012M029 = c(2, 0, 1), PCT012L030 = c(24, 9, 258), PCT012M030 = c(1, + 0, 1), PCT012L031 = c(15, 12, 260), PCT012M031 = c(0, 0, + 0), PCT012L032 = c(18, 11, 251), PCT012M032 = c(1, 2, 0), + PCT012L132 = c(28, 14, 238), PCT012M132 = c(0, 0, 2), PCT012L133 = c(16, + 15, 258), PCT012M133 = c(1, 1, 0), PCT012L134 = c(22, 8, + 274), PCT012M134 = c(0, 0, 0), PCT012L135 = c(33, 18, 290 + ), PCT012M135 = c(0, 1, 2), PCT012L136 = c(26, 9, 286), PCT012M136 = c(0, + 1, 1), PCT012K028 = c(4, 5, 13), PCT012N028 = c(0, 4, 8), + PCT012O028 = c(32, 17, 46), PCT012K029 = c(7, 3, 5), PCT012N029 = c(1, + 2, 9), PCT012O029 = c(28, 9, 56), PCT012K030 = c(1, 7, 3), + PCT012N030 = c(0, 0, 8), PCT012O030 = c(32, 18, 61), PCT012K031 = c(7, + 2, 5), PCT012N031 = c(0, 6, 5), PCT012O031 = c(24, 14, 53 + ), PCT012K032 = c(3, 1, 6), PCT012N032 = c(2, 1, 11), PCT012O032 = c(25, + 17, 43), PCT012K132 = c(8, 3, 7), PCT012N132 = c(1, 2, 6), + PCT012O132 = c(36, 19, 59), PCT012K133 = c(6, 3, 4), PCT012N133 = c(2, + 2, 6), PCT012O133 = c(24, 15, 64), PCT012K134 = c(4, 1, 7 + ), PCT012N134 = c(6, 2, 6), PCT012O134 = c(30, 22, 72), PCT012K135 = c(9, + 8, 5), PCT012N135 = c(2, 0, 12), PCT012O135 = c(29, 21, 65 + ), PCT012K136 = c(3, 3, 4), PCT012N136 = c(1, 1, 6), PCT012O136 = c(28, + 19, 82), PCT012I033 = c(629, 664, 2026), PCT012I034 = c(621, + 690, 1855), PCT012I035 = c(626, 662, 1912), PCT012I036 = c(513, + 622, 1778), PCT012I037 = c(574, 576, 1682), PCT012I137 = c(667, + 672, 2047), PCT012I138 = c(632, 619, 1907), PCT012I139 = c(667, + 645, 1909), PCT012I140 = c(641, 624, 1780), PCT012I141 = c(619, + 608, 1757), PCT012J033 = c(194, 166, 880), PCT012J034 = c(176, + 148, 831), PCT012J035 = c(178, 130, 788), PCT012J036 = c(140, + 125, 720), PCT012J037 = c(216, 112, 722), PCT012J137 = c(249, + 159, 998), PCT012J138 = c(253, 154, 877), PCT012J139 = c(239, + 167, 853), PCT012J140 = c(232, 120, 844), PCT012J141 = c(229, + 146, 884), PCT012H033 = c(86, 217, 505), PCT012H034 = c(69, + 185, 450), PCT012H035 = c(65, 177, 401), PCT012H036 = c(72, + 168, 445), PCT012H037 = c(74, 183, 383), PCT012H137 = c(79, + 172, 441), PCT012H138 = c(63, 145, 389), PCT012H139 = c(70, + 156, 398), PCT012H140 = c(71, 156, 398), PCT012H141 = c(62, + 148, 377), PCT012L033 = c(16, 8, 231), PCT012M033 = c(0, + 0, 1), PCT012L034 = c(16, 17, 264), PCT012M034 = c(0, 0, + 2), PCT012L035 = c(20, 13, 273), PCT012M035 = c(0, 0, 1), + PCT012L036 = c(20, 19, 257), PCT012M036 = c(0, 0, 0), PCT012L037 = c(19, + 14, 250), PCT012M037 = c(0, 0, 0), PCT012L137 = c(25, 11, + 257), PCT012M137 = c(0, 1, 0), PCT012L138 = c(20, 11, 230 + ), PCT012M138 = c(0, 1, 0), PCT012L139 = c(33, 10, 264), + PCT012M139 = c(2, 0, 0), PCT012L140 = c(30, 12, 235), PCT012M140 = c(1, + 1, 1), PCT012L141 = c(33, 22, 229), PCT012M141 = c(1, 0, + 0), PCT012K033 = c(6, 4, 3), PCT012N033 = c(1, 1, 6), PCT012O033 = c(26, + 14, 64), PCT012K034 = c(6, 4, 2), PCT012N034 = c(3, 2, 5), + PCT012O034 = c(16, 9, 49), PCT012K035 = c(6, 6, 2), PCT012N035 = c(2, + 2, 8), PCT012O035 = c(25, 17, 54), PCT012K036 = c(3, 2, 6 + ), PCT012N036 = c(1, 2, 3), PCT012O036 = c(15, 14, 43), PCT012K037 = c(6, + 2, 10), PCT012N037 = c(0, 0, 5), PCT012O037 = c(12, 15, 50 + ), PCT012K137 = c(5, 4, 6), PCT012N137 = c(1, 1, 7), PCT012O137 = c(22, + 23, 88), PCT012K138 = c(2, 2, 10), PCT012N138 = c(2, 1, 13 + ), PCT012O138 = c(19, 13, 70), PCT012K139 = c(8, 6, 5), PCT012N139 = c(2, + 3, 5), PCT012O139 = c(28, 12, 58), PCT012K140 = c(4, 2, 9 + ), PCT012N140 = c(0, 2, 6), PCT012O140 = c(32, 13, 51), PCT012K141 = c(4, + 5, 5), PCT012N141 = c(3, 2, 9), PCT012O141 = c(21, 13, 55 + ), PCT012I038 = c(568, 602, 1874), PCT012I039 = c(551, 636, + 1751), PCT012I040 = c(584, 682, 1930), PCT012I041 = c(651, + 780, 2031), PCT012I042 = c(694, 823, 2308), PCT012I142 = c(656, + 682, 1822), PCT012I143 = c(564, 690, 1823), PCT012I144 = c(645, + 745, 1862), PCT012I145 = c(748, 784, 2031), PCT012I146 = c(714, + 817, 2338), PCT012J038 = c(195, 131, 792), PCT012J039 = c(171, + 128, 733), PCT012J040 = c(188, 143, 795), PCT012J041 = c(220, + 134, 762), PCT012J042 = c(243, 152, 1001), PCT012J142 = c(232, + 143, 885), PCT012J143 = c(228, 123, 826), PCT012J144 = c(239, + 129, 930), PCT012J145 = c(250, 131, 965), PCT012J146 = c(303, + 173, 1093), PCT012H038 = c(72, 153, 402), PCT012H039 = c(66, + 151, 374), PCT012H040 = c(49, 193, 388), PCT012H041 = c(56, + 152, 369), PCT012H042 = c(54, 118, 333), PCT012H142 = c(66, + 142, 384), PCT012H143 = c(71, 136, 376), PCT012H144 = c(62, + 148, 335), PCT012H145 = c(63, 109, 354), PCT012H146 = c(70, + 122, 335), PCT012L038 = c(20, 13, 240), PCT012M038 = c(0, + 2, 0), PCT012L039 = c(18, 10, 259), PCT012M039 = c(0, 1, + 0), PCT012L040 = c(32, 19, 224), PCT012M040 = c(1, 0, 0), + PCT012L041 = c(29, 16, 228), PCT012M041 = c(0, 2, 1), PCT012L042 = c(23, + 11, 219), PCT012M042 = c(0, 0, 0), PCT012L142 = c(36, 18, + 228), PCT012M142 = c(0, 2, 3), PCT012L143 = c(29, 18, 204 + ), PCT012M143 = c(0, 0, 2), PCT012L144 = c(31, 21, 217), + PCT012M144 = c(0, 0, 0), PCT012L145 = c(31, 15, 247), PCT012M145 = c(1, + 1, 2), PCT012L146 = c(32, 18, 193), PCT012M146 = c(0, 2, + 1), PCT012K038 = c(8, 6, 12), PCT012N038 = c(4, 2, 10), PCT012O038 = c(15, + 12, 44), PCT012K039 = c(7, 4, 13), PCT012N039 = c(0, 1, 6 + ), PCT012O039 = c(21, 6, 47), PCT012K040 = c(4, 3, 3), PCT012N040 = c(0, + 3, 8), PCT012O040 = c(15, 10, 39), PCT012K041 = c(2, 4, 7 + ), PCT012N041 = c(2, 0, 8), PCT012O041 = c(19, 5, 44), PCT012K042 = c(7, + 5, 4), PCT012N042 = c(0, 2, 7), PCT012O042 = c(19, 21, 49 + ), PCT012K142 = c(6, 10, 5), PCT012N142 = c(2, 1, 6), PCT012O142 = c(30, + 11, 63), PCT012K143 = c(2, 2, 10), PCT012N143 = c(1, 0, 5 + ), PCT012O143 = c(21, 6, 48), PCT012K144 = c(10, 7, 8), PCT012N144 = c(1, + 3, 4), PCT012O144 = c(23, 13, 57), PCT012K145 = c(3, 3, 7 + ), PCT012N145 = c(3, 2, 4), PCT012O145 = c(18, 13, 47), PCT012K146 = c(8, + 8, 5), PCT012N146 = c(2, 4, 6), PCT012O146 = c(15, 12, 59 + ), PCT012I043 = c(729, 839, 2167), PCT012I044 = c(710, 771, + 2177), PCT012I045 = c(652, 803, 2141), PCT012I046 = c(747, + 893, 2367), PCT012I047 = c(765, 976, 2383), PCT012I147 = c(764, + 857, 2325), PCT012I148 = c(646, 841, 2257), PCT012I149 = c(700, + 862, 2272), PCT012I150 = c(758, 904, 2360), PCT012I151 = c(789, + 971, 2434), PCT012J043 = c(222, 158, 923), PCT012J044 = c(190, + 154, 834), PCT012J045 = c(242, 159, 835), PCT012J046 = c(225, + 173, 881), PCT012J047 = c(236, 181, 881), PCT012J147 = c(243, + 178, 1078), PCT012J148 = c(276, 160, 1002), PCT012J149 = c(258, + 178, 1144), PCT012J150 = c(270, 200, 1030), PCT012J151 = c(274, + 202, 1146), PCT012H043 = c(52, 149, 369), PCT012H044 = c(42, + 138, 341), PCT012H045 = c(54, 101, 315), PCT012H046 = c(43, + 93, 305), PCT012H047 = c(49, 84, 307), PCT012H147 = c(53, + 112, 349), PCT012H148 = c(59, 83, 302), PCT012H149 = c(66, + 77, 310), PCT012H150 = c(47, 82, 331), PCT012H151 = c(40, + 94, 289), PCT012L043 = c(28, 14, 218), PCT012M043 = c(1, + 0, 1), PCT012L044 = c(26, 17, 180), PCT012M044 = c(0, 0, + 0), PCT012L045 = c(27, 12, 145), PCT012M045 = c(0, 0, 0), + PCT012L046 = c(20, 19, 153), PCT012M046 = c(2, 0, 1), PCT012L047 = c(24, + 21, 162), PCT012M047 = c(0, 0, 0), PCT012L147 = c(44, 20, + 227), PCT012M147 = c(0, 2, 2), PCT012L148 = c(24, 16, 211 + ), PCT012M148 = c(0, 0, 0), PCT012L149 = c(30, 10, 138), + PCT012M149 = c(0, 0, 1), PCT012L150 = c(40, 18, 165), PCT012M150 = c(0, + 0, 1), PCT012L151 = c(31, 22, 147), PCT012M151 = c(2, 0, + 0), PCT012K043 = c(4, 1, 6), PCT012N043 = c(1, 4, 7), PCT012O043 = c(12, + 15, 60), PCT012K044 = c(6, 10, 10), PCT012N044 = c(0, 1, + 3), PCT012O044 = c(13, 10, 44), PCT012K045 = c(7, 3, 10), + PCT012N045 = c(1, 2, 2), PCT012O045 = c(13, 11, 34), PCT012K046 = c(5, + 5, 5), PCT012N046 = c(1, 0, 4), PCT012O046 = c(15, 8, 43), + PCT012K047 = c(5, 12, 3), PCT012N047 = c(1, 2, 3), PCT012O047 = c(13, + 12, 36), PCT012K147 = c(7, 3, 5), PCT012N147 = c(3, 1, 4), + PCT012O147 = c(27, 9, 54), PCT012K148 = c(6, 5, 5), PCT012N148 = c(0, + 0, 8), PCT012O148 = c(14, 8, 59), PCT012K149 = c(7, 3, 6), + PCT012N149 = c(2, 1, 5), PCT012O149 = c(22, 5, 57), PCT012K150 = c(10, + 6, 10), PCT012N150 = c(0, 2, 4), PCT012O150 = c(22, 11, 39 + ), PCT012K151 = c(14, 10, 6), PCT012N151 = c(1, 3, 7), PCT012O151 = c(14, + 12, 47), PCT012I048 = c(811, 1065, 2584), PCT012I049 = c(842, + 1062, 2617), PCT012I050 = c(811, 1067, 2680), PCT012I051 = c(781, + 1062, 2780), PCT012I052 = c(863, 1121, 2652), PCT012I152 = c(872, + 1023, 2626), PCT012I153 = c(841, 1073, 2685), PCT012I154 = c(870, + 1150, 2751), PCT012I155 = c(858, 1129, 2779), PCT012I156 = c(839, + 1261, 2890), PCT012J048 = c(219, 170, 929), PCT012J049 = c(271, + 198, 942), PCT012J050 = c(269, 173, 898), PCT012J051 = c(237, + 185, 923), PCT012J052 = c(254, 180, 820), PCT012J152 = c(290, + 202, 1094), PCT012J153 = c(274, 174, 1097), PCT012J154 = c(292, + 205, 1117), PCT012J155 = c(254, 177, 1087), PCT012J156 = c(275, + 200, 997), PCT012H048 = c(40, 74, 281), PCT012H049 = c(45, + 77, 297), PCT012H050 = c(52, 66, 250), PCT012H051 = c(42, + 76, 252), PCT012H052 = c(50, 60, 210), PCT012H152 = c(61, + 80, 287), PCT012H153 = c(42, 67, 226), PCT012H154 = c(41, + 62, 223), PCT012H155 = c(48, 49, 199), PCT012H156 = c(53, + 53, 255), PCT012L048 = c(26, 14, 166), PCT012M048 = c(0, + 0, 1), PCT012L049 = c(24, 11, 155), PCT012M049 = c(0, 0, + 1), PCT012L050 = c(28, 12, 156), PCT012M050 = c(0, 1, 0), + PCT012L051 = c(16, 17, 147), PCT012M051 = c(0, 1, 0), PCT012L052 = c(12, + 14, 122), PCT012M052 = c(0, 0, 1), PCT012L152 = c(41, 26, + 161), PCT012M152 = c(1, 0, 0), PCT012L153 = c(33, 19, 162 + ), PCT012M153 = c(3, 1, 3), PCT012L154 = c(37, 16, 172), + PCT012M154 = c(0, 0, 1), PCT012L155 = c(35, 20, 148), PCT012M155 = c(0, + 0, 0), PCT012L156 = c(31, 17, 156), PCT012M156 = c(0, 0, + 1), PCT012K048 = c(3, 12, 8), PCT012N048 = c(2, 1, 2), PCT012O048 = c(14, + 15, 31), PCT012K049 = c(9, 17, 13), PCT012N049 = c(2, 1, + 3), PCT012O049 = c(10, 10, 48), PCT012K050 = c(8, 8, 8), + PCT012N050 = c(1, 0, 6), PCT012O050 = c(13, 11, 45), PCT012K051 = c(10, + 17, 5), PCT012N051 = c(2, 0, 4), PCT012O051 = c(16, 18, 39 + ), PCT012K052 = c(7, 9, 8), PCT012N052 = c(2, 2, 3), PCT012O052 = c(11, + 12, 35), PCT012K152 = c(9, 10, 6), PCT012N152 = c(4, 1, 4 + ), PCT012O152 = c(17, 14, 68), PCT012K153 = c(6, 9, 9), PCT012N153 = c(1, + 5, 3), PCT012O153 = c(16, 17, 39), PCT012K154 = c(9, 6, 17 + ), PCT012N154 = c(4, 0, 3), PCT012O154 = c(18, 12, 51), PCT012K155 = c(3, + 8, 7), PCT012N155 = c(0, 0, 6), PCT012O155 = c(16, 14, 50 + ), PCT012K156 = c(6, 8, 9), PCT012N156 = c(3, 4, 5), PCT012O156 = c(22, + 19, 42), PCT012I053 = c(797, 1113, 2702), PCT012I054 = c(813, + 1157, 2809), PCT012I055 = c(759, 1121, 2835), PCT012I056 = c(735, + 1084, 2718), PCT012I057 = c(727, 1143, 2661), PCT012I157 = c(869, + 1143, 2958), PCT012I158 = c(807, 1185, 2836), PCT012I159 = c(829, + 1290, 2927), PCT012I160 = c(829, 1284, 2767), PCT012I161 = c(785, + 1242, 2803), PCT012J053 = c(239, 179, 851), PCT012J054 = c(204, + 162, 832), PCT012J055 = c(204, 201, 803), PCT012J056 = c(215, + 141, 777), PCT012J057 = c(201, 182, 749), PCT012J157 = c(256, + 222, 1080), PCT012J158 = c(258, 216, 946), PCT012J159 = c(275, + 199, 997), PCT012J160 = c(258, 173, 970), PCT012J161 = c(248, + 163, 890), PCT012H053 = c(35, 75, 211), PCT012H054 = c(41, + 53, 207), PCT012H055 = c(28, 39, 193), PCT012H056 = c(35, + 43, 189), PCT012H057 = c(33, 49, 157), PCT012H157 = c(33, + 52, 173), PCT012H158 = c(39, 43, 156), PCT012H159 = c(38, + 43, 173), PCT012H160 = c(38, 34, 157), PCT012H161 = c(45, + 40, 157), PCT012L053 = c(17, 18, 147), PCT012M053 = c(1, + 0, 1), PCT012L054 = c(18, 14, 119), PCT012M054 = c(1, 0, + 0), PCT012L055 = c(18, 15, 126), PCT012M055 = c(1, 2, 0), + PCT012L056 = c(13, 12, 113), PCT012M056 = c(1, 0, 1), PCT012L057 = c(15, + 5, 115), PCT012M057 = c(2, 1, 0), PCT012L157 = c(25, 15, + 140), PCT012M157 = c(1, 1, 2), PCT012L158 = c(34, 19, 121 + ), PCT012M158 = c(3, 1, 0), PCT012L159 = c(19, 13, 113), + PCT012M159 = c(1, 0, 0), PCT012L160 = c(28, 15, 143), PCT012M160 = c(3, + 1, 0), PCT012L161 = c(28, 16, 126), PCT012M161 = c(0, 0, + 0), PCT012K053 = c(10, 13, 9), PCT012N053 = c(1, 1, 4), PCT012O053 = c(10, + 7, 27), PCT012K054 = c(3, 13, 8), PCT012N054 = c(1, 0, 5), + PCT012O054 = c(20, 11, 30), PCT012K055 = c(8, 9, 7), PCT012N055 = c(1, + 0, 6), PCT012O055 = c(14, 13, 26), PCT012K056 = c(5, 13, + 9), PCT012N056 = c(0, 0, 6), PCT012O056 = c(10, 13, 29), + PCT012K057 = c(17, 5, 11), PCT012N057 = c(1, 0, 5), PCT012O057 = c(13, + 11, 30), PCT012K157 = c(12, 6, 15), PCT012N157 = c(1, 1, + 4), PCT012O157 = c(17, 15, 45), PCT012K158 = c(7, 6, 10), + PCT012N158 = c(0, 1, 6), PCT012O158 = c(17, 7, 34), PCT012K159 = c(9, + 9, 11), PCT012N159 = c(1, 2, 3), PCT012O159 = c(18, 11, 31 + ), PCT012K160 = c(7, 7, 10), PCT012N160 = c(3, 1, 4), PCT012O160 = c(16, + 13, 37), PCT012K161 = c(10, 3, 7), PCT012N161 = c(3, 2, 6 + ), PCT012O161 = c(12, 8, 43), PCT012I058 = c(673, 1133, 2554 + ), PCT012I059 = c(682, 1169, 2356), PCT012I060 = c(622, 1131, + 2476), PCT012I061 = c(628, 1114, 2265), PCT012I062 = c(593, + 1166, 2124), PCT012I162 = c(766, 1267, 2624), PCT012I163 = c(717, + 1287, 2546), PCT012I164 = c(716, 1360, 2612), PCT012I165 = c(711, + 1345, 2464), PCT012I166 = c(673, 1311, 2374), PCT012J058 = c(181, + 139, 694), PCT012J059 = c(188, 147, 602), PCT012J060 = c(187, + 162, 572), PCT012J061 = c(170, 116, 510), PCT012J062 = c(176, + 132, 591), PCT012J162 = c(247, 148, 881), PCT012J163 = c(241, + 171, 767), PCT012J164 = c(228, 165, 789), PCT012J165 = c(230, + 172, 761), PCT012J166 = c(199, 141, 708), PCT012H058 = c(36, + 40, 168), PCT012H059 = c(18, 47, 152), PCT012H060 = c(30, + 32, 141), PCT012H061 = c(21, 37, 105), PCT012H062 = c(24, + 24, 113), PCT012H162 = c(25, 32, 139), PCT012H163 = c(33, + 37, 142), PCT012H164 = c(18, 33, 117), PCT012H165 = c(12, + 26, 119), PCT012H166 = c(25, 33, 109), PCT012L058 = c(13, + 7, 117), PCT012M058 = c(0, 0, 2), PCT012L059 = c(7, 14, 98 + ), PCT012M059 = c(0, 0, 1), PCT012L060 = c(14, 9, 97), PCT012M060 = c(1, + 1, 0), PCT012L061 = c(7, 6, 88), PCT012M061 = c(0, 0, 0), + PCT012L062 = c(10, 7, 88), PCT012M062 = c(0, 0, 0), PCT012L162 = c(35, + 18, 131), PCT012M162 = c(0, 0, 2), PCT012L163 = c(26, 13, + 101), PCT012M163 = c(1, 0, 1), PCT012L164 = c(20, 13, 100 + ), PCT012M164 = c(0, 0, 0), PCT012L165 = c(31, 6, 100), PCT012M165 = c(1, + 0, 0), PCT012L166 = c(21, 10, 115), PCT012M166 = c(0, 0, + 0), PCT012K058 = c(9, 12, 9), PCT012N058 = c(0, 1, 1), PCT012O058 = c(11, + 7, 15), PCT012K059 = c(7, 6, 8), PCT012N059 = c(1, 0, 6), + PCT012O059 = c(10, 7, 24), PCT012K060 = c(7, 4, 7), PCT012N060 = c(0, + 0, 2), PCT012O060 = c(6, 9, 24), PCT012K061 = c(10, 9, 12 + ), PCT012N061 = c(0, 3, 4), PCT012O061 = c(8, 8, 19), PCT012K062 = c(2, + 12, 4), PCT012N062 = c(0, 1, 3), PCT012O062 = c(4, 6, 21), + PCT012K162 = c(7, 14, 6), PCT012N162 = c(1, 2, 1), PCT012O162 = c(9, + 9, 37), PCT012K163 = c(11, 5, 7), PCT012N163 = c(1, 0, 2), + PCT012O163 = c(14, 12, 34), PCT012K164 = c(7, 7, 11), PCT012N164 = c(0, + 2, 5), PCT012O164 = c(9, 7, 35), PCT012K165 = c(7, 6, 2), + PCT012N165 = c(0, 0, 2), PCT012O165 = c(12, 12, 30), PCT012K166 = c(12, + 7, 8), PCT012N166 = c(0, 1, 4), PCT012O166 = c(6, 8, 31), + PCT012I063 = c(645, 1264, 2040), PCT012I064 = c(628, 1287, + 2088), PCT012I167 = c(638, 1364, 2305), PCT012I168 = c(744, + 1466, 2204), PCT012J063 = c(158, 107, 593), PCT012J064 = c(155, + 119, 508), PCT012J167 = c(215, 137, 680), PCT012J168 = c(202, + 137, 685), PCT012H063 = c(19, 29, 108), PCT012H064 = c(23, + 30, 94), PCT012H167 = c(32, 29, 99), PCT012H168 = c(29, 23, + 101), PCT012L063 = c(14, 11, 88), PCT012M063 = c(0, 1, 0), + PCT012L064 = c(13, 9, 75), PCT012M064 = c(0, 0, 0), PCT012L167 = c(22, + 11, 98), PCT012M167 = c(0, 1, 0), PCT012L168 = c(27, 16, + 83), PCT012M168 = c(0, 2, 2), PCT012K063 = c(9, 7, 1), PCT012N063 = c(0, + 0, 1), PCT012O063 = c(6, 5, 22), PCT012K064 = c(4, 4, 3), + PCT012N064 = c(2, 1, 3), PCT012O064 = c(10, 13, 18), PCT012K167 = c(8, + 6, 4), PCT012N167 = c(1, 0, 4), PCT012O167 = c(6, 8, 24), + PCT012K168 = c(8, 8, 4), PCT012N168 = c(0, 0, 2), PCT012O168 = c(7, + 9, 22), PCT012I065 = c(698, 1424, 2125), PCT012I066 = c(716, + 1481, 2257), PCT012I067 = c(511, 1173, 1444), PCT012I169 = c(746, + 1674, 2317), PCT012I170 = c(770, 1643, 2315), PCT012I171 = c(592, + 1279, 1746), PCT012J065 = c(178, 112, 489), PCT012J066 = c(162, + 110, 460), PCT012J067 = c(108, 93, 372), PCT012J169 = c(181, + 147, 682), PCT012J170 = c(187, 105, 585), PCT012J171 = c(178, + 96, 492), PCT012H065 = c(24, 16, 85), PCT012H066 = c(16, + 29, 84), PCT012H067 = c(11, 26, 64), PCT012H169 = c(28, 21, + 100), PCT012H170 = c(26, 17, 90), PCT012H171 = c(25, 30, + 82), PCT012L065 = c(12, 6, 77), PCT012M065 = c(0, 0, 1), + PCT012L066 = c(13, 14, 69), PCT012M066 = c(0, 0, 0), PCT012L067 = c(4, + 9, 70), PCT012M067 = c(0, 0, 1), PCT012L169 = c(26, 14, 109 + ), PCT012M169 = c(0, 0, 1), PCT012L170 = c(24, 6, 97), PCT012M170 = c(0, + 1, 1), PCT012L171 = c(24, 8, 83), PCT012M171 = c(0, 0, 0), + PCT012K065 = c(9, 5, 7), PCT012N065 = c(0, 1, 2), PCT012O065 = c(14, + 6, 13), PCT012K066 = c(7, 13, 3), PCT012N066 = c(1, 0, 2), + PCT012O066 = c(8, 9, 19), PCT012K067 = c(3, 6, 4), PCT012N067 = c(3, + 0, 2), PCT012O067 = c(2, 6, 12), PCT012K169 = c(10, 8, 5), + PCT012N169 = c(0, 0, 1), PCT012O169 = c(9, 16, 19), PCT012K170 = c(3, + 3, 4), PCT012N170 = c(0, 0, 1), PCT012O170 = c(6, 10, 34), + PCT012K171 = c(6, 9, 4), PCT012N171 = c(1, 0, 1), PCT012O171 = c(7, + 6, 13), PCT012I068 = c(544, 1217, 1557), PCT012I069 = c(579, + 1275, 1549), PCT012I172 = c(594, 1321, 1687), PCT012I173 = c(653, + 1407, 1759), PCT012J068 = c(110, 75, 366), PCT012J069 = c(118, + 84, 347), PCT012J172 = c(152, 118, 456), PCT012J173 = c(141, + 90, 471), PCT012H068 = c(7, 24, 71), PCT012H069 = c(10, 18, + 40), PCT012H172 = c(16, 14, 81), PCT012H173 = c(15, 19, 59 + ), PCT012L068 = c(10, 5, 67), PCT012M068 = c(0, 0, 0), PCT012L069 = c(14, + 2, 49), PCT012M069 = c(0, 0, 0), PCT012L172 = c(30, 11, 80 + ), PCT012M172 = c(0, 0, 0), PCT012L173 = c(20, 17, 64), PCT012M173 = c(0, + 0, 0), PCT012K068 = c(5, 3, 3), PCT012N068 = c(1, 0, 3), + PCT012O068 = c(9, 4, 13), PCT012K069 = c(3, 9, 3), PCT012N069 = c(0, + 0, 2), PCT012O069 = c(9, 10, 17), PCT012K172 = c(0, 3, 3), + PCT012N172 = c(1, 1, 1), PCT012O172 = c(8, 4, 12), PCT012K173 = c(3, + 1, 3), PCT012N173 = c(0, 0, 2), PCT012O173 = c(7, 8, 7), + PCT012I070 = c(598, 1359, 1497), PCT012I071 = c(504, 1113, + 1327), PCT012I072 = c(405, 1060, 1252), PCT012I174 = c(611, + 1472, 1837), PCT012I175 = c(562, 1238, 1610), PCT012I176 = c(533, + 1186, 1429), PCT012J070 = c(136, 80, 345), PCT012J071 = c(105, + 92, 300), PCT012J072 = c(88, 83, 273), PCT012J174 = c(142, + 95, 468), PCT012J175 = c(130, 101, 378), PCT012J176 = c(122, + 72, 358), PCT012H070 = c(13, 17, 47), PCT012H071 = c(7, 17, + 55), PCT012H072 = c(18, 16, 48), PCT012H174 = c(14, 14, 56 + ), PCT012H175 = c(27, 18, 53), PCT012H176 = c(14, 17, 60), + PCT012L070 = c(7, 6, 79), PCT012M070 = c(1, 0, 0), PCT012L071 = c(10, + 11, 55), PCT012M071 = c(0, 0, 0), PCT012L072 = c(8, 9, 66 + ), PCT012M072 = c(0, 0, 0), PCT012L174 = c(20, 4, 49), PCT012M174 = c(0, + 0, 0), PCT012L175 = c(26, 10, 75), PCT012M175 = c(0, 0, 0 + ), PCT012L176 = c(19, 6, 54), PCT012M176 = c(0, 1, 0), PCT012K070 = c(4, + 5, 5), PCT012N070 = c(0, 1, 0), PCT012O070 = c(4, 10, 14), + PCT012K071 = c(4, 6, 3), PCT012N071 = c(0, 0, 1), PCT012O071 = c(3, + 5, 12), PCT012K072 = c(5, 8, 0), PCT012N072 = c(0, 0, 1), + PCT012O072 = c(3, 4, 9), PCT012K174 = c(3, 5, 5), PCT012N174 = c(1, + 0, 0), PCT012O174 = c(12, 8, 17), PCT012K175 = c(3, 6, 3), + PCT012N175 = c(0, 1, 1), PCT012O175 = c(10, 11, 15), PCT012K176 = c(4, + 6, 5), PCT012N176 = c(0, 1, 1), PCT012O176 = c(9, 7, 9), + PCT012I073 = c(430, 987, 1065), PCT012I074 = c(422, 1015, + 1134), PCT012I075 = c(399, 932, 1057), PCT012I076 = c(363, + 809, 904), PCT012I077 = c(309, 809, 918), PCT012I177 = c(488, + 1087, 1294), PCT012I178 = c(442, 1063, 1320), PCT012I179 = c(495, + 972, 1243), PCT012I180 = c(433, 896, 1164), PCT012I181 = c(427, + 865, 1146), PCT012J073 = c(106, 72, 272), PCT012J074 = c(97, + 47, 236), PCT012J075 = c(91, 60, 203), PCT012J076 = c(74, + 44, 214), PCT012J077 = c(77, 38, 188), PCT012J177 = c(110, + 82, 334), PCT012J178 = c(105, 81, 320), PCT012J179 = c(109, + 68, 312), PCT012J180 = c(92, 51, 267), PCT012J181 = c(79, + 46, 292), PCT012H073 = c(13, 19, 36), PCT012H074 = c(8, 16, + 41), PCT012H075 = c(12, 11, 45), PCT012H076 = c(16, 12, 35 + ), PCT012H077 = c(5, 13, 27), PCT012H177 = c(13, 14, 52), + PCT012H178 = c(17, 13, 44), PCT012H179 = c(11, 18, 44), PCT012H180 = c(18, + 8, 30), PCT012H181 = c(6, 13, 34), PCT012L073 = c(10, 12, + 46), PCT012M073 = c(0, 0, 0), PCT012L074 = c(6, 8, 54), PCT012M074 = c(0, + 0, 0), PCT012L075 = c(10, 5, 48), PCT012M075 = c(0, 0, 0), + PCT012L076 = c(8, 6, 36), PCT012M076 = c(0, 0, 0), PCT012L077 = c(6, + 3, 50), PCT012M077 = c(0, 0, 0), PCT012L177 = c(16, 7, 58 + ), PCT012M177 = c(0, 0, 0), PCT012L178 = c(10, 7, 45), PCT012M178 = c(0, + 0, 0), PCT012L179 = c(14, 3, 58), PCT012M179 = c(0, 0, 0), + PCT012L180 = c(11, 5, 46), PCT012M180 = c(0, 0, 0), PCT012L181 = c(10, + 6, 37), PCT012M181 = c(0, 0, 0), PCT012K073 = c(5, 4, 3), + PCT012N073 = c(1, 0, 0), PCT012O073 = c(4, 6, 9), PCT012K074 = c(4, + 3, 3), PCT012N074 = c(0, 0, 0), PCT012O074 = c(5, 6, 6), + PCT012K075 = c(5, 8, 2), PCT012N075 = c(1, 0, 0), PCT012O075 = c(4, + 8, 11), PCT012K076 = c(0, 3, 1), PCT012N076 = c(0, 0, 1), + PCT012O076 = c(6, 5, 8), PCT012K077 = c(2, 9, 3), PCT012N077 = c(1, + 0, 0), PCT012O077 = c(2, 7, 4), PCT012K177 = c(3, 4, 3), + PCT012N177 = c(3, 0, 3), PCT012O177 = c(5, 4, 6), PCT012K178 = c(4, + 5, 2), PCT012N178 = c(0, 0, 0), PCT012O178 = c(7, 4, 8), + PCT012K179 = c(5, 10, 5), PCT012N179 = c(2, 0, 1), PCT012O179 = c(3, + 5, 8), PCT012K180 = c(3, 6, 2), PCT012N180 = c(0, 1, 0), + PCT012O180 = c(3, 8, 5), PCT012K181 = c(5, 6, 0), PCT012N181 = c(0, + 0, 1), PCT012O181 = c(4, 3, 7), PCT012I078 = c(360, 781, + 916), PCT012I079 = c(320, 660, 816), PCT012I080 = c(290, + 660, 812), PCT012I081 = c(300, 640, 812), PCT012I082 = c(293, + 606, 791), PCT012I182 = c(407, 800, 1158), PCT012I183 = c(342, + 720, 1078), PCT012I184 = c(366, 733, 1114), PCT012I185 = c(349, + 702, 1083), PCT012I186 = c(314, 606, 1159), PCT012J078 = c(75, + 42, 196), PCT012J079 = c(62, 34, 151), PCT012J080 = c(50, + 48, 133), PCT012J081 = c(48, 30, 129), PCT012J082 = c(50, + 31, 111), PCT012J182 = c(92, 60, 260), PCT012J183 = c(81, + 48, 197), PCT012J184 = c(85, 54, 203), PCT012J185 = c(68, + 63, 197), PCT012J186 = c(66, 55, 207), PCT012H078 = c(10, + 5, 31), PCT012H079 = c(6, 11, 19), PCT012H080 = c(5, 10, + 17), PCT012H081 = c(2, 10, 25), PCT012H082 = c(4, 3, 22), + PCT012H182 = c(10, 9, 36), PCT012H183 = c(7, 15, 38), PCT012H184 = c(9, + 4, 28), PCT012H185 = c(8, 8, 21), PCT012H186 = c(10, 4, 19 + ), PCT012L078 = c(4, 5, 29), PCT012M078 = c(0, 1, 0), PCT012L079 = c(3, + 3, 26), PCT012M079 = c(0, 0, 0), PCT012L080 = c(2, 6, 31), + PCT012M080 = c(0, 0, 0), PCT012L081 = c(2, 2, 14), PCT012M081 = c(0, + 0, 0), PCT012L082 = c(3, 1, 22), PCT012M082 = c(0, 0, 0), + PCT012L182 = c(16, 5, 33), PCT012M182 = c(0, 0, 0), PCT012L183 = c(10, + 4, 31), PCT012M183 = c(1, 1, 0), PCT012L184 = c(11, 5, 17 + ), PCT012M184 = c(0, 0, 2), PCT012L185 = c(10, 3, 26), PCT012M185 = c(0, + 0, 0), PCT012L186 = c(9, 4, 20), PCT012M186 = c(0, 0, 0), + PCT012K078 = c(2, 2, 2), PCT012N078 = c(0, 0, 1), PCT012O078 = c(1, + 6, 2), PCT012K079 = c(0, 2, 1), PCT012N079 = c(1, 0, 0), + PCT012O079 = c(3, 6, 3), PCT012K080 = c(5, 2, 0), PCT012N080 = c(0, + 1, 2), PCT012O080 = c(1, 1, 5), PCT012K081 = c(4, 2, 3), + PCT012N081 = c(0, 0, 0), PCT012O081 = c(0, 1, 7), PCT012K082 = c(2, + 2, 1), PCT012N082 = c(0, 0, 0), PCT012O082 = c(0, 1, 3), + PCT012K182 = c(2, 3, 2), PCT012N182 = c(1, 1, 0), PCT012O182 = c(7, + 7, 10), PCT012K183 = c(1, 4, 2), PCT012N183 = c(0, 0, 0), + PCT012O183 = c(4, 1, 12), PCT012K184 = c(3, 7, 1), PCT012N184 = c(0, + 2, 1), PCT012O184 = c(4, 6, 4), PCT012K185 = c(2, 3, 1), + PCT012N185 = c(1, 0, 0), PCT012O185 = c(1, 4, 3), PCT012K186 = c(3, + 6, 2), PCT012N186 = c(0, 0, 0), PCT012O186 = c(6, 4, 7), + PCT012I083 = c(235, 538, 748), PCT012I084 = c(218, 491, 680 + ), PCT012I085 = c(187, 457, 660), PCT012I086 = c(147, 333, + 588), PCT012I087 = c(142, 304, 528), PCT012I187 = c(288, + 554, 1052), PCT012I188 = c(293, 527, 971), PCT012I189 = c(243, + 548, 951), PCT012I190 = c(225, 514, 935), PCT012I191 = c(210, + 397, 824), PCT012J083 = c(26, 38, 108), PCT012J084 = c(27, + 23, 90), PCT012J085 = c(39, 23, 86), PCT012J086 = c(31, 17, + 65), PCT012J087 = c(18, 19, 46), PCT012J187 = c(42, 49, 168 + ), PCT012J188 = c(55, 51, 181), PCT012J189 = c(29, 26, 147 + ), PCT012J190 = c(38, 44, 136), PCT012J191 = c(36, 41, 105 + ), PCT012H083 = c(5, 3, 24), PCT012H084 = c(5, 6, 23), PCT012H085 = c(2, + 3, 15), PCT012H086 = c(6, 3, 12), PCT012H087 = c(4, 7, 10 + ), PCT012H187 = c(10, 7, 18), PCT012H188 = c(9, 8, 15), PCT012H189 = c(4, + 7, 16), PCT012H190 = c(9, 4, 24), PCT012H191 = c(7, 3, 15 + ), PCT012L083 = c(3, 3, 10), PCT012M083 = c(1, 1, 0), PCT012L084 = c(3, + 2, 19), PCT012M084 = c(0, 2, 0), PCT012L085 = c(1, 2, 6), + PCT012M085 = c(0, 0, 0), PCT012L086 = c(1, 0, 6), PCT012M086 = c(0, + 0, 0), PCT012L087 = c(3, 1, 12), PCT012M087 = c(0, 0, 0), + PCT012L187 = c(7, 1, 18), PCT012M187 = c(0, 0, 0), PCT012L188 = c(5, + 2, 18), PCT012M188 = c(0, 0, 0), PCT012L189 = c(6, 0, 14), + PCT012M189 = c(1, 0, 0), PCT012L190 = c(4, 6, 7), PCT012M190 = c(0, + 0, 0), PCT012L191 = c(4, 1, 15), PCT012M191 = c(0, 0, 0), + PCT012K083 = c(3, 4, 3), PCT012N083 = c(0, 0, 1), PCT012O083 = c(2, + 2, 6), PCT012K084 = c(2, 1, 1), PCT012N084 = c(0, 0, 1), + PCT012O084 = c(3, 0, 9), PCT012K085 = c(2, 1, 0), PCT012N085 = c(0, + 0, 0), PCT012O085 = c(0, 3, 3), PCT012K086 = c(0, 0, 0), + PCT012N086 = c(0, 0, 0), PCT012O086 = c(1, 4, 4), PCT012K087 = c(1, + 2, 1), PCT012N087 = c(0, 0, 0), PCT012O087 = c(0, 0, 1), + PCT012K187 = c(3, 3, 1), PCT012N187 = c(0, 0, 0), PCT012O187 = c(1, + 2, 4), PCT012K188 = c(6, 3, 1), PCT012N188 = c(0, 0, 0), + PCT012O188 = c(2, 1, 4), PCT012K189 = c(4, 3, 0), PCT012N189 = c(1, + 0, 0), PCT012O189 = c(2, 3, 6), PCT012K190 = c(1, 0, 0), + PCT012N190 = c(0, 0, 1), PCT012O190 = c(1, 4, 5), PCT012K191 = c(1, + 0, 2), PCT012N191 = c(0, 0, 1), PCT012O191 = c(1, 1, 5), + PCT012I088 = c(122, 282, 483), PCT012I089 = c(102, 186, 452 + ), PCT012I090 = c(80, 161, 353), PCT012I091 = c(77, 154, + 321), PCT012I092 = c(63, 117, 259), PCT012I093 = c(51, 104, + 198), PCT012I094 = c(36, 61, 159), PCT012I095 = c(36, 54, + 110), PCT012I096 = c(17, 41, 80), PCT012I097 = c(9, 15, 59 + ), PCT012I098 = c(8, 16, 45), PCT012I099 = c(5, 10, 29), + PCT012I100 = c(7, 7, 18), PCT012I101 = c(3, 6, 14), PCT012I102 = c(2, + 1, 5), PCT012I103 = c(2, 6, 5), PCT012I104 = c(0, 1, 1), + PCT012I105 = c(0, 0, 0), PCT012I192 = c(213, 384, 806), PCT012I193 = c(215, + 380, 792), PCT012I194 = c(160, 336, 687), PCT012I195 = c(164, + 278, 606), PCT012I196 = c(135, 234, 518), PCT012I197 = c(93, + 189, 404), PCT012I198 = c(77, 181, 358), PCT012I199 = c(69, + 120, 291), PCT012I200 = c(57, 115, 224), PCT012I201 = c(42, + 89, 165), PCT012I202 = c(39, 54, 148), PCT012I203 = c(24, + 54, 116), PCT012I204 = c(27, 33, 75), PCT012I205 = c(13, + 25, 49), PCT012I206 = c(18, 14, 25), PCT012I207 = c(6, 32, + 53), PCT012I208 = c(0, 0, 1), PCT012I209 = c(0, 0, 0), PCT012J088 = c(23, + 29, 59), PCT012J089 = c(19, 19, 47), PCT012J090 = c(17, 17, + 32), PCT012J091 = c(12, 7, 34), PCT012J092 = c(7, 8, 24), + PCT012J093 = c(11, 7, 22), PCT012J094 = c(10, 3, 23), PCT012J095 = c(2, + 8, 13), PCT012J096 = c(3, 4, 16), PCT012J097 = c(0, 6, 3), + PCT012J098 = c(0, 2, 4), PCT012J099 = c(0, 0, 5), PCT012J100 = c(0, + 0, 4), PCT012J101 = c(2, 0, 2), PCT012J102 = c(0, 2, 0), + PCT012J103 = c(1, 2, 4), PCT012J104 = c(0, 0, 0), PCT012J105 = c(0, + 0, 0), PCT012J192 = c(42, 31, 96), PCT012J193 = c(31, 38, + 98), PCT012J194 = c(20, 25, 106), PCT012J195 = c(33, 18, + 57), PCT012J196 = c(19, 22, 69), PCT012J197 = c(21, 11, 54 + ), PCT012J198 = c(13, 12, 36), PCT012J199 = c(11, 5, 34), + PCT012J200 = c(13, 14, 28), PCT012J201 = c(5, 7, 23), PCT012J202 = c(10, + 5, 23), PCT012J203 = c(4, 4, 20), PCT012J204 = c(4, 6, 13 + ), PCT012J205 = c(3, 2, 7), PCT012J206 = c(2, 0, 9), PCT012J207 = c(4, + 4, 16), PCT012J208 = c(0, 0, 2), PCT012J209 = c(0, 0, 0), + PCT012H088 = c(3, 2, 15), PCT012H089 = c(5, 4, 5), PCT012H090 = c(3, + 3, 8), PCT012H091 = c(4, 0, 5), PCT012H092 = c(2, 0, 4), + PCT012H093 = c(1, 1, 4), PCT012H094 = c(0, 2, 3), PCT012H095 = c(1, + 1, 4), PCT012H096 = c(0, 0, 0), PCT012H097 = c(0, 0, 1), + PCT012H098 = c(0, 0, 0), PCT012H099 = c(0, 0, 1), PCT012H100 = c(0, + 0, 0), PCT012H101 = c(0, 0, 0), PCT012H102 = c(0, 0, 0), + PCT012H103 = c(0, 0, 1), PCT012H104 = c(0, 0, 0), PCT012H105 = c(0, + 1, 0), PCT012H192 = c(3, 3, 15), PCT012H193 = c(7, 4, 10), + PCT012H194 = c(1, 0, 11), PCT012H195 = c(2, 5, 14), PCT012H196 = c(1, + 6, 8), PCT012H197 = c(1, 2, 9), PCT012H198 = c(0, 2, 2), + PCT012H199 = c(4, 0, 4), PCT012H200 = c(2, 2, 2), PCT012H201 = c(1, + 1, 3), PCT012H202 = c(1, 0, 1), PCT012H203 = c(2, 1, 0), + PCT012H204 = c(1, 0, 1), PCT012H205 = c(2, 0, 1), PCT012H206 = c(2, + 1, 0), PCT012H207 = c(0, 0, 2), PCT012H208 = c(0, 0, 0), + PCT012H209 = c(0, 0, 0), PCT012L088 = c(1, 0, 9), PCT012M088 = c(0, + 0, 0), PCT012L089 = c(1, 1, 7), PCT012M089 = c(1, 0, 0), + PCT012L090 = c(1, 0, 5), PCT012M090 = c(0, 0, 0), PCT012L091 = c(1, + 0, 6), PCT012M091 = c(0, 0, 0), PCT012L092 = c(0, 1, 2), + PCT012M092 = c(0, 0, 0), PCT012L093 = c(0, 0, 3), PCT012M093 = c(0, + 0, 0), PCT012L094 = c(0, 0, 3), PCT012M094 = c(0, 0, 0), + PCT012L095 = c(0, 0, 1), PCT012M095 = c(0, 0, 0), PCT012L096 = c(0, + 0, 1), PCT012M096 = c(0, 0, 0), PCT012L097 = c(0, 0, 1), + PCT012M097 = c(0, 0, 0), PCT012L098 = c(0, 0, 1), PCT012M098 = c(0, + 0, 0), PCT012L099 = c(1, 0, 1), PCT012M099 = c(0, 0, 0), + PCT012L100 = c(0, 0, 0), PCT012M100 = c(0, 0, 0), PCT012L101 = c(0, + 0, 0), PCT012M101 = c(0, 0, 0), PCT012L102 = c(0, 0, 0), + PCT012M102 = c(0, 0, 0), PCT012L103 = c(0, 0, 0), PCT012M103 = c(0, + 0, 0), PCT012L104 = c(0, 0, 0), PCT012M104 = c(0, 0, 0), + PCT012L105 = c(0, 0, 0), PCT012M105 = c(0, 0, 0), PCT012L192 = c(2, + 0, 16), PCT012M192 = c(0, 0, 0), PCT012L193 = c(5, 0, 10), + PCT012M193 = c(0, 0, 0), PCT012L194 = c(4, 2, 2), PCT012M194 = c(0, + 0, 0), PCT012L195 = c(0, 1, 11), PCT012M195 = c(0, 0, 0), + PCT012L196 = c(2, 0, 7), PCT012M196 = c(0, 1, 1), PCT012L197 = c(0, + 0, 2), PCT012M197 = c(0, 0, 0), PCT012L198 = c(2, 0, 6), + PCT012M198 = c(0, 0, 0), PCT012L199 = c(1, 0, 1), PCT012M199 = c(0, + 0, 0), PCT012L200 = c(0, 0, 2), PCT012M200 = c(0, 0, 0), + PCT012L201 = c(0, 0, 0), PCT012M201 = c(0, 0, 0), PCT012L202 = c(0, + 0, 2), PCT012M202 = c(0, 0, 0), PCT012L203 = c(0, 0, 1), + PCT012M203 = c(0, 0, 0), PCT012L204 = c(0, 1, 0), PCT012M204 = c(0, + 0, 0), PCT012L205 = c(0, 0, 1), PCT012M205 = c(0, 0, 0), + PCT012L206 = c(0, 0, 0), PCT012M206 = c(0, 0, 0), PCT012L207 = c(0, + 0, 0), PCT012M207 = c(0, 0, 0), PCT012L208 = c(0, 0, 0), + PCT012M208 = c(0, 0, 0), PCT012L209 = c(0, 0, 0), PCT012M209 = c(0, + 0, 0), PCT012K088 = c(1, 0, 1), PCT012N088 = c(0, 0, 0), + PCT012O088 = c(0, 2, 4), PCT012K089 = c(0, 1, 0), PCT012N089 = c(0, + 0, 0), PCT012O089 = c(1, 0, 4), PCT012K090 = c(0, 0, 1), + PCT012N090 = c(0, 0, 0), PCT012O090 = c(2, 1, 0), PCT012K091 = c(1, + 0, 0), PCT012N091 = c(0, 0, 0), PCT012O091 = c(0, 1, 1), + PCT012K092 = c(0, 0, 0), PCT012N092 = c(0, 0, 0), PCT012O092 = c(0, + 0, 1), PCT012K093 = c(0, 0, 0), PCT012N093 = c(0, 0, 0), + PCT012O093 = c(1, 1, 0), PCT012K094 = c(0, 1, 0), PCT012N094 = c(0, + 0, 0), PCT012O094 = c(0, 1, 2), PCT012K095 = c(0, 0, 0), + PCT012N095 = c(0, 0, 0), PCT012O095 = c(0, 0, 1), PCT012K096 = c(0, + 1, 0), PCT012N096 = c(0, 0, 0), PCT012O096 = c(0, 0, 1), + PCT012K097 = c(1, 0, 0), PCT012N097 = c(0, 0, 0), PCT012O097 = c(0, + 0, 1), PCT012K098 = c(0, 0, 0), PCT012N098 = c(0, 0, 0), + PCT012O098 = c(1, 0, 0), PCT012K099 = c(0, 0, 0), PCT012N099 = c(0, + 0, 0), PCT012O099 = c(0, 0, 0), PCT012K100 = c(0, 0, 1), + PCT012N100 = c(0, 0, 0), PCT012O100 = c(0, 0, 0), PCT012K101 = c(0, + 0, 0), PCT012N101 = c(0, 0, 0), PCT012O101 = c(0, 0, 0), + PCT012K102 = c(0, 1, 0), PCT012N102 = c(0, 0, 0), PCT012O102 = c(0, + 0, 0), PCT012K103 = c(0, 0, 0), PCT012N103 = c(0, 0, 0), + PCT012O103 = c(0, 0, 0), PCT012K104 = c(0, 0, 0), PCT012N104 = c(0, + 0, 0), PCT012O104 = c(0, 0, 0), PCT012K105 = c(0, 0, 0), + PCT012N105 = c(0, 0, 0), PCT012O105 = c(0, 0, 0), PCT012K192 = c(3, + 1, 4), PCT012N192 = c(0, 0, 0), PCT012O192 = c(1, 2, 5), + PCT012K193 = c(1, 3, 1), PCT012N193 = c(0, 0, 0), PCT012O193 = c(1, + 2, 7), PCT012K194 = c(1, 1, 0), PCT012N194 = c(0, 0, 1), + PCT012O194 = c(0, 3, 6), PCT012K195 = c(0, 1, 0), PCT012N195 = c(0, + 0, 0), PCT012O195 = c(0, 1, 5), PCT012K196 = c(0, 2, 0), + PCT012N196 = c(0, 0, 0), PCT012O196 = c(0, 0, 7), PCT012K197 = c(0, + 0, 1), PCT012N197 = c(0, 0, 0), PCT012O197 = c(0, 2, 2), + PCT012K198 = c(0, 2, 3), PCT012N198 = c(1, 0, 0), PCT012O198 = c(0, + 0, 3), PCT012K199 = c(0, 1, 0), PCT012N199 = c(0, 1, 0), + PCT012O199 = c(0, 1, 0), PCT012K200 = c(0, 1, 0), PCT012N200 = c(0, + 0, 0), PCT012O200 = c(1, 1, 0), PCT012K201 = c(0, 0, 0), + PCT012N201 = c(0, 0, 0), PCT012O201 = c(1, 0, 0), PCT012K202 = c(0, + 0, 0), PCT012N202 = c(0, 0, 0), PCT012O202 = c(0, 1, 1), + PCT012K203 = c(0, 2, 1), PCT012N203 = c(0, 0, 0), PCT012O203 = c(0, + 1, 0), PCT012K204 = c(0, 0, 0), PCT012N204 = c(0, 0, 0), + PCT012O204 = c(0, 0, 0), PCT012K205 = c(0, 0, 0), PCT012N205 = c(0, + 0, 0), PCT012O205 = c(0, 0, 0), PCT012K206 = c(0, 0, 0), + PCT012N206 = c(0, 0, 0), PCT012O206 = c(0, 0, 1), PCT012K207 = c(0, + 0, 0), PCT012N207 = c(0, 0, 0), PCT012O207 = c(0, 0, 2), + PCT012K208 = c(0, 0, 0), PCT012N208 = c(0, 0, 0), PCT012O208 = c(0, + 0, 0), PCT012K209 = c(0, 0, 0), PCT012N209 = c(0, 0, 0), + PCT012O209 = c(0, 0, 0), r_whi = c(0.180469772578534, 0.253982943390052, + 0.565547284031414), r_bla = c(0.20243920720412, 0.131404525061301, + 0.666156267734578), r_his = c(0.127640977315251, 0.23154559484301, + 0.64081342784174), r_asi = c(0.117004133678974, 0.0690814825194423, + 0.813914383801583), r_oth = c(0.261609154774003, 0.205452215791101, + 0.532938629434896), r_1_whi = c(0.0103484947643979, 0.0105427846858639, + 0.0260587096422339), r_1_bla = c(0.0146427385936546, 0.0104292704864494, + 0.0473814393249885), r_1_his = c(0.0166072574807774, 0.0334740033596919, + 0.0732713292634627), r_1_asi = c(0.00662089259440902, 0.00378336719680516, + 0.0621102781475513), r_1_oth = c(0.0414439093359254, 0.0350815181372332, + 0.0870852295320991), r_2_whi = c(0.0107285531195462, 0.0110932727966841, + 0.028253844895288), r_2_bla = c(0.0154029831568352, 0.00988317932134788, + 0.0493570044222677), r_2_his = c(0.0140533453517434, 0.0259078679613772, + 0.0676718427773453), r_2_asi = c(0.00802213970433686, 0.00521964548448119, + 0.0575562250402859), r_2_oth = c(0.0378208810144479, 0.0298678920160827, + 0.0763486943843061), r_3_whi = c(0.0107984293193717, 0.0116386480148342, + 0.0309943553664921), r_3_bla = c(0.0160615048559283, 0.00978681029221231, + 0.0528851816556199), r_3_his = c(0.0112126302563472, 0.0185875636770872, + 0.058330260444408), r_3_asi = c(0.00693617319414279, 0.00483430252925103, + 0.0459258740278848), r_3_oth = c(0.033756019970839, 0.0235496840896037, + 0.0615914814651173), r_4_whi = c(0.00678651287085515, 0.00793350512652705, + 0.0207890215968586), r_4_bla = c(0.00990459466115579, 0.00679401655405767, + 0.0347892195179407), r_4_his = c(0.00701984403381544, 0.0101473620955737, + 0.0336242334849292), r_4_asi = c(0.00399355426329433, 0.00304771246409304, + 0.0267287886218735), r_4_oth = c(0.016082711085583, 0.013254981663942, + 0.0329607210710025), r_5_whi = c(0.00460160340314136, 0.00465614092495637, + 0.0187779504799302), r_5_bla = c(0.0113394224282854, 0.00433660631110064, + 0.0224914606332516), r_5_his = c(0.00520342524685541, 0.00641892353286625, + 0.0257986096884773), r_5_asi = c(0.00266236950886289, 0.00157640299866882, + 0.0248721362012191), r_5_oth = c(0.011973666769761, 0.00600892502098705, + 0.0212521539345204), r_6_whi = c(0.00225137707242583, 0.00229909740401396, + 0.00987981293630017), r_6_bla = c(0.005375250291784, 0.00189525757299954, + 0.0105845317000567), r_6_his = c(0.00243099657202169, 0.00357820843747012, + 0.0131792791685445), r_6_asi = c(0.00147130946542423, 0.000735654732712114, + 0.0136971904995446), r_6_oth = c(0.00428577740467459, 0.00207661379401758, + 0.0104272522423011), r_7_whi = c(0.00251043030104712, 0.00198721095113438, + 0.00958496945898778), r_7_bla = c(0.00455611354413166, 0.00177211936910409, + 0.00946558019509375), r_7_his = c(0.0025675694131465, 0.00422010079075675, + 0.012482757678808), r_7_asi = c(0.00112099768794227, 0.000560498843971134, + 0.0114201639459119), r_7_oth = c(0.00477179339901913, 0.00273936287721469, + 0.00901338753148058), r_8_whi = c(0.007042157504363, 0.00668255071989529, + 0.022169502617801), r_8_bla = c(0.00990459466115579, 0.00549303466072748, + 0.0279577261192192), r_8_his = c(0.00767539367121454, 0.0150912989442919, + 0.034867046339165), r_8_asi = c(0.0050795207734884, 0.00238212008687732, + 0.0347859595039585), r_8_oth = c(0.0112225511421376, 0.00600892502098705, + 0.0211637873900941), r_9_whi = c(0.0113455088350785, 0.0116420566099476, + 0.034428514943281), r_9_bla = c(0.0128920345643585, 0.00869462796200919, + 0.0447634140334722), r_9_his = c(0.0117042924843966, 0.0273691973614127, + 0.0615533794949536), r_9_asi = c(0.00760176557135851, 0.00413367897428712, + 0.0875429131927415), r_9_oth = c(0.0156850616356647, 0.0100296027923828, + 0.0325630716210843), r_10_whi = c(0.010547897578534, 0.0108768270069808, + 0.0317902623254799), r_10_bla = c(0.0112751764088617, 0.0076399224764699, + 0.0449561520917433), r_10_his = c(0.00971032900397427, 0.0233129839800057, + 0.0571830485789596), r_10_asi = c(0.00826735794857423, 0.00490436488474743, + 0.0874027884817488), r_10_oth = c(0.0124154994918924, 0.00865992135377546, + 0.0312375734546901), r_11_whi = c(0.0108648969240838, 0.0123408186082024, + 0.0336939626963351), r_11_bla = c(0.0121478515060338, 0.00742576907839085, + 0.0470173785482541), r_11_his = c(0.0085904317067508, 0.0194479725761735, + 0.0498490870105571), r_11_asi = c(0.00991382330273944, 0.00592026903944511, + 0.0794507111329083), r_11_oth = c(0.0118411169531215, 0.00790880572615208, + 0.028056377855344), r_12_whi = c(0.0123732002617801, 0.0148563618019197, + 0.0389994409904014), r_12_bla = c(0.0130419419430138, 0.00933173432129434, + 0.0522213061215749), r_12_his = c(0.0068969284768031, 0.0138348288059437, + 0.0439491402739651), r_12_asi = c(0.0104743221467106, 0.0059903313949415, + 0.0613746234148392), r_12_oth = c(0.0108690849644325, 0.00773207263729952, + 0.0258913975169001), r_13_whi = c(0.0142956479057592, 0.0187694289921466, + 0.0460910231239093), r_13_bla = c(0.014107355098457, 0.00997954835048345, + 0.0530243813643713), r_13_his = c(0.00647355266931618, 0.00906843665068764, + 0.0338700645989539), r_13_asi = c(0.0100539480137322, 0.00592026903944511, + 0.0544034190429482), r_13_oth = c(0.0107807184200062, 0.0114876507754164, + 0.0254937480669818), r_14_whi = c(0.0135491655759162, 0.0200459478621291, + 0.0477476003490401), r_14_bla = c(0.0126243428167596, 0.00984034864173207, + 0.0476223618978274), r_14_his = c(0.00498490870105571, 0.00643258081697874, + 0.0242143647314295), r_14_asi = c(0.00802213970433686, 0.00518461430673299, + 0.0443845022069642), r_14_oth = c(0.0109132682366456, 0.00888083771484116, + 0.021119604117881), r_15_whi = c(0.0115568417321117, 0.0209338868891798, + 0.0415763388961606), r_15_bla = c(0.0109593001466951, 0.00799327558330032, + 0.0368076152948357), r_15_his = c(0.0033050627552205, 0.00465713388235616, + 0.0178227557667882), r_15_asi = c(0.00655083023891263, 0.00364324248581237, + 0.0364674560358719), r_15_oth = c(0.00755533954844696, 0.0078204391817258, + 0.0165245438077144), r_16_whi = c(0.00452491001308901, 0.00917082515270506, + 0.0147200179973822), r_16_bla = c(0.00390829951494255, 0.00267691747598805, + 0.0132025569915731), r_16_his = c(0.00140670026358558, 0.00151595853648543, + 0.00549022821321752), r_16_asi = c(0.00266236950886289, 0.00178659006515799, + 0.0121207875008758), r_16_oth = c(0.00269517960500155, 0.00269517960500155, + 0.00477179339901913), r_17_whi = c(0.00687343204624782, 0.0147830770069808, + 0.020799247382199), r_17_bla = c(0.00532171194226424, 0.00354959257316015, + 0.0164898116520864), r_17_his = c(0.00177544693462258, 0.00189836249163491, + 0.0068969284768031), r_17_asi = c(0.00360821130806418, 0.00203180830939536, + 0.0178308694738317), r_17_oth = c(0.00393231122696947, 0.00432996067688773, + 0.00645075774311846), r_18_whi = c(0.00403918520942408, 0.0088964332460733, + 0.011166557591623), r_18_bla = c(0.00278934800997955, 0.00196485742737523, + 0.0087802893212408), r_18_his = c(0.000655549637399107, 0.0010242963084361, + 0.00342797831223283), r_18_asi = c(0.0025923071533665, 0.00122609122118686, + 0.00910810621453093), r_18_oth = c(0.00203243052180444, 0.00189988070516502, + 0.00304864578270667), r_19_whi = c(0.00547590804973822, 0.0126595222513089, + 0.0152568717277487), r_19_bla = c(0.00387082267027872, 0.0028000556798835, + 0.0113608377680933), r_19_his = c(0.00127012742246077, 0.00135207112713566, + 0.00435667363188156), r_19_asi = c(0.00318783717508583, 0.00164646535416521, + 0.013241785188818), r_19_oth = c(0.00287191269385411, 0.00371139486590377, + 0.00446251049352715), r_20_whi = c(0.00717168411867365, 0.016080047447644, + 0.0191648260253054), r_20_bla = c(0.00503260485485753, 0.00315340878671392, + 0.014123416603313), r_20_his = c(0.00162521680938529, 0.00187104792340995, + 0.00529902623564278), r_20_asi = c(0.00353814895256779, 0.00217193302038815, + 0.0167449029636376), r_20_oth = c(0.00384394468254319, 0.00508107630451111, + 0.00450669376574029), r_21_whi = c(0.00569405813699825, 0.011773287521815, + 0.0165981539048866), r_21_bla = c(0.00362454626248782, 0.00248953325266889, + 0.00955124155432536), r_21_his = c(0.000969667171986179, + 0.00107892544488603, 0.00349626473279524), r_21_asi = c(0.00248721362012191, + 0.00140124710992784, 0.00879282561479717), r_21_oth = c(0.00238589669950957, + 0.00326956214377237, 0.00331374541598551), r_22_whi = c(0.00372900305410122, + 0.0079471395069808, 0.0135270097076789), r_22_bla = c(0.00182565771862385, + 0.00177211936910409, 0.00606054116563695), r_22_his = c(0.000833094330861365, + 0.000696521489736551, 0.0023490528673468), r_22_asi = c(0.00136621593217964, + 0.000735654732712114, 0.00437889721852449), r_22_oth = c(0.00163478107188618, + 0.00163478107188618, 0.00265099633278841), r_23_whi = c(0.00336087478184991, + 0.00637407286212914, 0.0134792893760908), r_23_bla = c(0.00183101155357583, + 0.0017025195147284, 0.00526281975779251), r_23_his = c(0.000669206921511588, + 0.000559948648611737, 0.00183007607107251), r_23_asi = c(0.00077068591046031, + 0.000245218244237371, 0.00357318013031598), r_23_oth = c(0.000795298899836522, + 0.00172314761631246, 0.00300446251049353)), class = c("tbl_df", "tbl", "data.frame"), row.names = c(NA, -3L)) diff --git a/tests/testthat/test-census_geo_api.R b/tests/testthat/test-census_geo_api.R index d4d3534..8588bf6 100644 --- a/tests/testthat/test-census_geo_api.R +++ b/tests/testthat/test-census_geo_api.R @@ -2,7 +2,6 @@ skip_if_not(nzchar(Sys.getenv("CENSUS_API_KEY"))) test_that("snapshot", { # TODO: Test that sub-geographies sum to match pooled geographies (e.g. blocks sum to block groups, sum to tracts, sum to counties) - # TODO: Test that age/sex subsets sum to race totals # These snapshots were generated using the calculations in v2.0.0 # and verified that the calculations resulted in the same numbers for PR #120. @@ -23,3 +22,51 @@ test_that("snapshot", { style = "deparse" ) }) + +expect_subset_sums_equal_overall_total <- function(data) { + `%>%` <- dplyr::`%>%` + + sums <- data %>% + dplyr::select(-dplyr::starts_with("r_")) %>% + tidyr::pivot_longer(dplyr::starts_with("P")) %>% + dplyr::mutate( + subset = dplyr::case_when( + grepl("001", name) ~ "overall", + .default = "subset" + ), + name = sub("_.+", "", name) + ) %>% + dplyr::summarise( + value = sum(value), + .by = -value + ) %>% + dplyr::summarize( + are_equal = length(unique(value)) <= 1, + .by = c(-subset, -value) + ) + + expect_true(all(sums$are_equal)) +} + +test_that("sums", { + expect_subset_sums_equal_overall_total( + census_geo_api(state = "DE", geo = "county", year = "2020", sex = TRUE) + ) + expect_subset_sums_equal_overall_total( + census_geo_api(state = "DE", geo = "county", year = "2020", age = TRUE) + ) + expect_subset_sums_equal_overall_total( + census_geo_api(state = "DE", geo = "county", year = "2020", age = TRUE, sex = TRUE) + ) + expect_subset_sums_equal_overall_total( + census_geo_api(state = "DE", geo = "county", year = "2010", sex = TRUE) + ) + expect_subset_sums_equal_overall_total( + census_geo_api(state = "DE", geo = "county", year = "2010", age = TRUE) + ) + expect_subset_sums_equal_overall_total( + census_geo_api(state = "DE", geo = "county", year = "2010", age = TRUE, sex = TRUE) + ) +}) + +# TODO: Test that all variables sum to total population of geography diff --git a/tests/testthat/test-census_geo_api_names.R b/tests/testthat/test-census_geo_api_names.R index 2d7bd79..85d6c0e 100644 --- a/tests/testthat/test-census_geo_api_names.R +++ b/tests/testthat/test-census_geo_api_names.R @@ -5,10 +5,10 @@ test_that("census_geo_api_names() for 2020", { census_geo_api_names("2020"), list( r_whi = "P12I_001N", - r_bla = "P12B_001N", + r_bla = "P12J_001N", r_his = "P12H_001N", - r_asi = c("P12D_001N", "P12E_001N"), - r_oth = c("P12C_001N", "P12F_001N", "P12G_001N") + r_asi = c("P12L_001N", "P12M_001N"), + r_oth = c("P12K_001N", "P12N_001N", "P12O_001N") ) ) @@ -18,15 +18,15 @@ test_that("census_geo_api_names() for 2020", { census_geo_api_names("2020", age = FALSE, sex = FALSE), list( r_mal_whi = "P12I_002N", - r_mal_bla = "P12B_002N", + r_mal_bla = "P12J_002N", r_mal_his = "P12H_002N", - r_mal_asi = c("P12D_002N", "P12E_002N"), - r_mal_oth = c("P12C_002N", "P12F_002N", "P12G_002N"), + r_mal_asi = c("P12L_002N", "P12M_002N"), + r_mal_oth = c("P12K_002N", "P12N_002N", "P12O_002N"), r_fem_whi = "P12I_026N", - r_fem_bla = "P12B_026N", + r_fem_bla = "P12J_026N", r_fem_his = "P12H_026N", - r_fem_asi = c("P12D_026N", "P12E_026N"), - r_fem_oth = c("P12C_026N", "P12F_026N", "P12G_026N") + r_fem_asi = c("P12L_026N", "P12M_026N"), + r_fem_oth = c("P12K_026N", "P12N_026N", "P12O_026N") ) ) ) @@ -37,120 +37,120 @@ test_that("census_geo_api_names() for 2020", { census_geo_api_names("2020", age = FALSE, sex = FALSE), list( r_1_whi = c("P12I_003N", "P12I_027N"), - r_1_bla = c("P12B_003N", "P12B_027N"), + r_1_bla = c("P12J_003N", "P12J_027N"), r_1_his = c("P12H_003N", "P12H_027N"), - r_1_asi = c("P12D_003N", "P12E_003N", "P12D_027N", "P12E_027N"), - r_1_oth = c("P12C_003N", "P12F_003N", "P12G_003N", "P12C_027N", "P12F_027N", "P12G_027N"), + r_1_asi = c("P12L_003N", "P12M_003N", "P12L_027N", "P12M_027N"), + r_1_oth = c("P12K_003N", "P12N_003N", "P12O_003N", "P12K_027N", "P12N_027N", "P12O_027N"), r_2_whi = c("P12I_004N", "P12I_028N"), - r_2_bla = c("P12B_004N", "P12B_028N"), + r_2_bla = c("P12J_004N", "P12J_028N"), r_2_his = c("P12H_004N", "P12H_028N"), - r_2_asi = c("P12D_004N", "P12E_004N", "P12D_028N", "P12E_028N"), - r_2_oth = c("P12C_004N", "P12F_004N", "P12G_004N", "P12C_028N", "P12F_028N", "P12G_028N"), + r_2_asi = c("P12L_004N", "P12M_004N", "P12L_028N", "P12M_028N"), + r_2_oth = c("P12K_004N", "P12N_004N", "P12O_004N", "P12K_028N", "P12N_028N", "P12O_028N"), r_3_whi = c("P12I_005N", "P12I_029N"), - r_3_bla = c("P12B_005N", "P12B_029N"), + r_3_bla = c("P12J_005N", "P12J_029N"), r_3_his = c("P12H_005N", "P12H_029N"), - r_3_asi = c("P12D_005N", "P12E_005N", "P12D_029N", "P12E_029N"), - r_3_oth = c("P12C_005N", "P12F_005N", "P12G_005N", "P12C_029N", "P12F_029N", "P12G_029N"), + r_3_asi = c("P12L_005N", "P12M_005N", "P12L_029N", "P12M_029N"), + r_3_oth = c("P12K_005N", "P12N_005N", "P12O_005N", "P12K_029N", "P12N_029N", "P12O_029N"), r_4_whi = c("P12I_006N", "P12I_030N"), - r_4_bla = c("P12B_006N", "P12B_030N"), + r_4_bla = c("P12J_006N", "P12J_030N"), r_4_his = c("P12H_006N", "P12H_030N"), - r_4_asi = c("P12D_006N", "P12E_006N", "P12D_030N", "P12E_030N"), - r_4_oth = c("P12C_006N", "P12F_006N", "P12G_006N", "P12C_030N", "P12F_030N", "P12G_030N"), + r_4_asi = c("P12L_006N", "P12M_006N", "P12L_030N", "P12M_030N"), + r_4_oth = c("P12K_006N", "P12N_006N", "P12O_006N", "P12K_030N", "P12N_030N", "P12O_030N"), r_5_whi = c("P12I_007N", "P12I_031N"), - r_5_bla = c("P12B_007N", "P12B_031N"), + r_5_bla = c("P12J_007N", "P12J_031N"), r_5_his = c("P12H_007N", "P12H_031N"), - r_5_asi = c("P12D_007N", "P12E_007N", "P12D_031N", "P12E_031N"), - r_5_oth = c("P12C_007N", "P12F_007N", "P12G_007N", "P12C_031N", "P12F_031N", "P12G_031N"), + r_5_asi = c("P12L_007N", "P12M_007N", "P12L_031N", "P12M_031N"), + r_5_oth = c("P12K_007N", "P12N_007N", "P12O_007N", "P12K_031N", "P12N_031N", "P12O_031N"), r_6_whi = c("P12I_008N", "P12I_032N"), - r_6_bla = c("P12B_008N", "P12B_032N"), + r_6_bla = c("P12J_008N", "P12J_032N"), r_6_his = c("P12H_008N", "P12H_032N"), - r_6_asi = c("P12D_008N", "P12E_008N", "P12D_032N", "P12E_032N"), - r_6_oth = c("P12C_008N", "P12F_008N", "P12G_008N", "P12C_032N", "P12F_032N", "P12G_032N"), + r_6_asi = c("P12L_008N", "P12M_008N", "P12L_032N", "P12M_032N"), + r_6_oth = c("P12K_008N", "P12N_008N", "P12O_008N", "P12K_032N", "P12N_032N", "P12O_032N"), r_7_whi = c("P12I_009N", "P12I_033N"), - r_7_bla = c("P12B_009N", "P12B_033N"), + r_7_bla = c("P12J_009N", "P12J_033N"), r_7_his = c("P12H_009N", "P12H_033N"), - r_7_asi = c("P12D_009N", "P12E_009N", "P12D_033N", "P12E_033N"), - r_7_oth = c("P12C_009N", "P12F_009N", "P12G_009N", "P12C_033N", "P12F_033N", "P12G_033N"), + r_7_asi = c("P12L_009N", "P12M_009N", "P12L_033N", "P12M_033N"), + r_7_oth = c("P12K_009N", "P12N_009N", "P12O_009N", "P12K_033N", "P12N_033N", "P12O_033N"), r_8_whi = c("P12I_010N", "P12I_034N"), - r_8_bla = c("P12B_010N", "P12B_034N"), + r_8_bla = c("P12J_010N", "P12J_034N"), r_8_his = c("P12H_010N", "P12H_034N"), - r_8_asi = c("P12D_010N", "P12E_010N", "P12D_034N", "P12E_034N"), - r_8_oth = c("P12C_010N", "P12F_010N", "P12G_010N", "P12C_034N", "P12F_034N", "P12G_034N"), + r_8_asi = c("P12L_010N", "P12M_010N", "P12L_034N", "P12M_034N"), + r_8_oth = c("P12K_010N", "P12N_010N", "P12O_010N", "P12K_034N", "P12N_034N", "P12O_034N"), r_9_whi = c("P12I_011N", "P12I_035N"), - r_9_bla = c("P12B_011N", "P12B_035N"), + r_9_bla = c("P12J_011N", "P12J_035N"), r_9_his = c("P12H_011N", "P12H_035N"), - r_9_asi = c("P12D_011N", "P12E_011N", "P12D_035N", "P12E_035N"), - r_9_oth = c("P12C_011N", "P12F_011N", "P12G_011N", "P12C_035N", "P12F_035N", "P12G_035N"), + r_9_asi = c("P12L_011N", "P12M_011N", "P12L_035N", "P12M_035N"), + r_9_oth = c("P12K_011N", "P12N_011N", "P12O_011N", "P12K_035N", "P12N_035N", "P12O_035N"), r_10_whi = c("P12I_012N", "P12I_036N"), - r_10_bla = c("P12B_012N", "P12B_036N"), + r_10_bla = c("P12J_012N", "P12J_036N"), r_10_his = c("P12H_012N", "P12H_036N"), - r_10_asi = c("P12D_012N", "P12E_012N", "P12D_036N", "P12E_036N"), - r_10_oth = c("P12C_012N", "P12F_012N", "P12G_012N", "P12C_036N", "P12F_036N", "P12G_036N"), + r_10_asi = c("P12L_012N", "P12M_012N", "P12L_036N", "P12M_036N"), + r_10_oth = c("P12K_012N", "P12N_012N", "P12O_012N", "P12K_036N", "P12N_036N", "P12O_036N"), r_11_whi = c("P12I_013N", "P12I_037N"), - r_11_bla = c("P12B_013N", "P12B_037N"), + r_11_bla = c("P12J_013N", "P12J_037N"), r_11_his = c("P12H_013N", "P12H_037N"), - r_11_asi = c("P12D_013N", "P12E_013N", "P12D_037N", "P12E_037N"), - r_11_oth = c("P12C_013N", "P12F_013N", "P12G_013N", "P12C_037N", "P12F_037N", "P12G_037N"), + r_11_asi = c("P12L_013N", "P12M_013N", "P12L_037N", "P12M_037N"), + r_11_oth = c("P12K_013N", "P12N_013N", "P12O_013N", "P12K_037N", "P12N_037N", "P12O_037N"), r_12_whi = c("P12I_014N", "P12I_038N"), - r_12_bla = c("P12B_014N", "P12B_038N"), + r_12_bla = c("P12J_014N", "P12J_038N"), r_12_his = c("P12H_014N", "P12H_038N"), - r_12_asi = c("P12D_014N", "P12E_014N", "P12D_038N", "P12E_038N"), - r_12_oth = c("P12C_014N", "P12F_014N", "P12G_014N", "P12C_038N", "P12F_038N", "P12G_038N"), + r_12_asi = c("P12L_014N", "P12M_014N", "P12L_038N", "P12M_038N"), + r_12_oth = c("P12K_014N", "P12N_014N", "P12O_014N", "P12K_038N", "P12N_038N", "P12O_038N"), r_13_whi = c("P12I_015N", "P12I_039N"), - r_13_bla = c("P12B_015N", "P12B_039N"), + r_13_bla = c("P12J_015N", "P12J_039N"), r_13_his = c("P12H_015N", "P12H_039N"), - r_13_asi = c("P12D_015N", "P12E_015N", "P12D_039N", "P12E_039N"), - r_13_oth = c("P12C_015N", "P12F_015N", "P12G_015N", "P12C_039N", "P12F_039N", "P12G_039N"), + r_13_asi = c("P12L_015N", "P12M_015N", "P12L_039N", "P12M_039N"), + r_13_oth = c("P12K_015N", "P12N_015N", "P12O_015N", "P12K_039N", "P12N_039N", "P12O_039N"), r_14_whi = c("P12I_016N", "P12I_040N"), - r_14_bla = c("P12B_016N", "P12B_040N"), + r_14_bla = c("P12J_016N", "P12J_040N"), r_14_his = c("P12H_016N", "P12H_040N"), - r_14_asi = c("P12D_016N", "P12E_016N", "P12D_040N", "P12E_040N"), - r_14_oth = c("P12C_016N", "P12F_016N", "P12G_016N", "P12C_040N", "P12F_040N", "P12G_040N"), + r_14_asi = c("P12L_016N", "P12M_016N", "P12L_040N", "P12M_040N"), + r_14_oth = c("P12K_016N", "P12N_016N", "P12O_016N", "P12K_040N", "P12N_040N", "P12O_040N"), r_15_whi = c("P12I_017N", "P12I_041N"), - r_15_bla = c("P12B_017N", "P12B_041N"), + r_15_bla = c("P12J_017N", "P12J_041N"), r_15_his = c("P12H_017N", "P12H_041N"), - r_15_asi = c("P12D_017N", "P12E_017N", "P12D_041N", "P12E_041N"), - r_15_oth = c("P12C_017N", "P12F_017N", "P12G_017N", "P12C_041N", "P12F_041N", "P12G_041N"), + r_15_asi = c("P12L_017N", "P12M_017N", "P12L_041N", "P12M_041N"), + r_15_oth = c("P12K_017N", "P12N_017N", "P12O_017N", "P12K_041N", "P12N_041N", "P12O_041N"), r_16_whi = c("P12I_018N", "P12I_042N"), - r_16_bla = c("P12B_018N", "P12B_042N"), + r_16_bla = c("P12J_018N", "P12J_042N"), r_16_his = c("P12H_018N", "P12H_042N"), - r_16_asi = c("P12D_018N", "P12E_018N", "P12D_042N", "P12E_042N"), - r_16_oth = c("P12C_018N", "P12F_018N", "P12G_018N", "P12C_042N", "P12F_042N", "P12G_042N"), + r_16_asi = c("P12L_018N", "P12M_018N", "P12L_042N", "P12M_042N"), + r_16_oth = c("P12K_018N", "P12N_018N", "P12O_018N", "P12K_042N", "P12N_042N", "P12O_042N"), r_17_whi = c("P12I_019N", "P12I_043N"), - r_17_bla = c("P12B_019N", "P12B_043N"), + r_17_bla = c("P12J_019N", "P12J_043N"), r_17_his = c("P12H_019N", "P12H_043N"), - r_17_asi = c("P12D_019N", "P12E_019N", "P12D_043N", "P12E_043N"), - r_17_oth = c("P12C_019N", "P12F_019N", "P12G_019N", "P12C_043N", "P12F_043N", "P12G_043N"), + r_17_asi = c("P12L_019N", "P12M_019N", "P12L_043N", "P12M_043N"), + r_17_oth = c("P12K_019N", "P12N_019N", "P12O_019N", "P12K_043N", "P12N_043N", "P12O_043N"), r_18_whi = c("P12I_020N", "P12I_044N"), - r_18_bla = c("P12B_020N", "P12B_044N"), + r_18_bla = c("P12J_020N", "P12J_044N"), r_18_his = c("P12H_020N", "P12H_044N"), - r_18_asi = c("P12D_020N", "P12E_020N", "P12D_044N", "P12E_044N"), - r_18_oth = c("P12C_020N", "P12F_020N", "P12G_020N", "P12C_044N", "P12F_044N", "P12G_044N"), + r_18_asi = c("P12L_020N", "P12M_020N", "P12L_044N", "P12M_044N"), + r_18_oth = c("P12K_020N", "P12N_020N", "P12O_020N", "P12K_044N", "P12N_044N", "P12O_044N"), r_19_whi = c("P12I_021N", "P12I_045N"), - r_19_bla = c("P12B_021N", "P12B_045N"), + r_19_bla = c("P12J_021N", "P12J_045N"), r_19_his = c("P12H_021N", "P12H_045N"), - r_19_asi = c("P12D_021N", "P12E_021N", "P12D_045N", "P12E_045N"), - r_19_oth = c("P12C_021N", "P12F_021N", "P12G_021N", "P12C_045N", "P12F_045N", "P12G_045N"), + r_19_asi = c("P12L_021N", "P12M_021N", "P12L_045N", "P12M_045N"), + r_19_oth = c("P12K_021N", "P12N_021N", "P12O_021N", "P12K_045N", "P12N_045N", "P12O_045N"), r_20_whi = c("P12I_022N", "P12I_046N"), - r_20_bla = c("P12B_022N", "P12B_046N"), + r_20_bla = c("P12J_022N", "P12J_046N"), r_20_his = c("P12H_022N", "P12H_046N"), - r_20_asi = c("P12D_022N", "P12E_022N", "P12D_046N", "P12E_046N"), - r_20_oth = c("P12C_022N", "P12F_022N", "P12G_022N", "P12C_046N", "P12F_046N", "P12G_046N"), + r_20_asi = c("P12L_022N", "P12M_022N", "P12L_046N", "P12M_046N"), + r_20_oth = c("P12K_022N", "P12N_022N", "P12O_022N", "P12K_046N", "P12N_046N", "P12O_046N"), r_21_whi = c("P12I_023N", "P12I_047N"), - r_21_bla = c("P12B_023N", "P12B_047N"), + r_21_bla = c("P12J_023N", "P12J_047N"), r_21_his = c("P12H_023N", "P12H_047N"), - r_21_asi = c("P12D_023N", "P12E_023N", "P12D_047N", "P12E_047N"), - r_21_oth = c("P12C_023N", "P12F_023N", "P12G_023N", "P12C_047N", "P12F_047N", "P12G_047N"), + r_21_asi = c("P12L_023N", "P12M_023N", "P12L_047N", "P12M_047N"), + r_21_oth = c("P12K_023N", "P12N_023N", "P12O_023N", "P12K_047N", "P12N_047N", "P12O_047N"), r_22_whi = c("P12I_024N", "P12I_048N"), - r_22_bla = c("P12B_024N", "P12B_048N"), + r_22_bla = c("P12J_024N", "P12J_048N"), r_22_his = c("P12H_024N", "P12H_048N"), - r_22_asi = c("P12D_024N", "P12E_024N", "P12D_048N", "P12E_048N"), - r_22_oth = c("P12C_024N", "P12F_024N", "P12G_024N", "P12C_048N", "P12F_048N", "P12G_048N"), + r_22_asi = c("P12L_024N", "P12M_024N", "P12L_048N", "P12M_048N"), + r_22_oth = c("P12K_024N", "P12N_024N", "P12O_024N", "P12K_048N", "P12N_048N", "P12O_048N"), r_23_whi = c("P12I_025N", "P12I_049N"), - r_23_bla = c("P12B_025N", "P12B_049N"), + r_23_bla = c("P12J_025N", "P12J_049N"), r_23_his = c("P12H_025N", "P12H_049N"), - r_23_asi = c("P12D_025N", "P12E_025N", "P12D_049N", "P12E_049N"), - r_23_oth = c("P12C_025N", "P12F_025N", "P12G_025N", "P12C_049N", "P12F_049N", "P12G_049N") + r_23_asi = c("P12L_025N", "P12M_025N", "P12L_049N", "P12M_049N"), + r_23_oth = c("P12K_025N", "P12N_025N", "P12O_025N", "P12K_049N", "P12N_049N", "P12O_049N") ) ) ) @@ -161,235 +161,235 @@ test_that("census_geo_api_names() for 2020", { census_geo_api_names("2020", age = FALSE, sex = FALSE), list( r_mal_1_whi = "P12I_003N", - r_mal_1_bla = "P12B_003N", + r_mal_1_bla = "P12J_003N", r_mal_1_his = "P12H_003N", - r_mal_1_asi = c("P12D_003N", "P12E_003N"), - r_mal_1_oth = c("P12C_003N", "P12F_003N", "P12G_003N"), + r_mal_1_asi = c("P12L_003N", "P12M_003N"), + r_mal_1_oth = c("P12K_003N", "P12N_003N", "P12O_003N"), r_mal_2_whi = "P12I_004N", - r_mal_2_bla = "P12B_004N", + r_mal_2_bla = "P12J_004N", r_mal_2_his = "P12H_004N", - r_mal_2_asi = c("P12D_004N", "P12E_004N"), - r_mal_2_oth = c("P12C_004N", "P12F_004N", "P12G_004N"), + r_mal_2_asi = c("P12L_004N", "P12M_004N"), + r_mal_2_oth = c("P12K_004N", "P12N_004N", "P12O_004N"), r_mal_3_whi = "P12I_005N", - r_mal_3_bla = "P12B_005N", + r_mal_3_bla = "P12J_005N", r_mal_3_his = "P12H_005N", - r_mal_3_asi = c("P12D_005N", "P12E_005N"), - r_mal_3_oth = c("P12C_005N", "P12F_005N", "P12G_005N"), + r_mal_3_asi = c("P12L_005N", "P12M_005N"), + r_mal_3_oth = c("P12K_005N", "P12N_005N", "P12O_005N"), r_mal_4_whi = "P12I_006N", - r_mal_4_bla = "P12B_006N", + r_mal_4_bla = "P12J_006N", r_mal_4_his = "P12H_006N", - r_mal_4_asi = c("P12D_006N", "P12E_006N"), - r_mal_4_oth = c("P12C_006N", "P12F_006N", "P12G_006N"), + r_mal_4_asi = c("P12L_006N", "P12M_006N"), + r_mal_4_oth = c("P12K_006N", "P12N_006N", "P12O_006N"), r_mal_5_whi = "P12I_007N", - r_mal_5_bla = "P12B_007N", + r_mal_5_bla = "P12J_007N", r_mal_5_his = "P12H_007N", - r_mal_5_asi = c("P12D_007N", "P12E_007N"), - r_mal_5_oth = c("P12C_007N", "P12F_007N", "P12G_007N"), + r_mal_5_asi = c("P12L_007N", "P12M_007N"), + r_mal_5_oth = c("P12K_007N", "P12N_007N", "P12O_007N"), r_mal_6_whi = "P12I_008N", - r_mal_6_bla = "P12B_008N", + r_mal_6_bla = "P12J_008N", r_mal_6_his = "P12H_008N", - r_mal_6_asi = c("P12D_008N", "P12E_008N"), - r_mal_6_oth = c("P12C_008N", "P12F_008N", "P12G_008N"), + r_mal_6_asi = c("P12L_008N", "P12M_008N"), + r_mal_6_oth = c("P12K_008N", "P12N_008N", "P12O_008N"), r_mal_7_whi = "P12I_009N", - r_mal_7_bla = "P12B_009N", + r_mal_7_bla = "P12J_009N", r_mal_7_his = "P12H_009N", - r_mal_7_asi = c("P12D_009N", "P12E_009N"), - r_mal_7_oth = c("P12C_009N", "P12F_009N", "P12G_009N"), + r_mal_7_asi = c("P12L_009N", "P12M_009N"), + r_mal_7_oth = c("P12K_009N", "P12N_009N", "P12O_009N"), r_mal_8_whi = "P12I_010N", - r_mal_8_bla = "P12B_010N", + r_mal_8_bla = "P12J_010N", r_mal_8_his = "P12H_010N", - r_mal_8_asi = c("P12D_010N", "P12E_010N"), - r_mal_8_oth = c("P12C_010N", "P12F_010N", "P12G_010N"), + r_mal_8_asi = c("P12L_010N", "P12M_010N"), + r_mal_8_oth = c("P12K_010N", "P12N_010N", "P12O_010N"), r_mal_9_whi = "P12I_011N", - r_mal_9_bla = "P12B_011N", + r_mal_9_bla = "P12J_011N", r_mal_9_his = "P12H_011N", - r_mal_9_asi = c("P12D_011N", "P12E_011N"), - r_mal_9_oth = c("P12C_011N", "P12F_011N", "P12G_011N"), + r_mal_9_asi = c("P12L_011N", "P12M_011N"), + r_mal_9_oth = c("P12K_011N", "P12N_011N", "P12O_011N"), r_mal_10_whi = "P12I_012N", - r_mal_10_bla = "P12B_012N", + r_mal_10_bla = "P12J_012N", r_mal_10_his = "P12H_012N", - r_mal_10_asi = c("P12D_012N", "P12E_012N"), - r_mal_10_oth = c("P12C_012N", "P12F_012N", "P12G_012N"), + r_mal_10_asi = c("P12L_012N", "P12M_012N"), + r_mal_10_oth = c("P12K_012N", "P12N_012N", "P12O_012N"), r_mal_11_whi = "P12I_013N", - r_mal_11_bla = "P12B_013N", + r_mal_11_bla = "P12J_013N", r_mal_11_his = "P12H_013N", - r_mal_11_asi = c("P12D_013N", "P12E_013N"), - r_mal_11_oth = c("P12C_013N", "P12F_013N", "P12G_013N"), + r_mal_11_asi = c("P12L_013N", "P12M_013N"), + r_mal_11_oth = c("P12K_013N", "P12N_013N", "P12O_013N"), r_mal_12_whi = "P12I_014N", - r_mal_12_bla = "P12B_014N", + r_mal_12_bla = "P12J_014N", r_mal_12_his = "P12H_014N", - r_mal_12_asi = c("P12D_014N", "P12E_014N"), - r_mal_12_oth = c("P12C_014N", "P12F_014N", "P12G_014N"), + r_mal_12_asi = c("P12L_014N", "P12M_014N"), + r_mal_12_oth = c("P12K_014N", "P12N_014N", "P12O_014N"), r_mal_13_whi = "P12I_015N", - r_mal_13_bla = "P12B_015N", + r_mal_13_bla = "P12J_015N", r_mal_13_his = "P12H_015N", - r_mal_13_asi = c("P12D_015N", "P12E_015N"), - r_mal_13_oth = c("P12C_015N", "P12F_015N", "P12G_015N"), + r_mal_13_asi = c("P12L_015N", "P12M_015N"), + r_mal_13_oth = c("P12K_015N", "P12N_015N", "P12O_015N"), r_mal_14_whi = "P12I_016N", - r_mal_14_bla = "P12B_016N", + r_mal_14_bla = "P12J_016N", r_mal_14_his = "P12H_016N", - r_mal_14_asi = c("P12D_016N", "P12E_016N"), - r_mal_14_oth = c("P12C_016N", "P12F_016N", "P12G_016N"), + r_mal_14_asi = c("P12L_016N", "P12M_016N"), + r_mal_14_oth = c("P12K_016N", "P12N_016N", "P12O_016N"), r_mal_15_whi = "P12I_017N", - r_mal_15_bla = "P12B_017N", + r_mal_15_bla = "P12J_017N", r_mal_15_his = "P12H_017N", - r_mal_15_asi = c("P12D_017N", "P12E_017N"), - r_mal_15_oth = c("P12C_017N", "P12F_017N", "P12G_017N"), + r_mal_15_asi = c("P12L_017N", "P12M_017N"), + r_mal_15_oth = c("P12K_017N", "P12N_017N", "P12O_017N"), r_mal_16_whi = "P12I_018N", - r_mal_16_bla = "P12B_018N", + r_mal_16_bla = "P12J_018N", r_mal_16_his = "P12H_018N", - r_mal_16_asi = c("P12D_018N", "P12E_018N"), - r_mal_16_oth = c("P12C_018N", "P12F_018N", "P12G_018N"), + r_mal_16_asi = c("P12L_018N", "P12M_018N"), + r_mal_16_oth = c("P12K_018N", "P12N_018N", "P12O_018N"), r_mal_17_whi = "P12I_019N", - r_mal_17_bla = "P12B_019N", + r_mal_17_bla = "P12J_019N", r_mal_17_his = "P12H_019N", - r_mal_17_asi = c("P12D_019N", "P12E_019N"), - r_mal_17_oth = c("P12C_019N", "P12F_019N", "P12G_019N"), + r_mal_17_asi = c("P12L_019N", "P12M_019N"), + r_mal_17_oth = c("P12K_019N", "P12N_019N", "P12O_019N"), r_mal_18_whi = "P12I_020N", - r_mal_18_bla = "P12B_020N", + r_mal_18_bla = "P12J_020N", r_mal_18_his = "P12H_020N", - r_mal_18_asi = c("P12D_020N", "P12E_020N"), - r_mal_18_oth = c("P12C_020N", "P12F_020N", "P12G_020N"), + r_mal_18_asi = c("P12L_020N", "P12M_020N"), + r_mal_18_oth = c("P12K_020N", "P12N_020N", "P12O_020N"), r_mal_19_whi = "P12I_021N", - r_mal_19_bla = "P12B_021N", + r_mal_19_bla = "P12J_021N", r_mal_19_his = "P12H_021N", - r_mal_19_asi = c("P12D_021N", "P12E_021N"), - r_mal_19_oth = c("P12C_021N", "P12F_021N", "P12G_021N"), + r_mal_19_asi = c("P12L_021N", "P12M_021N"), + r_mal_19_oth = c("P12K_021N", "P12N_021N", "P12O_021N"), r_mal_20_whi = "P12I_022N", - r_mal_20_bla = "P12B_022N", + r_mal_20_bla = "P12J_022N", r_mal_20_his = "P12H_022N", - r_mal_20_asi = c("P12D_022N", "P12E_022N"), - r_mal_20_oth = c("P12C_022N", "P12F_022N", "P12G_022N"), + r_mal_20_asi = c("P12L_022N", "P12M_022N"), + r_mal_20_oth = c("P12K_022N", "P12N_022N", "P12O_022N"), r_mal_21_whi = "P12I_023N", - r_mal_21_bla = "P12B_023N", + r_mal_21_bla = "P12J_023N", r_mal_21_his = "P12H_023N", - r_mal_21_asi = c("P12D_023N", "P12E_023N"), - r_mal_21_oth = c("P12C_023N", "P12F_023N", "P12G_023N"), + r_mal_21_asi = c("P12L_023N", "P12M_023N"), + r_mal_21_oth = c("P12K_023N", "P12N_023N", "P12O_023N"), r_mal_22_whi = "P12I_024N", - r_mal_22_bla = "P12B_024N", + r_mal_22_bla = "P12J_024N", r_mal_22_his = "P12H_024N", - r_mal_22_asi = c("P12D_024N", "P12E_024N"), - r_mal_22_oth = c("P12C_024N", "P12F_024N", "P12G_024N"), + r_mal_22_asi = c("P12L_024N", "P12M_024N"), + r_mal_22_oth = c("P12K_024N", "P12N_024N", "P12O_024N"), r_mal_23_whi = "P12I_025N", - r_mal_23_bla = "P12B_025N", + r_mal_23_bla = "P12J_025N", r_mal_23_his = "P12H_025N", - r_mal_23_asi = c("P12D_025N", "P12E_025N"), - r_mal_23_oth = c("P12C_025N", "P12F_025N", "P12G_025N"), + r_mal_23_asi = c("P12L_025N", "P12M_025N"), + r_mal_23_oth = c("P12K_025N", "P12N_025N", "P12O_025N"), r_fem_1_whi = "P12I_027N", - r_fem_1_bla = "P12B_027N", + r_fem_1_bla = "P12J_027N", r_fem_1_his = "P12H_027N", - r_fem_1_asi = c("P12D_027N", "P12E_027N"), - r_fem_1_oth = c("P12C_027N", "P12F_027N", "P12G_027N"), + r_fem_1_asi = c("P12L_027N", "P12M_027N"), + r_fem_1_oth = c("P12K_027N", "P12N_027N", "P12O_027N"), r_fem_2_whi = "P12I_028N", - r_fem_2_bla = "P12B_028N", + r_fem_2_bla = "P12J_028N", r_fem_2_his = "P12H_028N", - r_fem_2_asi = c("P12D_028N", "P12E_028N"), - r_fem_2_oth = c("P12C_028N", "P12F_028N", "P12G_028N"), + r_fem_2_asi = c("P12L_028N", "P12M_028N"), + r_fem_2_oth = c("P12K_028N", "P12N_028N", "P12O_028N"), r_fem_3_whi = "P12I_029N", - r_fem_3_bla = "P12B_029N", + r_fem_3_bla = "P12J_029N", r_fem_3_his = "P12H_029N", - r_fem_3_asi = c("P12D_029N", "P12E_029N"), - r_fem_3_oth = c("P12C_029N", "P12F_029N", "P12G_029N"), + r_fem_3_asi = c("P12L_029N", "P12M_029N"), + r_fem_3_oth = c("P12K_029N", "P12N_029N", "P12O_029N"), r_fem_4_whi = "P12I_030N", - r_fem_4_bla = "P12B_030N", + r_fem_4_bla = "P12J_030N", r_fem_4_his = "P12H_030N", - r_fem_4_asi = c("P12D_030N", "P12E_030N"), - r_fem_4_oth = c("P12C_030N", "P12F_030N", "P12G_030N"), + r_fem_4_asi = c("P12L_030N", "P12M_030N"), + r_fem_4_oth = c("P12K_030N", "P12N_030N", "P12O_030N"), r_fem_5_whi = "P12I_031N", - r_fem_5_bla = "P12B_031N", + r_fem_5_bla = "P12J_031N", r_fem_5_his = "P12H_031N", - r_fem_5_asi = c("P12D_031N", "P12E_031N"), - r_fem_5_oth = c("P12C_031N", "P12F_031N", "P12G_031N"), + r_fem_5_asi = c("P12L_031N", "P12M_031N"), + r_fem_5_oth = c("P12K_031N", "P12N_031N", "P12O_031N"), r_fem_6_whi = "P12I_032N", - r_fem_6_bla = "P12B_032N", + r_fem_6_bla = "P12J_032N", r_fem_6_his = "P12H_032N", - r_fem_6_asi = c("P12D_032N", "P12E_032N"), - r_fem_6_oth = c("P12C_032N", "P12F_032N", "P12G_032N"), + r_fem_6_asi = c("P12L_032N", "P12M_032N"), + r_fem_6_oth = c("P12K_032N", "P12N_032N", "P12O_032N"), r_fem_7_whi = "P12I_033N", - r_fem_7_bla = "P12B_033N", + r_fem_7_bla = "P12J_033N", r_fem_7_his = "P12H_033N", - r_fem_7_asi = c("P12D_033N", "P12E_033N"), - r_fem_7_oth = c("P12C_033N", "P12F_033N", "P12G_033N"), + r_fem_7_asi = c("P12L_033N", "P12M_033N"), + r_fem_7_oth = c("P12K_033N", "P12N_033N", "P12O_033N"), r_fem_8_whi = "P12I_034N", - r_fem_8_bla = "P12B_034N", + r_fem_8_bla = "P12J_034N", r_fem_8_his = "P12H_034N", - r_fem_8_asi = c("P12D_034N", "P12E_034N"), - r_fem_8_oth = c("P12C_034N", "P12F_034N", "P12G_034N"), + r_fem_8_asi = c("P12L_034N", "P12M_034N"), + r_fem_8_oth = c("P12K_034N", "P12N_034N", "P12O_034N"), r_fem_9_whi = "P12I_035N", - r_fem_9_bla = "P12B_035N", + r_fem_9_bla = "P12J_035N", r_fem_9_his = "P12H_035N", - r_fem_9_asi = c("P12D_035N", "P12E_035N"), - r_fem_9_oth = c("P12C_035N", "P12F_035N", "P12G_035N"), + r_fem_9_asi = c("P12L_035N", "P12M_035N"), + r_fem_9_oth = c("P12K_035N", "P12N_035N", "P12O_035N"), r_fem_10_whi = "P12I_036N", - r_fem_10_bla = "P12B_036N", + r_fem_10_bla = "P12J_036N", r_fem_10_his = "P12H_036N", - r_fem_10_asi = c("P12D_036N", "P12E_036N"), - r_fem_10_oth = c("P12C_036N", "P12F_036N", "P12G_036N"), + r_fem_10_asi = c("P12L_036N", "P12M_036N"), + r_fem_10_oth = c("P12K_036N", "P12N_036N", "P12O_036N"), r_fem_11_whi = "P12I_037N", - r_fem_11_bla = "P12B_037N", + r_fem_11_bla = "P12J_037N", r_fem_11_his = "P12H_037N", - r_fem_11_asi = c("P12D_037N", "P12E_037N"), - r_fem_11_oth = c("P12C_037N", "P12F_037N", "P12G_037N"), + r_fem_11_asi = c("P12L_037N", "P12M_037N"), + r_fem_11_oth = c("P12K_037N", "P12N_037N", "P12O_037N"), r_fem_12_whi = "P12I_038N", - r_fem_12_bla = "P12B_038N", + r_fem_12_bla = "P12J_038N", r_fem_12_his = "P12H_038N", - r_fem_12_asi = c("P12D_038N", "P12E_038N"), - r_fem_12_oth = c("P12C_038N", "P12F_038N", "P12G_038N"), + r_fem_12_asi = c("P12L_038N", "P12M_038N"), + r_fem_12_oth = c("P12K_038N", "P12N_038N", "P12O_038N"), r_fem_13_whi = "P12I_039N", - r_fem_13_bla = "P12B_039N", + r_fem_13_bla = "P12J_039N", r_fem_13_his = "P12H_039N", - r_fem_13_asi = c("P12D_039N", "P12E_039N"), - r_fem_13_oth = c("P12C_039N", "P12F_039N", "P12G_039N"), + r_fem_13_asi = c("P12L_039N", "P12M_039N"), + r_fem_13_oth = c("P12K_039N", "P12N_039N", "P12O_039N"), r_fem_14_whi = "P12I_040N", - r_fem_14_bla = "P12B_040N", + r_fem_14_bla = "P12J_040N", r_fem_14_his = "P12H_040N", - r_fem_14_asi = c("P12D_040N", "P12E_040N"), - r_fem_14_oth = c("P12C_040N", "P12F_040N", "P12G_040N"), + r_fem_14_asi = c("P12L_040N", "P12M_040N"), + r_fem_14_oth = c("P12K_040N", "P12N_040N", "P12O_040N"), r_fem_15_whi = "P12I_041N", - r_fem_15_bla = "P12B_041N", + r_fem_15_bla = "P12J_041N", r_fem_15_his = "P12H_041N", - r_fem_15_asi = c("P12D_041N", "P12E_041N"), - r_fem_15_oth = c("P12C_041N", "P12F_041N", "P12G_041N"), + r_fem_15_asi = c("P12L_041N", "P12M_041N"), + r_fem_15_oth = c("P12K_041N", "P12N_041N", "P12O_041N"), r_fem_16_whi = "P12I_042N", - r_fem_16_bla = "P12B_042N", + r_fem_16_bla = "P12J_042N", r_fem_16_his = "P12H_042N", - r_fem_16_asi = c("P12D_042N", "P12E_042N"), - r_fem_16_oth = c("P12C_042N", "P12F_042N", "P12G_042N"), + r_fem_16_asi = c("P12L_042N", "P12M_042N"), + r_fem_16_oth = c("P12K_042N", "P12N_042N", "P12O_042N"), r_fem_17_whi = "P12I_043N", - r_fem_17_bla = "P12B_043N", + r_fem_17_bla = "P12J_043N", r_fem_17_his = "P12H_043N", - r_fem_17_asi = c("P12D_043N", "P12E_043N"), - r_fem_17_oth = c("P12C_043N", "P12F_043N", "P12G_043N"), + r_fem_17_asi = c("P12L_043N", "P12M_043N"), + r_fem_17_oth = c("P12K_043N", "P12N_043N", "P12O_043N"), r_fem_18_whi = "P12I_044N", - r_fem_18_bla = "P12B_044N", + r_fem_18_bla = "P12J_044N", r_fem_18_his = "P12H_044N", - r_fem_18_asi = c("P12D_044N", "P12E_044N"), - r_fem_18_oth = c("P12C_044N", "P12F_044N", "P12G_044N"), + r_fem_18_asi = c("P12L_044N", "P12M_044N"), + r_fem_18_oth = c("P12K_044N", "P12N_044N", "P12O_044N"), r_fem_19_whi = "P12I_045N", - r_fem_19_bla = "P12B_045N", + r_fem_19_bla = "P12J_045N", r_fem_19_his = "P12H_045N", - r_fem_19_asi = c("P12D_045N", "P12E_045N"), - r_fem_19_oth = c("P12C_045N", "P12F_045N", "P12G_045N"), + r_fem_19_asi = c("P12L_045N", "P12M_045N"), + r_fem_19_oth = c("P12K_045N", "P12N_045N", "P12O_045N"), r_fem_20_whi = "P12I_046N", - r_fem_20_bla = "P12B_046N", + r_fem_20_bla = "P12J_046N", r_fem_20_his = "P12H_046N", - r_fem_20_asi = c("P12D_046N", "P12E_046N"), - r_fem_20_oth = c("P12C_046N", "P12F_046N", "P12G_046N"), + r_fem_20_asi = c("P12L_046N", "P12M_046N"), + r_fem_20_oth = c("P12K_046N", "P12N_046N", "P12O_046N"), r_fem_21_whi = "P12I_047N", - r_fem_21_bla = "P12B_047N", + r_fem_21_bla = "P12J_047N", r_fem_21_his = "P12H_047N", - r_fem_21_asi = c("P12D_047N", "P12E_047N"), - r_fem_21_oth = c("P12C_047N", "P12F_047N", "P12G_047N"), + r_fem_21_asi = c("P12L_047N", "P12M_047N"), + r_fem_21_oth = c("P12K_047N", "P12N_047N", "P12O_047N"), r_fem_22_whi = "P12I_048N", - r_fem_22_bla = "P12B_048N", + r_fem_22_bla = "P12J_048N", r_fem_22_his = "P12H_048N", - r_fem_22_asi = c("P12D_048N", "P12E_048N"), - r_fem_22_oth = c("P12C_048N", "P12F_048N", "P12G_048N"), + r_fem_22_asi = c("P12L_048N", "P12M_048N"), + r_fem_22_oth = c("P12K_048N", "P12N_048N", "P12O_048N"), r_fem_23_whi = "P12I_049N", - r_fem_23_bla = "P12B_049N", + r_fem_23_bla = "P12J_049N", r_fem_23_his = "P12H_049N", - r_fem_23_asi = c("P12D_049N", "P12E_049N"), - r_fem_23_oth = c("P12C_049N", "P12F_049N", "P12G_049N") + r_fem_23_asi = c("P12L_049N", "P12M_049N"), + r_fem_23_oth = c("P12K_049N", "P12N_049N", "P12O_049N") ) ) ) diff --git a/tests/testthat/test-get_census_data.R b/tests/testthat/test-get_census_data.R index f43eece..598890e 100644 --- a/tests/testthat/test-get_census_data.R +++ b/tests/testthat/test-get_census_data.R @@ -67,5 +67,14 @@ if(Sys.getenv("CENSUS_API_KEY") != "") { )) expect_named(r$DC, c("state", "age", "sex", "year", "zcta")) expect_true(all(r$DC$zcta$state == "DC")) + + r <- suppressMessages(get_census_data( + key = NULL, + state = "DC", + census.geo = "zcta", + year = "2010" + )) + expect_named(r$DC, c("state", "age", "sex", "year", "zcta")) + expect_true(all(r$DC$zcta$state == "DC")) }) }