Unmet Need
using direct
method 2024
Sapta Hastho Ponco
2025-04-02
Load dan install packages yang akan digunakan
pacman::p_load(haven, dplyr, tidyr, survey, janitor, magrittr, ggplot2)1 Load Data
Data yang akan digunakan yaitu data Susenas KOR 2024 Maret, sampel individu
datakor <- read_sav("KOR24GAB_0208.sav")2 Variabel Analisis
Untuk meringankan beban kerja komputasi, maka diambil beberapa variabel saja.
Variabel yang digunakan antara lain:
kode provinsi (R101), kode kab/kota (R102), kode renumbering NKS (NURT), kode renumbering no urut ruta, strata (R105), klasifikasi kota desa,
dataku <- datakor %>%
select(R101, R102, WI1, WI2, PSU, SSU, STRATA, R105, R407, R1102, R1103, R1105, FWT, fpc1, fpc2)3 Filter Variabel
Seleksi berdasarkan variabel yang akan digunakan
Dengan kategori UNPK berdasarkan kondisi R1102 = 1 dan R1103 = 1 dan R1105 = 5 –> UNPK (atau diberi kategori 1)
library(dplyr)
indicator = dataku %>%
mutate(
kabu = R101*100+R102,
jml_pddk = 100,
theta =
dplyr::case_when(
((dataku$R1102 == 1 & dataku$R1103 == 1 & dataku$R1105 == 5))~ 100,))%>%
mutate(theta = ifelse(is.na(theta), 0, theta)) 4 Set Design Sampling
Gunakan penimbang yang diberikan dalam Susenas 2023 KOR Maret
susenas.design<- svydesign(id=~PSU, strata=~STRATA,
data = indicator, weights=~FWT, 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.86447585 Set Output Statistics
Set Output Statistics yang akan dikeluarkan
options(survey.adjust.domain.lonely=TRUE)
options(survey.lonely.psu="adjust")
#Provinsi
resultP <- svyby(
formula = ~theta,
denom = ~jml_pddk,
by = ~R101,
design = susenas.design,
deff = TRUE,
svyratio,
vartype = c("se", "ci", "ci", "cv", "cvpct", "var"))## Warning in matrix(1:(ns * reps), ncol = reps, nrow = ns, byrow = TRUE):
## non-empty data for zero-extent matrixresultP[is.na(resultP)] <- 0
resultP <- resultP %>%
mutate(
theta = round(`theta/jml_pddk` * 100, 2),
SE = round(`se.theta/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] <- 0options(survey.adjust.domain.lonely=TRUE)
options(survey.lonely.psu="adjust")
#Kabupaten/kota
resultK <- svyby(
formula = ~theta,
denom = ~jml_pddk,
by = ~kabu,
design = susenas.design,
deff = TRUE,
svyratio,
vartype = c("se", "ci", "ci", "cv", "cvpct", "var"))## Warning in matrix(1:(ns * reps), ncol = reps, nrow = ns, byrow = TRUE):
## non-empty data for zero-extent matrixresultK[is.na(resultK)] <- 0
resultK <- resultK %>%
mutate(
theta = round(`theta/jml_pddk` * 100, 2),
SE = round(`se.theta/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] <- 06 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):
Jika \(RSE \le 25\%\), estimasi bersifat presisi.
Jika \(25\% < RSE \le 50\%\), estimasi perlu dilakukan dengan hati-hati.
Jika \(RSE > 50\%\), estimasi dianggap sangat tidak presisi.
6.1 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 5.24 0.34 0.12 4.58 5.90 6.44 11.38
## 12 12 4.62 0.23 0.05 4.17 5.08 4.99 9.44
## 13 13 4.80 0.29 0.08 4.24 5.36 5.96 7.79
## 14 14 4.82 0.34 0.12 4.14 5.49 7.15 7.96
## 15 15 4.84 0.39 0.15 4.07 5.61 8.08 8.01
## 16 16 4.92 0.28 0.08 4.37 5.47 5.68 6.59
## 17 17 6.35 0.64 0.41 5.10 7.60 10.04 13.80
## 18 18 5.19 0.36 0.13 4.48 5.89 6.93 9.34
## 19 19 4.78 0.59 0.35 3.63 5.92 12.28 10.42
## 21 21 4.61 0.67 0.45 3.30 5.91 14.48 15.10
## 31 31 3.06 0.36 0.13 2.36 3.76 11.66 8.18
## 32 32 5.74 0.24 0.06 5.26 6.22 4.26 9.35
## 33 33 5.70 0.21 0.04 5.30 6.11 3.62 7.89
## 34 34 4.32 0.34 0.12 3.66 4.97 7.78 3.60
## 35 35 5.05 0.19 0.04 4.67 5.43 3.82 8.12
## 36 36 4.12 0.39 0.15 3.35 4.89 9.56 10.06
## 51 51 2.04 0.21 0.04 1.62 2.46 10.52 5.41
## 52 52 8.24 0.49 0.24 7.27 9.20 6.00 7.42
## 53 53 7.06 0.33 0.11 6.42 7.70 4.61 8.19
## 61 61 5.53 0.39 0.15 4.76 6.30 7.10 9.54
## 62 62 4.97 0.44 0.19 4.10 5.83 8.88 11.08
## 63 63 6.27 0.41 0.17 5.46 7.07 6.55 7.77
## 64 64 3.92 0.46 0.21 3.02 4.81 11.66 12.28
## 65 65 4.11 0.49 0.24 3.16 5.06 11.81 6.33
## 71 71 5.08 0.36 0.13 4.37 5.79 7.10 7.88
## 72 72 7.28 0.53 0.28 6.24 8.32 7.28 11.42
## 73 73 5.34 0.26 0.07 4.83 5.85 4.89 7.56
## 74 74 6.38 0.34 0.12 5.71 7.04 5.34 6.86
## 75 75 7.71 0.60 0.36 6.54 8.87 7.73 6.38
## 76 76 6.24 0.51 0.26 5.24 7.24 8.20 6.28
## 81 81 3.73 0.31 0.10 3.12 4.33 8.27 7.00
## 82 82 4.22 0.36 0.13 3.52 4.93 8.51 7.01
## 91 91 3.66 0.52 0.27 2.64 4.68 14.22 11.50
## 92 92 4.08 0.77 0.59 2.56 5.60 19.00 20.09
## 94 94 3.02 0.40 0.16 2.23 3.81 13.27 9.69
## 95 95 4.37 0.60 0.36 3.19 5.54 13.75 7.86
## 96 96 1.39 0.34 0.12 0.72 2.06 24.65 12.18
## 97 97 2.27 0.30 0.09 1.68 2.86 13.32 5.85Visualisasi
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",
x = "Provinsi",
y = "Nilai"
) +
scale_color_manual(name = "Legend", values = c("RSE" = "orange", "Theta" = "blue")) +
theme_minimal()## Warning: Using `size` aesthetic for lines was deprecated in ggplot2 3.4.0.
## ℹ Please use `linewidth` instead.
## This warning is displayed once every 8 hours.
## Call `lifecycle::last_lifecycle_warnings()` to see where this warning was
## generated.
6.2 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 5.36 1.09 1.19 3.21 7.50 20.44 5.07
## 1102 1102 4.24 0.82 0.67 2.62 5.85 19.42 3.88
## 1103 1103 2.10 0.48 0.23 1.15 3.05 23.02 2.53
## 1104 1104 12.68 2.85 8.12 7.09 18.27 22.50 14.48
## 1105 1105 8.02 2.13 4.54 3.85 12.20 26.56 15.24
## 1106 1106 8.32 2.06 4.24 4.28 12.35 24.76 11.27
## 1107 1107 8.18 2.19 4.80 3.89 12.47 26.78 12.85
## 1108 1108 1.28 0.40 0.16 0.51 2.05 30.87 3.30
## 1109 1109 4.42 1.10 1.21 2.27 6.57 24.86 7.02
## 1110 1110 3.90 1.01 1.02 1.91 5.88 26.01 6.72
## 1111 1111 4.74 1.13 1.28 2.52 6.96 23.87 7.10
## 1112 1112 5.74 1.83 3.35 2.15 9.33 31.90 13.89
## 1113 1113 5.56 1.58 2.50 2.47 8.65 28.33 8.88
## 1114 1114 5.57 1.64 2.69 2.36 8.78 29.42 11.36
## 1115 1115 6.20 1.65 2.72 2.97 9.43 26.56 8.65
## 1116 1116 6.11 1.38 1.90 3.41 8.81 22.52 6.26
## 1117 1117 6.29 1.08 1.17 4.17 8.40 17.16 3.81
## 1118 1118 2.84 1.00 1.00 0.89 4.79 35.08 6.34
## 1171 1171 4.57 0.74 0.55 3.12 6.03 16.25 2.91
## 1172 1172 3.36 0.86 0.74 1.67 5.05 25.71 4.11
## 1173 1173 6.72 1.48 2.19 3.83 9.61 21.96 7.22
## 1174 1174 3.39 0.86 0.74 1.70 5.07 25.39 5.01
## 1175 1175 4.24 0.88 0.77 2.51 5.97 20.85 4.39
## 1201 1201 9.31 1.66 2.76 6.06 12.56 17.81 8.10
## 1202 1202 2.90 0.54 0.29 1.85 3.95 18.52 2.71
## 1203 1203 5.27 1.41 1.99 2.51 8.03 26.71 8.70
## 1204 1204 3.65 0.69 0.48 2.30 5.00 18.89 3.33
## 1205 1205 4.85 1.21 1.46 2.47 7.23 25.04 7.83
## 1206 1206 2.21 0.50 0.25 1.24 3.19 22.49 2.13
## 1207 1207 1.79 0.48 0.23 0.84 2.73 27.05 3.12
## 1208 1208 5.21 1.10 1.21 3.05 7.36 21.13 6.83
## 1209 1209 8.56 1.48 2.19 5.65 11.47 17.33 7.59
## 1210 1210 8.60 2.23 4.97 4.22 12.97 25.97 13.32
## 1211 1211 3.48 0.87 0.76 1.77 5.19 25.08 5.20
## 1212 1212 4.42 0.73 0.53 2.99 5.86 16.56 4.61
## 1213 1213 6.16 1.23 1.51 3.76 8.57 19.92 7.23
## 1214 1214 2.10 0.55 0.30 1.02 3.18 26.22 3.58
## 1215 1215 2.91 0.66 0.44 1.63 4.20 22.50 3.15
## 1216 1216 6.96 1.57 2.46 3.88 10.05 22.61 7.92
## 1217 1217 6.47 1.32 1.74 3.88 9.05 20.41 5.74
## 1218 1218 7.18 1.83 3.35 3.60 10.76 25.42 13.05
## 1219 1219 7.26 1.38 1.90 4.55 9.96 19.01 6.90
## 1220 1220 4.47 0.94 0.88 2.63 6.30 21.01 4.44
## 1221 1221 5.34 1.11 1.23 3.16 7.52 20.81 5.24
## 1222 1222 5.41 1.00 1.00 3.45 7.36 18.47 4.71
## 1223 1223 5.93 1.56 2.43 2.86 8.99 26.38 10.22
## 1224 1224 6.71 1.69 2.86 3.40 10.02 25.17 10.62
## 1225 1225 12.42 2.26 5.11 7.98 16.85 18.23 11.39
## 1271 1271 1.75 0.44 0.19 0.90 2.61 24.91 2.40
## 1272 1272 2.74 0.69 0.48 1.39 4.09 25.17 3.63
## 1273 1273 4.84 1.19 1.42 2.52 7.17 24.52 6.58
## 1274 1274 3.07 0.69 0.48 1.72 4.43 22.49 3.00
## 1275 1275 1.87 0.46 0.21 0.97 2.77 24.63 4.19
## 1276 1276 4.64 1.04 1.08 2.60 6.69 22.44 4.67
## 1277 1277 1.35 0.29 0.08 0.79 1.92 21.20 1.35
## 1278 1278 10.20 2.33 5.43 5.62 14.78 22.89 14.46
## 1301 1301 8.79 1.94 3.76 4.99 12.60 22.09 8.65
## 1302 1302 3.59 0.70 0.49 2.22 4.96 19.49 3.87
## 1303 1303 4.52 1.00 1.00 2.57 6.47 22.04 5.46
## 1304 1304 4.93 1.19 1.42 2.61 7.26 24.02 6.59
## 1305 1305 4.88 1.01 1.02 2.90 6.86 20.72 5.23
## 1306 1306 4.14 1.09 1.19 2.01 6.28 26.27 7.49
## 1307 1307 3.25 0.64 0.41 1.99 4.51 19.81 3.65
## 1308 1308 5.83 1.38 1.90 3.12 8.54 23.73 7.57
## 1309 1309 4.97 1.14 1.30 2.73 7.20 23.00 6.20
## 1310 1310 10.97 1.92 3.69 7.20 14.73 17.52 7.75
## 1311 1311 9.20 1.45 2.10 6.35 12.05 15.81 5.95
## 1312 1312 2.79 1.06 1.12 0.71 4.87 38.08 10.38
## 1371 1371 4.40 0.90 0.81 2.64 6.15 20.37 5.98
## 1372 1372 3.19 0.73 0.53 1.77 4.62 22.77 3.57
## 1373 1373 3.92 1.02 1.04 1.92 5.93 26.07 5.64
## 1374 1374 3.41 0.90 0.81 1.64 5.18 26.51 4.86
## 1375 1375 6.62 1.27 1.61 4.12 9.12 19.25 5.43
## 1376 1376 4.07 0.86 0.74 2.38 5.76 21.22 3.99
## 1377 1377 3.82 0.78 0.61 2.30 5.35 20.31 3.61
## 1401 1401 1.35 0.48 0.23 0.41 2.29 35.65 3.70
## 1402 1402 7.26 1.37 1.88 4.58 9.94 18.85 6.95
## 1403 1403 8.48 1.53 2.34 5.48 11.49 18.07 8.41
## 1404 1404 2.11 0.75 0.56 0.64 3.59 35.62 6.64
## 1405 1405 4.21 0.92 0.85 2.41 6.01 21.81 4.94
## 1406 1406 3.78 0.87 0.76 2.08 5.47 22.91 5.94
## 1407 1407 3.93 1.16 1.35 1.66 6.19 29.43 9.31
## 1408 1408 3.25 1.07 1.14 1.16 5.35 32.84 9.76
## 1409 1409 7.60 1.59 2.53 4.49 10.71 20.86 9.92
## 1410 1410 8.44 1.52 2.31 5.47 11.41 17.96 6.16
## 1471 1471 4.55 0.75 0.56 3.07 6.02 16.56 4.11
## 1473 1473 2.60 0.49 0.24 1.64 3.56 18.82 2.27
## 1501 1501 8.01 1.53 2.34 5.02 11.00 19.07 6.12
## 1502 1502 5.60 1.56 2.43 2.54 8.66 27.89 10.51
## 1503 1503 4.76 1.11 1.23 2.59 6.94 23.33 5.85
## 1504 1504 5.29 1.80 3.24 1.76 8.82 34.06 14.57
## 1505 1505 3.71 0.96 0.92 1.82 5.60 25.97 6.31
## 1506 1506 5.79 1.21 1.46 3.42 8.17 20.90 5.30
## 1507 1507 3.09 0.79 0.62 1.55 4.64 25.45 4.49
## 1508 1508 4.89 1.07 1.14 2.80 6.99 21.83 5.40
## 1509 1509 7.04 1.38 1.90 4.33 9.75 19.64 6.83
## 1571 1571 2.96 1.03 1.06 0.93 4.99 34.93 9.74
## 1572 1572 5.15 0.95 0.90 3.29 7.01 18.46 3.34
## 1601 1601 5.09 0.99 0.98 3.14 7.04 19.54 4.99
## 1602 1602 2.98 0.52 0.27 1.95 4.00 17.55 2.57
## 1603 1603 7.81 1.42 2.02 5.02 10.59 18.20 8.06
## 1604 1604 4.93 0.87 0.76 3.21 6.64 17.76 4.06
## 1605 1605 3.35 0.69 0.48 2.00 4.71 20.58 2.77
## 1606 1606 4.89 1.19 1.42 2.56 7.22 24.31 7.45
## 1607 1607 6.07 1.13 1.28 3.85 8.28 18.62 6.34
## 1608 1608 4.30 1.45 2.10 1.46 7.15 33.73 9.89
## 1609 1609 4.57 0.96 0.92 2.68 6.45 21.06 5.14
## 1610 1610 8.52 1.57 2.46 5.44 11.60 18.42 7.59
## 1611 1611 6.35 1.13 1.28 4.13 8.57 17.81 3.98
## 1612 1612 7.44 2.61 6.81 2.32 12.56 35.14 19.16
## 1613 1613 7.75 1.87 3.50 4.10 11.41 24.07 8.56
## 1671 1671 2.90 0.58 0.34 1.76 4.04 20.05 4.02
## 1672 1672 5.79 1.65 2.72 2.55 9.03 28.55 10.43
## 1673 1673 4.84 1.25 1.56 2.39 7.28 25.78 6.64
## 1674 1674 5.52 1.08 1.17 3.40 7.65 19.64 4.76
## 1701 1701 7.64 1.98 3.92 3.75 11.52 25.98 11.27
## 1702 1702 5.29 1.09 1.19 3.15 7.43 20.64 4.80
## 1703 1703 6.47 1.43 2.04 3.67 9.28 22.09 7.20
## 1704 1704 4.19 0.78 0.61 2.66 5.72 18.65 3.03
## 1705 1705 6.72 1.88 3.53 3.04 10.41 27.96 10.82
## 1706 1706 6.40 1.27 1.61 3.90 8.90 19.91 5.82
## 1707 1707 5.52 1.06 1.12 3.44 7.61 19.26 3.96
## 1708 1708 4.36 0.93 0.86 2.54 6.18 21.34 3.85
## 1709 1709 6.11 1.54 2.37 3.08 9.13 25.30 7.86
## 1771 1771 7.97 2.62 6.86 2.84 13.11 32.88 22.60
## 1801 1801 9.64 2.60 6.76 4.54 14.74 27.00 14.66
## 1802 1802 5.36 0.97 0.94 3.47 7.26 18.02 4.75
## 1803 1803 3.88 1.46 2.13 1.03 6.74 37.54 15.77
## 1804 1804 4.11 1.08 1.17 2.00 6.21 26.21 8.09
## 1805 1805 5.77 1.01 1.02 3.78 7.76 17.56 5.44
## 1806 1806 5.37 1.04 1.08 3.34 7.40 19.31 5.44
## 1807 1807 8.85 1.81 3.28 5.29 12.40 20.50 9.47
## 1808 1808 5.39 1.00 1.00 3.43 7.35 18.55 4.59
## 1809 1809 5.51 1.71 2.92 2.17 8.85 30.94 12.67
## 1810 1810 6.51 1.31 1.72 3.95 9.07 20.07 6.81
## 1811 1811 9.38 1.73 2.99 5.99 12.76 18.41 6.33
## 1812 1812 4.63 1.26 1.59 2.15 7.11 27.31 7.16
## 1813 1813 3.04 1.07 1.14 0.94 5.13 35.21 7.87
## 1871 1871 3.24 0.85 0.72 1.56 4.92 26.38 6.86
## 1872 1872 2.51 0.59 0.35 1.36 3.67 23.49 2.85
## 1901 1901 5.59 2.15 4.62 1.38 9.80 38.44 18.65
## 1902 1902 4.07 0.64 0.41 2.81 5.33 15.81 1.99
## 1903 1903 4.75 0.99 0.98 2.82 6.69 20.72 4.23
## 1904 1904 5.31 1.14 1.30 3.08 7.53 21.43 5.09
## 1905 1905 3.49 0.89 0.79 1.74 5.24 25.55 4.77
## 1906 1906 4.76 1.31 1.72 2.20 7.32 27.47 6.88
## 1971 1971 4.92 1.14 1.30 2.68 7.16 23.25 5.54
## 2101 2101 4.39 1.03 1.06 2.38 6.40 23.39 5.44
## 2102 2102 2.82 0.61 0.37 1.63 4.01 21.53 2.85
## 2103 2103 4.33 0.76 0.58 2.83 5.82 17.62 2.69
## 2104 2104 2.71 0.77 0.59 1.20 4.21 28.36 3.97
## 2105 2105 2.22 0.50 0.25 1.24 3.20 22.59 2.34
## 2171 2171 5.14 0.97 0.94 3.23 7.05 18.96 5.92
## 2172 2172 3.47 0.80 0.64 1.90 5.03 23.04 3.87
## 3101 3101 5.20 1.27 1.61 2.72 7.68 24.36 6.19
## 3171 3171 2.02 0.41 0.17 1.22 2.82 20.22 3.16
## 3172 3172 3.19 0.86 0.74 1.50 4.88 27.04 9.41
## 3173 3173 2.63 0.59 0.35 1.47 3.79 22.44 3.31
## 3174 3174 2.87 0.75 0.56 1.40 4.35 26.22 7.51
## 3175 3175 4.61 0.92 0.85 2.81 6.42 19.95 6.69
## 3201 3201 5.08 0.92 0.85 3.29 6.88 18.04 7.90
## 3202 3202 6.29 1.01 1.02 4.31 8.27 16.05 6.04
## 3203 3203 7.12 0.91 0.83 5.33 8.91 12.82 4.44
## 3204 3204 5.55 0.90 0.81 3.79 7.30 16.19 6.20
## 3205 3205 8.41 1.08 1.17 6.30 10.52 12.82 5.87
## 3206 3206 6.95 1.05 1.10 4.89 9.00 15.08 5.19
## 3207 3207 11.50 2.02 4.08 7.54 15.46 17.57 12.18
## 3208 3208 7.94 1.27 1.61 5.46 10.43 15.95 6.00
## 3209 3209 9.37 1.76 3.10 5.93 12.81 18.73 12.57
## 3210 3210 7.70 1.13 1.28 5.49 9.91 14.65 5.01
## 3211 3211 2.73 0.52 0.27 1.71 3.75 19.11 2.76
## 3212 3212 9.80 1.66 2.76 6.54 13.05 16.95 9.29
## 3213 3213 4.51 1.09 1.19 2.37 6.64 24.16 7.77
## 3214 3214 6.00 1.08 1.17 3.87 8.12 18.09 5.88
## 3215 3215 4.98 1.14 1.30 2.74 7.22 22.93 9.62
## 3216 3216 4.73 1.23 1.51 2.33 7.13 25.93 13.26
## 3217 3217 5.90 0.78 0.61 4.37 7.44 13.28 3.84
## 3218 3218 9.50 1.48 2.19 6.59 12.40 15.61 5.53
## 3271 3271 4.75 0.93 0.86 2.92 6.58 19.63 5.54
## 3272 3272 5.25 0.94 0.88 3.40 7.09 17.91 4.14
## 3273 3273 2.47 0.42 0.18 1.64 3.29 17.01 2.61
## 3274 3274 4.44 1.04 1.08 2.41 6.47 23.38 5.09
## 3275 3275 3.86 0.79 0.62 2.31 5.41 20.50 6.60
## 3276 3276 1.82 0.44 0.19 0.95 2.69 24.45 4.12
## 3277 3277 2.79 0.48 0.23 1.86 3.73 17.08 2.20
## 3278 3278 9.28 2.31 5.34 4.76 13.80 24.84 15.90
## 3279 3279 3.29 0.65 0.42 2.02 4.56 19.70 2.46
## 3301 3301 5.13 1.03 1.06 3.10 7.15 20.16 7.51
## 3302 3302 4.78 1.30 1.69 2.23 7.34 27.21 13.10
## 3303 3303 9.66 1.21 1.46 7.28 12.04 12.56 5.23
## 3304 3304 6.66 1.38 1.90 3.95 9.37 20.78 8.97
## 3305 3305 13.28 2.22 4.93 8.93 17.63 16.73 13.30
## 3306 3306 6.50 0.99 0.98 4.57 8.43 15.17 4.04
## 3307 3307 8.48 1.15 1.32 6.22 10.74 13.61 4.84
## 3308 3308 7.07 1.03 1.06 5.05 9.09 14.56 4.94
## 3309 3309 3.65 0.64 0.41 2.40 4.90 17.46 3.35
## 3310 3310 5.16 0.82 0.67 3.55 6.77 15.94 4.10
## 3311 3311 3.85 0.80 0.64 2.29 5.41 20.67 4.61
## 3312 3312 4.07 0.70 0.49 2.70 5.43 17.10 3.65
## 3313 3313 2.59 0.57 0.32 1.48 3.71 21.89 3.69
## 3314 3314 5.40 1.70 2.89 2.06 8.74 31.57 15.20
## 3315 3315 3.44 0.64 0.41 2.18 4.69 18.64 4.02
## 3316 3316 5.24 0.71 0.50 3.84 6.64 13.63 2.66
## 3317 3317 4.38 0.80 0.64 2.81 5.95 18.29 3.90
## 3318 3318 3.32 0.95 0.90 1.47 5.18 28.46 7.77
## 3319 3319 5.26 1.83 3.35 1.68 8.85 34.77 18.90
## 3320 3320 7.26 1.84 3.39 3.66 10.86 25.29 14.84
## 3321 3321 5.40 1.15 1.32 3.15 7.65 21.29 7.57
## 3322 3322 3.44 0.75 0.56 1.97 4.90 21.74 5.04
## 3323 3323 5.56 1.02 1.04 3.55 7.56 18.38 5.74
## 3324 3324 5.51 0.97 0.94 3.60 7.42 17.67 5.26
## 3325 3325 4.96 0.83 0.69 3.33 6.58 16.69 4.14
## 3326 3326 7.39 1.09 1.19 5.26 9.53 14.71 5.50
## 3327 3327 8.41 1.07 1.14 6.32 10.50 12.67 4.87
## 3328 3328 7.39 1.19 1.42 5.06 9.72 16.08 6.59
## 3329 3329 6.67 0.77 0.59 5.16 8.18 11.54 3.33
## 3371 3371 3.39 1.02 1.04 1.39 5.38 30.01 5.34
## 3372 3372 5.23 1.09 1.19 3.08 7.37 20.92 5.53
## 3373 3373 4.48 1.03 1.06 2.45 6.50 23.04 4.68
## 3374 3374 2.57 0.61 0.37 1.37 3.76 23.82 5.16
## 3375 3375 6.42 1.42 2.02 3.64 9.20 22.09 7.32
## 3376 3376 6.53 1.69 2.86 3.21 9.85 25.93 9.44
## 3401 3401 2.76 0.54 0.29 1.70 3.82 19.58 2.70
## 3402 3402 4.64 0.72 0.52 3.23 6.06 15.53 3.37
## 3403 3403 3.60 0.52 0.27 2.59 4.62 14.34 2.17
## 3404 3404 4.39 0.68 0.46 3.05 5.72 15.56 3.26
## 3471 3471 6.12 1.05 1.10 4.05 8.18 17.21 4.02
## 3501 3501 6.13 1.58 2.50 3.02 9.23 25.87 10.79
## 3502 3502 7.67 1.00 1.00 5.71 9.62 13.01 3.73
## 3503 3503 4.53 0.74 0.55 3.08 5.99 16.41 3.35
## 3504 3504 8.70 1.79 3.20 5.19 12.22 20.62 11.20
## 3505 3505 6.25 0.83 0.69 4.62 7.88 13.30 3.36
## 3506 3506 7.36 1.46 2.13 4.50 10.22 19.82 10.37
## 3507 3507 4.48 0.84 0.71 2.83 6.13 18.80 6.27
## 3508 3508 3.27 0.63 0.40 2.05 4.50 19.13 3.55
## 3509 3509 4.26 0.76 0.58 2.77 5.76 17.91 5.50
## 3510 3510 8.28 1.09 1.19 6.14 10.41 13.16 5.16
## 3511 3511 4.28 0.84 0.71 2.63 5.93 19.72 4.26
## 3512 3512 5.87 1.18 1.39 3.55 8.19 20.15 5.91
## 3513 3513 4.84 0.97 0.94 2.93 6.75 20.10 6.08
## 3514 3514 3.51 0.79 0.62 1.96 5.07 22.60 6.26
## 3515 3515 2.61 0.52 0.27 1.59 3.64 20.06 3.97
## 3516 3516 4.06 0.94 0.88 2.21 5.90 23.20 6.69
## 3517 3517 4.55 0.83 0.69 2.93 6.17 18.18 4.72
## 3518 3518 4.94 0.87 0.76 3.23 6.64 17.60 4.48
## 3519 3519 5.15 0.95 0.90 3.29 7.01 18.41 4.44
## 3520 3520 5.37 1.21 1.46 3.01 7.73 22.46 7.48
## 3521 3521 5.92 1.24 1.54 3.49 8.34 20.94 6.52
## 3522 3522 4.34 1.16 1.35 2.07 6.61 26.69 10.12
## 3523 3523 4.42 0.78 0.61 2.89 5.94 17.60 4.36
## 3524 3524 2.44 0.54 0.29 1.39 3.50 22.05 3.71
## 3525 3525 2.22 0.56 0.31 1.11 3.32 25.39 4.57
## 3526 3526 4.76 1.08 1.17 2.65 6.87 22.61 7.15
## 3527 3527 9.96 2.37 5.62 5.31 14.61 23.80 17.34
## 3528 3528 5.94 1.43 2.04 3.15 8.74 24.02 10.39
## 3529 3529 4.91 1.64 2.69 1.70 8.13 33.38 15.07
## 3571 3571 4.55 1.24 1.54 2.11 6.99 27.33 7.57
## 3572 3572 3.75 0.72 0.52 2.34 5.15 19.12 2.60
## 3573 3573 5.01 1.00 1.00 3.06 6.97 19.92 5.59
## 3574 3574 5.53 1.42 2.02 2.75 8.31 25.64 7.92
## 3575 3575 3.93 0.84 0.71 2.27 5.58 21.50 3.58
## 3576 3576 3.15 0.65 0.42 1.87 4.43 20.70 2.52
## 3577 3577 2.82 0.74 0.55 1.36 4.28 26.33 3.35
## 3578 3578 6.02 0.98 0.96 4.11 7.93 16.22 6.92
## 3579 3579 4.39 1.10 1.21 2.23 6.55 25.05 5.81
## 3601 3601 6.40 1.10 1.21 4.24 8.56 17.23 6.20
## 3602 3602 5.10 0.83 0.69 3.46 6.73 16.36 4.29
## 3603 3603 5.03 1.12 1.25 2.83 7.24 22.35 10.69
## 3604 3604 5.10 0.92 0.85 3.31 6.90 17.96 5.99
## 3671 3671 1.81 0.37 0.14 1.08 2.55 20.55 2.81
## 3672 3672 2.70 0.64 0.41 1.44 3.95 23.68 3.75
## 3673 3673 4.04 0.63 0.40 2.81 5.28 15.60 2.96
## 3674 3674 2.40 0.66 0.44 1.11 3.70 27.55 6.01
## 5101 5101 2.49 0.67 0.45 1.18 3.80 26.77 3.73
## 5102 5102 2.30 0.53 0.28 1.27 3.33 22.85 3.30
## 5103 5103 1.53 0.48 0.23 0.60 2.47 31.07 4.48
## 5104 5104 0.81 0.28 0.08 0.27 1.36 34.23 2.90
## 5105 5105 2.82 1.30 1.69 0.27 5.37 46.07 13.01
## 5106 5106 2.45 0.55 0.30 1.38 3.53 22.34 2.93
## 5107 5107 3.17 0.92 0.85 1.37 4.97 28.95 6.58
## 5108 5108 1.99 0.61 0.37 0.80 3.19 30.51 4.97
## 5171 5171 2.14 0.54 0.29 1.07 3.21 25.45 4.63
## 5201 5201 7.66 1.28 1.64 5.15 10.17 16.72 5.83
## 5202 5202 12.01 1.45 2.10 9.16 14.86 12.11 5.19
## 5203 5203 6.35 1.02 1.04 4.36 8.34 15.98 4.76
## 5204 5204 10.99 2.17 4.71 6.73 15.25 19.79 11.37
## 5205 5205 7.15 1.40 1.96 4.40 9.90 19.62 6.45
## 5206 5206 10.40 1.67 2.79 7.13 13.66 16.02 7.09
## 5207 5207 4.12 1.17 1.37 1.83 6.41 28.31 6.85
## 5208 5208 8.18 1.44 2.07 5.36 10.99 17.59 5.10
## 5271 5271 3.36 0.72 0.52 1.96 4.77 21.27 3.73
## 5272 5272 9.38 1.62 2.62 6.20 12.55 17.27 6.24
## 5301 5301 1.62 0.44 0.19 0.77 2.47 26.91 2.88
## 5302 5302 12.89 2.18 4.75 8.62 17.15 16.90 10.40
## 5303 5303 4.58 0.94 0.88 2.73 6.42 20.61 4.70
## 5304 5304 7.75 1.24 1.54 5.32 10.19 16.01 5.68
## 5305 5305 9.27 1.82 3.31 5.69 12.84 19.67 8.59
## 5306 5306 5.12 1.12 1.25 2.93 7.31 21.87 5.66
## 5307 5307 8.17 1.90 3.61 4.45 11.90 23.24 11.08
## 5308 5308 8.69 1.75 3.06 5.25 12.13 20.20 7.88
## 5309 5309 5.82 1.03 1.06 3.81 7.84 17.64 4.21
## 5310 5310 11.39 1.60 2.56 8.25 14.53 14.06 6.57
## 5311 5311 4.21 0.65 0.42 2.94 5.48 15.41 2.37
## 5312 5312 12.11 2.32 5.38 7.56 16.67 19.18 10.98
## 5313 5313 4.92 0.86 0.74 3.23 6.61 17.52 3.88
## 5314 5314 7.27 1.48 2.19 4.37 10.17 20.33 6.93
## 5315 5315 7.77 2.02 4.08 3.80 11.74 26.06 12.10
## 5316 5316 2.69 0.63 0.40 1.46 3.93 23.42 3.71
## 5317 5317 4.98 1.40 1.96 2.24 7.73 28.11 12.01
## 5318 5318 8.06 1.67 2.79 4.78 11.35 20.77 8.97
## 5319 5319 9.05 1.77 3.13 5.58 12.53 19.60 8.13
## 5320 5320 10.52 2.36 5.57 5.89 15.14 22.44 12.87
## 5321 5321 6.39 1.54 2.37 3.37 9.42 24.14 7.82
## 5371 5371 4.47 0.81 0.66 2.88 6.05 18.12 3.80
## 6101 6101 7.56 1.38 1.90 4.85 10.27 18.28 7.21
## 6102 6102 5.09 1.95 3.80 1.26 8.91 38.36 16.59
## 6103 6103 7.93 1.44 2.07 5.11 10.75 18.13 6.70
## 6104 6104 9.62 2.10 4.41 5.50 13.74 21.84 12.02
## 6105 6105 6.72 1.76 3.10 3.26 10.17 26.28 12.01
## 6106 6106 5.16 1.31 1.72 2.58 7.73 25.47 8.29
## 6107 6107 3.06 0.64 0.41 1.80 4.32 20.99 3.30
## 6108 6108 2.78 0.86 0.74 1.09 4.48 31.09 5.40
## 6109 6109 3.39 1.07 1.14 1.29 5.49 31.62 6.60
## 6110 6110 5.92 1.21 1.46 3.55 8.30 20.43 5.31
## 6111 6111 7.93 1.76 3.10 4.48 11.38 22.20 8.78
## 6112 6112 4.51 1.23 1.51 2.09 6.93 27.33 9.45
## 6171 6171 4.72 1.18 1.39 2.42 7.03 24.88 8.43
## 6172 6172 4.12 1.02 1.04 2.12 6.11 24.74 6.24
## 6201 6201 1.90 0.64 0.41 0.64 3.16 33.78 4.43
## 6202 6202 3.89 0.91 0.83 2.11 5.67 23.40 4.77
## 6203 6203 7.93 1.49 2.22 5.02 10.84 18.72 6.61
## 6204 6204 5.72 1.42 2.02 2.93 8.50 24.87 6.86
## 6205 6205 6.72 1.97 3.88 2.86 10.58 29.30 10.97
## 6206 6206 4.84 1.48 2.19 1.94 7.73 30.58 8.75
## 6207 6207 3.24 0.84 0.71 1.60 4.89 25.88 4.11
## 6208 6208 4.94 1.44 2.07 2.13 7.75 29.06 7.97
## 6209 6209 5.12 1.35 1.82 2.47 7.77 26.43 6.94
## 6210 6210 5.57 2.21 4.88 1.24 9.90 39.70 16.70
## 6211 6211 6.46 2.21 4.88 2.12 10.79 34.28 15.90
## 6212 6212 5.30 1.30 1.69 2.75 7.85 24.55 6.08
## 6213 6213 2.04 0.64 0.41 0.79 3.29 31.26 3.91
## 6271 6271 5.81 2.31 5.34 1.29 10.33 39.70 21.18
## 6301 6301 5.45 1.22 1.49 3.07 7.83 22.29 6.10
## 6302 6302 9.20 2.01 4.04 5.25 13.15 21.89 10.81
## 6303 6303 7.35 1.20 1.44 5.00 9.71 16.35 5.28
## 6304 6304 9.53 1.84 3.39 5.93 13.14 19.30 8.17
## 6305 6305 8.36 2.68 7.18 3.10 13.62 32.11 16.71
## 6306 6306 7.55 1.68 2.82 4.26 10.84 22.24 7.78
## 6307 6307 8.13 1.51 2.28 5.17 11.08 18.54 6.25
## 6308 6308 8.95 2.94 8.64 3.18 14.71 32.89 20.16
## 6309 6309 5.33 0.97 0.94 3.44 7.23 18.15 3.80
## 6310 6310 4.53 1.01 1.02 2.54 6.51 22.35 5.54
## 6311 6311 7.78 1.95 3.80 3.97 11.60 24.99 9.75
## 6371 6371 1.62 0.34 0.12 0.94 2.29 21.29 1.78
## 6372 6372 5.64 1.31 1.72 3.07 8.20 23.20 6.10
## 6401 6401 2.70 0.74 0.55 1.25 4.15 27.41 4.01
## 6402 6402 3.65 1.22 1.49 1.25 6.04 33.48 8.50
## 6403 6403 5.27 1.54 2.37 2.26 8.28 29.18 12.47
## 6404 6404 2.84 0.93 0.86 1.02 4.66 32.71 6.80
## 6405 6405 3.80 0.94 0.88 1.96 5.65 24.71 5.12
## 6409 6409 4.33 1.03 1.06 2.32 6.34 23.71 5.08
## 6411 6411 1.26 0.64 0.41 0.01 2.50 50.57 6.34
## 6471 6471 3.28 0.83 0.69 1.66 4.91 25.26 5.36
## 6472 6472 4.26 1.01 1.02 2.27 6.24 23.77 7.44
## 6474 6474 3.43 1.11 1.23 1.25 5.61 32.49 7.71
## 6501 6501 3.86 0.85 0.72 2.19 5.54 22.12 4.25
## 6502 6502 3.40 0.70 0.49 2.02 4.78 20.70 3.11
## 6503 6503 5.43 1.19 1.42 3.11 7.75 21.83 5.42
## 6504 6504 3.72 0.96 0.92 1.83 5.60 25.86 5.47
## 6571 6571 4.69 0.98 0.96 2.77 6.61 20.90 5.07
## 7101 7101 7.32 1.68 2.82 4.02 10.62 22.97 9.35
## 7102 7102 4.33 0.85 0.72 2.66 6.01 19.68 3.70
## 7103 7103 6.18 2.65 7.02 0.99 11.37 42.86 21.47
## 7104 7104 2.89 0.73 0.53 1.47 4.31 25.11 3.62
## 7105 7105 10.11 1.85 3.42 6.47 13.74 18.34 6.54
## 7106 7106 4.91 1.34 1.80 2.29 7.54 27.29 6.92
## 7107 7107 6.35 1.56 2.43 3.28 9.41 24.62 7.77
## 7108 7108 2.08 0.45 0.20 1.20 2.96 21.54 1.83
## 7109 7109 5.63 1.55 2.40 2.60 8.67 27.49 8.47
## 7110 7110 8.89 1.47 2.16 6.01 11.76 16.50 5.46
## 7111 7111 7.06 2.36 5.57 2.44 11.69 33.40 16.82
## 7171 7171 3.09 0.71 0.50 1.69 4.49 23.11 3.77
## 7172 7172 3.61 0.58 0.34 2.48 4.74 15.96 1.96
## 7173 7173 3.82 0.94 0.88 1.98 5.65 24.58 4.27
## 7174 7174 2.71 0.66 0.44 1.41 4.01 24.51 3.35
## 7201 7201 4.72 0.88 0.77 3.00 6.44 18.64 3.08
## 7202 7202 5.99 1.49 2.22 3.07 8.91 24.84 9.05
## 7203 7203 4.54 1.93 3.72 0.76 8.32 42.43 16.43
## 7204 7204 8.94 1.26 1.59 6.47 11.40 14.08 3.87
## 7205 7205 8.74 2.73 7.45 3.39 14.08 31.22 21.03
## 7206 7206 7.85 1.28 1.64 5.34 10.35 16.28 4.80
## 7207 7207 10.28 2.83 8.01 4.74 15.82 27.49 20.17
## 7208 7208 8.54 1.73 2.99 5.15 11.94 20.27 10.22
## 7209 7209 5.85 1.40 1.96 3.10 8.60 24.00 7.39
## 7210 7210 10.52 2.07 4.28 6.45 14.58 19.70 10.17
## 7211 7211 4.05 1.33 1.77 1.45 6.66 32.76 8.59
## 7212 7212 3.11 0.70 0.49 1.74 4.47 22.42 2.84
## 7271 7271 5.34 0.90 0.81 3.58 7.10 16.81 3.43
## 7301 7301 5.78 0.97 0.94 3.89 7.68 16.73 3.45
## 7302 7302 7.09 1.77 3.13 3.62 10.56 24.98 10.40
## 7303 7303 15.02 3.01 9.06 9.12 20.93 20.06 14.17
## 7304 7304 3.82 1.28 1.64 1.32 6.32 33.40 10.59
## 7305 7305 4.45 0.90 0.81 2.68 6.21 20.21 4.44
## 7306 7306 6.71 1.22 1.49 4.33 9.10 18.11 6.34
## 7307 7307 8.08 2.06 4.24 4.05 12.11 25.45 12.42
## 7308 7308 6.19 1.13 1.28 3.98 8.39 18.21 5.52
## 7309 7309 3.89 1.14 1.30 1.65 6.13 29.42 8.72
## 7310 7310 6.65 1.82 3.31 3.08 10.22 27.37 10.99
## 7311 7311 3.91 0.50 0.25 2.92 4.89 12.88 1.98
## 7312 7312 3.02 0.54 0.29 1.95 4.08 18.02 2.08
## 7313 7313 5.35 1.30 1.69 2.80 7.89 24.28 7.68
## 7314 7314 1.78 0.46 0.21 0.87 2.68 25.98 2.92
## 7315 7315 8.10 1.60 2.56 4.95 11.24 19.80 8.05
## 7316 7316 5.74 1.08 1.17 3.61 7.86 18.90 4.59
## 7317 7317 6.18 0.89 0.79 4.44 7.93 14.43 3.63
## 7318 7318 5.42 1.13 1.28 3.21 7.62 20.79 5.10
## 7322 7322 7.33 1.50 2.25 4.39 10.28 20.49 7.27
## 7325 7325 5.79 1.25 1.56 3.33 8.24 21.66 6.71
## 7326 7326 2.56 0.57 0.32 1.45 3.67 22.19 2.96
## 7371 7371 3.50 0.65 0.42 2.23 4.76 18.45 4.12
## 7372 7372 7.35 1.65 2.72 4.13 10.58 22.39 8.48
## 7373 7373 5.97 1.43 2.04 3.18 8.77 23.87 8.01
## 7401 7401 6.45 0.89 0.79 4.70 8.20 13.87 3.06
## 7402 7402 6.34 0.94 0.88 4.51 8.17 14.76 2.92
## 7403 7403 6.25 1.27 1.61 3.75 8.75 20.37 5.67
## 7404 7404 3.17 0.62 0.38 1.97 4.38 19.38 2.72
## 7405 7405 6.72 1.22 1.49 4.32 9.12 18.22 4.95
## 7406 7406 6.22 1.45 2.10 3.38 9.06 23.31 7.07
## 7407 7407 6.05 1.22 1.49 3.66 8.45 20.20 5.24
## 7408 7408 7.90 1.30 1.69 5.36 10.44 16.39 4.69
## 7409 7409 15.61 3.17 10.05 9.39 21.83 20.32 16.16
## 7410 7410 5.75 1.01 1.02 3.78 7.73 17.49 4.29
## 7411 7411 8.56 1.93 3.72 4.78 12.34 22.52 9.17
## 7412 7412 10.21 2.08 4.33 6.13 14.28 20.38 8.84
## 7413 7413 11.67 1.75 3.06 8.23 15.11 15.03 5.54
## 7414 7414 2.89 0.67 0.45 1.57 4.22 23.33 3.36
## 7415 7415 8.43 1.47 2.16 5.55 11.32 17.45 6.58
## 7471 7471 5.96 1.16 1.35 3.68 8.24 19.54 5.55
## 7472 7472 3.71 0.73 0.53 2.28 5.14 19.72 3.29
## 7501 7501 9.52 1.20 1.44 7.16 11.87 12.63 3.49
## 7502 7502 6.15 1.11 1.23 3.98 8.32 17.99 5.08
## 7503 7503 8.57 1.47 2.16 5.69 11.45 17.16 5.82
## 7504 7504 7.77 1.53 2.34 4.77 10.78 19.71 7.48
## 7505 7505 7.60 1.53 2.34 4.60 10.60 20.13 6.34
## 7571 7571 8.24 1.72 2.96 4.87 11.62 20.89 8.00
## 7601 7601 7.50 1.06 1.12 5.42 9.57 14.12 3.79
## 7602 7602 4.63 0.86 0.74 2.94 6.31 18.60 4.78
## 7603 7603 3.82 0.72 0.52 2.41 5.24 18.91 2.80
## 7604 7604 7.52 1.23 1.51 5.11 9.93 16.34 5.85
## 7605 7605 10.69 2.16 4.67 6.47 14.92 20.16 10.52
## 7606 7606 3.79 0.75 0.56 2.33 5.25 19.69 3.09
## 8101 8101 1.51 0.45 0.20 0.64 2.39 29.56 2.90
## 8102 8102 3.38 1.07 1.14 1.27 5.48 31.80 8.51
## 8103 8103 3.14 0.75 0.56 1.66 4.61 23.97 4.84
## 8104 8104 4.26 1.14 1.30 2.03 6.49 26.68 6.96
## 8105 8105 1.38 0.32 0.10 0.75 2.00 23.12 1.84
## 8106 8106 6.31 1.18 1.39 4.00 8.62 18.69 5.88
## 8107 8107 4.28 1.04 1.08 2.24 6.33 24.36 6.59
## 8108 8108 3.11 0.80 0.64 1.53 4.68 25.89 5.11
## 8109 8109 5.79 1.23 1.51 3.37 8.20 21.30 7.31
## 8171 8171 3.91 0.71 0.50 2.52 5.30 18.11 3.21
## 8172 8172 3.21 0.65 0.42 1.94 4.49 20.28 3.25
## 8201 8201 6.16 2.09 4.37 2.07 10.26 33.91 15.37
## 8202 8202 3.16 0.89 0.79 1.41 4.90 28.22 5.26
## 8203 8203 3.62 0.89 0.79 1.86 5.37 24.75 5.09
## 8204 8204 4.58 0.77 0.59 3.07 6.10 16.86 3.16
## 8205 8205 4.41 0.77 0.59 2.90 5.92 17.47 3.35
## 8206 8206 2.68 0.68 0.46 1.35 4.01 25.37 3.65
## 8207 8207 4.77 1.26 1.59 2.31 7.23 26.33 8.24
## 8208 8208 4.83 0.84 0.71 3.18 6.47 17.41 3.16
## 8271 8271 4.28 1.05 1.10 2.23 6.33 24.45 6.17
## 8272 8272 2.65 0.64 0.41 1.40 3.91 24.16 3.58
## 9101 9101 1.85 0.52 0.27 0.83 2.88 28.30 3.23
## 9102 9102 3.96 1.07 1.14 1.86 6.06 27.01 7.26
## 9103 9103 1.39 0.53 0.28 0.35 2.43 38.04 4.05
## 9104 9104 3.19 0.82 0.67 1.58 4.81 25.74 5.25
## 9105 9105 5.47 1.32 1.74 2.87 8.06 24.23 7.38
## 9111 9111 2.03 0.53 0.28 0.99 3.07 26.17 2.88
## 9112 9112 2.80 0.83 0.69 1.17 4.44 29.72 5.04
## 9201 9201 3.44 0.65 0.42 2.17 4.71 18.85 3.13
## 9202 9202 4.39 1.05 1.10 2.33 6.45 23.93 5.79
## 9203 9203 2.39 0.77 0.59 0.88 3.90 32.14 7.84
## 9204 9204 1.78 0.40 0.16 0.99 2.57 22.60 1.42
## 9205 9205 1.94 0.67 0.45 0.63 3.26 34.55 3.48
## 9271 9271 4.86 1.40 1.96 2.11 7.61 28.86 10.75
## 9403 9403 3.65 0.83 0.69 2.01 5.28 22.90 4.26
## 9408 9408 1.52 0.57 0.32 0.41 2.63 37.42 4.74
## 9409 9409 4.92 1.50 2.25 1.97 7.86 30.54 10.95
## 9419 9419 1.26 0.43 0.18 0.42 2.10 34.05 2.72
## 9420 9420 3.23 0.83 0.69 1.61 4.86 25.65 4.05
## 9426 9426 4.57 0.90 0.81 2.80 6.33 19.70 4.02
## 9427 9427 2.15 0.81 0.66 0.56 3.74 37.79 7.47
## 9428 9428 1.00 0.33 0.11 0.36 1.64 32.84 1.08
## 9471 9471 2.57 0.71 0.50 1.18 3.97 27.69 4.32
## 9501 9501 5.37 0.87 0.76 3.66 7.07 16.18 3.22
## 9502 9502 4.98 2.74 7.51 0.00 10.35 54.94 33.22
## 9503 9503 4.74 1.10 1.21 2.58 6.90 23.23 6.46
## 9504 9504 1.31 0.49 0.24 0.35 2.28 37.58 4.64
## 9601 9601 0.93 0.33 0.11 0.28 1.58 35.64 2.33
## 9602 9602 0.99 0.38 0.14 0.25 1.74 38.35 2.71
## 9603 9603 0.44 0.32 0.10 0.00 1.06 72.56 3.99
## 9604 9604 3.16 0.81 0.66 1.58 4.74 25.56 4.39
## 9605 9605 0.35 0.27 0.07 0.00 0.87 75.36 3.26
## 9606 9606 1.12 0.73 0.53 0.00 2.55 65.61 9.02
## 9607 9607 0.14 0.14 0.02 0.00 0.40 98.99 2.36
## 9608 9608 3.55 2.24 5.02 0.00 7.94 63.01 21.52
## 9701 9701 2.20 0.57 0.32 1.09 3.31 25.79 2.52
## 9702 9702 2.89 0.87 0.76 1.19 4.60 30.11 4.71
## 9703 9703 0.60 0.42 0.18 0.00 1.42 69.64 4.06
## 9704 9704 1.73 0.77 0.59 0.22 3.24 44.43 6.86
## 9705 9705 2.99 0.55 0.30 1.92 4.07 18.26 1.98
## 9706 9706 2.02 0.65 0.42 0.75 3.29 32.18 2.72
## 9707 9707 2.58 0.98 0.96 0.67 4.50 37.82 8.42
## 9708 9708 5.12 1.30 1.69 2.57 7.67 25.37 7.09Visualisasi
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_UNMET_PROV_2024.csv")
write.csv(resultK,"RSE_DIRECT_UNMET_KAKO_2024.csv")
rm(list=ls()) # Clear environmentDirektorat Statistik Kesejahteraan Rakyat, BPS, saptahas@bps.go.id