pacman::p_load(haven, dplyr, tidyr, survey, janitor, magrittr, ggplot2)

Indikator Rata-rata Lama Sekolah (Mean Years of Schooling, MYS) adalah rata-rata jumlah tahun yang ditempuh oleh penduduk berumur 15 tahun ke atas untuk menempuh semua jenjang pendidikan yang pernah dijalani. Untuk mereka yang tamat SD diperhitungkan lama sekolah selama 6 tahun, tamat SMP diperhitungkan lama sekolah selama 9 tahun, tamat SM diperhitungkan lama sekolah selama 12 tahun tanpa memperhitungkan apakah pernah tinggal kelas atau tidak (BPS, 2023).

df <- filter(datakor, R407 %in% (15:98))

1 Manipulasi working variable

Untuk meringankan beban kerja dari proses penghitungan (running data), perlu dilakukan filter variabel yang digunakan.

Terkait indikator ini, variabel yang digunakan antara lain: kode provinsi, kode kab/kota, kode renumbering NKS, kode renumbering no urut ruta, klasifikasi kota desa, umur (tahun), variabel-variabel pendidikan, dan penimbang.

2 Calculating indicator

Statistik dari indikator MYS adalah rata-rata jumlah tahun yang telah diselesaikan oleh penduduk usia 15 tahun ke atas dalam menjalani pendidikan formal.

df <- df %>%
  select(R101, R102, WI1, WI2, R105, R407, R610:R614, FWT) %>%
  mutate(kabu = R101*100 + R102) %>%
  rename(prov="R101", kab="R102") %>%
  mutate(psu=prov*10000000 + kab*100000 + WI1, 
         denom=100, 
         strata= prov*10000000 + kab*100000 + R105,
         nasional = 99)
df <- df %>%
  mutate(
  JENJANG = case_when(
    R612 >= 1 & R612 <= 5 ~ 1,
    R612 >= 6 & R612 <= 10 ~ 2,
    R612 >= 11 & R612 <= 17 ~ 3,
    R612 == 18 ~ 4,
    R612 == 19 ~ 5,
    R612 == 20 ~ 6,
    R612 == 21 ~ 7,
    R612 == 23 ~ 8,
    R612 == 24 ~ 9,
    R612 == 22 ~ 10),
  IJASAH = case_when(
    R614 == 25 ~ 1,
    R614 >= 1 & R614 <= 5 ~ 2,
    R614 >= 6 & R614 <= 10 ~ 3,
    R614 >= 11 & R614 <= 17 ~ 4,
    R614 == 18 ~ 5,
    R614 == 19 ~ 6,
    R614 == 20 ~ 7,
    R614 == 21 ~ 8,
    R614 == 23 ~ 9,
    R614 == 24 ~ 10,
    R614 == 22 ~ 11),
  TAMT = case_when(
    IJASAH == 1 ~ 0,
    IJASAH == 2 ~ 6,
    IJASAH == 3 ~ 9,
    IJASAH == 4 ~ 12,
    IJASAH == 5 ~ 14,
    IJASAH == 6 ~ 15,
    IJASAH == 7 ~ 16,
    IJASAH == 8 ~ 17,
    IJASAH == 9 ~ 19,
    IJASAH == 10 ~ 22,
    IJASAH == 11 ~ 18)
  )

Kebutuhan dari indikator ini diantaranya adalah:

  1. Menggambarkan kompetensi SDM

    Secara umum, indikator menunjukkan tingkat keterampilan dan kompetensi penduduk suatu negara, yang dapat dilihat sebagai proksi dari aspek kuantitatif dan kualitatif dari stok sumber daya manusia. Nilai yang relatif tinggi menunjukkan besarnya proporsi penduduk dewasa menurut tingkat pendidikan tertinggi yang dicapai atau diselesaikan, dan mencerminkan sistem pendidikan yang berkinerja baik (UNESCO, 2009).

  2. Evaluasi Kebijakan Pendidikan

    MYS digunakan untuk merefleksikan kemampuan masyarakat dalam mengakses pendidikan, khususnya pendidikan formal (BPS, 2022).

  3. Analisis Ketimpangan Pendidikan

    Indikator ini digunakan untuk menganalisis ketimpangan pendidikan antarwilayah atau antargender.

3 Set Design Sampling

susenas.design<- svydesign(id=~psu, strata=~strata, 
                            data = df, weights=~FWT)
summary(susenas.design$prob)
##      Min.   1st Qu.    Median      Mean   3rd Qu.      Max. 
## 0.0002617 0.0032893 0.0073179 0.0227792 0.0185233 0.4888525

4 Set Output Statistics

options(survey.adjust.domain.lonely=TRUE)
options(survey.lonely.psu="adjust")

#Nasional
resultN <- svyby(
  formula = ~MYS15, 
  denom = ~denom, 
  by = ~nasional, 
  design = susenas.design, 
  deff = TRUE, 
  svyratio, 
  vartype = c("se", "ci", "ci", "cv", "cvpct", "var"))

resultN[is.na(resultN)] <- 0

resultN <- resultN %>%
  mutate(
    theta = round(`MYS15/denom` * 100, 2),
    SE = round(`se.MYS15/denom` * 100, 2),
    VAR = round(SE^2, 2),
    CI_LOWER = round(`ci_l` * 100, 2),
    CI_UPPER = round(`ci_u` * 100, 2),
    RSE = round(`cv%`, 2),
    DEFF = round(DEff, 2))

resultN$CI_LOWER[resultN$CI_LOWER < 0] <- 0
options(survey.adjust.domain.lonely=TRUE)
options(survey.lonely.psu="adjust")

#Provinsi
resultP <- svyby(
  formula = ~MYS15, 
  denom = ~denom, 
  by = ~prov, 
  design = susenas.design, 
  deff = TRUE, 
  svyratio, 
  vartype = c("se", "ci", "ci", "cv", "cvpct", "var"))

resultP[is.na(resultP)] <- 0

resultP <- resultP %>%
  mutate(
    theta = round(`MYS15/denom` * 100, 2),
    SE = round(`se.MYS15/denom` * 100, 2),
    VAR = round(SE^2, 2),
    CI_LOWER = round(`ci_l` * 100, 2),
    CI_UPPER = round(`ci_u` * 100, 2),
    RSE = round(`cv%`, 2),
    DEFF = round(DEff, 2))

resultP$CI_LOWER[resultP$CI_LOWER < 0] <- 0
options(survey.adjust.domain.lonely=TRUE)
options(survey.lonely.psu="adjust")

#Kabupaten/kota
resultK <- svyby(
  formula = ~MYS15, 
  denom = ~denom, 
  by = ~kabu, 
  design = susenas.design, 
  deff = TRUE, 
  svyratio, 
  vartype = c("se", "ci", "ci", "cv", "cvpct", "var"))

resultK[is.na(resultK)] <- 0

resultK <- resultK %>%
  mutate(
    theta = round(`MYS15/denom` * 100, 2),
    SE = round(`se.MYS15/denom` * 100, 2),
    VAR = round(SE^2, 2),
    CI_LOWER = round(`ci_l` * 100, 2),
    CI_UPPER = round(`ci_u` * 100, 2),
    RSE = round(`cv%`, 2),
    DEFF = round(DEff, 2))

resultK$CI_LOWER[resultK$CI_LOWER < 0] <- 0

5 Level Estimate

Kualitas presisi hasil estimasi suatu survei bisa diamati dari nilai RSE yang dihasilkan. Kesalahan sampling dari beberapa estimasi harus digunakan secara hati-hati. Untuk estimasi yang berdasarkan jumlah kasus yang kecil, kesalahan relatif cenderung besar.

Batasan nilai RSE (Australian Bureau Statistics):

  1. Jika \(RSE \le 25\%\), estimasi bersifat presisi.

  2. Jika \(25\% < RSE \le 50\%\), estimasi perlu dilakukan dengan hati-hati.

  3. Jika \(RSE > 50\%\), estimasi dianggap sangat tidak presisi.

5.1 Estimation by national

outputN = resultN[,c(1,10:16)]
names(outputN) = c("Nasional","Estimasi","SE","VAR","CI LOWER","CI UPPER","RSE","DEFF")
outputN[order(outputN$Nasional, decreasing = FALSE), ]
##    Nasional Estimasi   SE VAR CI LOWER CI UPPER  RSE DEFF
## 99       99     9.13 0.01   0      9.1     9.16 0.15 8.37

5.2 Estimation by province

outputP = resultP[,c(1,10:16)]
names(outputP) = c("Prov","Estimasi","SE","VAR","CI LOWER","CI UPPER","RSE","DEFF")
outputP[order(outputP$Prov, decreasing = FALSE), ]
##    Prov Estimasi   SE  VAR CI LOWER CI UPPER  RSE  DEFF
## 11   11     9.89 0.04 0.00     9.80     9.98 0.45  3.70
## 12   12    10.08 0.05 0.00     9.98    10.17 0.48  7.58
## 13   13     9.60 0.06 0.00     9.47     9.72 0.66  5.94
## 14   14     9.61 0.06 0.00     9.48     9.73 0.67  5.09
## 15   15     9.17 0.06 0.00     9.04     9.30 0.71  3.96
## 16   16     8.90 0.06 0.00     8.78     9.03 0.72  6.09
## 17   17     9.36 0.08 0.01     9.20     9.51 0.83  4.16
## 18   18     8.72 0.05 0.00     8.61     8.82 0.63  4.19
## 19   19     8.66 0.10 0.01     8.47     8.86 1.14  4.71
## 21   21    10.52 0.10 0.01    10.33    10.71 0.92  7.35
## 31   31    11.42 0.07 0.00    11.28    11.56 0.63  5.38
## 32   32     9.17 0.04 0.00     9.08     9.25 0.48  7.10
## 33   33     8.44 0.03 0.00     8.38     8.50 0.37  3.86
## 34   34    10.17 0.09 0.01    10.00    10.34 0.85  3.75
## 35   35     8.53 0.04 0.00     8.46     8.60 0.43  5.02
## 36   36     9.48 0.09 0.01     9.30     9.66 0.95  9.17
## 51   51     9.75 0.07 0.00     9.61     9.89 0.74  4.55
## 52   52     8.39 0.08 0.01     8.24     8.54 0.90  3.79
## 53   53     8.31 0.05 0.00     8.21     8.41 0.62  4.39
## 61   61     8.17 0.07 0.00     8.04     8.30 0.81  4.65
## 62   62     9.07 0.08 0.01     8.92     9.22 0.85  6.62
## 63   63     8.95 0.07 0.00     8.81     9.09 0.82  5.73
## 64   64    10.17 0.08 0.01    10.02    10.32 0.77  5.93
## 65   65     9.54 0.11 0.01     9.32     9.76 1.18  4.85
## 71   71     9.94 0.06 0.00     9.81    10.06 0.64  5.72
## 72   72     9.22 0.08 0.01     9.06     9.37 0.86  6.78
## 73   73     9.13 0.07 0.00     9.00     9.27 0.75  8.27
## 74   74     9.62 0.08 0.01     9.47     9.78 0.82  6.93
## 75   75     8.48 0.10 0.01     8.28     8.69 1.21  4.32
## 76   76     8.48 0.11 0.01     8.28     8.69 1.24  4.16
## 81   81    10.38 0.07 0.00    10.24    10.52 0.70  6.10
## 82   82     9.61 0.07 0.00     9.48     9.75 0.74  4.22
## 91   91    10.15 0.07 0.00    10.01    10.28 0.69  4.76
## 94   94     7.35 0.09 0.01     7.17     7.52 1.23 10.08

Visualisasi

ggplot(data = resultP, aes(x = prov)) +
  geom_line(aes(y = RSE, color = "RSE"), size = 1) +                   # Garis untuk RSE
  geom_line(aes(y = theta, color = "Theta"), size = 1) +               # Garis untuk Theta
  geom_hline(yintercept = 25, color = "red", linetype = "dashed") +    # Garis horizontal pada y = 25
  labs(
    title = "Provinsi vs Nilai RSE dan Theta",
    x = "Provinsi",
    y = "Nilai"
  ) +
  scale_color_manual(name = "Legend", values = c("RSE" = "orange", "Theta" = "blue")) +
  theme_minimal()

5.3 Estimation by district

outputK = resultK[,c(1,10:16)]
names(outputK) = c("Kako","Estimasi","SE","VAR","CI LOWER","CI UPPER","RSE","DEFF")
(outputK[order(outputK$Kako, decreasing = FALSE), ])
##      Kako Estimasi   SE  VAR CI LOWER CI UPPER   RSE  DEFF
## 1101 1101    10.05 0.19 0.04     9.69    10.42  1.84  2.63
## 1102 1102     9.07 0.21 0.04     8.65     9.48  2.33  3.00
## 1103 1103     8.77 0.23 0.05     8.32     9.21  2.60  3.76
## 1104 1104    10.28 0.13 0.02    10.03    10.53  1.24  1.45
## 1105 1105     9.01 0.12 0.01     8.76     9.25  1.37  1.83
## 1106 1106     9.90 0.18 0.03     9.54    10.26  1.86  2.80
## 1107 1107    10.30 0.23 0.05     9.85    10.75  2.22  4.10
## 1108 1108    10.52 0.20 0.04    10.13    10.91  1.90  4.06
## 1109 1109     9.42 0.19 0.04     9.03     9.80  2.07  3.62
## 1110 1110     9.64 0.17 0.03     9.32     9.97  1.72  2.36
## 1111 1111     9.39 0.15 0.02     9.09     9.69  1.63  2.43
## 1112 1112     9.20 0.18 0.03     8.84     9.55  1.96  2.47
## 1113 1113     8.63 0.27 0.07     8.11     9.15  3.07  3.60
## 1114 1114     9.75 0.18 0.03     9.39    10.10  1.86  2.56
## 1115 1115     9.35 0.20 0.04     8.95     9.75  2.17  3.05
## 1116 1116     8.99 0.20 0.04     8.59     9.38  2.24  2.67
## 1117 1117    10.29 0.20 0.04     9.89    10.69  1.98  3.47
## 1118 1118     9.85 0.17 0.03     9.51    10.19  1.77  2.27
## 1171 1171    12.65 0.15 0.02    12.36    12.94  1.17  2.83
## 1172 1172    11.01 0.19 0.04    10.64    11.37  1.69  2.77
## 1173 1173    11.29 0.21 0.04    10.88    11.70  1.86  3.81
## 1174 1174    11.29 0.20 0.04    10.89    11.69  1.80  4.21
## 1175 1175     9.00 0.20 0.04     8.62     9.39  2.17  2.40
## 1201 1201     7.71 0.31 0.10     7.11     8.31  3.96  5.35
## 1202 1202     9.32 0.18 0.03     8.96     9.67  1.94  3.50
## 1203 1203     9.88 0.23 0.05     9.44    10.33  2.30  5.40
## 1204 1204     9.11 0.22 0.05     8.68     9.55  2.41  4.14
## 1205 1205    10.20 0.20 0.04     9.81    10.59  1.96  3.98
## 1206 1206    10.58 0.19 0.04    10.20    10.96  1.84  3.47
## 1207 1207     9.79 0.21 0.04     9.38    10.19  2.11  4.65
## 1208 1208     9.24 0.16 0.03     8.93     9.55  1.71  2.89
## 1209 1209     9.90 0.18 0.03     9.56    10.25  1.77  3.81
## 1210 1210     9.96 0.16 0.03     9.65    10.28  1.62  2.37
## 1211 1211    10.16 0.15 0.02     9.86    10.46  1.51  2.67
## 1212 1212    10.46 0.18 0.03    10.10    10.81  1.74  5.83
## 1213 1213     9.22 0.15 0.02     8.92     9.53  1.68  2.86
## 1214 1214     7.51 0.21 0.04     7.10     7.93  2.82  2.33
## 1215 1215    10.06 0.20 0.04     9.67    10.44  1.95  3.09
## 1216 1216     9.93 0.23 0.05     9.48    10.39  2.34  3.72
## 1217 1217     9.42 0.24 0.06     8.96     9.88  2.51  3.70
## 1218 1218     9.30 0.17 0.03     8.98     9.63  1.77  2.88
## 1219 1219     9.18 0.20 0.04     8.79     9.58  2.20  3.96
## 1220 1220     9.80 0.22 0.05     9.37    10.23  2.25  5.79
## 1221 1221     9.82 0.19 0.04     9.45    10.19  1.91  3.32
## 1222 1222     9.17 0.26 0.07     8.67     9.67  2.78  6.36
## 1223 1223     9.42 0.16 0.03     9.10     9.75  1.75  2.47
## 1224 1224     7.80 0.22 0.05     7.36     8.24  2.87  3.07
## 1225 1225     8.03 0.23 0.05     7.58     8.47  2.83  2.60
## 1271 1271    10.35 0.23 0.05     9.91    10.79  2.19  5.26
## 1272 1272     9.89 0.27 0.07     9.36    10.42  2.73  6.35
## 1273 1273    11.59 0.14 0.02    11.32    11.86  1.19  2.43
## 1274 1274    10.91 0.20 0.04    10.52    11.30  1.82  4.18
## 1275 1275    11.59 0.16 0.03    11.28    11.90  1.37  5.74
## 1276 1276    11.12 0.17 0.03    10.79    11.46  1.53  3.19
## 1277 1277    10.98 0.18 0.03    10.62    11.33  1.66  3.21
## 1278 1278     9.15 0.24 0.06     8.68     9.62  2.62  3.60
## 1301 1301     8.36 0.29 0.08     7.80     8.93  3.45  4.87
## 1302 1302     8.94 0.19 0.04     8.57     9.32  2.13  3.33
## 1303 1303     8.14 0.28 0.08     7.60     8.68  3.40  5.27
## 1304 1304     9.13 0.23 0.05     8.68     9.59  2.54  4.08
## 1305 1305     9.30 0.27 0.07     8.77     9.82  2.90  5.52
## 1306 1306     8.90 0.21 0.04     8.49     9.31  2.36  3.59
## 1307 1307     9.44 0.21 0.04     9.04     9.85  2.17  4.07
## 1308 1308     8.53 0.28 0.08     7.97     9.08  3.33  6.26
## 1309 1309     8.36 0.24 0.06     7.89     8.84  2.90  4.34
## 1310 1310     9.13 0.22 0.05     8.70     9.56  2.40  3.60
## 1311 1311     9.01 0.17 0.03     8.68     9.35  1.91  2.42
## 1312 1312     9.51 0.22 0.05     9.09     9.94  2.27  4.30
## 1371 1371    11.27 0.19 0.04    10.90    11.65  1.70  5.22
## 1372 1372    11.12 0.20 0.04    10.73    11.50  1.75  3.42
## 1373 1373    10.24 0.24 0.06     9.77    10.71  2.36  4.29
## 1374 1374    11.71 0.16 0.03    11.39    12.03  1.38  2.68
## 1375 1375    11.37 0.19 0.04    11.00    11.74  1.67  3.42
## 1376 1376    10.87 0.17 0.03    10.53    11.21  1.60  2.83
## 1377 1377    10.85 0.23 0.05    10.39    11.30  2.15  4.93
## 1401 1401     9.17 0.22 0.05     8.75     9.60  2.37  3.93
## 1402 1402     8.72 0.22 0.05     8.28     9.15  2.56  4.37
## 1403 1403     7.83 0.26 0.07     7.32     8.34  3.35  6.92
## 1404 1404     9.09 0.17 0.03     8.76     9.42  1.83  2.67
## 1405 1405    10.06 0.19 0.04     9.70    10.42  1.84  3.50
## 1406 1406     9.58 0.14 0.02     9.31     9.85  1.46  2.72
## 1407 1407     9.24 0.16 0.03     8.92     9.57  1.77  2.94
## 1408 1408     9.68 0.19 0.04     9.32    10.04  1.91  3.74
## 1409 1409     9.02 0.19 0.04     8.66     9.38  2.06  3.67
## 1410 1410     8.49 0.22 0.05     8.06     8.91  2.56  3.44
## 1471 1471    11.73 0.18 0.03    11.38    12.08  1.52  5.44
## 1473 1473    10.18 0.21 0.04     9.78    10.59  2.04  4.50
## 1501 1501     8.66 0.21 0.04     8.25     9.06  2.38  2.68
## 1502 1502     8.81 0.30 0.09     8.22     9.39  3.40  8.40
## 1503 1503     8.52 0.24 0.06     8.04     8.99  2.87  4.87
## 1504 1504     8.54 0.19 0.04     8.16     8.91  2.23  3.11
## 1505 1505     9.17 0.19 0.04     8.80     9.55  2.10  3.99
## 1506 1506     8.21 0.18 0.03     7.87     8.56  2.16  2.63
## 1507 1507     8.73 0.19 0.04     8.36     9.10  2.16  3.39
## 1508 1508     8.40 0.20 0.04     8.01     8.80  2.41  3.40
## 1509 1509     8.79 0.16 0.03     8.48     9.10  1.81  2.45
## 1571 1571    11.38 0.13 0.02    11.12    11.64  1.15  2.34
## 1572 1572    10.35 0.23 0.05     9.89    10.80  2.24  3.01
## 1601 1601     9.19 0.18 0.03     8.84     9.54  1.95  3.10
## 1602 1602     7.57 0.19 0.04     7.20     7.94  2.50  4.00
## 1603 1603     8.53 0.19 0.04     8.16     8.90  2.22  3.99
## 1604 1604     8.89 0.17 0.03     8.55     9.23  1.96  3.03
## 1605 1605     8.00 0.20 0.04     7.61     8.39  2.47  3.29
## 1606 1606     8.29 0.20 0.04     7.89     8.68  2.43  3.81
## 1607 1607     7.69 0.27 0.07     7.17     8.21  3.47  6.96
## 1608 1608     8.53 0.16 0.03     8.21     8.85  1.91  2.74
## 1609 1609     8.43 0.18 0.03     8.08     8.78  2.13  3.25
## 1610 1610     8.58 0.22 0.05     8.15     9.01  2.57  5.24
## 1611 1611     7.91 0.18 0.03     7.55     8.26  2.28  2.39
## 1612 1612     7.72 0.27 0.07     7.19     8.26  3.55  4.22
## 1613 1613     8.22 0.24 0.06     7.75     8.68  2.88  4.06
## 1671 1671    11.01 0.15 0.02    10.71    11.31  1.39  3.66
## 1672 1672    10.48 0.19 0.04    10.10    10.86  1.84  3.33
## 1673 1673     9.63 0.17 0.03     9.30     9.96  1.74  2.19
## 1674 1674     9.56 0.24 0.06     9.08    10.04  2.56  4.82
## 1701 1701     9.55 0.19 0.04     9.18     9.91  1.96  2.66
## 1702 1702     9.27 0.20 0.04     8.87     9.67  2.21  3.19
## 1703 1703     8.49 0.20 0.04     8.11     8.88  2.34  2.79
## 1704 1704     8.74 0.17 0.03     8.41     9.07  1.92  2.00
## 1705 1705     8.52 0.18 0.03     8.15     8.88  2.17  2.71
## 1706 1706     8.96 0.18 0.03     8.60     9.31  2.03  2.65
## 1707 1707     8.68 0.25 0.06     8.20     9.17  2.85  3.93
## 1708 1708     8.67 0.17 0.03     8.33     9.01  2.00  1.93
## 1709 1709     8.05 0.35 0.12     7.36     8.74  4.36  6.42
## 1771 1771    11.71 0.21 0.04    11.30    12.12  1.78  4.59
## 1801 1801     8.70 0.24 0.06     8.22     9.17  2.78  5.10
## 1802 1802     7.62 0.16 0.03     7.31     7.93  2.06  2.47
## 1803 1803     8.32 0.19 0.04     7.95     8.70  2.31  4.11
## 1804 1804     8.52 0.14 0.02     8.25     8.79  1.61  2.11
## 1805 1805     8.26 0.15 0.02     7.95     8.56  1.87  2.86
## 1806 1806     8.74 0.16 0.03     8.44     9.05  1.78  3.02
## 1807 1807     8.11 0.16 0.03     7.80     8.41  1.94  2.45
## 1808 1808     8.03 0.28 0.08     7.48     8.57  3.49  6.49
## 1809 1809     8.49 0.20 0.04     8.11     8.87  2.31  3.68
## 1810 1810     8.73 0.16 0.03     8.42     9.04  1.79  2.43
## 1811 1811     7.43 0.24 0.06     6.96     7.91  3.27  4.27
## 1812 1812     8.22 0.22 0.05     7.78     8.65  2.70  4.04
## 1813 1813     8.87 0.18 0.03     8.53     9.22  1.97  3.44
## 1871 1871    10.99 0.22 0.05    10.56    11.41  1.97  5.84
## 1872 1872    10.85 0.25 0.06    10.37    11.34  2.28  4.80
## 1901 1901     8.87 0.21 0.04     8.45     9.28  2.38  3.47
## 1902 1902     9.21 0.21 0.04     8.79     9.63  2.32  3.63
## 1903 1903     7.82 0.28 0.08     7.28     8.36  3.52  4.81
## 1904 1904     7.45 0.26 0.07     6.93     7.97  3.56  4.67
## 1905 1905     7.60 0.25 0.06     7.12     8.08  3.23  4.28
## 1906 1906     9.13 0.23 0.05     8.68     9.59  2.54  4.56
## 1971 1971    10.54 0.27 0.07    10.01    11.07  2.55  5.50
## 2101 2101     9.18 0.21 0.04     8.76     9.59  2.30  4.64
## 2102 2102     9.72 0.22 0.05     9.29    10.15  2.26  4.59
## 2103 2103     9.43 0.21 0.04     9.01     9.84  2.25  2.62
## 2104 2104     7.91 0.27 0.07     7.38     8.45  3.45  4.21
## 2105 2105     8.45 0.24 0.06     7.99     8.92  2.79  3.54
## 2171 2171    11.01 0.13 0.02    10.75    11.28  1.22  3.70
## 2172 2172    10.89 0.21 0.04    10.48    11.29  1.91  3.88
## 3101 3101     9.18 0.24 0.06     8.72     9.65  2.60  4.34
## 3171 3171    11.72 0.15 0.02    11.42    12.02  1.30  4.58
## 3172 3172    11.80 0.12 0.01    11.56    12.04  1.03  3.44
## 3173 3173    11.46 0.16 0.03    11.14    11.78  1.43  4.20
## 3174 3174    11.21 0.18 0.03    10.86    11.55  1.58  6.22
## 3175 3175    10.78 0.16 0.03    10.47    11.09  1.47  4.78
## 3201 3201     8.86 0.18 0.03     8.51     9.22  2.02  6.45
## 3202 3202     7.84 0.16 0.03     7.53     8.15  2.01  4.53
## 3203 3203     7.38 0.15 0.02     7.09     7.67  2.01  3.96
## 3204 3204     9.36 0.18 0.03     9.00     9.72  1.96  6.64
## 3205 3205     8.25 0.15 0.02     7.95     8.55  1.84  4.91
## 3206 3206     8.31 0.19 0.04     7.94     8.68  2.28  6.68
## 3207 3207     8.45 0.16 0.03     8.13     8.76  1.88  4.58
## 3208 3208     8.23 0.24 0.06     7.75     8.70  2.93  7.14
## 3209 3209     8.30 0.16 0.03     8.00     8.61  1.88  3.36
## 3210 3210     7.96 0.14 0.02     7.69     8.23  1.72  2.77
## 3211 3211     9.01 0.18 0.03     8.66     9.35  1.95  4.91
## 3212 3212     7.51 0.19 0.04     7.14     7.89  2.56  3.53
## 3213 3213     8.00 0.15 0.02     7.70     8.29  1.86  2.67
## 3214 3214     8.36 0.21 0.04     7.95     8.76  2.48  5.26
## 3215 3215     8.50 0.23 0.05     8.05     8.94  2.67  7.80
## 3216 3216     9.87 0.19 0.04     9.49    10.25  1.97  7.39
## 3217 3217     8.41 0.16 0.03     8.10     8.72  1.89  5.26
## 3218 3218     8.37 0.14 0.02     8.10     8.64  1.64  2.32
## 3271 3271    10.63 0.23 0.05    10.17    11.08  2.18  7.29
## 3272 3272    10.45 0.16 0.03    10.14    10.77  1.54  2.96
## 3273 3273    11.09 0.19 0.04    10.73    11.45  1.67  6.95
## 3274 3274    10.47 0.26 0.07     9.96    10.98  2.49  6.42
## 3275 3275    11.75 0.18 0.03    11.40    12.10  1.52  6.94
## 3276 3276    11.55 0.17 0.03    11.22    11.89  1.49  5.60
## 3277 3277    11.34 0.19 0.04    10.97    11.71  1.67  5.71
## 3278 3278     9.73 0.19 0.04     9.36    10.10  1.96  4.70
## 3279 3279     8.99 0.21 0.04     8.57     9.40  2.37  3.85
## 3301 3301     7.85 0.14 0.02     7.57     8.14  1.84  2.98
## 3302 3302     8.32 0.16 0.03     8.01     8.64  1.93  3.47
## 3303 3303     7.78 0.16 0.03     7.47     8.09  2.01  3.39
## 3304 3304     7.37 0.16 0.03     7.05     7.68  2.19  3.42
## 3305 3305     8.11 0.16 0.03     7.79     8.43  2.02  3.52
## 3306 3306     8.77 0.16 0.03     8.46     9.08  1.81  2.65
## 3307 3307     7.20 0.18 0.03     6.85     7.56  2.50  4.93
## 3308 3308     8.11 0.16 0.03     7.79     8.42  2.00  3.68
## 3309 3309     8.48 0.15 0.02     8.18     8.78  1.81  2.92
## 3310 3310     9.54 0.15 0.02     9.25     9.82  1.53  2.58
## 3311 3311    10.04 0.17 0.03     9.71    10.37  1.69  3.06
## 3312 3312     7.99 0.15 0.02     7.69     8.29  1.90  3.12
## 3313 3313     9.34 0.18 0.03     8.99     9.69  1.93  4.18
## 3314 3314     8.27 0.17 0.03     7.94     8.60  2.03  3.10
## 3315 3315     7.80 0.12 0.01     7.55     8.04  1.59  2.54
## 3316 3316     7.58 0.18 0.03     7.24     7.93  2.33  2.95
## 3317 3317     8.23 0.17 0.03     7.90     8.56  2.04  3.43
## 3318 3318     8.25 0.16 0.03     7.93     8.57  1.98  3.03
## 3319 3319     9.60 0.17 0.03     9.26     9.94  1.81  3.65
## 3320 3320     8.69 0.16 0.03     8.37     9.01  1.88  3.40
## 3321 3321     8.84 0.15 0.02     8.55     9.13  1.66  2.76
## 3322 3322     8.58 0.19 0.04     8.21     8.94  2.16  3.93
## 3323 3323     7.89 0.13 0.02     7.64     8.13  1.60  2.27
## 3324 3324     8.28 0.16 0.03     7.96     8.59  1.94  3.06
## 3325 3325     7.60 0.17 0.03     7.27     7.93  2.20  3.69
## 3326 3326     7.97 0.19 0.04     7.61     8.34  2.33  4.80
## 3327 3327     7.17 0.17 0.03     6.83     7.51  2.42  3.97
## 3328 3328     7.89 0.17 0.03     7.55     8.23  2.19  3.67
## 3329 3329     7.09 0.18 0.03     6.73     7.44  2.58  4.37
## 3371 3371    11.11 0.20 0.04    10.71    11.50  1.80  3.72
## 3372 3372    11.02 0.14 0.02    10.75    11.29  1.24  2.21
## 3373 3373    11.40 0.20 0.04    11.01    11.80  1.76  3.65
## 3374 3374    10.87 0.19 0.04    10.50    11.24  1.74  5.20
## 3375 3375     9.49 0.17 0.03     9.15     9.84  1.84  3.28
## 3376 3376    10.02 0.18 0.03     9.67    10.36  1.75  2.65
## 3401 3401     9.43 0.16 0.03     9.12     9.74  1.67  2.77
## 3402 3402    10.15 0.22 0.05     9.72    10.57  2.13  4.87
## 3403 3403     7.70 0.16 0.03     7.40     8.01  2.04  2.89
## 3404 3404    11.22 0.15 0.02    10.93    11.51  1.32  2.79
## 3471 3471    12.23 0.18 0.03    11.87    12.59  1.50  3.70
## 3501 3501     8.19 0.19 0.04     7.81     8.57  2.35  4.95
## 3502 3502     8.01 0.19 0.04     7.63     8.38  2.39  3.38
## 3503 3503     8.03 0.16 0.03     7.73     8.34  1.94  2.92
## 3504 3504     8.92 0.17 0.03     8.59     9.25  1.88  3.47
## 3505 3505     8.04 0.18 0.03     7.69     8.40  2.27  3.72
## 3506 3506     8.62 0.14 0.02     8.34     8.89  1.64  2.76
## 3507 3507     8.21 0.18 0.03     7.85     8.57  2.24  5.11
## 3508 3508     7.57 0.18 0.03     7.22     7.91  2.32  3.62
## 3509 3509     6.97 0.21 0.04     6.55     7.39  3.08  5.36
## 3510 3510     8.17 0.17 0.03     7.84     8.50  2.05  3.18
## 3511 3511     6.95 0.24 0.06     6.48     7.42  3.43  4.20
## 3512 3512     7.44 0.23 0.05     6.98     7.89  3.11  4.24
## 3513 3513     6.91 0.23 0.05     6.46     7.35  3.30  5.05
## 3514 3514     7.84 0.18 0.03     7.49     8.20  2.32  4.36
## 3515 3515    10.85 0.15 0.02    10.56    11.13  1.35  4.10
## 3516 3516     9.40 0.16 0.03     9.09     9.70  1.65  3.10
## 3517 3517     9.06 0.15 0.02     8.77     9.35  1.62  2.86
## 3518 3518     8.58 0.17 0.03     8.25     8.91  1.94  2.89
## 3519 3519     8.20 0.16 0.03     7.88     8.52  1.99  2.68
## 3520 3520     8.80 0.17 0.03     8.47     9.13  1.90  3.02
## 3521 3521     8.20 0.22 0.05     7.77     8.64  2.68  4.58
## 3522 3522     7.88 0.17 0.03     7.54     8.21  2.17  3.42
## 3523 3523     7.87 0.14 0.02     7.59     8.15  1.83  2.57
## 3524 3524     8.67 0.14 0.02     8.40     8.94  1.58  2.31
## 3525 3525    10.30 0.18 0.03     9.95    10.65  1.73  4.81
## 3526 3526     6.57 0.22 0.05     6.14     7.01  3.41  5.41
## 3527 3527     5.96 0.17 0.03     5.63     6.29  2.82  2.88
## 3528 3528     7.93 0.24 0.06     7.46     8.39  3.01  4.89
## 3529 3529     6.34 0.22 0.05     5.92     6.77  3.42  3.63
## 3571 3571    11.00 0.18 0.03    10.66    11.35  1.60  3.12
## 3572 3572    10.81 0.23 0.05    10.36    11.26  2.13  3.76
## 3573 3573    11.19 0.29 0.08    10.62    11.76  2.60  8.68
## 3574 3574     9.81 0.33 0.11     9.17    10.44  3.32  7.40
## 3575 3575     9.93 0.28 0.08     9.38    10.48  2.81  5.73
## 3576 3576    11.22 0.18 0.03    10.87    11.58  1.61  3.14
## 3577 3577    11.63 0.21 0.04    11.23    12.04  1.78  4.14
## 3578 3578    10.82 0.21 0.04    10.41    11.23  1.92  7.22
## 3579 3579    10.11 0.22 0.05     9.68    10.55  2.17  4.57
## 3601 3601     7.57 0.20 0.04     7.17     7.96  2.67  6.40
## 3602 3602     7.20 0.16 0.03     6.90     7.51  2.17  3.98
## 3603 3603     9.33 0.21 0.04     8.91     9.75  2.29  9.00
## 3604 3604     8.38 0.17 0.03     8.05     8.72  2.06  4.82
## 3671 3671    10.97 0.21 0.04    10.56    11.37  1.89  7.62
## 3672 3672    10.44 0.21 0.04    10.04    10.85  1.98  6.31
## 3673 3673     8.89 0.23 0.05     8.44     9.33  2.54  6.35
## 3674 3674    11.46 0.21 0.04    11.05    11.86  1.81  7.07
## 5101 5101     8.97 0.19 0.04     8.59     9.34  2.12  3.18
## 5102 5102     9.56 0.22 0.05     9.12     9.99  2.31  5.26
## 5103 5103    11.21 0.16 0.03    10.89    11.52  1.43  3.64
## 5104 5104    10.20 0.21 0.04     9.79    10.62  2.08  5.07
## 5105 5105     8.86 0.23 0.05     8.41     9.31  2.59  3.80
## 5106 5106     8.07 0.25 0.06     7.58     8.57  3.14  4.97
## 5107 5107     7.00 0.21 0.04     6.59     7.41  3.01  3.36
## 5108 5108     8.07 0.15 0.02     7.77     8.37  1.91  2.46
## 5171 5171    11.51 0.16 0.03    11.20    11.82  1.38  4.43
## 5201 5201     7.72 0.21 0.04     7.32     8.12  2.66  3.00
## 5202 5202     7.37 0.18 0.03     7.03     7.72  2.40  2.27
## 5203 5203     7.87 0.20 0.04     7.47     8.27  2.59  3.34
## 5204 5204     9.09 0.19 0.04     8.72     9.46  2.08  3.03
## 5205 5205     9.43 0.18 0.03     9.07     9.79  1.93  2.18
## 5206 5206     9.15 0.17 0.03     8.82     9.48  1.85  2.04
## 5207 5207     9.24 0.21 0.04     8.83     9.65  2.28  2.80
## 5208 5208     7.20 0.23 0.05     6.76     7.65  3.13  2.72
## 5271 5271     9.88 0.26 0.07     9.38    10.39  2.60  5.24
## 5272 5272    10.85 0.28 0.08    10.29    11.40  2.61  5.47
## 5301 5301     7.54 0.24 0.06     7.06     8.02  3.23  4.45
## 5302 5302     8.22 0.23 0.05     7.78     8.67  2.78  3.93
## 5303 5303     8.07 0.21 0.04     7.65     8.49  2.64  3.77
## 5304 5304     7.52 0.21 0.04     7.11     7.93  2.77  3.88
## 5305 5305     8.50 0.21 0.04     8.08     8.92  2.51  3.35
## 5306 5306     7.62 0.18 0.03     7.26     7.98  2.42  2.59
## 5307 5307     8.55 0.22 0.05     8.13     8.97  2.52  3.73
## 5308 5308     8.08 0.23 0.05     7.62     8.54  2.89  4.07
## 5309 5309     8.30 0.24 0.06     7.83     8.77  2.89  4.88
## 5310 5310     7.42 0.21 0.04     7.00     7.83  2.87  3.61
## 5311 5311     8.65 0.30 0.09     8.07     9.24  3.45  6.00
## 5312 5312     8.97 0.18 0.03     8.61     9.33  2.06  3.21
## 5313 5313     8.12 0.23 0.05     7.67     8.56  2.81  4.41
## 5314 5314     8.20 0.19 0.04     7.82     8.57  2.31  2.46
## 5315 5315     8.37 0.22 0.05     7.94     8.80  2.62  3.70
## 5316 5316     7.82 0.31 0.10     7.21     8.43  3.98  6.86
## 5317 5317     7.03 0.21 0.04     6.62     7.43  2.93  3.29
## 5318 5318     8.46 0.17 0.03     8.12     8.79  2.00  2.65
## 5319 5319     8.30 0.23 0.05     7.85     8.75  2.75  4.74
## 5320 5320     7.60 0.27 0.07     7.07     8.13  3.55  4.84
## 5321 5321     7.43 0.19 0.04     7.06     7.81  2.59  2.31
## 5371 5371    11.30 0.19 0.04    10.92    11.67  1.69  4.49
## 6101 6101     6.97 0.16 0.03     6.65     7.28  2.33  2.23
## 6102 6102     8.27 0.20 0.04     7.88     8.65  2.39  2.45
## 6103 6103     7.91 0.18 0.03     7.56     8.25  2.24  2.57
## 6104 6104     7.81 0.32 0.10     7.18     8.44  4.12  7.89
## 6105 6105     7.86 0.26 0.07     7.35     8.37  3.31  6.35
## 6106 6106     8.06 0.20 0.04     7.66     8.45  2.49  4.19
## 6107 6107     8.13 0.24 0.06     7.65     8.60  2.99  4.44
## 6108 6108     8.37 0.27 0.07     7.83     8.90  3.28  6.19
## 6109 6109     7.59 0.32 0.10     6.97     8.21  4.15  7.00
## 6110 6110     7.96 0.34 0.12     7.30     8.62  4.24  6.65
## 6111 6111     6.94 0.23 0.05     6.48     7.39  3.37  3.26
## 6112 6112     7.75 0.20 0.04     7.35     8.15  2.63  3.44
## 6171 6171    10.30 0.22 0.05     9.87    10.74  2.17  4.68
## 6172 6172     8.46 0.28 0.08     7.92     9.00  3.26  5.39
## 6201 6201     9.03 0.21 0.04     8.62     9.44  2.31  3.50
## 6202 6202     8.14 0.25 0.06     7.65     8.63  3.05  5.10
## 6203 6203     8.53 0.22 0.05     8.10     8.97  2.60  4.68
## 6204 6204     9.53 0.20 0.04     9.13     9.92  2.12  3.68
## 6205 6205     8.98 0.22 0.05     8.54     9.42  2.50  4.23
## 6206 6206     8.67 0.23 0.05     8.22     9.12  2.66  4.31
## 6207 6207     8.98 0.23 0.05     8.52     9.44  2.60  4.36
## 6208 6208     8.14 0.23 0.05     7.69     8.59  2.82  4.28
## 6209 6209     9.21 0.22 0.05     8.77     9.65  2.44  5.69
## 6210 6210     8.59 0.22 0.05     8.16     9.01  2.53  4.25
## 6211 6211     9.49 0.21 0.04     9.07     9.91  2.26  4.52
## 6212 6212    10.03 0.29 0.08     9.46    10.59  2.87  7.17
## 6213 6213     8.42 0.23 0.05     7.96     8.87  2.75  4.67
## 6271 6271    11.58 0.31 0.10    10.97    12.19  2.68  8.81
## 6301 6301     8.28 0.18 0.03     7.93     8.62  2.14  2.91
## 6302 6302     7.97 0.23 0.05     7.52     8.41  2.86  4.49
## 6303 6303     8.44 0.25 0.06     7.96     8.93  2.95  5.49
## 6304 6304     8.20 0.25 0.06     7.72     8.68  2.99  4.75
## 6305 6305     8.44 0.25 0.06     7.94     8.94  3.01  5.00
## 6306 6306     8.22 0.23 0.05     7.76     8.68  2.85  4.16
## 6307 6307     8.63 0.21 0.04     8.23     9.04  2.40  4.23
## 6308 6308     8.20 0.24 0.06     7.73     8.67  2.91  4.15
## 6309 6309     9.41 0.21 0.04     9.00     9.82  2.20  4.30
## 6310 6310     8.76 0.22 0.05     8.33     9.19  2.50  4.33
## 6311 6311     8.71 0.18 0.03     8.36     9.06  2.06  2.56
## 6371 6371    10.42 0.24 0.06     9.96    10.89  2.27  6.02
## 6372 6372    10.92 0.29 0.08    10.34    11.49  2.68  7.30
## 6401 6401     9.33 0.22 0.05     8.89     9.76  2.39  4.46
## 6402 6402     9.16 0.20 0.04     8.76     9.56  2.21  2.98
## 6403 6403     9.39 0.18 0.03     9.03     9.75  1.94  3.75
## 6404 6404     9.52 0.22 0.05     9.08     9.96  2.36  4.70
## 6405 6405     9.66 0.27 0.07     9.13    10.19  2.78  5.79
## 6409 6409     8.85 0.19 0.04     8.48     9.23  2.17  2.79
## 6411 6411     9.15 0.19 0.04     8.78     9.51  2.04  2.97
## 6471 6471    10.79 0.16 0.03    10.47    11.11  1.51  3.74
## 6472 6472    11.38 0.17 0.03    11.04    11.72  1.52  4.47
## 6474 6474    11.02 0.31 0.10    10.42    11.63  2.79  8.47
## 6501 6501     9.45 0.24 0.06     8.97     9.92  2.55  3.87
## 6502 6502     9.14 0.28 0.08     8.59     9.70  3.09  5.55
## 6503 6503     9.19 0.28 0.08     8.64     9.73  3.03  4.68
## 6504 6504     8.64 0.22 0.05     8.22     9.07  2.52  3.61
## 6571 6571    10.48 0.19 0.04    10.11    10.86  1.82  4.00
## 7101 7101     8.66 0.17 0.03     8.33     9.00  1.97  3.40
## 7102 7102    10.38 0.19 0.04    10.00    10.76  1.86  3.99
## 7103 7103     8.91 0.23 0.05     8.45     9.37  2.62  4.10
## 7104 7104    10.02 0.19 0.04     9.64    10.40  1.92  3.15
## 7105 7105     9.62 0.19 0.04     9.25     9.98  1.95  3.38
## 7106 7106    10.44 0.21 0.04    10.03    10.85  2.00  4.06
## 7107 7107     8.69 0.22 0.05     8.27     9.11  2.47  3.80
## 7108 7108     9.51 0.18 0.03     9.15     9.87  1.93  3.24
## 7109 7109     9.10 0.18 0.03     8.76     9.45  1.93  2.75
## 7110 7110     8.35 0.20 0.04     7.95     8.74  2.43  3.95
## 7111 7111     8.56 0.19 0.04     8.19     8.93  2.20  3.43
## 7171 7171    11.13 0.22 0.05    10.70    11.56  1.98  6.11
## 7172 7172    10.02 0.21 0.04     9.61    10.44  2.10  4.59
## 7173 7173    11.01 0.20 0.04    10.61    11.41  1.85  4.19
## 7174 7174     9.99 0.22 0.05     9.56    10.41  2.17  5.03
## 7201 7201     8.93 0.27 0.07     8.39     9.47  3.08  5.82
## 7202 7202     8.94 0.24 0.06     8.48     9.41  2.65  5.18
## 7203 7203     9.07 0.23 0.05     8.62     9.53  2.56  4.70
## 7204 7204     9.58 0.18 0.03     9.21     9.94  1.93  3.14
## 7205 7205     8.11 0.26 0.07     7.61     8.61  3.16  6.69
## 7206 7206     9.06 0.29 0.08     8.50     9.62  3.16  6.29
## 7207 7207     9.27 0.19 0.04     8.89     9.64  2.06  3.35
## 7208 7208     8.43 0.22 0.05     7.99     8.86  2.63  5.01
## 7209 7209     8.85 0.18 0.03     8.48     9.21  2.09  2.98
## 7210 7210     9.24 0.25 0.06     8.75     9.74  2.75  5.78
## 7211 7211     8.66 0.24 0.06     8.19     9.12  2.76  4.32
## 7212 7212     9.15 0.21 0.04     8.75     9.56  2.26  3.81
## 7271 7271    11.41 0.25 0.06    10.91    11.90  2.22  5.82
## 7301 7301     8.69 0.34 0.12     8.02     9.37  3.97  5.79
## 7302 7302     8.64 0.27 0.07     8.12     9.16  3.09  4.87
## 7303 7303     7.73 0.28 0.08     7.18     8.28  3.63  3.93
## 7304 7304     7.66 0.28 0.08     7.10     8.22  3.70  5.74
## 7305 7305     8.15 0.23 0.05     7.69     8.60  2.85  3.52
## 7306 7306     8.82 0.22 0.05     8.38     9.26  2.53  4.11
## 7307 7307     8.05 0.27 0.07     7.52     8.58  3.35  4.70
## 7308 7308     8.31 0.29 0.08     7.75     8.88  3.48  5.76
## 7309 7309     9.01 0.26 0.07     8.50     9.51  2.84  6.29
## 7310 7310     8.88 0.32 0.10     8.25     9.52  3.65  6.73
## 7311 7311     7.96 0.20 0.04     7.58     8.34  2.46  3.77
## 7312 7312     8.66 0.37 0.14     7.94     9.38  4.24  8.70
## 7313 7313     8.02 0.22 0.05     7.59     8.46  2.77  3.60
## 7314 7314     8.53 0.23 0.05     8.08     8.99  2.73  4.31
## 7315 7315     8.90 0.26 0.07     8.40     9.41  2.89  5.00
## 7316 7316     9.08 0.23 0.05     8.63     9.53  2.54  3.23
## 7317 7317     9.21 0.20 0.04     8.82     9.60  2.18  3.94
## 7318 7318     8.90 0.26 0.07     8.38     9.42  2.97  4.35
## 7322 7322     8.58 0.18 0.03     8.21     8.94  2.15  2.99
## 7325 7325     8.92 0.19 0.04     8.54     9.30  2.17  2.94
## 7326 7326     8.93 0.19 0.04     8.56     9.31  2.12  2.29
## 7371 7371    11.42 0.26 0.07    10.92    11.93  2.26 10.44
## 7372 7372    10.79 0.23 0.05    10.34    11.24  2.14  4.13
## 7373 7373    11.22 0.23 0.05    10.77    11.67  2.05  4.96
## 7401 7401     8.97 0.32 0.10     8.34     9.60  3.59  7.18
## 7402 7402     8.71 0.30 0.09     8.12     9.31  3.48  5.18
## 7403 7403     9.82 0.21 0.04     9.41    10.23  2.12  3.68
## 7404 7404     9.83 0.22 0.05     9.41    10.25  2.19  4.14
## 7405 7405     8.94 0.23 0.05     8.50     9.39  2.54  3.99
## 7406 7406     8.31 0.33 0.11     7.67     8.95  3.94  5.83
## 7407 7407     8.76 0.34 0.12     8.10     9.42  3.84  6.52
## 7408 7408     9.20 0.31 0.10     8.60     9.80  3.35  6.33
## 7409 7409     9.22 0.23 0.05     8.76     9.67  2.52  3.30
## 7410 7410     9.78 0.22 0.05     9.35    10.20  2.23  4.21
## 7411 7411     8.78 0.36 0.13     8.07     9.49  4.11  8.51
## 7412 7412     9.61 0.27 0.07     9.07    10.14  2.84  4.89
## 7413 7413     8.71 0.29 0.08     8.13     9.28  3.36  4.54
## 7414 7414     7.34 0.22 0.05     6.91     7.77  2.97  2.30
## 7415 7415     8.12 0.23 0.05     7.67     8.57  2.83  3.75
## 7471 7471    12.03 0.26 0.07    11.53    12.54  2.14  6.16
## 7472 7472    11.09 0.23 0.05    10.64    11.53  2.05  4.01
## 7501 7501     7.76 0.26 0.07     7.25     8.28  3.37  4.50
## 7502 7502     8.22 0.23 0.05     7.78     8.67  2.76  4.13
## 7503 7503     8.00 0.23 0.05     7.55     8.44  2.83  3.34
## 7504 7504     8.54 0.17 0.03     8.20     8.89  2.05  2.40
## 7505 7505     7.49 0.20 0.04     7.10     7.88  2.66  2.32
## 7571 7571    10.23 0.23 0.05     9.78    10.69  2.27  3.88
## 7601 7601     9.67 0.26 0.07     9.15    10.19  2.74  4.78
## 7602 7602     8.21 0.24 0.06     7.75     8.68  2.87  4.17
## 7603 7603     8.57 0.28 0.08     8.02     9.12  3.28  3.82
## 7604 7604     8.48 0.22 0.05     8.06     8.91  2.56  3.27
## 7605 7605     8.28 0.20 0.04     7.90     8.67  2.38  2.47
## 7606 7606     8.04 0.19 0.04     7.66     8.42  2.40  2.28
## 8101 8101    10.22 0.19 0.04     9.86    10.59  1.82  4.46
## 8102 8102     9.92 0.15 0.02     9.62    10.22  1.55  3.08
## 8103 8103    10.22 0.21 0.04     9.81    10.64  2.08  5.41
## 8104 8104     9.36 0.18 0.03     9.02     9.71  1.90  2.79
## 8105 8105     9.30 0.19 0.04     8.94     9.67  2.01  3.85
## 8106 8106     9.68 0.22 0.05     9.26    10.11  2.22  4.68
## 8107 8107     9.12 0.16 0.03     8.81     9.43  1.74  2.44
## 8108 8108     9.44 0.21 0.04     9.03     9.86  2.26  4.90
## 8109 8109     8.79 0.34 0.12     8.13     9.45  3.82  8.16
## 8171 8171    11.81 0.14 0.02    11.53    12.09  1.20  2.97
## 8172 8172    10.78 0.17 0.03    10.44    11.12  1.60  3.20
## 8201 8201     8.82 0.22 0.05     8.38     9.26  2.53  4.31
## 8202 8202     9.49 0.31 0.10     8.88    10.10  3.28  6.61
## 8203 8203     9.47 0.18 0.03     9.12     9.82  1.88  3.17
## 8204 8204     8.66 0.15 0.02     8.37     8.95  1.70  1.82
## 8205 8205     9.07 0.21 0.04     8.66     9.48  2.32  4.32
## 8206 8206     9.08 0.25 0.06     8.58     9.57  2.77  5.15
## 8207 8207     8.19 0.25 0.06     7.69     8.68  3.10  4.99
## 8208 8208     8.72 0.21 0.04     8.32     9.13  2.37  3.59
## 8271 8271    11.84 0.13 0.02    11.59    12.09  1.08  2.22
## 8272 8272    10.19 0.20 0.04     9.80    10.59  1.97  3.77
## 9101 9101    10.58 0.18 0.03    10.23    10.94  1.71  2.86
## 9102 9102     9.93 0.21 0.04     9.53    10.34  2.08  4.05
## 9103 9103     8.85 0.42 0.18     8.03     9.67  4.71 11.03
## 9104 9104    10.26 0.17 0.03     9.93    10.59  1.62  2.89
## 9105 9105    10.07 0.19 0.04     9.70    10.44  1.87  2.67
## 9106 9106     8.87 0.29 0.08     8.30     9.44  3.26  6.65
## 9107 9107     9.49 0.21 0.04     9.07     9.90  2.23  3.57
## 9108 9108     8.73 0.16 0.03     8.42     9.05  1.83  2.54
## 9109 9109     8.75 0.33 0.11     8.11     9.39  3.75  3.65
## 9110 9110    10.20 0.27 0.07     9.67    10.73  2.66  3.50
## 9111 9111     8.86 0.36 0.13     8.14     9.57  4.10  8.12
## 9112 9112     7.61 0.29 0.08     7.03     8.18  3.85  4.25
## 9171 9171    11.37 0.15 0.02    11.07    11.66  1.32  2.92
## 9401 9401    10.07 0.22 0.05     9.65    10.50  2.15  3.40
## 9402 9402     6.74 0.49 0.24     5.77     7.70  7.33 10.18
## 9403 9403    10.31 0.21 0.04     9.90    10.72  2.04  2.97
## 9404 9404    10.50 0.19 0.04    10.12    10.87  1.82  2.89
## 9408 9408     9.78 0.30 0.09     9.19    10.37  3.09  6.33
## 9409 9409    10.35 0.19 0.04     9.97    10.73  1.88  3.14
## 9410 9410     4.32 0.50 0.25     3.33     5.31 11.67 11.46
## 9411 9411     4.53 0.61 0.37     3.33     5.72 13.50 12.79
## 9412 9412    10.42 0.24 0.06     9.96    10.89  2.27  4.26
## 9413 9413     9.49 0.33 0.11     8.84    10.15  3.50  6.92
## 9414 9414     7.65 0.32 0.10     7.01     8.29  4.25  5.90
## 9415 9415     7.15 0.42 0.18     6.32     7.97  5.90  7.93
## 9416 9416     5.00 0.33 0.11     4.36     5.65  6.59  8.94
## 9417 9417     4.73 0.33 0.11     4.07     5.38  7.04  5.49
## 9418 9418     3.00 0.45 0.20     2.11     3.89 15.10 11.32
## 9419 9419    10.01 0.30 0.09     9.42    10.59  2.97  5.88
## 9420 9420     8.79 0.27 0.07     8.27     9.31  3.02  4.25
## 9426 9426     9.22 0.28 0.08     8.68     9.76  2.99  5.28
## 9427 9427     9.64 0.26 0.07     9.14    10.14  2.65  4.46
## 9428 9428     9.23 0.32 0.10     8.61     9.85  3.42  4.81
## 9429 9429     5.11 0.32 0.10     4.48     5.74  6.31  6.53
## 9430 9430     5.24 0.35 0.12     4.56     5.92  6.62  4.83
## 9431 9431     5.64 0.42 0.18     4.82     6.46  7.39  6.49
## 9432 9432     6.95 0.35 0.12     6.26     7.65  5.10  7.49
## 9433 9433     1.37 0.19 0.04     1.00     1.74 13.86  5.03
## 9434 9434     4.71 0.42 0.18     3.89     5.53  8.85 10.83
## 9435 9435     2.78 0.35 0.12     2.10     3.47 12.56  7.97
## 9436 9436     3.54 0.38 0.14     2.79     4.29 10.84  8.14
## 9471 9471    11.84 0.24 0.06    11.37    12.31  2.01  4.54

Visualisasi

ggplot(data = resultK, aes(x =kabu)) +
  geom_line(aes(y = RSE, color = "RSE"), size = 1) +  
  geom_line(aes(y = theta, color = "Theta"), size = 1) +
  geom_hline(yintercept = 25, color = "red", linetype = "dashed") +
  labs(
    title = "Kab/kota vs Nilai RSE dan Theta",
    x = "Kab/kota",
    y = "Nilai"
  ) +
  scale_color_manual(name = "Legend", values = c("RSE" = "magenta", "Theta" = "purple")) +
  theme_minimal()

Eksport to CSV

write.csv(resultP,"RSE_DIRECT_MYS15_NAS_2023.csv")
write.csv(resultP,"RSE_DIRECT_MYS15_PROV_2023.csv")
write.csv(resultK,"RSE_DIRECT_MYS15_KAKO_2023.csv")
rm(list=ls())  # Clear environment

6 References

[1] BPS. 2023. Statistik Pendidikan 2023. Jakarta: Badan Pusat Statistik.

[2] BPS. 2022. Indeks Pembangunan Manusia 2021. Jakarta: Badan Pusat Statistik.

[3] Confidentiality and relative standard error, Australian Bureau Statistics. https://www.abs.gov.au/statistics/microdata-tablebuilder/tablebuilder/confidentiality-and-relative-standard-error

[4] UNESCO. 2009. Education Indicators Technical guidelines.


Direktorat Statistik Kesejahteraan Rakyat, BPS,