Load Packages

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

Load Dataset

Dilakukan simulasi pada indikator keluhan kesehatan dari data Susenas Maret 2024.

datakor <- read_sav("KOR24GAB_0208.sav")

1 Manipulasi working variable

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

dataku <- datakor %>%
  select(R101, R102, WI1, WI2, PSU, SSU, STRATA, R105, R407, R1102, FWT, fpc1, fpc2)

2 Calculating indicator

Persentase penduduk yang mengalami keluhan kesehatan.

Perbandingan antara jumlah penduduk yang mempunyai keluhan kesehatan dalam sebulan terakhir terhadap jumlah penduduk, dinyatakan dalam satuan persen (%).

Konsep definisi:

Keadaan seseorang yang mempunyai gangguan kesehatan atau kejiwaan, baik karena gangguan/penyakit yang sering dialami, seperti: panas, batuk, pilek, diare, sakit kepala, maupun karena penyakit akut, penyakit kronis (meskipun selama sebulan terakhir tidak mempunyai keluhan), karena kecelakaan, kriminalitas atau keluhan kesehatan lainnya. Keluhan yang dimaksud adalah keluhan fisik maupun psikis yang dialami dalam 1 (satu) bulan terakhir, yang berakhir 1 (satu) hari sebelum pencacahan (BPS, 2023).

Formula:

\(P_{Keluhan}=\frac{Jumlah\,penduduk\,yang\,mengalami\,keluhan\,kesehatan}{Jumlah\,penduduk}\times 100\%\)

indicator <- dataku %>%
  mutate(
    kabu = R101 * 100 + R102,
    jml_pddk = 100,
    keluh = case_when(is.na(R1102) ~ 0, R1102 == 1 ~ 100, R1102 == 5 ~ 0))

3 Penghitungan dengan R

3.1 Set Design Sampling

Desain sampling yang digunakan yaitu Stratified 2 - level Sampling design.

susenas.design<- svydesign(id=~PSU+SSU, strata=~STRATA, 
                            data = indicator, weights=~FWT,fpc = ~fpc1+fpc2, nest=TRUE)
summary(susenas.design$prob)
##      Min.   1st Qu.    Median      Mean   3rd Qu.      Max. 
## 0.0002383 0.0034118 0.0076019 0.0285624 0.0196124 0.8644758

Set Output Statistics

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

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

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

resultP <- resultP %>%
  mutate(
    theta = round(`keluh/jml_pddk` * 100, 2),
    SE = round(`se.keluh/jml_pddk` * 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 = ~keluh, 
  denom = ~jml_pddk, 
  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(`keluh/jml_pddk` * 100, 2),
    SE = round(`se.keluh/jml_pddk` * 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

3.2 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.

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    25.58 0.65 0.42    24.31    26.85  2.54 10.99
## 12   12    23.99 0.59 0.35    22.85    25.14  2.44 14.71
## 13   13    30.10 0.72 0.52    28.69    31.51  2.39 10.66
## 14   14    23.60 0.81 0.66    22.02    25.18  3.42 11.10
## 15   15    23.73 0.84 0.71    22.08    25.37  3.54  9.42
## 16   16    29.89 0.90 0.81    28.11    31.66  3.03 15.39
## 17   17    28.27 1.12 1.25    26.08    30.46  3.96 12.47
## 18   18    29.54 0.85 0.72    27.87    31.22  2.89 12.46
## 19   19    28.80 1.04 1.08    26.76    30.84  3.61  7.27
## 21   21    26.26 2.23 4.97    21.89    30.62  8.48 38.22
## 31   31    24.07 1.05 1.10    22.02    26.13  4.35 11.41
## 32   32    28.88 0.55 0.30    27.81    29.96  1.90 12.40
## 33   33    31.86 0.44 0.19    31.00    32.73  1.38  8.90
## 34   34    28.82 0.96 0.92    26.93    30.70  3.33  5.94
## 35   35    28.75 0.44 0.19    27.90    29.61  1.52  9.73
## 36   36    25.77 1.01 1.02    23.79    27.75  3.93 13.71
## 51   51    19.32 0.75 0.56    17.86    20.79  3.86  8.38
## 52   52    38.88 1.05 1.10    36.82    40.94  2.70 10.67
## 53   53    29.56 0.73 0.53    28.13    31.00  2.48 13.11
## 61   61    24.37 0.80 0.64    22.80    25.95  3.30 11.35
## 62   62    21.37 0.86 0.74    19.69    23.06  4.02 11.82
## 63   63    29.00 0.84 0.71    27.36    30.65  2.89  9.27
## 64   64    20.01 0.92 0.85    18.21    21.81  4.58 11.64
## 65   65    19.93 1.19 1.42    17.60    22.25  5.95  9.34
## 71   71    23.30 0.88 0.77    21.58    25.01  3.76 12.51
## 72   72    22.41 0.88 0.77    20.67    24.14  3.95 12.35
## 73   73    23.43 0.61 0.37    22.23    24.63  2.62 11.78
## 74   74    20.87 0.73 0.53    19.43    22.30  3.51 11.49
## 75   75    33.80 1.27 1.61    31.31    36.29  3.76  9.26
## 76   76    23.33 1.26 1.59    20.85    25.80  5.42 12.52
## 81   81    14.38 0.78 0.61    12.85    15.92  5.45 13.19
## 82   82    13.22 0.68 0.46    11.88    14.55  5.16  8.89
## 91   91    14.18 1.08 1.17    12.06    16.29  7.61 14.35
## 92   92    17.43 1.22 1.49    15.04    19.81  6.99 13.50
## 94   94    14.43 0.92 0.85    12.63    16.23  6.36 12.05
## 95   95    21.24 1.35 1.82    18.60    23.88  6.34  9.88
## 96   96    10.72 1.20 1.44     8.38    13.07 11.16 21.29
## 97   97     8.88 0.64 0.41     7.63    10.13  7.17  7.12

Visualisasi

ggplot(data = resultP, aes(x = R101)) +
  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 - R",
    x = "Provinsi",
    y = "Nilai"
  ) +
  scale_color_manual(name = " ", values = c("RSE" = "orange", "Theta" = "blue")) +
  theme_minimal()

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    18.72 2.63  6.92    13.55    23.88 14.07  9.78
## 1102 1102    15.60 1.78  3.17    12.10    19.09 11.44  5.62
## 1103 1103    13.04 1.82  3.31     9.48    16.60 13.93  6.47
## 1104 1104    24.44 3.17 10.05    18.23    30.65 12.96 10.71
## 1105 1105    29.18 3.48 12.11    22.36    35.99 11.92 14.50
## 1106 1106    31.16 3.51 12.32    24.27    38.04 11.27 11.66
## 1107 1107    27.15 2.53  6.40    22.20    32.11  9.31  6.50
## 1108 1108    24.67 2.12  4.49    20.53    28.82  8.58  6.44
## 1109 1109    31.53 2.36  5.57    26.91    36.15  7.48  6.33
## 1110 1110    33.73 2.29  5.24    29.25    38.22  6.78  5.74
## 1111 1111    20.27 2.63  6.92    15.12    25.42 12.96 10.70
## 1112 1112    30.50 2.34  5.48    25.91    35.09  7.68  5.80
## 1113 1113    19.62 2.63  6.92    14.46    24.77 13.40  8.24
## 1114 1114    22.14 2.24  5.02    17.75    26.54 10.14  6.50
## 1115 1115    23.91 3.36 11.29    17.32    30.51 14.07 11.53
## 1116 1116    39.05 2.62  6.86    33.93    44.18  6.70  5.45
## 1117 1117    23.63 2.83  8.01    18.09    29.17 11.96  8.53
## 1118 1118    25.78 2.61  6.81    20.67    30.89 10.12  6.27
## 1171 1171    26.86 2.57  6.60    21.83    31.89  9.55  7.71
## 1172 1172    23.94 2.60  6.76    18.84    29.03 10.85  6.63
## 1173 1173    29.03 3.54 12.53    22.08    35.97 12.21 12.67
## 1174 1174    21.38 2.01  4.04    17.44    25.32  9.40  5.32
## 1175 1175    16.87 1.76  3.10    13.42    20.32 10.44  5.04
## 1201 1201    31.11 3.08  9.49    25.09    37.14  9.88 10.99
## 1202 1202    16.91 1.89  3.57    13.21    20.61 11.16  6.70
## 1203 1203    23.48 3.09  9.55    17.43    29.54 13.16 11.64
## 1204 1204    21.56 1.65  2.72    18.34    24.79  7.64  3.96
## 1205 1205    13.51 2.31  5.34     8.99    18.03 17.08 11.16
## 1206 1206    16.36 1.94  3.76    12.55    20.17 11.88  5.15
## 1207 1207    14.09 2.10  4.41     9.96    18.21 14.94  8.58
## 1208 1208    27.87 2.40  5.76    23.18    32.57  8.60  7.96
## 1209 1209    24.83 2.30  5.29    20.33    29.33  9.25  7.64
## 1210 1210    64.94 3.61 13.03    57.87    72.00  5.55 11.99
## 1211 1211    19.54 2.43  5.90    14.78    24.30 12.43  8.61
## 1212 1212    20.69 2.11  4.45    16.57    24.82 10.18  9.81
## 1213 1213    29.61 2.57  6.60    24.57    34.66  8.69  8.81
## 1214 1214    12.53 2.12  4.49     8.38    16.68 16.91  9.95
## 1215 1215    16.07 1.97  3.88    12.21    19.93 12.25  5.97
## 1216 1216    32.78 2.66  7.08    27.56    37.99  8.11  6.64
## 1217 1217    28.95 2.61  6.81    23.83    34.07  9.03  6.62
## 1218 1218    25.02 2.88  8.29    19.37    30.67 11.52 11.55
## 1219 1219    30.86 2.72  7.40    25.53    36.19  8.81  8.45
## 1220 1220    17.18 2.34  5.48    12.59    21.78 13.64  8.32
## 1221 1221    32.42 4.10 16.81    24.38    40.46 12.65 16.46
## 1222 1222    19.07 2.24  5.02    14.68    23.47 11.76  7.87
## 1223 1223    24.89 3.67 13.47    17.71    32.08 14.73 16.76
## 1224 1224    29.87 3.94 15.52    22.15    37.58 13.18 17.26
## 1225 1225    33.66 3.38 11.42    27.04    40.28 10.03 12.35
## 1271 1271    27.10 1.62  2.62    23.92    30.29  5.99  2.90
## 1272 1272    25.41 2.30  5.29    20.90    29.92  9.05  5.67
## 1273 1273    24.33 2.80  7.84    18.84    29.82 11.52  9.17
## 1274 1274    25.79 3.11  9.67    19.70    31.88 12.04  9.44
## 1275 1275    22.70 1.89  3.57    18.99    26.40  8.33  7.38
## 1276 1276    25.67 2.58  6.66    20.60    30.73 10.07  6.67
## 1277 1277    32.08 3.00  9.00    26.19    37.97  9.36  9.08
## 1278 1278    35.84 3.53 12.46    28.92    42.75  9.85 13.16
## 1301 1301    30.23 4.69 22.00    21.03    39.43 15.53 19.21
## 1302 1302    26.69 2.22  4.93    22.34    31.03  8.31  6.88
## 1303 1303    24.77 2.33  5.43    20.21    29.33  9.39  6.90
## 1304 1304    29.10 2.96  8.76    23.29    34.91 10.18  9.36
## 1305 1305    41.40 2.73  7.45    36.05    46.76  6.60  7.32
## 1306 1306    37.77 2.78  7.73    32.32    43.22  7.36  8.26
## 1307 1307    28.48 2.24  5.02    24.08    32.87  7.87  6.82
## 1308 1308    30.09 3.42 11.70    23.39    36.79 11.36 12.06
## 1309 1309    28.14 2.98  8.88    22.31    33.97 10.57  9.83
## 1310 1310    31.77 3.40 11.56    25.11    38.42 10.69 10.90
## 1311 1311    32.29 2.60  6.76    27.20    37.39  8.05  7.27
## 1312 1312    24.84 3.10  9.61    18.75    30.92 12.50 12.87
## 1371 1371    30.16 2.08  4.33    26.08    34.23  6.89  6.42
## 1372 1372    29.32 2.34  5.48    24.73    33.92  7.99  5.53
## 1373 1373    26.74 1.94  3.76    22.94    30.55  7.26  3.91
## 1374 1374    31.01 2.86  8.18    25.40    36.62  9.24  7.52
## 1375 1375    23.99 2.18  4.75    19.73    28.26  9.07  5.36
## 1376 1376    29.77 2.30  5.29    25.27    34.28  7.72  5.28
## 1377 1377    39.38 2.77  7.67    33.95    44.81  7.04  7.08
## 1401 1401    12.94 1.53  2.34     9.95    15.94 11.80  4.40
## 1402 1402    27.08 2.80  7.84    21.60    32.57 10.34  9.92
## 1403 1403    34.72 2.95  8.70    28.93    40.50  8.50 10.66
## 1404 1404    12.31 1.79  3.20     8.81    15.81 14.52  7.18
## 1405 1405    24.82 3.40 11.56    18.16    31.48 13.69 14.60
## 1406 1406    19.03 2.06  4.24    14.99    23.08 10.85  7.98
## 1407 1407    18.47 2.40  5.76    13.76    23.18 13.01 10.08
## 1408 1408    13.10 1.78  3.17     9.62    16.58 13.56  7.46
## 1409 1409    35.52 3.17 10.05    29.31    41.72  8.91 12.13
## 1410 1410    22.60 2.50  6.25    17.71    27.50 11.06  7.40
## 1471 1471    28.87 2.44  5.95    24.10    33.65  8.45  9.11
## 1473 1473    19.48 2.02  4.08    15.52    23.43 10.36  6.24
## 1501 1501    37.13 3.18 10.11    30.89    43.37  8.57  8.40
## 1502 1502    23.49 2.78  7.73    18.04    28.95 11.84  9.81
## 1503 1503    23.78 2.36  5.57    19.16    28.39  9.91  6.58
## 1504 1504    20.72 3.14  9.86    14.57    26.87 15.14 13.48
## 1505 1505    19.57 2.71  7.34    14.26    24.88 13.85 11.34
## 1506 1506    21.76 2.77  7.67    16.33    27.20 12.75  8.92
## 1507 1507    13.28 1.79  3.20     9.77    16.80 13.50  6.06
## 1508 1508    27.04 2.98  8.88    21.21    32.88 11.01  9.90
## 1509 1509    30.71 2.88  8.29    25.08    36.35  9.36  9.09
## 1571 1571    23.58 2.10  4.41    19.48    27.69  8.88  6.38
## 1572 1572    21.43 2.04  4.16    17.44    25.43  9.50  4.44
## 1601 1601    39.83 3.09  9.55    33.76    45.89  7.77  9.72
## 1602 1602    28.98 2.95  8.70    23.19    34.77 10.20 11.54
## 1603 1603    30.29 2.98  8.88    24.45    36.12  9.82 12.05
## 1604 1604    34.87 3.05  9.30    28.89    40.84  8.74 10.16
## 1605 1605    25.54 3.01  9.06    19.63    31.44 11.80  9.00
## 1606 1606    22.07 3.09  9.55    16.01    28.13 14.01 13.62
## 1607 1607    21.43 2.92  8.53    15.70    27.15 13.63 14.33
## 1608 1608    23.55 3.56 12.67    16.56    30.53 15.13 13.63
## 1609 1609    45.41 3.66 13.40    38.24    52.58  8.05 13.07
## 1610 1610    42.38 2.58  6.66    37.33    47.43  6.08  6.53
## 1611 1611    34.83 3.70 13.69    27.58    42.09 10.63 11.16
## 1612 1612    37.42 4.23 17.89    29.13    45.71 11.31 14.75
## 1613 1613    24.14 3.10  9.61    18.07    30.21 12.83  9.21
## 1671 1671    27.83 2.86  8.18    22.23    33.43 10.26 13.61
## 1672 1672    23.68 2.67  7.13    18.45    28.91 11.26  8.19
## 1673 1673    27.44 3.48 12.11    20.63    34.26 12.67 11.93
## 1674 1674    25.53 2.66  7.08    20.32    30.74 10.42  7.86
## 1701 1701    33.63 3.61 13.03    26.55    40.72 10.75 11.82
## 1702 1702    24.69 2.77  7.67    19.25    30.12 11.24  8.35
## 1703 1703    28.20 3.08  9.49    22.16    34.24 10.93  9.99
## 1704 1704    26.13 2.55  6.50    21.13    31.13  9.77  6.72
## 1705 1705    35.59 4.46 19.89    26.85    44.32 12.52 16.64
## 1706 1706    21.33 2.66  7.08    16.12    26.54 12.47  9.05
## 1707 1707    18.27 1.49  2.22    15.35    21.19  8.15  2.71
## 1708 1708    34.23 3.60 12.96    27.17    41.28 10.52 10.67
## 1709 1709    30.97 3.28 10.76    24.55    37.40 10.59  9.50
## 1771 1771    28.72 3.33 11.09    22.19    35.25 11.60 13.09
## 1801 1801    29.04 4.84 23.43    19.55    38.52 16.67 21.43
## 1802 1802    30.41 2.71  7.34    25.09    35.73  8.92  8.98
## 1803 1803    23.74 2.64  6.97    18.56    28.92 11.13 10.69
## 1804 1804    26.09 2.61  6.81    20.96    31.21 10.02  9.75
## 1805 1805    34.33 2.46  6.05    29.50    39.15  7.17  7.75
## 1806 1806    21.68 2.17  4.71    17.43    25.93 10.00  7.13
## 1807 1807    45.76 3.90 15.21    38.11    53.40  8.52 14.22
## 1808 1808    32.22 3.31 10.96    25.73    38.72 10.28 11.76
## 1809 1809    35.07 4.68 21.90    25.89    44.24 13.35 21.84
## 1810 1810    34.88 3.26 10.63    28.49    41.27  9.35 11.36
## 1811 1811    31.44 2.46  6.05    26.62    36.27  7.83  5.07
## 1812 1812    33.04 4.05 16.40    25.10    40.97 12.26 14.67
## 1813 1813    21.00 2.67  7.13    15.76    26.25 12.73  8.74
## 1871 1871    24.98 2.58  6.66    19.92    30.03 10.32 10.45
## 1872 1872    27.32 2.05  4.20    23.31    31.34  7.50  4.23
## 1901 1901    30.79 2.79  7.78    25.33    36.26  9.06  7.78
## 1902 1902    26.84 2.64  6.97    21.66    32.01  9.84  6.65
## 1903 1903    30.13 2.34  5.48    25.54    34.72  7.78  5.15
## 1904 1904    26.97 2.62  6.86    21.83    32.11  9.72  6.91
## 1905 1905    22.21 2.36  5.57    17.58    26.84 10.63  6.52
## 1906 1906    31.66 2.47  6.10    26.82    36.49  7.80  5.14
## 1971 1971    32.49 2.83  8.01    26.95    38.03  8.70  7.21
## 2101 2101    20.52 2.67  7.13    15.28    25.76 13.03  9.49
## 2102 2102    18.01 1.95  3.80    14.20    21.83 10.82  5.43
## 2103 2103    17.53 1.99  3.96    13.64    21.42 11.32  5.23
## 2104 2104    16.53 2.20  4.84    12.21    20.85 13.34  6.25
## 2105 2105    15.14 1.96  3.84    11.29    18.99 12.97  6.06
## 2171 2171    28.96 3.26 10.63    22.58    35.35 11.24 15.68
## 2172 2172    27.20 2.67  7.13    21.97    32.43  9.81  7.30
## 3101 3101    35.06 2.74  7.51    29.70    40.42  7.80  6.25
## 3171 3171    19.79 1.72  2.96    16.41    23.16  8.70  7.00
## 3172 3172    18.73 2.05  4.20    14.71    22.76 10.97 10.82
## 3173 3173    26.85 3.39 11.49    20.21    33.49 12.61 14.21
## 3174 3174    28.58 2.51  6.30    23.67    33.50  8.78 11.38
## 3175 3175    29.88 2.09  4.37    25.77    33.98  7.01  7.28
## 3201 3201    28.81 2.39  5.71    24.14    33.49  8.28 12.58
## 3202 3202    31.36 2.34  5.48    26.77    35.95  7.47  8.90
## 3203 3203    37.80 2.23  4.97    33.42    42.17  5.91  7.47
## 3204 3204    25.80 1.91  3.65    22.07    29.54  7.39  7.65
## 3205 3205    37.88 2.40  5.76    33.17    42.59  6.34  9.55
## 3206 3206    33.57 2.78  7.73    28.13    39.01  8.27 10.56
## 3207 3207    43.22 3.05  9.30    37.25    49.20  7.05 11.50
## 3208 3208    38.23 2.76  7.62    32.83    43.63  7.21  8.79
## 3209 3209    39.75 2.80  7.84    34.26    45.24  7.05 11.36
## 3210 3210    36.83 2.88  8.29    31.18    42.48  7.83  9.99
## 3211 3211    28.72 2.05  4.20    24.71    32.74  7.14  5.53
## 3212 3212    32.66 3.42 11.70    25.96    39.36 10.46 15.81
## 3213 3213    19.58 2.07  4.28    15.53    23.64 10.56  7.66
## 3214 3214    32.75 2.52  6.35    27.80    37.70  7.71  8.16
## 3215 3215    24.30 2.51  6.30    19.38    29.22 10.33 11.96
## 3216 3216    20.06 2.13  4.54    15.89    24.24 10.62 11.26
## 3217 3217    30.66 2.06  4.24    26.61    34.70  6.74  6.96
## 3218 3218    47.56 2.98  8.88    41.71    53.41  6.27  7.72
## 3271 3271    27.44 2.58  6.66    22.39    32.49  9.39  9.60
## 3272 3272    29.09 2.15  4.62    24.88    33.30  7.38  5.21
## 3273 3273    20.95 2.15  4.62    16.73    25.17 10.28  9.99
## 3274 3274    29.67 3.22 10.37    23.35    35.98 10.86  9.98
## 3275 3275    23.46 1.87  3.50    19.79    27.13  7.98  7.64
## 3276 3276    19.36 1.61  2.59    16.20    22.53  8.34  6.22
## 3277 3277    24.29 1.98  3.92    20.40    28.18  8.17  5.63
## 3278 3278    30.61 2.64  6.97    25.44    35.79  8.63  8.27
## 3279 3279    25.14 2.48  6.15    20.29    29.99  9.85  6.07
## 3301 3301    28.04 2.13  4.54    23.87    32.21  7.59  7.68
## 3302 3302    28.49 1.95  3.80    24.67    32.32  6.85  6.59
## 3303 3303    38.29 1.98  3.92    34.41    42.16  5.16  5.12
## 3304 3304    28.93 2.37  5.62    24.28    33.58  8.20  7.98
## 3305 3305    55.41 3.13  9.80    49.28    61.53  5.64 12.28
## 3306 3306    34.56 2.67  7.13    29.33    39.78  7.72  7.93
## 3307 3307    38.66 2.68  7.18    33.41    43.91  6.93  8.54
## 3308 3308    34.38 2.10  4.41    30.25    38.50  6.12  6.00
## 3309 3309    30.16 2.49  6.20    25.27    35.04  8.27  8.57
## 3310 3310    34.60 2.23  4.97    30.23    38.97  6.44  6.51
## 3311 3311    25.13 2.31  5.34    20.59    29.66  9.21  7.65
## 3312 3312    23.19 1.80  3.24    19.66    26.72  7.76  5.36
## 3313 3313    20.62 2.13  4.54    16.46    24.79 10.31  7.98
## 3314 3314    26.55 2.54  6.45    21.57    31.53  9.57  8.86
## 3315 3315    31.67 2.25  5.06    27.25    36.08  7.11  7.60
## 3316 3316    36.68 3.11  9.67    30.59    42.77  8.47 10.78
## 3317 3317    32.65 2.80  7.84    27.16    38.15  8.58  9.09
## 3318 3318    28.06 3.03  9.18    22.12    34.00 10.80 12.70
## 3319 3319    32.38 3.05  9.30    26.41    38.35  9.41 11.93
## 3320 3320    28.17 2.66  7.08    22.97    33.38  9.43 10.33
## 3321 3321    30.49 2.64  6.97    25.32    35.66  8.65  9.61
## 3322 3322    28.02 2.09  4.37    23.92    32.13  7.47  6.51
## 3323 3323    31.23 2.60  6.76    26.14    36.33  8.32  9.09
## 3324 3324    31.26 2.38  5.66    26.60    35.93  7.61  7.62
## 3325 3325    28.67 2.92  8.53    22.94    34.40 10.19 11.92
## 3326 3326    33.02 2.57  6.60    27.98    38.06  7.80  9.54
## 3327 3327    40.67 2.34  5.48    36.07    45.26  5.77  7.53
## 3328 3328    32.99 1.86  3.46    29.34    36.63  5.63  4.98
## 3329 3329    39.75 2.41  5.81    35.02    44.47  6.07  8.50
## 3371 3371    25.19 3.14  9.86    19.03    31.35 12.47  8.87
## 3372 3372    32.92 2.33  5.43    28.36    37.47  7.07  5.61
## 3373 3373    29.40 2.27  5.15    24.95    33.85  7.72  4.67
## 3374 3374    23.60 2.03  4.12    19.62    27.59  8.61  7.91
## 3375 3375    39.86 2.69  7.24    34.58    45.14  6.76  6.62
## 3376 3376    25.03 2.23  4.97    20.67    29.40  8.90  5.31
## 3401 3401    30.10 2.29  5.24    25.61    34.59  7.61  6.21
## 3402 3402    31.76 1.94  3.76    27.94    35.57  6.12  5.01
## 3403 3403    23.44 1.68  2.82    20.14    26.74  7.18  4.45
## 3404 3404    28.54 2.01  4.04    24.60    32.48  7.04  5.81
## 3471 3471    30.65 2.31  5.34    26.12    35.17  7.53  5.23
## 3501 3501    36.73 2.95  8.70    30.95    42.52  8.04  9.27
## 3502 3502    39.36 2.57  6.60    34.33    44.39  6.52  7.34
## 3503 3503    28.47 2.10  4.41    24.36    32.59  7.38  5.67
## 3504 3504    30.81 3.09  9.55    24.77    36.86 10.02 12.35
## 3505 3505    36.61 2.99  8.94    30.74    42.47  8.17 10.98
## 3506 3506    30.57 2.81  7.90    25.06    36.08  9.20 12.38
## 3507 3507    28.04 1.95  3.80    24.22    31.86  6.95  7.12
## 3508 3508    24.62 2.31  5.34    20.09    29.14  9.38  8.24
## 3509 3509    30.10 2.16  4.67    25.85    34.34  7.19  8.58
## 3510 3510    32.95 2.36  5.57    28.32    37.58  7.16  8.32
## 3511 3511    35.76 2.79  7.78    30.29    41.23  7.81  8.31
## 3512 3512    26.96 2.56  6.55    21.93    31.98  9.51  7.80
## 3513 3513    32.04 2.86  8.18    26.43    37.65  8.93 11.12
## 3514 3514    20.38 2.04  4.16    16.38    24.37 10.01  8.63
## 3515 3515    23.15 1.86  3.46    19.50    26.80  8.05  7.16
## 3516 3516    31.59 2.85  8.12    26.00    37.18  9.03 11.06
## 3517 3517    36.72 2.10  4.41    32.61    40.82  5.71  5.66
## 3518 3518    34.19 2.12  4.49    30.04    38.33  6.19  5.55
## 3519 3519    26.02 2.50  6.25    21.13    30.91  9.59  7.81
## 3520 3520    31.07 2.44  5.95    26.28    35.87  7.87  7.29
## 3521 3521    29.09 2.59  6.71    24.01    34.18  8.91  7.71
## 3522 3522    21.28 1.86  3.46    17.65    24.92  8.72  6.43
## 3523 3523    34.29 2.64  6.97    29.12    39.47  7.70  9.41
## 3524 3524    27.13 1.89  3.57    23.42    30.84  6.97  5.51
## 3525 3525    25.67 2.19  4.80    21.37    29.96  8.54  7.88
## 3526 3526    15.88 1.76  3.10    12.42    19.33 11.11  6.52
## 3527 3527    34.92 3.11  9.67    28.82    41.02  8.92 11.80
## 3528 3528    28.33 3.01  9.06    22.44    34.22 10.61 12.68
## 3529 3529    22.80 2.66  7.08    17.60    28.01 11.65 10.50
## 3571 3571    24.60 1.75  3.06    21.16    28.04  7.13  3.52
## 3572 3572    23.55 2.15  4.62    19.33    27.77  9.14  4.69
## 3573 3573    24.26 1.95  3.80    20.43    28.09  8.05  5.54
## 3574 3574    29.45 2.83  8.01    23.91    34.98  9.59  7.91
## 3575 3575    34.97 2.91  8.47    29.27    40.67  8.32  7.06
## 3576 3576    27.96 2.47  6.10    23.12    32.80  8.84  5.48
## 3577 3577    33.83 2.96  8.76    28.03    39.62  8.74  6.50
## 3578 3578    27.09 1.93  3.72    23.31    30.87  7.12  7.74
## 3579 3579    28.06 2.95  8.70    22.27    33.85 10.52  8.70
## 3601 3601    25.13 2.16  4.67    20.90    29.36  8.60  7.58
## 3602 3602    26.05 2.42  5.86    21.31    30.80  9.29  9.07
## 3603 3603    28.38 2.57  6.60    23.34    33.43  9.07 13.16
## 3604 3604    29.33 2.97  8.82    23.51    35.15 10.13 14.70
## 3671 3671    24.30 2.02  4.08    20.35    28.25  8.30  7.96
## 3672 3672    26.14 1.96  3.84    22.30    29.99  7.51  4.81
## 3673 3673    31.46 2.23  4.97    27.08    35.84  7.10  6.68
## 3674 3674    16.77 2.00  4.00    12.85    20.69 11.93  9.23
## 5101 5101    15.02 1.82  3.31    11.46    18.59 12.11  5.28
## 5102 5102    28.15 2.66  7.08    22.93    33.36  9.45  9.41
## 5103 5103    17.38 1.94  3.76    13.58    21.17 11.15  7.79
## 5104 5104    17.44 1.55  2.40    14.40    20.47  8.88  5.04
## 5105 5105    10.47 2.11  4.45     6.35    14.60 20.10  9.97
## 5106 5106    24.32 2.00  4.00    20.40    28.23  8.21  5.06
## 5107 5107    29.06 2.79  7.78    23.60    34.53  9.60  9.05
## 5108 5108    20.18 2.25  5.06    15.77    24.58 11.15  8.25
## 5171 5171    14.97 1.64  2.69    11.75    18.19 10.96  6.92
## 5201 5201    38.22 2.76  7.62    32.82    43.62  7.21  8.09
## 5202 5202    45.29 2.90  8.41    39.61    50.97  6.40  8.79
## 5203 5203    38.23 2.61  6.81    33.11    43.34  6.83  7.94
## 5204 5204    35.73 3.08  9.49    29.69    41.77  8.63  9.74
## 5205 5205    34.92 2.82  7.95    29.38    40.45  8.09  7.63
## 5206 5206    40.09 2.75  7.56    34.69    45.49  6.87  7.51
## 5207 5207    26.71 2.71  7.34    21.40    32.02 10.15  7.46
## 5208 5208    40.35 3.46 11.97    33.57    47.12  8.57  9.20
## 5271 5271    31.92 2.89  8.35    26.25    37.59  9.06  9.11
## 5272 5272    51.35 4.01 16.08    43.49    59.21  7.81 13.01
## 5301 5301    12.66 1.63  2.66     9.46    15.85 12.87  5.80
## 5302 5302    35.42 3.11  9.67    29.32    41.53  8.79 10.44
## 5303 5303    25.65 2.89  8.35    19.99    31.32 11.27 10.11
## 5304 5304    29.59 3.04  9.24    23.64    35.54 10.26 11.66
## 5305 5305    29.05 3.46 11.97    22.26    35.84 11.93 12.66
## 5306 5306    23.47 3.43 11.76    16.73    30.20 14.64 14.39
## 5307 5307    31.00 2.82  7.95    25.47    36.53  9.10  8.58
## 5308 5308    38.11 3.40 11.56    31.44    44.78  8.93  9.97
## 5309 5309    32.07 2.89  8.35    26.42    37.73  9.00  8.36
## 5310 5310    55.26 3.30 10.89    48.78    61.73  5.98 11.41
## 5311 5311    25.35 3.58 12.82    18.34    32.36 14.11 15.34
## 5312 5312    34.23 3.52 12.39    27.33    41.13 10.28 11.92
## 5313 5313    23.06 2.40  5.76    18.36    27.75 10.39  7.91
## 5314 5314    21.49 2.65  7.02    16.29    26.68 12.33  8.89
## 5315 5315    23.30 3.14  9.86    17.14    29.47 13.49 11.69
## 5316 5316    26.79 3.51 12.32    19.90    33.68 13.12 15.41
## 5317 5317    31.88 3.80 14.44    24.43    39.33 11.93 19.29
## 5318 5318    27.09 2.91  8.47    21.39    32.80 10.74 10.16
## 5319 5319    25.84 2.34  5.48    21.25    30.42  9.05  6.07
## 5320 5320    42.89 3.32 11.02    36.38    49.39  7.74  9.78
## 5321 5321    20.94 3.50 12.25    14.08    27.80 16.71 14.55
## 5371 5371    30.92 3.47 12.04    24.13    37.71 11.21 13.92
## 6101 6101    27.85 2.50  6.25    22.95    32.75  8.98  8.22
## 6102 6102    25.39 3.41 11.63    18.71    32.08 13.43 12.91
## 6103 6103    29.62 3.50 12.25    22.76    36.48 11.82 13.90
## 6104 6104    27.72 2.71  7.34    22.41    33.04  9.78  8.68
## 6105 6105    24.90 3.03  9.18    18.96    30.84 12.17 11.85
## 6106 6106    23.10 3.94 15.52    15.37    30.82 17.06 20.54
## 6107 6107    17.99 1.78  3.17    14.50    21.47  9.88  5.09
## 6108 6108    17.78 3.52 12.39    10.89    24.68 19.78 16.52
## 6109 6109    17.80 4.34 18.84     9.29    26.31 24.39 24.24
## 6110 6110    27.28 2.40  5.76    22.57    31.99  8.81  5.88
## 6111 6111    34.07 3.64 13.25    26.93    41.21 10.70 12.23
## 6112 6112    20.39 2.15  4.62    16.18    24.60 10.53  7.60
## 6171 6171    25.69 1.75  3.06    22.27    29.12  6.80  4.39
## 6172 6172    28.52 2.48  6.15    23.66    33.38  8.69  7.15
## 6201 6201    15.37 2.05  4.20    11.34    19.39 13.36  6.51
## 6202 6202    20.30 2.46  6.05    15.48    25.13 12.11  8.04
## 6203 6203    28.13 2.55  6.50    23.13    33.14  9.07  7.06
## 6204 6204    19.09 2.82  7.95    13.56    24.62 14.77  9.42
## 6205 6205    19.22 2.72  7.40    13.89    24.54 14.13  8.42
## 6206 6206    15.34 2.25  5.06    10.93    19.75 14.68  7.18
## 6207 6207    26.61 3.24 10.50    20.26    32.96 12.18  9.85
## 6208 6208    17.43 3.31 10.96    10.94    23.93 19.01 13.87
## 6209 6209    19.38 3.02  9.12    13.46    25.31 15.59 10.76
## 6210 6210    19.50 3.30 10.89    13.03    25.98 16.94 12.49
## 6211 6211    22.31 3.69 13.62    15.08    29.55 16.55 15.42
## 6212 6212    26.85 2.36  5.57    22.23    31.48  8.79  5.11
## 6213 6213     7.40 1.38  1.90     4.70    10.10 18.63  5.32
## 6271 6271    31.00 3.70 13.69    23.74    38.25 11.94 13.96
## 6301 6301    27.72 3.09  9.55    21.66    33.77 11.14 10.14
## 6302 6302    27.38 2.75  7.56    21.98    32.78 10.06  8.49
## 6303 6303    35.83 2.68  7.18    30.58    41.07  7.47  7.74
## 6304 6304    33.33 3.41 11.63    26.65    40.01 10.22 10.88
## 6305 6305    27.24 3.19 10.18    20.98    33.49 11.72  9.14
## 6306 6306    46.43 2.96  8.76    40.62    52.23  6.38  6.80
## 6307 6307    29.98 2.67  7.13    24.75    35.21  8.90  6.97
## 6308 6308    34.09 3.75 14.06    26.74    41.44 11.00 11.87
## 6309 6309    26.87 2.89  8.35    21.21    32.53 10.74  8.67
## 6310 6310    22.22 2.53  6.40    17.26    27.18 11.40  8.68
## 6311 6311    40.36 3.75 14.06    33.01    47.72  9.29 10.81
## 6371 6371    19.26 2.18  4.75    14.98    23.53 11.32  7.30
## 6372 6372    26.13 3.05  9.30    20.14    32.11 11.68  9.16
## 6401 6401    14.49 2.56  6.55     9.48    19.50 17.65 10.16
## 6402 6402    17.68 2.60  6.76    12.58    22.77 14.70  9.29
## 6403 6403    18.61 2.68  7.18    13.36    23.86 14.40 12.49
## 6404 6404    16.60 2.57  6.60    11.56    21.63 15.48 10.36
## 6405 6405    24.44 3.32 11.02    17.93    30.95 13.59 12.65
## 6409 6409    23.67 2.34  5.48    19.08    28.25  9.89  6.06
## 6411 6411    17.48 3.89 15.13     9.85    25.12 22.28 20.49
## 6471 6471    22.16 1.78  3.17    18.67    25.65  8.04  4.55
## 6472 6472    22.61 2.06  4.24    18.57    26.65  9.12  7.20
## 6474 6474    16.14 1.86  3.46    12.49    19.78 11.53  5.26
## 6501 6501    20.08 2.28  5.20    15.62    24.54 11.34  6.98
## 6502 6502    20.78 2.58  6.66    15.72    25.84 12.43  8.34
## 6503 6503    27.24 3.36 11.29    20.65    33.83 12.34 11.29
## 6504 6504    16.52 2.04  4.16    12.51    20.52 12.37  6.40
## 6571 6571    21.30 2.37  5.62    16.65    25.94 11.13  7.90
## 7101 7101    30.93 4.62 21.34    21.88    39.98 14.93 22.39
## 7102 7102    29.44 3.55 12.60    22.48    36.41 12.06 12.80
## 7103 7103    25.84 2.75  7.56    20.45    31.23 10.65  7.01
## 7104 7104    13.52 2.15  4.62     9.31    17.74 15.90  7.63
## 7105 7105    36.74 2.59  6.71    31.67    41.81  7.04  4.98
## 7106 7106    23.05 3.73 13.91    15.74    30.36 16.17 14.08
## 7107 7107    21.10 2.29  5.24    16.60    25.60 10.87  5.98
## 7108 7108    12.17 2.08  4.33     8.10    16.24 17.08  7.51
## 7109 7109    26.89 3.25 10.56    20.52    33.26 12.09 10.08
## 7110 7110    28.28 2.59  6.71    23.20    33.35  9.16  6.80
## 7111 7111    22.46 2.68  7.18    17.21    27.70 11.91  8.16
## 7171 7171    12.36 1.31  1.72     9.79    14.93 10.60  3.51
## 7172 7172    19.52 1.87  3.50    15.86    23.19  9.58  4.58
## 7173 7173    22.95 2.18  4.75    18.67    27.23  9.51  4.81
## 7174 7174    20.13 1.87  3.50    16.46    23.80  9.31  4.37
## 7201 7201    13.77 1.82  3.31    10.21    17.32 13.19  4.97
## 7202 7202    23.36 2.82  7.95    17.83    28.88 12.07 10.23
## 7203 7203    10.46 2.46  6.05     5.64    15.28 23.51 12.39
## 7204 7204    24.25 2.59  6.71    19.18    29.32 10.67  7.25
## 7205 7205    25.12 4.03 16.24    17.21    33.02 16.05 19.48
## 7206 7206    24.87 2.31  5.34    20.35    29.40  9.28  6.07
## 7207 7207    31.44 3.07  9.42    25.42    37.46  9.77 10.20
## 7208 7208    22.41 2.59  6.71    17.34    27.48 11.54 10.24
## 7209 7209    18.00 2.04  4.16    14.00    22.01 11.35  5.84
## 7210 7210    28.85 2.89  8.35    23.18    34.52 10.03  9.10
## 7211 7211     8.94 1.73  2.99     5.56    12.33 19.32  6.95
## 7212 7212    12.98 1.79  3.20     9.48    16.48 13.76  4.98
## 7271 7271    22.83 2.97  8.82    17.01    28.65 13.02 10.77
## 7301 7301    19.06 1.77  3.13    15.59    22.53  9.29  4.08
## 7302 7302    23.14 2.59  6.71    18.06    28.21 11.19  8.22
## 7303 7303    54.16 4.35 18.92    45.63    62.70  8.04 15.21
## 7304 7304    21.19 2.33  5.43    16.62    25.76 11.00  7.77
## 7305 7305    29.17 2.62  6.86    24.04    34.31  8.98  7.76
## 7306 7306    29.90 3.32 11.02    23.40    36.40 11.09 14.10
## 7307 7307    24.87 2.58  6.66    19.81    29.93 10.38  7.78
## 7308 7308    23.19 2.38  5.66    18.54    27.85 10.24  7.99
## 7309 7309    23.57 2.10  4.41    19.46    27.69  8.91  6.10
## 7310 7310    24.85 3.73 13.91    17.54    32.16 15.00 15.32
## 7311 7311    13.83 1.41  1.99    11.06    16.59 10.21  4.90
## 7312 7312    17.98 2.48  6.15    13.13    22.84 13.77  8.55
## 7313 7313    25.90 2.98  8.88    20.06    31.74 11.50 10.66
## 7314 7314    28.64 3.64 13.25    21.51    35.77 12.71 15.53
## 7315 7315    35.74 3.82 14.59    28.26    43.23 10.69 14.80
## 7316 7316    17.48 1.53  2.34    14.48    20.49  8.76  3.43
## 7317 7317    23.69 2.28  5.20    19.21    28.16  9.65  7.65
## 7318 7318    22.49 2.31  5.34    17.96    27.02 10.28  6.32
## 7322 7322    25.25 2.70  7.29    19.97    30.54 10.68  8.43
## 7325 7325    22.15 2.15  4.62    17.94    26.37  9.70  6.24
## 7326 7326    18.81 2.18  4.75    14.53    23.09 11.62  7.14
## 7371 7371    17.82 1.85  3.42    14.19    21.44 10.38  7.80
## 7372 7372    27.12 2.85  8.12    21.54    32.70 10.50  8.75
## 7373 7373    23.99 2.39  5.71    19.31    28.67  9.95  6.92
## 7401 7401    23.32 2.85  8.12    17.73    28.92 12.24 10.51
## 7402 7402    25.07 2.88  8.29    19.43    30.71 11.48  8.73
## 7403 7403    18.56 2.53  6.40    13.59    23.52 13.65  8.71
## 7404 7404    14.49 1.83  3.35    10.91    18.08 12.61  5.95
## 7405 7405    18.33 2.11  4.45    14.20    22.46 11.49  6.13
## 7406 7406    22.05 3.16  9.99    15.85    28.25 14.34 11.42
## 7407 7407    13.35 1.60  2.56    10.21    16.49 12.00  4.42
## 7408 7408    28.67 3.19 10.18    22.41    34.93 11.14 10.15
## 7409 7409    39.08 3.31 10.96    32.60    45.55  8.46  9.71
## 7410 7410    16.28 2.05  4.20    12.27    20.30 12.59  7.08
## 7411 7411    21.88 3.17 10.05    15.67    28.09 14.48 11.34
## 7412 7412    18.60 2.02  4.08    14.64    22.57 10.88  5.07
## 7413 7413    35.08 3.45 11.90    28.32    41.84  9.83  9.70
## 7414 7414    17.61 2.14  4.58    13.42    21.80 12.13  6.52
## 7415 7415    20.44 2.33  5.43    15.87    25.00 11.40  7.83
## 7471 7471    20.70 2.74  7.51    15.33    26.06 13.22 10.46
## 7472 7472    18.40 2.05  4.20    14.38    22.41 11.12  6.12
## 7501 7501    37.67 2.82  7.95    32.14    43.21  7.49  7.05
## 7502 7502    28.18 2.79  7.78    22.71    33.66  9.91  9.23
## 7503 7503    38.42 3.21 10.30    32.14    44.70  8.34  9.15
## 7504 7504    30.66 2.61  6.81    25.53    35.78  8.53  7.35
## 7505 7505    36.58 3.89 15.13    28.95    44.21 10.64 12.42
## 7571 7571    37.52 2.56  6.55    32.50    42.55  6.83  5.72
## 7601 7601    21.08 2.53  6.40    16.12    26.04 12.01  9.04
## 7602 7602    22.39 2.90  8.41    16.72    28.07 12.93 13.73
## 7603 7603    22.93 2.62  6.86    17.79    28.07 11.44  7.67
## 7604 7604    23.37 2.70  7.29    18.08    28.66 11.56 10.98
## 7605 7605    30.14 3.12  9.73    24.03    36.26 10.35 10.00
## 7606 7606    20.37 2.52  6.35    15.43    25.31 12.37  7.93
## 8101 8101     9.04 1.42  2.02     6.26    11.82 15.67  5.28
## 8102 8102    20.60 2.86  8.18    14.99    26.21 13.90 12.08
## 8103 8103    14.71 2.08  4.33    10.63    18.80 14.15  8.98
## 8104 8104    14.15 2.45  6.00     9.35    18.95 17.30 10.85
## 8105 8105     5.31 0.66  0.44     4.02     6.60 12.42  2.13
## 8106 8106    14.43 1.63  2.66    11.23    17.62 11.30  5.38
## 8107 8107    13.04 1.48  2.19    10.15    15.93 11.31  4.76
## 8108 8108    13.42 1.87  3.50     9.75    17.09 13.95  7.17
## 8109 8109    15.14 1.74  3.03    11.72    18.56 11.53  6.22
## 8171 8171    15.87 2.06  4.24    11.83    19.91 12.99  7.66
## 8172 8172    16.55 1.95  3.80    12.73    20.37 11.78  6.55
## 8201 8201    14.27 2.25  5.06     9.85    18.69 15.80  8.45
## 8202 8202    17.78 3.57 12.74    10.78    24.77 20.08 17.65
## 8203 8203    13.11 2.20  4.84     8.79    17.42 16.80  9.43
## 8204 8204    12.29 1.35  1.82     9.65    14.94 10.98  3.91
## 8205 8205    11.97 1.33  1.77     9.36    14.59 11.13  4.02
## 8206 8206    13.48 2.19  4.80     9.18    17.77 16.26  8.50
## 8207 8207    14.30 1.87  3.50    10.63    17.97 13.09  6.78
## 8208 8208    11.38 1.57  2.46     8.31    14.46 13.77  5.00
## 8271 8271    14.76 2.36  5.57    10.13    19.39 16.01 10.24
## 8272 8272    10.21 1.35  1.82     7.57    12.85 13.20  4.46
## 9101 9101     9.56 1.36  1.85     6.89    12.24 14.26  4.60
## 9102 9102    15.72 1.66  2.76    12.46    18.98 10.58  5.04
## 9103 9103    11.98 2.47  6.10     7.14    16.82 20.61 11.47
## 9104 9104    14.63 1.96  3.84    10.79    18.47 13.41  7.39
## 9105 9105    16.32 2.68  7.18    11.08    21.57 16.40 11.39
## 9111 9111    14.52 2.22  4.93    10.18    18.86 15.26  8.03
## 9112 9112    12.39 2.04  4.16     8.39    16.38 16.44  7.56
## 9201 9201    16.75 1.67  2.79    13.48    20.03  9.98  4.96
## 9202 9202    17.67 1.88  3.53    13.98    21.36 10.66  5.37
## 9203 9203    25.15 2.86  8.18    19.55    30.75 11.36 13.45
## 9204 9204    18.48 3.30 10.89    12.01    24.94 17.85 11.11
## 9205 9205     8.62 1.47  2.16     5.75    11.50 17.02  4.02
## 9271 9271    16.42 2.08  4.33    12.35    20.50 12.65  7.94
## 9403 9403    13.48 1.89  3.57     9.77    17.18 14.02  6.57
## 9408 9408    12.86 2.14  4.58     8.66    17.07 16.67  8.99
## 9409 9409    18.66 2.52  6.35    13.72    23.59 13.51  9.49
## 9419 9419    10.31 1.47  2.16     7.42    13.20 14.30  4.32
## 9420 9420    16.30 1.99  3.96    12.41    20.20 12.18  5.33
## 9426 9426    14.14 1.65  2.72    10.91    17.37 11.66  4.85
## 9427 9427     8.67 1.22  1.49     6.27    11.06 14.10  4.50
## 9428 9428     4.30 1.05  1.10     2.25     6.35 24.34  2.64
## 9471 9471    14.74 1.94  3.76    10.93    18.54 13.18  6.41
## 9501 9501    23.00 2.13  4.54    18.82    27.17  9.26  5.56
## 9502 9502    27.60 3.65 13.32    20.45    34.74 13.21 13.96
## 9503 9503    27.80 3.62 13.10    20.71    34.90 13.02 15.71
## 9504 9504     6.46 1.03  1.06     4.44     8.48 15.93  4.32
## 9601 9601     8.17 1.60  2.56     5.04    11.30 19.56  6.64
## 9602 9602    24.11 6.02 36.24    12.31    35.91 24.97 36.43
## 9603 9603     2.28 0.89  0.79     0.54     4.03 38.93  6.10
## 9604 9604    17.23 2.56  6.55    12.21    22.26 14.86  9.47
## 9605 9605     5.48 2.15  4.62     1.28     9.69 39.14 14.43
## 9606 9606     3.71 1.05  1.10     1.64     5.77 28.40  5.76
## 9607 9607     1.39 0.52  0.27     0.37     2.41 37.43  3.46
## 9608 9608    20.22 5.67 32.15     9.10    31.34 28.06 29.36
## 9701 9701     7.74 1.12  1.25     5.54     9.95 14.50  2.97
## 9702 9702     8.78 1.90  3.61     5.06    12.49 21.60  7.83
## 9703 9703     1.31 0.49  0.24     0.35     2.27 37.40  2.57
## 9704 9704     9.37 1.21  1.46     6.99    11.75 12.95  3.42
## 9705 9705    14.14 0.96  0.92    12.25    16.02  6.80  1.47
## 9706 9706    10.64 2.77  7.67     5.22    16.07 26.00 10.24
## 9707 9707    14.55 2.16  4.67    10.33    18.78 14.81  8.30
## 9708 9708    10.25 1.67  2.79     6.98    13.51 16.26  6.16

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 - R",
    x = "Kab/kota",
    y = "Nilai"
  ) +
  scale_color_manual(name = " ", values = c("RSE" = "magenta", "Theta" = "purple")) +
  theme_minimal()

Eksport to CSV

write.csv(resultP,"RSE_DIRECT_KELUH_PROV_2024.csv")
write.csv(resultK,"RSE_DIRECT_KELUH_KAKO_2024.csv")
rm(list=ls())  # Clear environment

4 Penghitungan dengan SPSS

4.1 Set Design Sampling

Desain sampling untuk estimasi indikator yang dibangun menggunakan SPSS memiliki keterbatasan, terutama pada jenis desain sampling tertentu. Hal ini disebabkan oleh fakta bahwa SPSS adalah sebuah paket perangkat lunak statistik yang terbatas pada fitur-fitur yang tersedia di dalam menu-menunya, sehingga pengembangan lebih lanjut menjadi sangat terbatas.

Perintah CSPLAN dalam SPSS digunakan untuk mendefinisikan dan menerapkan rencana desain sampling yang kompleks dalam analisis data. Dengan melibatkan FPC pada penghitungan, penulisan sintaks CSPLAN SPSS adalah sebagai berikut:

* Analysis Preparation Wizard.
CSPLAN ANALYSIS  
/PLAN FILE='D:\RSE SPSS\RSE.csaplan'
/PLANVARS ANALYSISWEIGHT=FWT         
/PRINT PLAN  
/DESIGN STRATA=STRATA CLUSTER=PSU
/ESTIMATOR TYPE=EQUAL_WOR  
/POPSIZE VARIABLE=fpc1  
/DESIGN CLUSTER=SSU
/ESTIMATOR TYPE=EQUAL_WOR  
/POPSIZE VARIABLE=fpc2.

Dari CSPLAN di atas tampak:

  • variabel STRATA sebagai variabel stratifikasi yang berarti data dibagi menjadi strata yang berbeda. Variabel STRATA sudah disediakan dalam dataset.

  • variabel PSU ( Primary Sampling Unit) sebagai variabel klaster untuk level pertama dari desain sampling. Variabel PSU sudah disediakan dalam dataset.

  • TYPE=EQUAL_WOR menunjukkan estimator yang digunakan adalah tanpa pengembalian ( Without Replacement), dan semua unit sampel dianggap memiliki peluang yang sama untuk dipilih.

  • VARIABLE=fpc1 menetapkan variabel fpc1 sebagai ukuran populasi untuk strata yang berkaitan dengan level pertama dari desain sampling.

  • variabel SSU ( Secondary Sampling Unit) sebagai variabel klaster untuk level kedua dari desain sampling. Variabel SSU sudah disediakan dalam dataset.

  • VARIABLE=fpc2 menetapkan variabel fpc2 sebagai ukuran populasi untuk strata yang terkait dengan level kedua dari desain sampling.

4.2 Level Estimate

Estimation by province

Hasil estimasi dengan SPSS menurut provinsi ditunjukkan:

dataP <- read_excel("HEALTH_KELUH PROV 2024.xls")
dataP[, 2:7] <- round(dataP[, 2:7], 2)
## # A tibble: 38 × 8
##     PROV EST   SE    `CI LOWER` `CI UPPER` CV    DEFF   POPSIZE
##    <dbl> <chr> <chr> <chr>      <chr>      <chr> <chr>    <dbl>
##  1    11 25.58 0.65  24.31      26.85      0.03  5.46   5652120
##  2    12 23.99 0.59  22.85      25.14      0.02  12.59 15368583
##  3    13 30.10 0.72  28.69      31.51      0.02  6.18   5786360
##  4    14 23.60 0.81  22.02      25.18      0.03  11.60  7381474
##  5    15 23.73 0.84  22.08      25.37      0.04  6.36   3739406
##  6    16 29.89 0.90  28.11      31.66      0.03  15.25  8975590
##  7    17 28.27 1.12  26.08      30.46      0.04  5.58   2075023
##  8    18 29.54 0.85  27.87      31.22      0.03  13.45  8806627
##  9    19 28.80 1.04  26.76      30.84      0.04  3.53   1537917
## 10    21 26.26 2.23  21.89      30.62      0.08  28.75  2575911
## 11    31 24.07 1.05  22.02      26.13      0.04  28.25 10820116
## 12    32 28.88 0.55  27.81      29.96      0.02  33.01 51568693
## 13    33 31.86 0.44  31.00      32.73      0.01  13.75 35379588
## 14    34 28.82 0.96  26.93      30.70      0.03  8.07   4113253
## 15    35 28.75 0.44  27.90      29.61      0.02  16.43 40660605
## 16    36 25.77 1.01  23.79      27.75      0.04  31.59 13557010
## 17    51 19.32 0.75  17.86      20.79      0.04  7.17   4609533
## 18    52 38.88 1.05  36.82      40.94      0.03  11.11  5493860
## 19    53 29.56 0.73  28.13      31.00      0.02  6.51   5789683
## 20    61 24.37 0.80  22.80      25.95      0.03  8.14   5318285
## 21    62 21.37 0.86  19.69      23.06      0.04  5.40   2818749
## 22    63 29.00 0.84  27.36      30.65      0.03  6.65   4457578
## 23    64 20.01 0.92  18.21      21.81      0.05  8.76   3827332
## 24    65 19.93 1.19  17.60      22.25      0.06  2.91    756927
## 25    71 23.30 0.88  21.58      25.01      0.04  4.81   2577009
## 26    72 22.41 0.88  20.67      24.14      0.04  6.33   3226963
## 27    73 23.43 0.61  22.23      24.63      0.03  8.35   9134428
## 28    74 20.87 0.73  19.43      22.30      0.04  4.03   2852319
## 29    75 33.80 1.27  31.31      36.29      0.04  3.85   1221372
## 30    76 23.33 1.26  20.85      25.80      0.05  5.62   1446583
## 31    81 14.38 0.78  12.85      15.92      0.05  4.03   1854801
## 32    82 13.22 0.68  11.88      14.55      0.05  2.32   1314102
## 33    91 14.18 1.08  12.06      16.29      0.08  2.12    508533
## 34    92 17.43 1.22  15.04      19.81      0.07  2.53    564043
## 35    94 14.43 0.92  12.63      16.23      0.06  2.63    885767
## 36    95 21.24 1.35  18.60      23.88      0.06  2.50    528706
## 37    96 10.72 1.20  8.38       13.07      0.11  6.75   1036419
## 38    97 8.88  0.64  7.63       10.13      0.07  2.42   1108541

Visualisasi

ggplot(data = dataP, aes(x = PROV)) +
  geom_line(aes(y = CV*100, color = "RSE", group = 1), size = 1) +      
  geom_line(aes(y = EST, color = "Theta", group = 1), size = 1) +       
  geom_hline(yintercept = 25, color = "red", linetype = "dashed") +     
  labs(
    title = "Provinsi vs Nilai RSE dan Theta - SPSS",
    x = "Provinsi",
    y = "Nilai"
  ) +
  scale_color_manual(name = " ", values = c("RSE" = "orange", "Theta" = "blue")) +
  theme_minimal()

Estimation by district

Hasil estimasi menurut kab/kota ditunjukkan sebagai berikut:

dataK <- read_excel("HEALTH_KELUH KAKO 2024.xls")
dataK[, 2:7] <- round(dataK[, 2:7], 2)
## # A tibble: 514 × 8
##      KAKO EST   SE    `CI LOWER` `CI UPPER` CV    DEFF  POPSIZE
##     <dbl> <chr> <chr> <chr>      <chr>      <chr> <chr>   <dbl>
##   1  1101 18.72 2.63  13.55      23.88      0.14  1.97    99426
##   2  1102 15.60 1.78  12.10      19.09      0.11  1.37   130276
##   3  1103 13.04 1.82  9.48       16.60      0.14  3.22   253955
##   4  1104 24.44 3.17  18.23      30.65      0.13  5.39   227500
##   5  1105 29.18 3.48  22.36      35.99      0.12  11.70  458974
##   6  1106 31.16 3.51  24.27      38.04      0.11  5.59   223163
##   7  1107 27.15 2.53  22.20      32.11      0.09  3.10   220313
##   8  1108 24.67 2.12  20.53      28.82      0.09  4.68   446424
##   9  1109 31.53 2.36  26.91      36.15      0.07  5.24   467276
##  10  1110 33.73 2.29  29.25      38.22      0.07  5.05   495113
##  11  1111 20.27 2.63  15.12      25.42      0.13  12.06  648511
##  12  1112 30.50 2.34  25.91      35.09      0.08  1.80   159547
##  13  1113 19.62 2.63  14.46      24.77      0.13  1.87    97994
##  14  1114 22.14 2.24  17.75      26.54      0.10  3.95   310171
##  15  1115 23.91 3.36  17.32      30.51      0.14  4.73   174713
##  16  1116 39.05 2.62  33.93      44.18      0.07  1.25    99921
##  17  1117 23.63 2.83  18.09      29.17      0.12  3.01   156023
##  18  1118 25.78 2.61  20.67      30.89      0.10  2.62   169321
##  19  1171 26.86 2.57  21.83      31.89      0.10  4.19   287217
##  20  1172 23.94 2.60  18.84      29.03      0.11  0.57    35478
##  21  1173 29.03 3.54  22.08      35.97      0.12  4.98   187458
##  22  1174 21.38 2.01  17.44      25.32      0.09  2.28   218190
##  23  1175 16.87 1.76  13.42      20.32      0.10  0.82    85156
##  24  1201 31.11 3.08  25.09      37.14      0.10  2.81   145999
##  25  1202 16.91 1.89  13.21      20.61      0.11  5.15   466484
##  26  1203 23.48 3.09  17.43      29.54      0.13  6.65   287561
##  27  1204 21.56 1.65  18.34      24.79      0.08  2.92   417125
##  28  1205 13.51 2.31  8.99       18.03      0.17  6.17   310674
##  29  1206 16.36 1.94  12.55      20.17      0.12  2.24   186348
##  30  1207 14.09 2.10  9.96       18.21      0.15  8.63   541503
##  31  1208 27.87 2.40  23.18      32.57      0.09  9.41   755818
##  32  1209 24.83 2.30  20.33      29.33      0.09  10.86  882045
##  33  1210 64.94 3.61  57.87      72.00      0.06  7.14   287131
##  34  1211 19.54 2.43  14.78      24.30      0.12  7.42   453850
##  35  1212 20.69 2.11  16.57      24.82      0.10  28.86 2450944
##  36  1213 29.61 2.57  24.57      34.66      0.09  14.76 1067648
##  37  1214 12.53 2.12  8.38       16.68      0.17  5.96   333734
##  38  1215 16.07 1.97  12.21      19.93      0.12  2.49   198723
##  39  1216 32.78 2.66  27.56      37.99      0.08  0.76    54227
##  40  1217 28.95 2.61  23.83      34.07      0.09  1.85   127694
##  41  1218 25.02 2.88  19.37      30.67      0.12  11.90  616922
##  42  1219 30.86 2.72  25.53      36.19      0.09  6.57   435459
##  43  1220 17.18 2.34  12.59      21.78      0.14  5.10   303350
##  44  1221 32.42 4.10  24.38      40.46      0.13  10.62  317322
##  45  1222 19.07 2.24  14.68      23.47      0.12  5.36   377497
##  46  1223 24.89 3.67  17.70      32.08      0.15  11.86  378629
##  47  1224 29.87 3.94  22.15      37.58      0.13  4.58   141992
##  48  1225 33.66 3.38  27.04      40.28      0.10  1.84    82826
##  49  1271 27.10 1.62  23.92      30.29      0.06  0.51    87715
##  50  1272 25.41 2.30  20.90      29.92      0.09  2.26   185832
##  51  1273 24.33 2.80  18.84      29.82      0.12  4.90   263559
##  52  1274 25.79 3.11  19.70      31.88      0.12  3.82   174052
##  53  1275 22.70 1.89  18.99      26.40      0.08  20.90 2354446
##  54  1276 25.67 2.58  20.60      30.73      0.10  4.45   291942
##  55  1277 32.08 3.00  26.19      37.97      0.09  4.31   238962
##  56  1278 35.84 3.53  28.92      42.75      0.10  3.55   150570
##  57  1301 30.23 4.69  21.03      39.43      0.16  4.65   102199
##  58  1302 26.69 2.22  22.34      31.03      0.08  5.26   481360
##  59  1303 24.77 2.33  20.21      29.33      0.09  4.88   385506
##  60  1304 29.10 2.96  23.29      34.91      0.10  4.81   259518
##  61  1305 41.40 2.73  36.05      46.76      0.07  4.70   350522
##  62  1306 37.77 2.78  32.32      43.22      0.07  6.11   426846
##  63  1307 28.48 2.24  24.08      32.87      0.08  5.48   509672
##  64  1308 30.09 3.42  23.39      36.79      0.11  9.70   400832
##  65  1309 28.14 2.98  22.31      33.98      0.11  5.66   296719
##  66  1310 31.77 3.40  25.11      38.42      0.11  4.35   187880
##  67  1311 32.29 2.60  27.19      37.39      0.08  3.86   286551
##  68  1312 24.84 3.10  18.75      30.92      0.12  11.11  494395
##  69  1371 30.16 2.08  26.08      34.23      0.07  9.11  1019532
##  70  1372 29.32 2.34  24.73      33.92      0.08  0.90    78202
##  71  1373 26.74 1.94  22.94      30.55      0.07  0.55    65393
##  72  1374 31.01 2.86  25.39      36.63      0.09  0.96    57641
##  73  1375 23.99 2.18  19.73      28.26      0.09  1.61   142636
##  74  1376 29.77 2.30  25.27      34.28      0.08  1.62   146993
##  75  1377 39.38 2.77  33.95      44.81      0.07  1.32    93963
##  76  1401 12.94 1.53  9.95       15.94      0.12  2.92   324023
##  77  1402 27.08 2.80  21.60      32.57      0.10  7.90   457023
##  78  1403 34.72 2.95  28.93      40.50      0.08  12.13  725358
##  79  1404 12.31 1.79  8.81       15.81      0.15  7.46   579154
##  80  1405 24.82 3.40  18.16      31.48      0.14  14.16  525525
##  81  1406 19.03 2.06  14.99      23.08      0.11  11.09  920605
##  82  1407 18.47 2.40  13.76      23.18      0.13  13.18  788577
##  83  1408 13.10 1.78  9.62       16.58      0.14  6.95   574962
##  84  1409 35.52 3.17  29.31      41.72      0.09  14.55  763165
##  85  1410 22.60 2.50  17.71      27.50      0.11  2.98   191497
##  86  1471 28.87 2.44  24.09      33.65      0.08  15.29 1212318
##  87  1473 19.48 2.02  15.52      23.43      0.10  3.61   319267
##  88  1501 37.13 3.18  30.89      43.37      0.09  4.59   242910
##  89  1502 23.49 2.78  18.04      28.95      0.12  7.47   398645
##  90  1503 23.78 2.36  19.16      28.39      0.10  4.24   317622
##  91  1504 20.72 3.14  14.57      26.87      0.15  7.18   274958
##  92  1505 19.57 2.71  14.26      24.88      0.14  9.75   479461
##  93  1506 21.76 2.77  16.33      27.20      0.13  4.24   215592
##  94  1507 13.28 1.79  9.77       16.80      0.14  4.22   347103
##  95  1508 27.04 2.98  21.21      32.88      0.11  7.19   367132
##  96  1509 30.71 2.88  25.08      36.35      0.09  6.70   396075
##  97  1571 23.58 2.10  19.48      27.69      0.09  6.48   611120
##  98  1572 21.43 2.04  17.44      25.43      0.10  0.95    88788
##  99  1601 39.83 3.09  33.76      45.89      0.08  6.76   389038
## 100  1602 28.98 2.95  23.19      34.77      0.10  16.26  880355
## 101  1603 30.29 2.98  24.45      36.12      0.10  12.39  678350
## 102  1604 34.87 3.05  28.89      40.84      0.09  7.59   425588
## 103  1605 25.54 3.01  19.63      31.44      0.12  8.89   427444
## 104  1606 22.07 3.09  16.01      28.13      0.14  16.71  690433
## 105  1607 21.43 2.92  15.70      27.15      0.14  20.01  906160
## 106  1608 23.55 3.56  16.56      30.53      0.15  11.77  383322
## 107  1609 45.41 3.66  38.24      52.58      0.08  16.57  705125
## 108  1610 42.38 2.58  37.33      47.43      0.06  5.35   451114
## 109  1611 34.83 3.70  27.58      42.09      0.11  6.93   263568
## 110  1612 37.42 4.23  29.13      45.71      0.11  6.71   201487
## 111  1613 24.14 3.10  18.07      30.21      0.13  4.63   202982
## 112  1671 27.83 2.86  22.23      33.43      0.10  31.43 1776384
## 113  1672 23.68 2.67  18.45      28.91      0.11  3.45   201469
## 114  1673 27.44 3.48  20.63      34.26      0.13  3.82   144488
## 115  1674 25.53 2.66  20.32      30.74      0.10  4.02   248283
## 116  1701 33.63 3.61  26.55      40.72      0.11  4.74   186023
## 117  1702 24.69 2.77  19.25      30.12      0.11  5.07   281047
## 118  1703 28.20 3.08  22.16      34.24      0.11  6.44   315079
## 119  1704 26.13 2.55  21.13      31.13      0.10  1.87   127020
## 120  1705 35.59 4.46  26.85      44.32      0.13  7.68   203628
## 121  1706 21.33 2.66  16.12      26.54      0.12  3.56   194086
## 122  1707 18.27 1.49  15.35      21.19      0.08  0.77   119756
## 123  1708 34.23 3.60  27.17      41.28      0.11  3.64   145273
## 124  1709 30.97 3.28  24.55      37.40      0.11  2.59   118172
## 125  1771 28.72 3.33  22.19      35.25      0.12  9.09   384939
## 126  1801 29.04 4.84  19.55      38.52      0.17  15.49  312857
## 127  1802 30.41 2.71  25.09      35.73      0.09  9.47   624609
## 128  1803 23.74 2.64  18.56      28.92      0.11  17.70 1053039
## 129  1804 26.09 2.61  20.96      31.21      0.10  16.70 1081438
## 130  1805 34.33 2.46  29.50      39.15      0.07  15.51 1324972
## 131  1806 21.68 2.17  17.43      25.93      0.10  7.55   626000
## 132  1807 45.76 3.90  38.11      53.40      0.09  12.51  468686
## 133  1808 32.22 3.31  25.73      38.72      0.10  10.40  475235
## 134  1809 35.07 4.68  25.89      44.24      0.13  19.44  463865
## 135  1810 34.88 3.26  28.49      41.27      0.09  8.44   413963
## 136  1811 31.44 2.46  26.62      36.27      0.08  2.50   204464
## 137  1812 33.04 4.05  25.10      40.97      0.12  9.10   281777
## 138  1813 21.00 2.67  15.76      26.25      0.13  3.00   159922
## 139  1871 24.98 2.58  19.92      30.03      0.10  17.59 1137850
## 140  1872 27.32 2.05  23.31      31.34      0.07  1.64   177950
## 141  1901 30.79 2.79  25.33      36.26      0.09  5.54   348555
## 142  1902 26.84 2.64  21.66      32.01      0.10  3.03   196086
## 143  1903 30.13 2.34  25.54      34.72      0.08  2.51   220451
## 144  1904 26.97 2.62  21.83      32.11      0.10  3.08   202582
## 145  1905 22.21 2.36  17.58      26.84      0.11  3.05   216798
## 146  1906 31.66 2.47  26.82      36.49      0.08  1.64   133701
## 147  1971 32.49 2.83  26.95      38.03      0.09  3.48   219744
## 148  2101 20.52 2.67  15.28      25.76      0.13  4.63   242189
## 149  2102 18.01 1.95  14.20      21.83      0.11  1.89   168677
## 150  2103 17.53 1.99  13.64      21.42      0.11  0.98    82798
## 151  2104 16.53 2.20  12.21      20.85      0.13  1.39    90358
## 152  2105 15.14 1.96  11.29      18.99      0.13  0.58    44661
## 153  2171 28.96 3.26  22.58      35.35      0.11  38.69 1722927
## 154  2172 27.20 2.67  21.97      32.43      0.10  3.51   224301
## 155  3101 35.06 2.74  29.69      40.42      0.08  0.38    26768
## 156  3171 19.79 1.72  16.41      23.16      0.09  18.98 2333842
## 157  3172 18.73 2.05  14.71      22.76      0.11  35.88 2971674
## 158  3173 26.85 3.39  20.21      33.49      0.13  23.27  914918
## 159  3174 28.58 2.51  23.67      33.50      0.09  36.29 2701771
## 160  3175 29.88 2.09  25.77      33.98      0.07  17.07 1871143
## 161  3201 28.81 2.39  24.14      33.49      0.08  76.62 6335412
## 162  3202 31.36 2.34  26.77      35.95      0.07  28.42 2560066
## 163  3203 37.80 2.23  33.42      42.17      0.06  21.80 2360231
## 164  3204 25.80 1.91  22.07      29.54      0.07  32.03 3873970
## 165  3205 37.88 2.40  33.17      42.59      0.06  28.63 2679580
## 166  3206 33.57 2.78  28.13      39.01      0.08  27.47 1825834
## 167  3207 43.22 3.05  37.25      49.20      0.07  20.27 1229132
## 168  3208 38.23 2.76  32.83      43.63      0.07  15.55 1109831
## 169  3209 39.75 2.80  34.26      45.24      0.07  31.93 2235596
## 170  3210 36.83 2.88  31.18      42.48      0.08  19.35 1243280
## 171  3211 28.72 2.05  24.71      32.74      0.07  10.70 1196056
## 172  3212 32.66 3.42  25.96      39.36      0.10  41.13 1778745
## 173  3213 19.58 2.07  15.53      23.64      0.11  19.03 1607575
## 174  3214 32.75 2.52  27.80      37.70      0.08  12.21  968614
## 175  3215 24.30 2.51  19.38      29.22      0.10  35.65 2390062
## 176  3216 20.06 2.13  15.89      24.24      0.11  52.56 4261112
## 177  3217 30.66 2.06  26.61      34.70      0.07  14.96 1712644
## 178  3218 47.56 2.98  41.71      53.41      0.06  6.38   410654
## 179  3271 27.44 2.58  22.39      32.49      0.09  16.44 1132548
## 180  3272 29.09 2.15  24.88      33.30      0.07  3.26   334447
## 181  3273 20.95 2.15  16.73      25.17      0.10  31.85 2609884
## 182  3274 29.67 3.22  23.35      35.98      0.11  7.00   323333
## 183  3275 23.46 1.87  19.79      27.13      0.08  27.35 3215031
## 184  3276 19.36 1.61  16.20      22.53      0.08  19.47 2677666
## 185  3277 24.29 1.98  20.40      28.18      0.08  5.74   615197
## 186  3278 30.61 2.64  25.44      35.79      0.09  9.89   691414
## 187  3279 25.14 2.48  20.29      29.99      0.10  2.71   190779
## 188  3301 28.04 2.13  23.87      32.21      0.08  17.03 1741630
## 189  3302 28.49 1.95  24.67      32.32      0.07  14.17 1738812
## 190  3303 38.29 1.98  34.41      42.16      0.05  6.93   963996
## 191  3304 28.93 2.37  24.28      33.58      0.08  11.16  935798
## 192  3305 55.41 3.13  49.28      61.53      0.06  20.62 1196334
## 193  3306 34.56 2.67  29.33      39.78      0.08  9.83   717488
## 194  3307 38.66 2.68  33.41      43.92      0.07  10.49  795330
## 195  3308 34.38 2.10  30.25      38.50      0.06  11.33 1327563
## 196  3309 30.16 2.49  25.27      35.04      0.08  12.81  996316
## 197  3310 34.60 2.23  30.23      38.97      0.06  11.27 1177782
## 198  3311 25.13 2.31  20.59      29.66      0.09  11.31  912591
## 199  3312 23.19 1.80  19.66      26.72      0.08  7.59   958129
## 200  3313 20.62 2.13  16.46      24.79      0.10  10.96  911144
## 201  3314 26.55 2.54  21.57      31.53      0.10  12.86  891604
## 202  3315 31.67 2.25  27.25      36.08      0.07  14.18 1390457
## 203  3316 36.68 3.11  30.59      42.77      0.08  15.72  868433
## 204  3317 32.65 2.80  27.16      38.15      0.09  10.15  652499
## 205  3318 28.06 3.03  22.12      34.00      0.11  25.24 1273559
## 206  3319 32.38 3.05  26.41      38.35      0.09  16.76  908215
## 207  3320 28.17 2.66  22.96      33.38      0.09  20.15 1327012
## 208  3321 30.49 2.64  25.32      35.66      0.09  17.18 1201479
## 209  3322 28.02 2.09  23.92      32.13      0.07  10.46 1104647
## 210  3323 31.23 2.60  26.14      36.33      0.08  10.86  793099
## 211  3324 31.26 2.38  26.59      35.93      0.08  11.39  992005
## 212  3325 28.67 2.92  22.94      34.40      0.10  14.34  788749
## 213  3326 33.02 2.57  27.98      38.07      0.08  11.92  913423
## 214  3327 40.67 2.34  36.07      45.26      0.06  12.91 1301423
## 215  3328 32.99 1.86  29.34      36.63      0.06  9.78  1438682
## 216  3329 39.75 2.41  35.02      44.47      0.06  19.21 1816307
## 217  3371 25.19 3.14  19.03      31.35      0.12  2.78   122027
## 218  3372 32.92 2.33  28.36      37.47      0.07  5.56   520966
## 219  3373 29.40 2.27  24.95      33.85      0.08  2.21   204174
## 220  3374 23.60 2.03  19.62      27.59      0.09  19.28 1931626
## 221  3375 39.86 2.69  34.58      45.14      0.07  4.16   315484
## 222  3376 25.03 2.23  20.67      29.40      0.09  2.89   250805
## 223  3401 30.10 2.29  25.61      34.59      0.08  4.98   457699
## 224  3402 31.76 1.94  27.94      35.57      0.06  8.28  1089070
## 225  3403 23.44 1.68  20.14      26.74      0.07  5.46   793053
## 226  3404 28.54 2.01  24.60      32.48      0.07  11.33 1313634
## 227  3471 30.65 2.31  26.12      35.17      0.08  5.02   459797
## 228  3501 36.73 2.95  30.95      42.52      0.08  9.13   558451
## 229  3502 39.36 2.57  34.33      44.39      0.07  10.58  879148
## 230  3503 28.47 2.10  24.36      32.59      0.07  6.63   702340
## 231  3504 30.81 3.09  24.76      36.86      0.10  20.58 1057038
## 232  3505 36.61 2.99  30.74      42.47      0.08  19.73 1174881
## 233  3506 30.57 2.81  25.06      36.08      0.09  25.97 1601435
## 234  3507 28.04 1.95  24.22      31.86      0.07  21.91 2673149
## 235  3508 24.62 2.31  20.09      29.14      0.09  13.17 1052577
## 236  3509 30.10 2.16  25.85      34.34      0.07  24.20 2494548
## 237  3510 32.95 2.36  28.32      37.58      0.07  17.91 1631182
## 238  3511 35.76 2.79  30.29      41.23      0.08  11.68  790478
## 239  3512 26.96 2.56  21.93      31.98      0.10  10.13  696579
## 240  3513 32.04 2.86  26.43      37.65      0.09  19.62 1198556
## 241  3514 20.38 2.04  16.38      24.37      0.10  18.74 1677757
## 242  3515 23.15 1.86  19.50      26.80      0.08  20.50 2414283
## 243  3516 31.59 2.85  26.00      37.18      0.09  19.01 1160398
## 244  3517 36.72 2.10  32.61      40.82      0.06  10.59 1286116
## 245  3518 34.19 2.12  30.04      38.33      0.06  9.24  1065640
## 246  3519 26.02 2.50  21.13      30.91      0.10  9.69   687887
## 247  3520 31.07 2.44  26.28      35.87      0.08  7.76   638373
## 248  3521 29.09 2.59  24.01      34.18      0.09  11.98  843265
## 249  3522 21.28 1.86  17.65      24.92      0.09  11.28 1260189
## 250  3523 34.29 2.64  29.12      39.47      0.08  16.05 1192337
## 251  3524 27.13 1.89  23.42      30.84      0.07  9.53  1208650
## 252  3525 25.67 2.19  21.37      29.96      0.09  15.15 1380081
## 253  3526 15.88 1.76  12.42      19.33      0.11  10.39 1023410
## 254  3527 34.92 3.11  28.82      41.02      0.09  19.12 1029037
## 255  3528 28.33 3.01  22.44      34.22      0.11  17.84  920720
## 256  3529 22.80 2.66  17.60      28.01      0.12  19.31 1105012
## 257  3571 24.60 1.75  21.16      28.04      0.07  2.13   295682
## 258  3572 23.55 2.15  19.33      27.77      0.09  1.64   146199
## 259  3573 24.26 1.95  20.43      28.09      0.08  8.06   891157
## 260  3574 29.45 2.83  23.91      34.98      0.10  4.13   246613
## 261  3575 34.97 2.91  29.27      40.67      0.08  3.35   206634
## 262  3576 27.96 2.47  23.12      32.80      0.09  1.76   133420
## 263  3577 33.83 2.96  28.03      39.62      0.09  3.04   178848
## 264  3578 27.09 1.93  23.31      30.87      0.07  24.17 2942845
## 265  3579 28.06 2.95  22.27      33.85      0.11  4.06   215690
## 266  3601 25.13 2.16  20.90      29.36      0.09  13.34 1235309
## 267  3602 26.05 2.42  21.31      30.80      0.09  17.54 1323472
## 268  3603 28.38 2.57  23.34      33.43      0.09  57.72 4065537
## 269  3604 29.33 2.97  23.51      35.15      0.10  28.31 1527436
## 270  3671 24.30 2.02  20.35      28.25      0.08  22.86 2372958
## 271  3672 26.14 1.96  22.30      29.99      0.08  4.01   462145
## 272  3673 31.46 2.23  27.08      35.84      0.07  7.42   736041
## 273  3674 16.77 2.00  12.85      20.69      0.12  22.91 1834112
## 274  5101 15.02 1.82  11.46      18.59      0.12  3.23   286224
## 275  5102 28.15 2.66  22.93      33.36      0.09  6.97   457786
## 276  5103 17.38 1.94  13.58      21.17      0.11  8.63   757511
## 277  5104 17.44 1.55  14.40      20.47      0.09  3.89   535763
## 278  5105 10.47 2.11  6.35       14.60      0.20  3.76   182631
## 279  5106 24.32 2.00  20.40      28.23      0.08  2.20   232846
## 280  5107 29.06 2.79  23.60      34.53      0.10  6.99   425476
## 281  5108 20.18 2.25  15.77      24.59      0.11  9.28   678110
## 282  5171 14.97 1.64  11.75      18.19      0.11  9.71  1053186
## 283  5201 38.22 2.76  32.82      43.62      0.07  10.70  763556
## 284  5202 45.29 2.90  39.61      50.97      0.06  14.96 1013580
## 285  5203 38.23 2.61  33.11      43.34      0.07  16.03 1275192
## 286  5204 35.73 3.08  29.69      41.77      0.09  8.83   489640
## 287  5205 34.92 2.82  29.38      40.45      0.08  4.23   276639
## 288  5206 40.09 2.75  34.69      45.49      0.07  7.22   525378
## 289  5207 26.71 2.71  21.40      32.02      0.10  2.82   172807
## 290  5208 40.35 3.46  33.57      47.12      0.09  5.09   235452
## 291  5271 31.92 2.89  26.25      37.59      0.09  9.17   546815
## 292  5272 51.35 4.01  43.49      59.21      0.08  5.46   194801
## 293  5301 12.66 1.63  9.46       15.85      0.13  1.45   138472
## 294  5302 35.42 3.11  29.32      41.53      0.09  5.02   272057
## 295  5303 25.65 2.89  19.99      31.32      0.11  8.07   422755
## 296  5304 29.59 3.04  23.64      35.54      0.10  9.16   475667
## 297  5305 29.05 3.46  22.26      35.84      0.12  6.68   263458
## 298  5306 23.47 3.43  16.73      30.20      0.15  6.72   235028
## 299  5307 31.00 2.82  25.47      36.53      0.09  3.42   210782
## 300  5308 38.11 3.40  31.44      44.78      0.09  3.33   155726
## 301  5309 32.07 2.89  26.41      37.73      0.09  4.41   264473
## 302  5310 55.26 3.30  48.78      61.73      0.06  6.27   325896
## 303  5311 25.35 3.58  18.34      32.36      0.14  8.14   276287
## 304  5312 34.23 3.52  27.33      41.13      0.10  4.13   172461
## 305  5313 23.06 2.40  18.36      27.75      0.10  5.07   359770
## 306  5314 21.49 2.65  16.29      26.68      0.12  3.72   205168
## 307  5315 23.30 3.14  17.14      29.47      0.13  7.29   302765
## 308  5316 26.79 3.51  19.90      33.68      0.13  2.12    77199
## 309  5317 31.88 3.80  24.43      39.33      0.12  10.88  375128
## 310  5318 27.09 2.91  21.39      32.80      0.11  2.83   151890
## 311  5319 25.84 2.34  21.25      30.42      0.09  3.72   299568
## 312  5320 42.89 3.32  36.38      49.39      0.08  2.16   110135
## 313  5321 20.94 3.50  14.08      27.80      0.17  6.61   205048
## 314  5371 30.92 3.47  24.13      37.71      0.11  12.00  489950
## 315  6101 27.85 2.50  22.95      32.75      0.09  7.45   549016
## 316  6102 25.39 3.41  18.71      32.08      0.13  7.29   272485
## 317  6103 29.62 3.50  22.76      36.48      0.12  10.08  393925
## 318  6104 27.72 2.71  22.41      33.04      0.10  4.37   273537
## 319  6105 24.90 3.03  18.96      30.84      0.12  10.54  492788
## 320  6106 23.10 3.94  15.37      30.83      0.17  20.87  547851
## 321  6107 17.99 1.78  14.50      21.47      0.10  4.09   437955
## 322  6108 17.78 3.52  10.89      24.68      0.20  10.34  280510
## 323  6109 17.80 4.34  9.29       26.31      0.24  11.57  206289
## 324  6110 27.28 2.40  22.57      31.99      0.09  2.79   219987
## 325  6111 34.07 3.64  26.93      41.22      0.11  3.08   119523
## 326  6112 20.39 2.15  16.18      24.60      0.11  7.48   605260
## 327  6171 25.69 1.75  22.27      29.12      0.07  4.74   680335
## 328  6172 28.52 2.48  23.66      33.38      0.09  3.13   238824
## 329  6201 15.37 2.05  11.34      19.39      0.13  4.62   327077
## 330  6202 20.30 2.46  15.48      25.13      0.12  7.68   471672
## 331  6203 28.13 2.55  23.13      33.14      0.09  5.21   370764
## 332  6204 19.09 2.82  13.56      24.62      0.15  3.19   142403
## 333  6205 19.22 2.72  13.89      24.54      0.14  2.78   134128
## 334  6206 15.34 2.25  10.93      19.75      0.15  1.19    69984
## 335  6207 26.61 3.24  20.26      32.96      0.12  2.01    85729
## 336  6208 17.43 3.31  10.94      23.93      0.19  7.54   226980
## 337  6209 19.38 3.02  13.46      25.31      0.16  4.62   181731
## 338  6210 19.50 3.30  13.03      25.98      0.17  3.92   129610
## 339  6211 22.31 3.69  15.08      29.55      0.17  4.14   120814
## 340  6212 26.85 2.36  22.23      31.48      0.09  1.62   131383
## 341  6213 7.40  1.38  4.70       10.10      0.19  1.48   122432
## 342  6271 31.00 3.70  23.74      38.25      0.12  8.48   304042
## 343  6301 27.72 3.09  21.66      33.78      0.11  7.43   357940
## 344  6302 27.38 2.75  21.98      32.78      0.10  5.96   358669
## 345  6303 35.83 2.68  30.58      41.07      0.07  8.29   611483
## 346  6304 33.33 3.41  26.65      40.01      0.10  7.32   321807
## 347  6305 27.24 3.19  20.98      33.49      0.12  4.41   196833
## 348  6306 46.43 2.96  40.62      52.23      0.06  3.70   240686
## 349  6307 29.98 2.67  24.75      35.21      0.09  4.08   275694
## 350  6308 34.09 3.75  26.74      41.44      0.11  6.60   242037
## 351  6309 26.87 2.89  21.21      32.53      0.11  4.90   265169
## 352  6310 22.22 2.53  17.26      27.18      0.11  6.94   429800
## 353  6311 40.36 3.75  33.01      47.72      0.09  3.49   137039
## 354  6371 19.26 2.18  14.98      23.53      0.11  9.70   728565
## 355  6372 26.13 3.05  20.14      32.11      0.12  6.14   291856
## 356  6401 14.49 2.56  9.48       19.50      0.18  6.73   292948
## 357  6402 17.68 2.60  12.58      22.77      0.15  3.06   151598
## 358  6403 18.61 2.68  13.36      23.86      0.14  16.74  810045
## 359  6404 16.60 2.57  11.56      21.63      0.15  8.75   421350
## 360  6405 24.44 3.32  17.93      30.95      0.14  6.33   243275
## 361  6409 23.67 2.34  19.08      28.25      0.10  2.13   161099
## 362  6411 17.48 3.89  9.85       25.12      0.22  1.23    26956
## 363  6471 22.16 1.78  18.67      25.65      0.08  5.22   651807
## 364  6472 22.61 2.06  18.57      26.65      0.09  9.29   878196
## 365  6474 16.14 1.86  12.49      19.78      0.12  2.12   190058
## 366  6501 20.08 2.28  15.62      24.54      0.11  1.31    92789
## 367  6502 20.78 2.58  15.72      25.84      0.12  2.40   136172
## 368  6503 27.24 3.36  20.65      33.83      0.12  0.83    33377
## 369  6504 16.52 2.04  12.51      20.52      0.12  2.90   219759
## 370  6571 21.30 2.37  16.65      25.94      0.11  4.01   274830
## 371  7101 30.93 4.62  21.88      39.98      0.15  11.34  260940
## 372  7102 29.44 3.55  22.48      36.41      0.12  9.23   349073
## 373  7103 25.84 2.75  20.45      31.23      0.11  2.28   132711
## 374  7104 13.52 2.15  9.31       17.74      0.16  1.63    94812
## 375  7105 36.74 2.59  31.67      41.81      0.07  2.68   213329
## 376  7106 23.05 3.73  15.74      30.36      0.16  7.03   206057
## 377  7107 21.10 2.29  16.60      25.60      0.11  1.15    83696
## 378  7108 12.17 2.08  8.10       16.24      0.17  1.19    67675
## 379  7109 26.89 3.25  20.52      33.26      0.12  2.51   107435
## 380  7110 28.28 2.59  23.20      33.35      0.09  1.00    69542
## 381  7111 22.46 2.68  17.21      27.70      0.12  1.35    75215
## 382  7171 12.36 1.31  9.79       14.93      0.11  3.00   434811
## 383  7172 19.52 1.87  15.86      23.19      0.10  2.24   231236
## 384  7173 22.95 2.18  18.67      27.23      0.10  1.32   112788
## 385  7174 20.13 1.87  16.46      23.80      0.09  1.31   137689
## 386  7201 13.77 1.82  10.21      17.32      0.13  1.46   120640
## 387  7202 23.36 2.82  17.83      28.88      0.12  7.73   399533
## 388  7203 10.46 2.46  5.64       15.28      0.24  3.65   129784
## 389  7204 24.25 2.59  19.18      29.32      0.11  4.42   278467
## 390  7205 25.12 4.03  17.21      33.02      0.16  11.73  311671
## 391  7206 24.87 2.31  20.35      29.40      0.09  3.03   244010
## 392  7207 31.44 3.07  25.42      37.46      0.10  3.36   176209
## 393  7208 22.41 2.59  17.34      27.48      0.12  8.80   525137
## 394  7209 18.00 2.04  14.00      22.01      0.11  1.97   159661
## 395  7210 28.85 2.89  23.18      34.52      0.10  4.39   247029
## 396  7211 8.94  1.73  5.56       12.33      0.19  1.29    80946
## 397  7212 12.98 1.79  9.48       16.48      0.14  1.72   140061
## 398  7271 22.83 2.97  17.01      28.65      0.13  9.03   413815
## 399  7301 19.06 1.77  15.59      22.54      0.09  1.23   138469
## 400  7302 23.14 2.59  18.06      28.21      0.11  6.99   425938
## 401  7303 54.16 4.35  45.63      62.70      0.08  6.36   191317
## 402  7304 21.19 2.33  16.62      25.76      0.11  5.26   371332
## 403  7305 29.17 2.62  24.04      34.31      0.09  4.39   303334
## 404  7306 29.90 3.32  23.40      36.40      0.11  18.41  805302
## 405  7307 24.87 2.58  19.81      29.93      0.10  3.86   248151
## 406  7308 23.19 2.38  18.54      27.85      0.10  5.03   364713
## 407  7309 23.57 2.10  19.46      27.69      0.09  3.69   346013
## 408  7310 24.85 3.73  17.54      32.16      0.15  5.71   176219
## 409  7311 13.83 1.41  11.06      16.59      0.10  5.55   762184
## 410  7312 17.98 2.48  13.13      22.84      0.14  4.17   230367
## 411  7313 25.90 2.98  20.06      31.74      0.11  8.26   410706
## 412  7314 28.64 3.64  21.51      35.77      0.13  8.70   308145
## 413  7315 35.74 3.82  28.26      43.23      0.11  10.70  386562
## 414  7316 17.48 1.53  14.48      20.49      0.09  1.52   213979
## 415  7317 23.69 2.28  19.21      28.16      0.10  4.75   378004
## 416  7318 22.49 2.31  17.96      27.03      0.10  3.26   243703
## 417  7322 25.25 2.70  19.97      30.54      0.11  5.44   324294
## 418  7325 22.15 2.15  17.94      26.37      0.10  3.69   316106
## 419  7326 18.81 2.18  14.53      23.09      0.12  3.30   242012
## 420  7371 17.82 1.85  14.19      21.44      0.10  16.29 1600645
## 421  7372 27.12 2.85  21.54      32.70      0.10  2.69   150576
## 422  7373 23.99 2.39  19.31      28.67      0.10  2.67   196357
## 423  7401 23.32 2.85  17.73      28.92      0.12  2.04   102994
## 424  7402 25.07 2.88  19.43      30.71      0.11  4.43   230808
## 425  7403 18.56 2.53  13.59      23.52      0.14  4.74   256449
## 426  7404 14.49 1.83  10.91      18.08      0.13  3.34   284403
## 427  7405 18.33 2.11  14.20      22.46      0.11  4.27   331033
## 428  7406 22.05 3.16  15.85      28.25      0.14  4.94   195077
## 429  7407 13.35 1.60  10.21      16.49      0.12  1.00   103984
## 430  7408 28.67 3.19  22.41      34.93      0.11  3.55   163316
## 431  7409 39.08 3.31  32.60      45.55      0.08  1.30    65191
## 432  7410 16.28 2.05  12.27      20.30      0.13  0.91    67916
## 433  7411 21.88 3.17  15.67      28.09      0.14  3.79   148309
## 434  7412 18.60 2.02  14.64      22.57      0.11  0.41    34471
## 435  7413 35.08 3.45  28.32      41.84      0.10  1.93    84890
## 436  7414 17.61 2.14  13.42      21.80      0.12  1.29    93899
## 437  7415 20.44 2.33  15.87      25.00      0.11  1.17    80549
## 438  7471 20.70 2.74  15.33      26.06      0.13  8.54   429880
## 439  7472 18.40 2.05  14.38      22.41      0.11  2.18   179150
## 440  7501 37.67 2.82  32.14      43.21      0.07  2.61   176819
## 441  7502 28.18 2.79  22.71      33.66      0.10  6.32   376901
## 442  7503 38.42 3.21  32.14      44.70      0.08  3.20   169035
## 443  7504 30.66 2.61  25.53      35.78      0.09  2.24   159938
## 444  7505 36.58 3.89  28.95      44.21      0.11  3.20   112601
## 445  7571 37.52 2.56  32.50      42.55      0.07  2.76   226078
## 446  7601 21.08 2.53  16.12      26.05      0.12  3.00   179010
## 447  7602 22.39 2.90  16.72      28.07      0.13  9.45   449550
## 448  7603 22.93 2.62  17.79      28.07      0.11  2.82   166208
## 449  7604 23.37 2.70  18.08      28.66      0.12  5.61   316426
## 450  7605 30.14 3.12  24.03      36.26      0.10  3.84   190805
## 451  7606 20.37 2.52  15.43      25.31      0.12  2.47   144584
## 452  8101 9.04  1.42  6.26       11.82      0.16  1.20   113067
## 453  8102 20.60 2.86  14.99      26.21      0.14  2.21   101176
## 454  8103 14.71 2.08  10.63      18.80      0.14  5.68   377069
## 455  8104 14.15 2.45  9.35       18.95      0.17  3.34   155550
## 456  8105 5.31  0.66  4.02       6.60       0.12  0.36    95951
## 457  8106 14.43 1.63  11.23      17.62      0.11  1.63   173804
## 458  8107 13.04 1.48  10.15      15.93      0.11  0.97   115774
## 459  8108 13.42 1.87  9.75       17.09      0.14  0.97    73830
## 460  8109 15.14 1.74  11.72      18.56      0.12  0.67    64430
## 461  8171 15.87 2.06  11.83      19.91      0.13  6.98   503334
## 462  8172 16.55 1.95  12.73      20.37      0.12  0.97    80816
## 463  8201 14.27 2.25  9.85       18.69      0.16  2.20   121645
## 464  8202 17.78 3.57  10.78      24.77      0.20  2.28    60113
## 465  8203 13.11 2.20  8.79       17.42      0.17  1.94   104440
## 466  8204 12.29 1.35  9.65       14.94      0.11  1.79   242888
## 467  8205 11.97 1.33  9.36       14.59      0.11  1.48   201133
## 468  8206 13.48 2.19  9.18       17.77      0.16  1.85   103010
## 469  8207 14.30 1.87  10.63      17.97      0.13  0.90    72115
## 470  8208 11.38 1.57  8.31       14.46      0.14  0.59    55587
## 471  8271 14.76 2.36  10.13      19.39      0.16  4.85   250692
## 472  8272 10.21 1.35  7.57       12.85      0.13  0.88   102479
## 473  9101 9.56  1.36  6.89       12.24      0.14  0.79    84670
## 474  9102 15.72 1.66  12.46      18.98      0.11  0.63    69825
## 475  9103 11.98 2.47  7.14       16.82      0.21  0.91    36180
## 476  9104 14.63 1.96  10.78      18.47      0.13  0.96    71622
## 477  9105 16.32 2.68  11.08      21.57      0.16  4.12   180391
## 478  9111 14.52 2.22  10.18      18.86      0.15  0.49    28466
## 479  9112 12.39 2.04  8.39       16.38      0.16  0.62    37379
## 480  9201 16.75 1.67  13.48      20.03      0.10  0.45    51358
## 481  9202 17.67 1.88  13.98      21.36      0.11  1.06    99948
## 482  9203 25.15 2.86  19.55      30.75      0.11  1.00    52715
## 483  9204 18.48 3.30  12.01      24.94      0.18  1.44    45766
## 484  9205 8.62  1.47  5.75       11.50      0.17  0.20    16784
## 485  9271 16.42 2.08  12.35      20.50      0.13  4.07   297472
## 486  9403 13.48 1.89  9.77       17.18      0.14  1.75   131084
## 487  9408 12.86 2.14  8.66       17.07      0.17  1.81   101572
## 488  9409 18.66 2.52  13.72      23.59      0.14  2.79   153174
## 489  9419 10.31 1.47  7.42       13.20      0.14  0.43    41772
## 490  9420 16.30 1.99  12.41      20.20      0.12  0.76    60132
## 491  9426 14.14 1.65  10.91      17.37      0.12  0.34    34503
## 492  9427 8.67  1.22  6.27       11.06      0.14  0.18    21882
## 493  9428 4.30  1.05  2.25       6.35       0.24  0.30    26054
## 494  9471 14.74 1.94  10.93      18.54      0.13  4.13   315594
## 495  9501 23.00 2.13  18.82      27.17      0.09  2.71   242745
## 496  9502 27.60 3.65  20.45      34.74      0.13  2.12    73053
## 497  9503 27.80 3.62  20.71      34.90      0.13  3.01   105774
## 498  9504 6.46  1.03  4.44       8.48       0.16  0.82   107134
## 499  9601 8.17  1.60  5.04       11.30      0.20  3.36   226309
## 500  9602 24.11 6.02  12.31      35.91      0.25  8.74   101298
## 501  9603 2.28  0.89  0.54       4.03       0.39  1.22    79077
## 502  9604 17.23 2.56  12.21      22.26      0.15  3.01   150099
## 503  9605 5.48  2.15  1.28       9.69       0.39  6.72   173546
## 504  9606 3.71  1.05  1.64       5.77       0.28  0.73    53732
## 505  9607 1.39  0.52  0.37       2.41       0.37  1.00   116660
## 506  9608 20.22 5.67  9.10       31.34      0.28  11.79  135698
## 507  9701 7.74  1.12  5.54       9.95       0.15  0.92   119680
## 508  9702 8.78  1.90  5.06       12.49      0.22  4.14   211757
## 509  9703 1.31  0.49  0.35       2.27       0.37  1.71   211298
## 510  9704 9.37  1.21  6.99       11.75      0.13  1.18   156112
## 511  9705 14.14 0.96  12.25      16.02      0.07  0.19    57255
## 512  9706 10.64 2.77  5.22       16.07      0.26  2.50    71415
## 513  9707 14.55 2.16  10.33      18.78      0.15  3.29   202023
## 514  9708 10.25 1.67  6.98       13.51      0.16  1.04    79001

Visualisasi

ggplot(data = dataK, aes(x = KAKO)) +
  geom_line(aes(y = CV*100, color = "RSE", group = 1), size = 1) +      
  geom_line(aes(y = EST, color = "Theta", group = 1), size = 1) +       
  geom_hline(yintercept = 25, color = "red", linetype = "dashed") +     
  labs(
    title = "Provinsi vs Nilai RSE dan Theta - SPSS",
    x = "Provinsi",
    y = "Nilai"
  ) +
  scale_color_manual(name = " ", values = c("RSE" = "magenta", "Theta" = "purple")) +
  theme_minimal()

5 Kesimpulan

Berdasarkan hasil tampak bahwa nilai SE dan besaran lainnya (kecuali DEFF = Desain Efek) baik penghitungan dengan R maupun SPSS sudah memberikan hasil yang sama. Oleh karena itu, jika ingin menghitung dari salah satu software (R atau SPSS), maka ikuti sintaks desain sampling yang diberikan.

6 References

[1] Badan Pusat Statistik. (2023). Konsep dan Definisi Susenas Maret 2023. Jakarta: BPS.

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


Direktorat Statistik Kesejahteraan Rakyat, BPS,