Direct Estimation
Improved Drinking
Water 2024
using Finite Population Correction (FPC)
Sapta Hastho Ponco
2025-05-05
Load Packages
pacman::p_load(readxl, haven, dplyr, tidyr, survey, janitor, magrittr, ggplot2)Load Dataset
datakor <- read_sav("KOR24GAB_0208.sav")Indikator air minum layak ditujukan untuk pertanyaan-pertanyaan terkait rumah tangga. Oleh karena itu dalam penghitungan perlu dilakukan penyeleksian unit analisis.
ruta <- filter(datakor , R403==1)1 Manipulasi working variable
Untuk meringankan beban kerja, perlu dilakukan filter variabel yang digunakan.
dataku <- ruta %>%
select(R101, R102, WI1, WI2, PSU, SSU, STRATA, R105, R403, R407, R1801:R1817, R301, FWT, fpc1, fpc2)2 Calculating indicator
Sumber air minum layak adalah jika rumah tangga menggunakan sumber air minum utama berupa ledeng, perpipaan, perpipaan eceran, kran halaman, hidran umum, air terlindungi, dan penampungan air hujan. Air terlindungi mencakup sumur bor/ pompa, sumur terlindungi, dan mata air terlindung. Bagi rumah tangga yang menggunakan sumber air minum berupa air kemasan atau air isi ulang, maka rumah tangga dikategorikan memiliki akses air minum layak sumber air untuk mandi/cuci berasal dari ledeng, sumur bor/pompa, sumur terlindung, mata air terlindung, dan air hujan (Metadata SDGs 2024).
indicator <- dataku %>%
mutate(kabu = R101*100 + R102) %>%
rename(prov="R101", kab="R102") %>%
mutate(psu=prov*10000000 + kab*100000 + WI1,
strata= prov*10000000 + kab*100000 + R105,
denom=100,
strata=strata,
indic1_layak = dplyr::case_when(dataku$R1809A<=3 & dataku$R1809B==1 & dataku$R1809C<=2~100,
(dataku$R105 ==2 & dataku$R1809A<=3 & dataku$R1809B==1 & dataku$R1809C==4~100)),
indic1_layak=case_when(indic1_layak==100~100,is.na(indic1_layak)~0)) %>%
mutate(drinkw=car::recode(dataku$R1810A,"3:5=100;7=100;10=100;else=0")) %>%
mutate(bottle=car::recode(dataku$R1810A,"1:2=100;else=0")) %>%
mutate(bathing=car::recode(dataku$R1814A,"3:5=100;7=100;10=100;else=0")) %>%
mutate(indic2_layak = dplyr::case_when(drinkw==100 | (bottle==100 & bathing==100)~100)) %>%
mutate(indic2_layak = ifelse(is.na(indic2_layak), 0, indic2_layak)) %>%
mutate(luas_cap=R1804/R301) %>%
mutate(indic3_layak = dplyr::case_when(luas_cap>=7.20~100,luas_cap<7.20~0)) %>%
mutate(roof=car::recode(dataku$R1806,"1:3=100;6=100;else=0")) %>%
mutate(floor=car::recode(dataku$R1808,"1:6=100;else=0")) %>%
mutate(wall=car::recode(dataku$R1807,"1:3=100;5=100;else=0")) %>%
mutate(indic4_layak = dplyr::case_when((roof==100 & floor==100 & wall==100)~100)) %>%
mutate(indic4_layak = ifelse(is.na(indic4_layak), 0, indic4_layak))
indicator <- indicator %>%
mutate(drink = dplyr::case_when((indic2_layak==100)~100)) %>%
mutate(drink = ifelse(is.na(drink), 0, drink))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.0037058 0.0084963 0.0337114 0.0236375 0.8644758Set Output Statistics
options(survey.adjust.domain.lonely=TRUE)
options(survey.lonely.psu="adjust")
#Provinsi
resultP <- svyby(
formula = ~drink,
denom = ~denom,
by = ~prov,
design = susenas.design,
deff = TRUE,
svyratio,
vartype = c("se", "ci", "ci", "cv", "cvpct", "var"))
resultP[is.na(resultP)] <- 0
resultP <- resultP %>%
mutate(
theta = round(`drink/denom` * 100, 2),
SE = round(`se.drink/denom` * 100, 2),
VAR = round(SE^2, 2),
CI_LOWER = round(`ci_l` * 100, 2),
CI_UPPER = round(`ci_u` * 100, 2),
RSE = round(`cv%`, 2),
DEFF = round(DEff, 2))
resultP$CI_LOWER[resultP$CI_LOWER < 0] <- 0options(survey.adjust.domain.lonely=TRUE)
options(survey.lonely.psu="adjust")
#Kabupaten/kota
resultK <- svyby(
formula = ~drink,
denom = ~denom,
by = ~kabu,
design = susenas.design,
deff = TRUE,
svyratio,
vartype = c("se", "ci", "ci", "cv", "cvpct", "var"))
resultK[is.na(resultK)] <- 0
resultK <- resultK %>%
mutate(
theta = round(`drink/denom` * 100, 2),
SE = round(`se.drink/denom` * 100, 2),
VAR = round(SE^2, 2),
CI_LOWER = round(`ci_l` * 100, 2),
CI_UPPER = round(`ci_u` * 100, 2),
RSE = round(`cv%`, 2),
DEFF = round(DEff, 2))
resultK$CI_LOWER[resultK$CI_LOWER < 0] <- 03.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):
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.
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 90.08 0.62 0.38 88.87 91.30 0.69 5.85
## 12 12 92.94 0.38 0.14 92.20 93.68 0.41 4.65
## 13 13 86.85 0.78 0.61 85.32 88.38 0.90 6.28
## 14 14 92.24 0.65 0.42 90.97 93.51 0.70 4.93
## 15 15 82.16 1.00 1.00 80.19 84.13 1.22 4.82
## 16 16 87.23 0.85 0.72 85.57 88.89 0.97 7.19
## 17 17 72.10 1.32 1.74 69.51 74.68 1.83 5.05
## 18 18 84.51 0.89 0.79 82.77 86.25 1.05 6.33
## 19 19 82.97 1.33 1.77 80.37 85.57 1.60 5.03
## 21 21 94.71 0.53 0.28 93.68 95.74 0.55 2.35
## 31 31 99.96 0.03 0.00 99.91 100.01 0.03 0.89
## 32 32 94.90 0.32 0.10 94.26 95.53 0.34 5.58
## 33 33 95.43 0.28 0.08 94.87 95.98 0.30 5.47
## 34 34 96.91 0.49 0.24 95.96 97.86 0.50 3.20
## 35 35 96.93 0.21 0.04 96.50 97.35 0.22 5.04
## 36 36 94.01 0.53 0.28 92.98 95.04 0.56 3.59
## 51 51 98.32 0.31 0.10 97.73 98.92 0.31 3.59
## 52 52 96.13 0.53 0.28 95.09 97.18 0.56 5.37
## 53 53 88.55 0.70 0.49 87.19 89.92 0.79 6.07
## 61 61 82.13 0.98 0.96 80.21 84.05 1.19 5.86
## 62 62 78.71 1.20 1.44 76.37 81.06 1.52 6.92
## 63 63 77.34 1.10 1.21 75.19 79.50 1.42 5.78
## 64 64 89.11 1.03 1.06 87.10 91.13 1.15 6.68
## 65 65 89.42 1.15 1.32 87.16 91.67 1.28 3.84
## 71 71 94.46 0.48 0.23 93.52 95.40 0.51 3.80
## 72 72 87.86 0.92 0.85 86.06 89.65 1.04 5.97
## 73 73 92.23 0.48 0.23 91.29 93.17 0.52 4.92
## 74 74 95.31 0.51 0.26 94.30 96.32 0.54 5.53
## 75 75 96.13 0.65 0.42 94.86 97.39 0.67 3.90
## 76 76 80.14 1.50 2.25 77.19 83.08 1.87 5.08
## 81 81 93.32 0.74 0.55 91.87 94.77 0.79 5.43
## 82 82 88.92 0.96 0.92 87.03 90.80 1.08 5.26
## 91 91 82.38 1.35 1.82 79.74 85.02 1.63 4.61
## 92 92 81.68 1.67 2.79 78.41 84.95 2.04 5.66
## 94 94 85.84 1.10 1.21 83.68 88.00 1.28 4.55
## 95 95 71.90 2.09 4.37 67.80 76.00 2.91 4.71
## 96 96 81.53 1.48 2.19 78.63 84.42 1.81 6.23
## 97 97 30.64 1.90 3.61 26.91 34.38 6.21 6.79Visualisasi
ggplot(data = resultP, aes(x = prov)) +
geom_line(aes(y = RSE, color = "RSE"), size = 1) + # Garis untuk RSE
geom_line(aes(y = theta, color = "Theta"), size = 1) + # Garis untuk Theta
geom_hline(yintercept = 25, color = "red", linetype = "dashed") + # Garis horizontal pada y = 25
labs(
title = "Provinsi vs Nilai RSE dan Theta - 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 82.79 2.72 7.40 77.47 88.11 3.28 2.91
## 1102 1102 79.60 3.36 11.29 73.02 86.19 4.22 3.88
## 1103 1103 85.87 3.62 13.10 78.77 92.96 4.22 6.32
## 1104 1104 91.40 3.10 9.61 85.32 97.48 3.39 7.13
## 1105 1105 84.43 3.27 10.69 78.02 90.84 3.87 5.39
## 1106 1106 88.27 3.59 12.89 81.23 95.31 4.07 7.16
## 1107 1107 96.52 1.53 2.34 93.51 99.53 1.59 4.03
## 1108 1108 90.68 2.74 7.51 85.30 96.05 3.02 5.87
## 1109 1109 87.57 3.07 9.42 81.55 93.58 3.51 6.09
## 1110 1110 90.06 2.37 5.62 85.41 94.70 2.63 4.14
## 1111 1111 94.56 1.35 1.82 91.91 97.20 1.43 2.48
## 1112 1112 95.12 1.35 1.82 92.47 97.76 1.42 2.16
## 1113 1113 80.93 3.10 9.61 74.86 87.01 3.83 3.49
## 1114 1114 89.02 2.59 6.71 83.95 94.09 2.91 4.26
## 1115 1115 90.75 2.93 8.58 85.01 96.49 3.23 5.68
## 1116 1116 87.15 3.48 12.11 80.32 93.98 4.00 5.96
## 1117 1117 84.98 3.84 14.75 77.44 92.51 4.52 6.43
## 1118 1118 94.54 2.11 4.45 90.40 98.67 2.23 4.77
## 1171 1171 99.05 0.62 0.38 97.84 100.26 0.62 2.52
## 1172 1172 89.92 2.88 8.29 84.28 95.56 3.20 4.80
## 1173 1173 98.06 0.93 0.86 96.24 99.89 0.95 2.69
## 1174 1174 97.41 1.09 1.19 95.26 99.55 1.12 2.71
## 1175 1175 57.25 4.51 20.34 48.40 66.10 7.89 4.71
## 1201 1201 48.88 5.24 27.46 38.60 59.16 10.73 6.05
## 1202 1202 84.72 3.04 9.24 78.77 90.67 3.58 5.03
## 1203 1203 74.84 4.85 23.52 65.33 84.35 6.48 7.68
## 1204 1204 82.98 2.41 5.81 78.26 87.70 2.90 2.53
## 1205 1205 92.91 2.47 6.10 88.06 97.75 2.66 6.16
## 1206 1206 96.01 1.60 2.56 92.87 99.14 1.67 3.68
## 1207 1207 93.42 2.23 4.97 89.05 97.79 2.39 5.33
## 1208 1208 94.60 1.62 2.62 91.42 97.78 1.72 4.04
## 1209 1209 99.09 0.50 0.25 98.10 100.08 0.51 2.34
## 1210 1210 97.62 0.89 0.79 95.87 99.38 0.92 2.15
## 1211 1211 96.61 1.57 2.46 93.53 99.69 1.63 5.26
## 1212 1212 96.37 1.20 1.44 94.02 98.72 1.24 4.12
## 1213 1213 92.49 1.92 3.69 88.73 96.25 2.07 4.40
## 1214 1214 69.38 5.12 26.21 59.34 79.41 7.38 7.58
## 1215 1215 97.73 0.86 0.74 96.05 99.42 0.88 1.84
## 1216 1216 68.97 3.38 11.42 62.36 75.59 4.89 2.82
## 1217 1217 80.38 3.37 11.36 73.78 86.97 4.19 4.01
## 1218 1218 98.43 0.63 0.40 97.21 99.66 0.64 1.91
## 1219 1219 97.11 1.75 3.06 93.67 100.55 1.81 7.29
## 1220 1220 82.23 4.32 18.66 73.76 90.70 5.26 7.85
## 1221 1221 89.96 2.18 4.75 85.69 94.23 2.42 3.06
## 1222 1222 90.05 2.58 6.66 84.99 95.11 2.87 4.91
## 1223 1223 91.50 2.52 6.35 86.57 96.43 2.75 5.06
## 1224 1224 71.72 4.60 21.16 62.72 80.73 6.41 5.83
## 1225 1225 70.95 3.92 15.37 63.27 78.63 5.52 4.23
## 1271 1271 96.86 1.92 3.69 93.09 100.63 1.99 6.87
## 1272 1272 91.27 2.60 6.76 86.17 96.37 2.85 4.73
## 1273 1273 99.74 0.15 0.02 99.45 100.03 0.15 0.51
## 1274 1274 99.53 0.27 0.07 99.00 100.06 0.27 0.86
## 1275 1275 99.84 0.15 0.02 99.55 100.12 0.15 1.37
## 1276 1276 98.88 0.61 0.37 97.68 100.08 0.62 1.93
## 1277 1277 64.99 4.71 22.18 55.76 74.21 7.24 5.60
## 1278 1278 80.17 3.43 11.76 73.45 86.88 4.27 4.13
## 1301 1301 70.23 6.01 36.12 58.45 82.02 8.56 9.43
## 1302 1302 90.00 1.94 3.76 86.19 93.80 2.16 2.97
## 1303 1303 81.84 3.13 9.80 75.71 87.97 3.82 4.37
## 1304 1304 72.02 3.84 14.75 64.49 79.54 5.33 4.26
## 1305 1305 92.74 2.25 5.06 88.33 97.14 2.42 4.92
## 1306 1306 91.77 1.69 2.86 88.47 95.07 1.84 2.48
## 1307 1307 89.68 1.99 3.96 85.77 93.58 2.22 3.20
## 1308 1308 68.64 7.08 50.13 54.76 82.52 10.32 15.35
## 1309 1309 86.06 2.48 6.15 81.20 90.92 2.88 3.18
## 1310 1310 79.09 3.31 10.96 72.61 85.57 4.18 3.67
## 1311 1311 70.00 3.70 13.69 62.75 77.26 5.29 4.03
## 1312 1312 83.07 3.01 9.06 77.17 88.96 3.62 4.26
## 1371 1371 97.65 0.70 0.49 96.27 99.02 0.72 1.75
## 1372 1372 97.04 1.33 1.77 94.44 99.64 1.37 3.39
## 1373 1373 89.62 2.27 5.15 85.18 94.06 2.53 3.04
## 1374 1374 98.20 0.73 0.53 96.77 99.63 0.74 1.61
## 1375 1375 99.45 0.38 0.14 98.70 100.20 0.38 1.48
## 1376 1376 99.44 0.41 0.17 98.63 100.24 0.41 1.67
## 1377 1377 95.62 1.28 1.64 93.12 98.12 1.33 2.18
## 1401 1401 86.73 3.04 9.24 80.77 92.69 3.51 5.01
## 1402 1402 77.63 3.33 11.09 71.10 84.17 4.29 4.22
## 1403 1403 94.83 1.51 2.28 91.86 97.80 1.60 3.64
## 1404 1404 83.26 4.03 16.24 75.36 91.15 4.84 7.68
## 1405 1405 91.79 1.96 3.84 87.95 95.64 2.14 3.36
## 1406 1406 96.00 1.59 2.53 92.88 99.12 1.66 5.20
## 1407 1407 93.45 1.96 3.84 89.61 97.29 2.10 4.63
## 1408 1408 95.72 1.41 1.99 92.95 98.49 1.48 3.52
## 1409 1409 86.82 3.03 9.18 80.88 92.76 3.49 5.91
## 1410 1410 93.44 2.24 5.02 89.05 97.83 2.40 4.71
## 1471 1471 99.65 0.33 0.11 99.00 100.29 0.33 2.47
## 1473 1473 94.82 1.56 2.43 91.75 97.88 1.65 3.07
## 1501 1501 76.34 4.94 24.40 66.65 86.03 6.48 8.30
## 1502 1502 60.74 4.65 21.62 51.63 69.85 7.65 6.00
## 1503 1503 76.41 3.35 11.22 69.85 82.97 4.38 3.87
## 1504 1504 84.80 3.13 9.80 78.67 90.94 3.69 4.74
## 1505 1505 78.46 3.01 9.06 72.56 84.36 3.84 3.54
## 1506 1506 87.12 2.55 6.50 82.12 92.13 2.93 3.62
## 1507 1507 93.84 2.16 4.67 89.61 98.06 2.30 4.97
## 1508 1508 80.71 3.70 13.69 73.46 87.97 4.59 5.84
## 1509 1509 81.93 3.00 9.00 76.06 87.81 3.66 4.03
## 1571 1571 96.06 1.10 1.21 93.91 98.21 1.14 2.24
## 1572 1572 89.50 3.46 11.97 82.73 96.27 3.86 6.97
## 1601 1601 82.55 2.97 8.82 76.72 88.38 3.60 4.07
## 1602 1602 78.41 3.70 13.69 71.17 85.65 4.71 6.31
## 1603 1603 91.18 2.08 4.33 87.11 95.26 2.28 4.21
## 1604 1604 84.86 2.83 8.01 79.31 90.42 3.34 4.39
## 1605 1605 91.39 2.55 6.50 86.40 96.39 2.79 4.71
## 1606 1606 86.11 2.93 8.58 80.36 91.85 3.40 5.11
## 1607 1607 81.42 3.31 10.96 74.93 87.91 4.07 5.66
## 1608 1608 82.60 3.47 12.04 75.81 89.40 4.20 4.79
## 1609 1609 96.19 1.16 1.35 93.92 98.47 1.21 2.62
## 1610 1610 73.69 4.75 22.56 64.38 83.00 6.44 7.71
## 1611 1611 69.33 5.03 25.30 59.48 79.18 7.25 6.57
## 1612 1612 94.09 2.28 5.20 89.62 98.57 2.43 5.12
## 1613 1613 88.41 2.53 6.40 83.44 93.37 2.87 3.41
## 1671 1671 95.00 2.68 7.18 89.75 100.26 2.82 12.94
## 1672 1672 94.94 1.62 2.62 91.76 98.12 1.71 3.15
## 1673 1673 89.81 3.04 9.24 83.85 95.77 3.39 5.64
## 1674 1674 97.18 0.89 0.79 95.43 98.93 0.92 1.66
## 1701 1701 68.26 4.95 24.50 58.56 77.96 7.25 6.29
## 1702 1702 55.60 4.79 22.94 46.21 64.98 8.61 5.79
## 1703 1703 80.40 3.02 9.12 74.48 86.31 3.76 3.59
## 1704 1704 79.79 2.65 7.02 74.60 84.98 3.32 2.42
## 1705 1705 70.46 4.51 20.34 61.62 79.29 6.40 5.60
## 1706 1706 75.82 3.94 15.52 68.10 83.54 5.20 4.93
## 1707 1707 57.41 4.31 18.58 48.95 65.86 7.51 4.25
## 1708 1708 61.20 5.68 32.26 50.07 72.34 9.28 7.56
## 1709 1709 50.44 4.96 24.60 40.73 60.16 9.82 5.39
## 1771 1771 92.21 2.14 4.58 88.02 96.41 2.32 4.20
## 1801 1801 64.37 5.35 28.62 53.89 74.86 8.31 7.14
## 1802 1802 68.48 4.06 16.48 60.52 76.44 5.93 5.75
## 1803 1803 95.52 1.10 1.21 93.36 97.67 1.15 2.37
## 1804 1804 87.01 2.96 8.76 81.20 92.82 3.41 6.77
## 1805 1805 84.59 3.31 10.96 78.11 91.07 3.91 7.31
## 1806 1806 69.31 4.38 19.18 60.73 77.89 6.32 6.68
## 1807 1807 71.40 3.78 14.29 64.00 78.80 5.29 4.91
## 1808 1808 88.55 2.92 8.53 82.82 94.27 3.30 5.98
## 1809 1809 89.80 2.51 6.30 84.88 94.72 2.80 4.56
## 1810 1810 93.33 1.66 2.76 90.08 96.59 1.78 3.12
## 1811 1811 81.81 3.88 15.05 74.22 89.41 4.74 5.80
## 1812 1812 63.16 6.38 40.70 50.65 75.68 10.11 10.75
## 1813 1813 73.77 3.61 13.03 66.69 80.85 4.90 3.61
## 1871 1871 99.09 0.51 0.26 98.10 100.09 0.51 2.27
## 1872 1872 96.25 1.72 2.96 92.88 99.63 1.79 4.65
## 1901 1901 77.35 3.97 15.76 69.56 85.14 5.14 5.62
## 1902 1902 83.36 2.96 8.76 77.56 89.17 3.55 3.62
## 1903 1903 84.94 3.43 11.76 78.22 91.65 4.03 5.27
## 1904 1904 88.25 2.30 5.29 83.75 92.76 2.60 2.93
## 1905 1905 78.30 3.71 13.76 71.04 85.56 4.73 4.56
## 1906 1906 70.94 3.88 15.05 63.33 78.54 5.47 4.04
## 1971 1971 96.47 1.15 1.32 94.21 98.73 1.19 2.24
## 2101 2101 85.89 3.39 11.49 79.24 92.54 3.95 5.94
## 2102 2102 84.57 3.15 9.92 78.40 90.75 3.72 4.44
## 2103 2103 91.53 2.34 5.48 86.94 96.12 2.56 3.97
## 2104 2104 85.48 3.55 12.60 78.52 92.43 4.15 5.66
## 2105 2105 24.65 4.61 21.25 15.61 33.68 18.71 6.23
## 2171 2171 99.40 0.38 0.14 98.65 100.14 0.38 2.06
## 2172 2172 91.96 1.93 3.72 88.18 95.74 2.10 2.85
## 3101 3101 92.93 0.76 0.58 91.43 94.43 0.82 0.49
## 3171 3171 100.00 0.00 0.00 100.00 100.00 0.00 0.00
## 3172 3172 100.00 0.00 0.00 100.00 100.00 0.00 0.00
## 3173 3173 100.00 0.00 0.00 100.00 100.00 0.00 0.00
## 3174 3174 100.00 0.00 0.00 100.00 100.00 0.00 0.00
## 3175 3175 99.85 0.15 0.02 99.57 100.14 0.15 1.52
## 3201 3201 93.93 1.56 2.43 90.87 96.99 1.66 5.57
## 3202 3202 89.30 1.87 3.50 85.64 92.97 2.09 3.96
## 3203 3203 91.62 1.63 2.66 88.41 94.82 1.78 3.77
## 3204 3204 96.28 1.14 1.30 94.04 98.52 1.19 4.28
## 3205 3205 85.37 2.19 4.80 81.09 89.65 2.56 4.32
## 3206 3206 86.78 2.34 5.48 82.20 91.37 2.70 4.77
## 3207 3207 96.42 0.90 0.81 94.67 98.18 0.93 2.42
## 3208 3208 97.92 0.98 0.96 96.00 99.85 1.00 3.98
## 3209 3209 97.12 0.87 0.76 95.42 98.83 0.89 2.91
## 3210 3210 94.38 1.67 2.79 91.12 97.65 1.77 4.77
## 3211 3211 96.23 1.26 1.59 93.76 98.70 1.31 3.95
## 3212 3212 97.39 0.92 0.85 95.59 99.20 0.94 3.47
## 3213 3213 98.69 0.44 0.19 97.81 99.56 0.45 1.51
## 3214 3214 91.02 2.26 5.11 86.59 95.44 2.48 5.19
## 3215 3215 97.23 0.72 0.52 95.82 98.63 0.74 2.15
## 3216 3216 96.76 1.09 1.19 94.63 98.89 1.12 4.40
## 3217 3217 89.25 2.34 5.48 84.67 93.83 2.62 5.92
## 3218 3218 88.36 2.64 6.97 83.20 93.53 2.98 5.02
## 3271 3271 99.11 0.35 0.12 98.43 99.79 0.35 1.13
## 3272 3272 99.46 0.35 0.12 98.78 100.14 0.35 1.48
## 3273 3273 97.62 1.04 1.08 95.58 99.67 1.07 5.32
## 3274 3274 98.23 0.81 0.66 96.64 99.82 0.83 2.50
## 3275 3275 98.89 0.84 0.71 97.24 100.54 0.85 7.34
## 3276 3276 99.85 0.11 0.01 99.63 100.07 0.11 0.85
## 3277 3277 99.60 0.30 0.09 99.02 100.19 0.30 1.79
## 3278 3278 93.79 2.26 5.11 89.36 98.23 2.41 6.52
## 3279 3279 96.97 1.18 1.39 94.67 99.28 1.21 2.72
## 3301 3301 93.33 1.46 2.13 90.47 96.19 1.56 3.69
## 3302 3302 94.33 1.53 2.34 91.34 97.33 1.62 4.55
## 3303 3303 88.24 2.00 4.00 84.33 92.15 2.26 3.23
## 3304 3304 91.75 1.99 3.96 87.85 95.65 2.17 4.38
## 3305 3305 87.29 2.47 6.10 82.45 92.12 2.83 5.00
## 3306 3306 90.04 2.54 6.45 85.06 95.02 2.82 5.70
## 3307 3307 95.24 1.46 2.13 92.37 98.11 1.54 4.11
## 3308 3308 93.27 1.85 3.42 89.65 96.90 1.98 4.97
## 3309 3309 94.83 1.95 3.80 91.01 98.64 2.05 6.71
## 3310 3310 99.48 0.35 0.12 98.79 100.17 0.35 2.19
## 3311 3311 99.15 0.41 0.17 98.36 99.95 0.41 1.62
## 3312 3312 96.96 1.19 1.42 94.62 99.30 1.23 4.22
## 3313 3313 99.64 0.21 0.04 99.24 100.05 0.21 1.02
## 3314 3314 98.45 0.60 0.36 97.27 99.62 0.61 2.05
## 3315 3315 87.09 2.29 5.24 82.61 91.58 2.63 4.84
## 3316 3316 97.02 1.15 1.32 94.78 99.27 1.18 3.83
## 3317 3317 95.86 1.67 2.79 92.57 99.14 1.75 5.59
## 3318 3318 98.32 1.06 1.12 96.23 100.40 1.08 6.23
## 3319 3319 99.60 0.28 0.08 99.05 100.14 0.28 1.60
## 3320 3320 89.46 2.69 7.24 84.18 94.74 3.01 7.02
## 3321 3321 95.99 1.91 3.65 92.25 99.72 1.98 8.56
## 3322 3322 99.60 0.20 0.04 99.22 99.99 0.20 0.85
## 3323 3323 97.83 1.00 1.00 95.86 99.80 1.03 3.97
## 3324 3324 97.65 1.36 1.85 94.99 100.32 1.39 7.07
## 3325 3325 96.46 1.40 1.96 93.71 99.21 1.45 4.78
## 3326 3326 96.91 1.39 1.93 94.17 99.64 1.44 5.44
## 3327 3327 94.38 2.08 4.33 90.32 98.45 2.20 7.41
## 3328 3328 98.28 0.92 0.85 96.48 100.07 0.93 4.76
## 3329 3329 94.66 1.41 1.99 91.88 97.43 1.49 4.12
## 3371 3371 99.93 0.07 0.00 99.78 100.07 0.07 0.40
## 3372 3372 99.72 0.28 0.08 99.18 100.27 0.28 2.07
## 3373 3373 98.98 0.67 0.45 97.67 100.28 0.67 2.52
## 3374 3374 99.87 0.13 0.02 99.63 100.12 0.13 1.29
## 3375 3375 99.27 0.37 0.14 98.55 99.99 0.37 1.23
## 3376 3376 99.66 0.34 0.12 99.01 100.32 0.34 2.08
## 3401 3401 92.93 1.68 2.82 89.63 96.22 1.81 3.07
## 3402 3402 97.96 0.76 0.58 96.47 99.44 0.77 2.47
## 3403 3403 98.47 0.47 0.22 97.54 99.40 0.48 1.23
## 3404 3404 95.37 1.19 1.42 93.05 97.70 1.25 3.03
## 3471 3471 99.94 0.06 0.00 99.83 100.06 0.06 0.43
## 3501 3501 84.35 2.93 8.58 78.60 90.09 3.47 4.89
## 3502 3502 95.57 1.58 2.50 92.48 98.65 1.65 4.86
## 3503 3503 81.47 2.96 8.76 75.67 87.26 3.63 4.82
## 3504 3504 98.62 0.77 0.59 97.11 100.13 0.78 3.75
## 3505 3505 98.66 0.79 0.62 97.12 100.20 0.80 4.26
## 3506 3506 96.72 1.45 2.10 93.88 99.57 1.50 6.94
## 3507 3507 97.94 0.85 0.72 96.27 99.61 0.87 4.23
## 3508 3508 97.31 0.87 0.76 95.59 99.02 0.90 2.54
## 3509 3509 95.30 1.02 1.04 93.30 97.30 1.07 2.82
## 3510 3510 97.06 0.95 0.90 95.20 98.93 0.98 3.47
## 3511 3511 96.54 1.38 1.90 93.83 99.25 1.43 4.77
## 3512 3512 95.17 2.14 4.58 90.97 99.36 2.25 8.58
## 3513 3513 96.90 1.24 1.54 94.47 99.33 1.28 4.66
## 3514 3514 98.62 0.60 0.36 97.44 99.80 0.61 2.74
## 3515 3515 98.07 0.85 0.72 96.41 99.73 0.86 4.11
## 3516 3516 97.46 1.45 2.10 94.62 100.31 1.49 7.42
## 3517 3517 98.76 0.64 0.41 97.52 100.01 0.64 3.01
## 3518 3518 99.71 0.23 0.05 99.25 100.16 0.23 1.55
## 3519 3519 97.55 1.10 1.21 95.40 99.70 1.12 4.17
## 3520 3520 98.56 0.66 0.44 97.27 99.85 0.67 2.38
## 3521 3521 97.07 1.13 1.28 94.86 99.28 1.16 3.85
## 3522 3522 96.13 1.68 2.82 92.84 99.42 1.75 7.30
## 3523 3523 96.96 1.39 1.93 94.24 99.68 1.43 5.96
## 3524 3524 84.60 3.45 11.90 77.84 91.37 4.08 8.34
## 3525 3525 97.65 0.92 0.85 95.85 99.45 0.94 3.34
## 3526 3526 98.00 0.90 0.81 96.23 99.76 0.92 3.48
## 3527 3527 97.66 1.04 1.08 95.62 99.70 1.06 3.93
## 3528 3528 96.98 0.92 0.85 95.19 98.77 0.94 2.40
## 3529 3529 99.41 0.23 0.05 98.96 99.85 0.23 0.80
## 3571 3571 99.73 0.20 0.04 99.34 100.12 0.20 0.99
## 3572 3572 99.60 0.34 0.12 98.94 100.27 0.34 1.63
## 3573 3573 99.48 0.34 0.12 98.82 100.14 0.34 1.82
## 3574 3574 100.00 0.00 0.00 100.00 100.00 0.00 0.00
## 3575 3575 99.94 0.06 0.00 99.83 100.06 0.06 0.33
## 3576 3576 99.77 0.15 0.02 99.47 100.07 0.15 0.58
## 3577 3577 100.00 0.00 0.00 100.00 100.00 0.00 0.00
## 3578 3578 100.00 0.00 0.00 100.00 100.00 0.00 0.00
## 3579 3579 98.27 1.06 1.12 96.19 100.36 1.08 3.89
## 3601 3601 78.51 2.78 7.73 73.06 83.97 3.55 4.18
## 3602 3602 77.92 3.31 10.96 71.43 84.41 4.25 5.79
## 3603 3603 97.79 0.84 0.71 96.15 99.42 0.85 3.75
## 3604 3604 94.75 1.35 1.82 92.10 97.39 1.43 3.34
## 3671 3671 98.04 0.59 0.35 96.89 99.20 0.60 1.95
## 3672 3672 99.61 0.28 0.08 99.07 100.16 0.28 1.34
## 3673 3673 99.27 0.31 0.10 98.65 99.88 0.32 1.01
## 3674 3674 98.37 0.69 0.48 97.01 99.72 0.70 2.81
## 5101 5101 95.14 1.55 2.40 92.11 98.17 1.63 3.23
## 5102 5102 98.38 0.70 0.49 97.01 99.74 0.71 2.14
## 5103 5103 100.00 0.00 0.00 100.00 100.00 0.00 0.00
## 5104 5104 99.55 0.37 0.14 98.82 100.27 0.37 2.13
## 5105 5105 99.02 0.94 0.88 97.19 100.86 0.94 5.30
## 5106 5106 96.07 1.75 3.06 92.64 99.49 1.82 4.96
## 5107 5107 99.50 0.35 0.12 98.82 100.18 0.35 1.72
## 5108 5108 94.36 1.64 2.69 91.14 97.59 1.74 4.02
## 5171 5171 99.96 0.04 0.00 99.89 100.04 0.04 0.33
## 5201 5201 94.92 1.65 2.72 91.69 98.15 1.74 4.45
## 5202 5202 94.66 1.66 2.76 91.40 97.91 1.75 4.69
## 5203 5203 96.17 1.20 1.44 93.82 98.52 1.25 3.55
## 5204 5204 95.75 1.34 1.80 93.12 98.38 1.40 3.14
## 5205 5205 98.85 0.65 0.42 97.57 100.12 0.66 2.31
## 5206 5206 97.13 1.17 1.37 94.84 99.42 1.20 3.46
## 5207 5207 94.73 1.53 2.34 91.73 97.74 1.62 2.61
## 5208 5208 95.02 2.03 4.12 91.04 99.00 2.14 5.08
## 5271 5271 98.77 0.69 0.48 97.43 100.12 0.70 2.73
## 5272 5272 99.57 0.30 0.09 98.98 100.15 0.30 1.17
## 5301 5301 84.00 3.71 13.76 76.73 91.27 4.41 5.72
## 5302 5302 73.31 4.08 16.65 65.31 81.32 5.57 4.86
## 5303 5303 90.16 2.20 4.84 85.85 94.47 2.44 3.33
## 5304 5304 75.48 3.38 11.42 68.85 82.11 4.48 4.34
## 5305 5305 92.51 2.89 8.35 86.85 98.17 3.12 6.60
## 5306 5306 86.36 3.07 9.42 80.34 92.38 3.56 4.60
## 5307 5307 87.28 2.71 7.34 81.97 92.58 3.10 3.67
## 5308 5308 96.72 1.99 3.96 92.82 100.62 2.06 6.96
## 5309 5309 99.80 0.14 0.02 99.52 100.07 0.14 0.57
## 5310 5310 93.95 2.41 5.81 89.23 98.68 2.56 6.38
## 5311 5311 97.64 0.82 0.67 96.04 99.24 0.84 1.66
## 5312 5312 97.20 1.39 1.93 94.48 99.93 1.43 3.96
## 5313 5313 95.54 2.05 4.20 91.52 99.56 2.15 6.13
## 5314 5314 82.48 3.69 13.62 75.26 89.71 4.47 5.09
## 5315 5315 89.72 3.55 12.60 82.76 96.69 3.96 7.58
## 5316 5316 59.59 4.15 17.22 51.47 67.72 6.96 3.97
## 5317 5317 86.12 3.81 14.52 78.66 93.58 4.42 6.94
## 5318 5318 93.46 2.23 4.97 89.09 97.84 2.39 4.55
## 5319 5319 88.18 3.17 10.05 81.98 94.39 3.59 5.31
## 5320 5320 72.95 4.04 16.32 65.04 80.86 5.53 4.54
## 5321 5321 90.98 2.12 4.49 86.83 95.13 2.33 2.98
## 5371 5371 93.68 3.56 12.67 86.70 100.66 3.80 14.07
## 6101 6101 95.38 1.93 3.72 91.59 99.17 2.03 5.96
## 6102 6102 86.41 4.43 19.62 77.72 95.09 5.13 9.60
## 6103 6103 71.64 4.79 22.94 62.25 81.02 6.68 6.92
## 6104 6104 90.19 2.38 5.66 85.52 94.86 2.64 4.00
## 6105 6105 82.66 3.52 12.39 75.75 89.57 4.26 6.09
## 6106 6106 69.49 4.63 21.44 60.42 78.56 6.66 7.15
## 6107 6107 73.21 3.97 15.76 65.44 80.99 5.42 5.72
## 6108 6108 67.25 4.65 21.62 58.13 76.37 6.92 5.64
## 6109 6109 60.12 6.08 36.97 48.21 72.04 10.11 8.86
## 6110 6110 51.24 5.05 25.50 41.35 61.14 9.85 5.95
## 6111 6111 98.26 0.73 0.53 96.83 99.69 0.74 1.70
## 6112 6112 92.84 1.90 3.61 89.12 96.56 2.04 3.80
## 6171 6171 98.54 0.59 0.35 97.38 99.69 0.60 1.80
## 6172 6172 91.96 1.89 3.57 88.26 95.67 2.06 2.78
## 6201 6201 88.90 2.65 7.02 83.71 94.09 2.98 4.42
## 6202 6202 83.59 3.68 13.54 76.38 90.80 4.40 6.52
## 6203 6203 67.14 4.43 19.62 58.46 75.81 6.59 5.85
## 6204 6204 59.88 5.74 32.95 48.64 71.13 9.58 7.52
## 6205 6205 75.34 5.01 25.10 65.51 85.17 6.65 7.55
## 6206 6206 89.48 2.37 5.62 84.84 94.13 2.65 3.32
## 6207 6207 70.34 4.61 21.25 61.31 79.38 6.55 5.64
## 6208 6208 70.99 5.37 28.84 60.47 81.52 7.56 7.75
## 6209 6209 83.07 4.26 18.15 74.73 91.41 5.12 7.32
## 6210 6210 85.05 3.39 11.49 78.40 91.70 3.99 5.04
## 6211 6211 64.49 5.01 25.10 54.67 74.32 7.77 6.06
## 6212 6212 72.46 3.69 13.62 65.22 79.70 5.10 3.81
## 6213 6213 60.70 4.87 23.72 51.16 70.24 8.02 5.52
## 6271 6271 99.46 0.47 0.22 98.54 100.38 0.47 2.50
## 6301 6301 72.54 4.15 17.22 64.40 80.67 5.72 5.73
## 6302 6302 71.50 3.96 15.68 63.75 79.26 5.53 5.10
## 6303 6303 68.89 4.17 17.39 60.72 77.07 6.05 6.10
## 6304 6304 57.63 5.17 26.73 47.49 67.77 8.98 7.27
## 6305 6305 83.79 3.66 13.40 76.62 90.96 4.37 5.64
## 6306 6306 69.89 5.16 26.63 59.77 80.00 7.38 7.90
## 6307 6307 86.12 3.01 9.06 80.23 92.01 3.49 5.03
## 6308 6308 85.99 3.31 10.96 79.50 92.47 3.85 5.30
## 6309 6309 77.71 4.18 17.47 69.51 85.91 5.38 6.31
## 6310 6310 64.37 3.97 15.76 56.60 72.15 6.16 4.52
## 6311 6311 84.75 3.31 10.96 78.26 91.24 3.91 4.72
## 6371 6371 99.53 0.25 0.06 99.04 100.02 0.25 1.01
## 6372 6372 76.06 4.72 22.28 66.80 85.32 6.21 7.02
## 6401 6401 72.31 4.82 23.23 62.86 81.77 6.67 6.65
## 6402 6402 85.92 2.87 8.24 80.30 91.53 3.34 3.78
## 6403 6403 86.76 3.26 10.63 80.36 93.16 3.76 6.67
## 6404 6404 73.58 4.70 22.09 64.38 82.79 6.38 6.99
## 6405 6405 89.17 2.87 8.24 83.55 94.80 3.22 4.76
## 6409 6409 84.25 3.07 9.42 78.24 90.26 3.64 3.94
## 6411 6411 52.45 4.57 20.88 43.49 61.41 8.72 4.49
## 6471 6471 99.50 0.25 0.06 99.01 99.99 0.25 0.87
## 6472 6472 97.50 0.88 0.77 95.78 99.22 0.90 2.49
## 6474 6474 99.73 0.20 0.04 99.33 100.13 0.20 0.89
## 6501 6501 72.60 4.70 22.09 63.38 81.81 6.48 6.13
## 6502 6502 73.67 4.10 16.81 65.64 81.70 5.56 4.82
## 6503 6503 94.40 1.40 1.96 91.65 97.15 1.49 2.00
## 6504 6504 93.13 1.90 3.61 89.40 96.86 2.04 3.19
## 6571 6571 99.23 0.63 0.40 98.00 100.47 0.64 2.98
## 7101 7101 90.16 1.88 3.53 86.47 93.85 2.09 2.46
## 7102 7102 91.50 1.79 3.20 87.99 95.00 1.96 2.72
## 7103 7103 91.96 2.15 4.62 87.75 96.18 2.34 3.62
## 7104 7104 92.61 1.78 3.17 89.11 96.11 1.93 2.61
## 7105 7105 96.16 1.67 2.79 92.88 99.44 1.74 4.35
## 7106 7106 97.29 1.27 1.61 94.80 99.78 1.30 3.51
## 7107 7107 91.14 1.83 3.35 87.56 94.72 2.01 2.28
## 7108 7108 97.34 0.98 0.96 95.42 99.26 1.01 2.09
## 7109 7109 93.41 1.47 2.16 90.53 96.28 1.57 1.92
## 7110 7110 89.09 2.71 7.34 83.78 94.40 3.04 4.21
## 7111 7111 94.10 1.56 2.43 91.04 97.16 1.66 2.41
## 7171 7171 97.28 1.24 1.54 94.86 99.70 1.27 4.05
## 7172 7172 95.31 1.53 2.34 92.31 98.31 1.61 3.00
## 7173 7173 98.21 0.69 0.48 96.86 99.56 0.70 1.53
## 7174 7174 97.37 1.38 1.90 94.67 100.08 1.42 4.12
## 7201 7201 95.78 1.63 2.66 92.59 98.97 1.70 3.58
## 7202 7202 90.98 3.00 9.00 85.10 96.85 3.29 7.25
## 7203 7203 88.92 3.32 11.02 82.41 95.43 3.74 6.12
## 7204 7204 95.18 1.99 3.96 91.28 99.09 2.09 4.97
## 7205 7205 89.41 2.71 7.34 84.10 94.72 3.03 4.76
## 7206 7206 73.05 4.26 18.15 64.71 81.39 5.83 5.16
## 7207 7207 80.61 3.47 12.04 73.81 87.41 4.30 4.29
## 7208 7208 81.69 3.43 11.76 74.96 88.42 4.20 5.60
## 7209 7209 95.08 1.64 2.69 91.86 98.30 1.73 3.21
## 7210 7210 78.34 4.19 17.56 70.12 86.56 5.35 6.48
## 7211 7211 85.95 3.90 15.21 78.31 93.59 4.54 6.77
## 7212 7212 84.58 2.91 8.47 78.87 90.29 3.44 3.53
## 7271 7271 99.25 0.62 0.38 98.03 100.47 0.63 2.96
## 7301 7301 95.23 2.42 5.86 90.50 99.97 2.54 7.14
## 7302 7302 86.77 2.73 7.45 81.42 92.12 3.15 4.30
## 7303 7303 97.02 1.08 1.17 94.90 99.14 1.12 2.33
## 7304 7304 97.43 1.16 1.35 95.16 99.69 1.19 3.55
## 7305 7305 98.34 0.70 0.49 96.97 99.72 0.71 1.89
## 7306 7306 94.88 1.92 3.69 91.11 98.65 2.03 5.71
## 7307 7307 90.69 2.04 4.16 86.70 94.68 2.25 2.81
## 7308 7308 91.34 2.75 7.56 85.95 96.74 3.01 6.36
## 7309 7309 95.07 1.40 1.96 92.33 97.81 1.47 2.77
## 7310 7310 98.26 0.71 0.50 96.88 99.65 0.72 1.68
## 7311 7311 91.80 2.25 5.06 87.40 96.20 2.45 5.31
## 7312 7312 95.09 3.31 10.96 88.60 101.58 3.48 14.67
## 7313 7313 91.39 3.25 10.56 85.01 97.77 3.56 8.91
## 7314 7314 96.78 1.41 1.99 94.03 99.54 1.45 4.22
## 7315 7315 90.47 2.85 8.12 84.89 96.05 3.15 6.23
## 7316 7316 83.34 3.83 14.67 75.84 90.85 4.59 6.19
## 7317 7317 82.05 3.13 9.80 75.91 88.18 3.82 4.41
## 7318 7318 66.50 4.41 19.45 57.86 75.15 6.63 5.11
## 7322 7322 88.12 1.97 3.88 84.27 91.98 2.23 2.30
## 7325 7325 88.27 3.21 10.30 81.98 94.56 3.64 6.21
## 7326 7326 70.34 4.02 16.16 62.46 78.22 5.72 4.46
## 7371 7371 98.57 0.53 0.28 97.53 99.61 0.54 1.81
## 7372 7372 99.53 0.26 0.07 99.02 100.04 0.26 0.80
## 7373 7373 99.26 0.38 0.14 98.52 100.00 0.38 1.06
## 7401 7401 97.21 1.08 1.17 95.09 99.33 1.11 2.34
## 7402 7402 97.66 1.06 1.12 95.58 99.75 1.09 2.70
## 7403 7403 86.53 3.42 11.70 79.83 93.22 3.95 5.60
## 7404 7404 93.67 1.99 3.96 89.76 97.57 2.13 3.84
## 7405 7405 91.82 2.33 5.43 87.24 96.39 2.54 4.13
## 7406 7406 96.88 1.67 2.79 93.61 100.15 1.72 5.02
## 7407 7407 99.19 0.52 0.27 98.17 100.20 0.52 1.83
## 7408 7408 94.08 1.80 3.24 90.56 97.60 1.91 3.17
## 7409 7409 91.92 1.93 3.72 88.13 95.70 2.10 2.81
## 7410 7410 85.87 2.56 6.55 80.86 90.88 2.98 3.07
## 7411 7411 96.35 1.73 2.99 92.96 99.75 1.80 4.64
## 7412 7412 96.42 0.95 0.90 94.56 98.27 0.98 1.41
## 7413 7413 95.52 1.66 2.76 92.27 98.77 1.73 3.49
## 7414 7414 97.67 1.27 1.61 95.17 100.17 1.31 3.87
## 7415 7415 99.46 0.35 0.12 98.77 100.15 0.35 1.27
## 7471 7471 99.70 0.27 0.07 99.17 100.23 0.27 1.39
## 7472 7472 100.00 0.00 0.00 99.99 100.00 0.00 0.03
## 7501 7501 95.37 1.15 1.32 93.12 97.63 1.21 1.67
## 7502 7502 96.20 1.61 2.59 93.05 99.35 1.67 4.68
## 7503 7503 96.20 1.13 1.28 93.99 98.41 1.17 1.93
## 7504 7504 96.05 1.92 3.69 92.28 99.83 2.00 5.65
## 7505 7505 90.71 2.58 6.66 85.65 95.76 2.85 4.34
## 7571 7571 99.37 0.42 0.18 98.55 100.18 0.42 1.58
## 7601 7601 98.50 0.79 0.62 96.94 100.05 0.80 2.31
## 7602 7602 83.60 3.36 11.29 77.01 90.18 4.02 5.88
## 7603 7603 36.70 5.48 30.03 25.95 47.45 14.94 7.18
## 7604 7604 84.36 2.99 8.94 78.50 90.22 3.54 4.50
## 7605 7605 87.28 2.37 5.62 82.64 91.92 2.71 2.81
## 7606 7606 81.48 4.10 16.81 73.44 89.53 5.04 6.09
## 8101 8101 92.89 2.46 6.05 88.08 97.71 2.65 4.97
## 8102 8102 93.48 2.49 6.20 88.60 98.36 2.66 5.72
## 8103 8103 93.40 2.54 6.45 88.42 98.39 2.72 6.85
## 8104 8104 91.55 2.42 5.86 86.81 96.28 2.64 4.13
## 8105 8105 83.85 3.28 10.76 77.42 90.28 3.91 4.32
## 8106 8106 93.58 1.92 3.69 89.82 97.35 2.05 3.42
## 8107 8107 80.48 3.85 14.82 72.93 88.04 4.79 5.17
## 8108 8108 96.14 1.84 3.39 92.54 99.74 1.91 5.02
## 8109 8109 84.81 3.45 11.90 78.04 91.57 4.07 5.10
## 8171 8171 98.22 0.82 0.67 96.61 99.83 0.84 2.36
## 8172 8172 95.56 1.28 1.64 93.04 98.07 1.34 2.13
## 8201 8201 93.86 2.03 4.12 89.89 97.84 2.16 3.85
## 8202 8202 89.34 2.97 8.82 83.51 95.16 3.33 4.77
## 8203 8203 85.20 3.32 11.02 78.70 91.70 3.89 4.86
## 8204 8204 84.34 3.56 12.67 77.36 91.32 4.22 5.42
## 8205 8205 89.50 2.94 8.64 83.74 95.27 3.28 5.26
## 8206 8206 71.17 4.43 19.62 62.49 79.85 6.22 5.21
## 8207 8207 87.62 2.48 6.15 82.77 92.48 2.83 3.15
## 8208 8208 85.77 3.19 10.18 79.51 92.03 3.72 4.69
## 8271 8271 97.91 0.93 0.86 96.10 99.72 0.95 2.65
## 8272 8272 93.42 2.02 4.08 89.46 97.39 2.17 4.02
## 9101 9101 95.43 1.54 2.37 92.41 98.44 1.61 2.97
## 9102 9102 85.15 2.81 7.90 79.64 90.66 3.30 3.38
## 9103 9103 30.49 4.19 17.56 22.27 38.71 13.75 3.29
## 9104 9104 89.90 2.53 6.40 84.93 94.87 2.82 3.97
## 9105 9105 84.40 2.97 8.82 78.58 90.23 3.52 3.91
## 9111 9111 83.12 2.83 8.01 77.58 88.66 3.40 2.94
## 9112 9112 72.50 4.06 16.48 64.55 80.46 5.60 4.77
## 9201 9201 88.48 2.04 4.16 84.49 92.47 2.30 2.33
## 9202 9202 88.89 2.60 6.76 83.80 93.98 2.92 3.77
## 9203 9203 57.76 5.17 26.73 47.62 67.90 8.96 6.04
## 9204 9204 45.05 4.80 23.04 35.66 54.45 10.64 4.19
## 9205 9205 57.96 7.91 62.57 42.45 73.47 13.65 8.83
## 9271 9271 88.99 2.67 7.13 83.76 94.22 3.00 4.51
## 9403 9403 81.58 3.62 13.10 74.49 88.67 4.44 4.96
## 9408 9408 76.77 4.42 19.54 68.11 85.43 5.76 6.22
## 9409 9409 87.17 3.21 10.30 80.88 93.46 3.68 5.07
## 9419 9419 90.28 2.95 8.70 84.51 96.06 3.26 5.24
## 9420 9420 60.69 3.47 12.04 53.88 67.50 5.72 2.82
## 9426 9426 93.47 1.81 3.28 89.92 97.03 1.94 2.96
## 9427 9427 68.84 2.78 7.73 63.40 74.28 4.04 1.80
## 9428 9428 62.23 7.22 52.13 48.09 76.38 11.60 5.74
## 9471 9471 96.86 1.07 1.14 94.75 98.97 1.11 2.12
## 9501 9501 75.52 3.63 13.18 68.40 82.63 4.81 4.04
## 9502 9502 67.26 3.89 15.13 59.64 74.89 5.78 3.69
## 9503 9503 38.55 5.13 26.32 28.51 48.60 13.29 5.96
## 9504 9504 98.53 0.63 0.40 97.29 99.77 0.64 1.49
## 9601 9601 95.88 1.63 2.66 92.67 99.08 1.70 3.59
## 9602 9602 70.66 4.64 21.53 61.56 79.76 6.57 5.61
## 9603 9603 89.60 3.10 9.61 83.52 95.68 3.46 5.60
## 9604 9604 87.49 3.14 9.86 81.33 93.64 3.59 4.97
## 9605 9605 90.45 4.00 16.00 82.62 98.28 4.42 10.09
## 9606 9606 90.15 2.80 7.84 84.67 95.63 3.10 4.84
## 9607 9607 30.09 6.78 45.97 16.80 43.37 22.52 11.73
## 9608 9608 88.90 3.69 13.62 81.66 96.14 4.16 7.12
## 9701 9701 16.73 3.99 15.92 8.91 24.55 23.84 4.58
## 9702 9702 47.45 5.86 34.34 35.97 58.93 12.34 7.81
## 9703 9703 1.08 0.53 0.28 0.04 2.13 49.21 1.14
## 9704 9704 18.37 4.54 20.61 9.48 27.25 24.69 7.64
## 9705 9705 1.83 0.64 0.41 0.58 3.09 34.99 1.25
## 9706 9706 44.58 4.33 18.75 36.10 53.05 9.70 2.68
## 9707 9707 46.74 6.49 42.12 34.03 59.46 13.88 9.83
## 9708 9708 78.31 3.99 15.92 70.49 86.14 5.10 5.16Visualisasi
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_AIR_LAYAK_PROV_2024.csv")
write.csv(resultK,"RSE_DIRECT_AIR_LAYAK_KAKO_2024.csv")
rm(list=ls()) # Clear environment4 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("HOME_AIR_LAYAK 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 90.08 0.62 88.87 91.30 0.01 2.79 1369082.
## 2 12 92.94 0.38 92.20 93.68 0.00 3.83 3707469.
## 3 13 86.85 0.78 85.32 88.38 0.01 3.53 1394062.
## 4 14 92.24 0.65 90.97 93.51 0.01 5.01 1790826.
## 5 15 82.16 1.00 80.19 84.13 0.01 3.10 948359.
## 6 16 87.23 0.85 85.57 88.89 0.01 6.75 2205567.
## 7 17 72.10 1.32 69.51 74.68 0.02 2.19 531866.
## 8 18 84.51 0.89 82.77 86.25 0.01 6.52 2281298.
## 9 19 82.97 1.33 80.37 85.57 0.02 2.36 399637.
## 10 21 94.71 0.53 93.68 95.74 0.01 1.81 690501.
## 11 31 99.96 0.03 99.91 100.01 0.00 2.09 2880360.
## 12 32 94.90 0.32 94.26 95.53 0.00 14.15 13924466.
## 13 33 95.43 0.28 94.87 95.98 0.00 8.21 9512380.
## 14 34 96.91 0.49 95.96 97.86 0.01 4.61 1238591.
## 15 35 96.93 0.21 96.50 97.35 0.00 8.24 11258540.
## 16 36 94.01 0.53 92.98 95.04 0.01 8.02 3450320.
## 17 51 98.32 0.31 97.73 98.92 0.00 3.25 1212690.
## 18 52 96.13 0.53 95.09 97.18 0.01 5.60 1537009.
## 19 53 88.55 0.70 87.19 89.92 0.01 2.90 1282950.
## 20 61 82.13 0.98 80.21 84.05 0.01 3.98 1280298.
## 21 62 78.71 1.20 76.37 81.06 0.02 3.07 757265.
## 22 63 77.34 1.10 75.19 79.50 0.01 3.96 1207351.
## 23 64 89.11 1.03 87.10 91.13 0.01 4.98 964156.
## 24 65 89.42 1.15 87.16 91.67 0.01 1.17 176652.
## 25 71 94.46 0.48 93.52 95.40 0.01 1.39 670640.
## 26 72 87.86 0.92 86.06 89.65 0.01 2.90 779354.
## 27 73 92.23 0.48 91.29 93.17 0.01 3.26 2149860.
## 28 74 95.31 0.51 94.30 96.32 0.01 1.85 660674.
## 29 75 96.13 0.65 94.86 97.39 0.01 1.56 293537.
## 30 76 80.14 1.50 77.19 83.08 0.02 2.21 328206.
## 31 81 93.32 0.74 91.87 94.77 0.01 1.65 396054.
## 32 82 88.92 0.96 87.03 90.80 0.01 1.28 287599.
## 33 91 82.38 1.35 79.74 85.02 0.02 0.69 117078.
## 34 92 81.68 1.67 78.41 84.95 0.02 1.12 126648.
## 35 94 85.84 1.10 83.68 88.00 0.01 1.00 212299.
## 36 95 71.90 2.09 67.80 76.00 0.03 1.28 124940.
## 37 96 81.53 1.48 78.63 84.42 0.02 1.92 279213.
## 38 97 30.64 1.90 26.91 34.38 0.06 2.34 288746.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("HOME_AIR_LAYAK 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 82.79 2.72 77.47 88.12 0.03 0.59 23829.
## 2 1102 79.60 3.36 73.02 86.19 0.04 0.98 29563.
## 3 1103 85.87 3.62 78.77 92.96 0.04 3.03 59057.
## 4 1104 91.40 3.10 85.32 97.48 0.03 3.26 56147.
## 5 1105 84.43 3.27 78.02 90.84 0.04 4.09 105925.
## 6 1106 88.27 3.59 81.23 95.31 0.04 3.42 57832.
## 7 1107 96.52 1.53 93.51 99.53 0.02 1.86 55906.
## 8 1108 90.68 2.74 85.30 96.05 0.03 4.37 103486.
## 9 1109 87.57 3.07 81.55 93.58 0.04 4.80 116781.
## 10 1110 90.06 2.37 85.41 94.70 0.03 3.45 116007.
## 11 1111 94.56 1.35 91.91 97.20 0.01 2.59 154677.
## 12 1112 95.12 1.35 92.47 97.76 0.01 0.68 36478.
## 13 1113 80.93 3.10 74.86 87.01 0.04 0.76 25750.
## 14 1114 89.02 2.59 83.95 94.09 0.03 2.43 74684.
## 15 1115 90.75 2.93 85.01 96.49 0.03 2.23 45954.
## 16 1116 87.15 3.48 80.32 93.98 0.04 1.41 27372.
## 17 1117 84.98 3.84 77.44 92.51 0.05 2.27 41370.
## 18 1118 94.54 2.11 90.40 98.67 0.02 1.81 44344.
## 19 1171 99.05 0.62 97.84 100.26 0.01 1.36 70449.
## 20 1172 89.92 2.88 84.28 95.56 0.03 0.40 9284.
## 21 1173 98.06 0.93 96.24 99.89 0.01 0.97 44597.
## 22 1174 97.41 1.09 95.26 99.55 0.01 1.14 50828.
## 23 1175 57.25 4.51 48.40 66.10 0.08 0.74 18764.
## 24 1201 48.88 5.24 38.60 59.16 0.11 1.53 29237.
## 25 1202 84.72 3.04 78.77 90.67 0.04 3.73 110358.
## 26 1203 74.84 4.85 65.33 84.35 0.06 4.24 71439.
## 27 1204 82.98 2.41 78.26 87.70 0.03 1.75 89893.
## 28 1205 92.91 2.47 88.06 97.75 0.03 3.26 74153.
## 29 1206 96.01 1.60 92.87 99.14 0.02 1.51 47688.
## 30 1207 93.42 2.23 89.05 97.79 0.02 4.90 127888.
## 31 1208 94.60 1.62 91.42 97.78 0.02 4.45 181490.
## 32 1209 99.09 0.50 98.10 100.08 0.01 3.08 228838.
## 33 1210 97.62 0.89 95.87 99.38 0.01 1.17 71310.
## 34 1211 96.61 1.57 93.53 99.69 0.02 4.48 125082.
## 35 1212 96.37 1.20 94.02 98.72 0.01 11.51 589305.
## 36 1213 92.49 1.92 88.73 96.25 0.02 6.86 272607.
## 37 1214 69.38 5.12 59.34 79.41 0.07 4.26 72720.
## 38 1215 97.73 0.86 96.05 99.42 0.01 0.75 47254.
## 39 1216 68.97 3.38 62.36 75.59 0.05 0.32 12498.
## 40 1217 80.38 3.37 73.78 86.97 0.04 1.09 31836.
## 41 1218 98.43 0.63 97.21 99.66 0.01 1.86 154731.
## 42 1219 97.11 1.75 93.67 100.55 0.02 5.31 101789.
## 43 1220 82.23 4.32 73.76 90.70 0.05 4.64 76466.
## 44 1221 89.96 2.18 85.69 94.23 0.02 1.86 74875.
## 45 1222 90.05 2.58 84.99 95.11 0.03 3.23 91338.
## 46 1223 91.50 2.52 86.57 96.43 0.03 3.47 89792.
## 47 1224 71.72 4.60 62.72 80.73 0.06 1.48 29934.
## 48 1225 70.95 3.92 63.27 78.63 0.06 0.59 16814.
## 49 1271 96.86 1.92 93.09 100.63 0.02 1.10 19118.
## 50 1272 91.27 2.60 86.17 96.37 0.03 1.70 42123.
## 51 1273 99.74 0.15 99.45 100.03 0.00 0.25 63697.
## 52 1274 99.53 0.27 99.00 100.06 0.00 0.32 43082.
## 53 1275 99.84 0.15 99.55 100.12 0.00 3.51 560951.
## 54 1276 98.88 0.61 97.68 100.08 0.01 1.18 73833.
## 55 1277 64.99 4.71 55.76 74.21 0.07 2.51 54218.
## 56 1278 80.17 3.43 73.45 86.88 0.04 1.09 31110.
## 57 1301 70.23 6.01 58.45 82.02 0.09 2.08 25312.
## 58 1302 90.00 1.94 86.19 93.80 0.02 2.22 111544.
## 59 1303 81.84 3.13 75.71 87.97 0.04 2.87 91738.
## 60 1304 72.02 3.84 64.49 79.54 0.05 2.17 62545.
## 61 1305 92.74 2.25 88.33 97.14 0.02 3.18 89284.
## 62 1306 91.77 1.69 88.47 95.07 0.02 1.73 96929.
## 63 1307 89.68 1.99 85.77 93.58 0.02 2.55 125569.
## 64 1308 68.64 7.08 54.76 82.52 0.10 11.64 105318.
## 65 1309 86.06 2.48 81.20 90.92 0.03 1.71 70475.
## 66 1310 79.09 3.31 72.61 85.57 0.04 1.44 45933.
## 67 1311 70.00 3.70 62.75 77.26 0.05 2.15 69541.
## 68 1312 83.07 3.01 77.17 88.96 0.04 3.66 119991.
## 69 1371 97.65 0.70 96.27 99.02 0.01 2.44 240836.
## 70 1372 97.04 1.33 94.44 99.64 0.01 0.54 18650.
## 71 1373 89.62 2.27 85.18 94.06 0.03 0.42 16171.
## 72 1374 98.20 0.73 96.77 99.63 0.01 0.19 13556.
## 73 1375 99.45 0.38 98.70 100.20 0.00 0.44 34882.
## 74 1376 99.44 0.41 98.63 100.24 0.00 0.51 35547.
## 75 1377 95.62 1.28 93.11 98.12 0.01 0.37 20240.
## 76 1401 86.73 3.04 80.77 92.69 0.04 3.16 82857.
## 77 1402 77.63 3.33 71.10 84.17 0.04 3.38 111147.
## 78 1403 94.83 1.51 91.86 97.80 0.02 4.01 180483.
## 79 1404 83.26 4.03 75.36 91.15 0.05 8.16 147691.
## 80 1405 91.79 1.96 87.95 95.64 0.02 3.12 128867.
## 81 1406 96.00 1.59 92.88 99.12 0.02 7.00 223277.
## 82 1407 93.45 1.96 89.61 97.29 0.02 5.78 194404.
## 83 1408 95.72 1.41 92.95 98.49 0.01 3.14 135878.
## 84 1409 86.82 3.03 80.88 92.76 0.03 6.70 175747.
## 85 1410 93.44 2.24 89.05 97.83 0.02 1.73 44546.
## 86 1471 99.65 0.33 99.00 100.29 0.00 4.21 290848.
## 87 1473 94.82 1.56 91.75 97.88 0.02 1.77 75081.
## 88 1501 76.34 4.94 66.65 86.03 0.06 4.65 72449.
## 89 1502 60.74 4.65 51.63 69.85 0.08 4.40 102410.
## 90 1503 76.41 3.35 69.85 82.97 0.04 2.36 79982.
## 91 1504 84.80 3.13 78.67 90.94 0.04 2.46 68345.
## 92 1505 78.46 3.01 72.56 84.36 0.04 3.05 119688.
## 93 1506 87.12 2.55 82.11 92.13 0.03 1.60 57840.
## 94 1507 93.84 2.16 89.61 98.06 0.02 3.32 87102.
## 95 1508 80.71 3.70 73.46 87.97 0.05 3.96 94801.
## 96 1509 81.93 3.00 76.06 87.81 0.04 2.79 96887.
## 97 1571 96.06 1.10 93.91 98.21 0.01 2.20 145098.
## 98 1572 89.50 3.46 82.73 96.27 0.04 1.43 23758.
## 99 1601 82.55 2.97 76.72 88.38 0.04 2.73 93743.
## 100 1602 78.41 3.70 71.17 85.65 0.05 8.65 225866.
## 101 1603 91.18 2.08 87.11 95.26 0.02 4.15 162582.
## 102 1604 84.86 2.83 79.31 90.42 0.03 3.20 108071.
## 103 1605 91.39 2.55 86.40 96.39 0.03 4.34 110917.
## 104 1606 86.11 2.93 80.36 91.85 0.03 5.85 171573.
## 105 1607 81.42 3.31 74.93 87.91 0.04 7.71 224120.
## 106 1608 82.60 3.47 75.81 89.40 0.04 4.09 102956.
## 107 1609 96.19 1.16 93.92 98.47 0.01 3.15 180710.
## 108 1610 73.69 4.75 64.38 83.00 0.06 6.04 109357.
## 109 1611 69.33 5.03 59.48 79.18 0.07 3.70 65704.
## 110 1612 94.09 2.28 89.62 98.57 0.02 2.19 49083.
## 111 1613 88.41 2.53 83.44 93.37 0.03 1.61 54083.
## 112 1671 95.00 2.68 89.75 100.26 0.03 28.79 400691.
## 113 1672 94.94 1.62 91.76 98.12 0.02 1.30 50169.
## 114 1673 89.81 3.04 83.85 95.77 0.03 1.75 36449.
## 115 1674 97.18 0.89 95.43 98.93 0.01 0.82 59492.
## 116 1701 68.26 4.95 58.56 77.96 0.07 2.42 45103.
## 117 1702 55.60 4.79 46.21 64.98 0.09 3.36 76188.
## 118 1703 80.40 3.02 74.48 86.31 0.04 2.26 82243.
## 119 1704 79.79 2.65 74.60 84.98 0.03 0.65 31655.
## 120 1705 70.46 4.51 61.62 79.29 0.06 2.41 52101.
## 121 1706 75.82 3.94 68.10 83.54 0.05 1.92 47840.
## 122 1707 57.41 4.31 48.95 65.86 0.08 1.15 31735.
## 123 1708 61.20 5.68 50.07 72.34 0.09 2.46 38185.
## 124 1709 50.44 4.96 40.73 60.16 0.10 1.40 30108.
## 125 1771 92.21 2.14 88.02 96.41 0.02 2.93 96708.
## 126 1801 64.37 5.35 53.89 74.86 0.08 4.86 81985.
## 127 1802 68.48 4.06 60.52 76.45 0.06 5.69 156759.
## 128 1803 95.52 1.10 93.36 97.67 0.01 3.70 276765.
## 129 1804 87.01 2.96 81.20 92.82 0.03 10.70 290341.
## 130 1805 84.59 3.31 78.11 91.07 0.04 13.98 351053.
## 131 1806 69.31 4.38 60.73 77.89 0.06 6.62 154690.
## 132 1807 71.40 3.78 64.00 78.80 0.05 4.09 123364.
## 133 1808 88.55 2.92 82.82 94.27 0.03 5.13 128604.
## 134 1809 89.80 2.51 84.88 94.73 0.03 3.97 121334.
## 135 1810 93.33 1.66 90.08 96.59 0.02 2.22 105453.
## 136 1811 81.81 3.88 74.22 89.41 0.05 2.69 56169.
## 137 1812 63.16 6.38 50.65 75.68 0.10 6.40 76933.
## 138 1813 73.77 3.61 66.69 80.85 0.05 1.32 41379.
## 139 1871 99.09 0.51 98.10 100.09 0.01 3.71 272181.
## 140 1872 96.25 1.72 92.88 99.63 0.02 1.73 44288.
## 141 1901 77.35 3.97 69.56 85.14 0.05 3.89 90875.
## 142 1902 83.36 2.96 77.56 89.17 0.04 1.54 51378.
## 143 1903 84.94 3.43 78.22 91.65 0.04 2.50 57385.
## 144 1904 88.25 2.30 83.75 92.76 0.03 1.26 51937.
## 145 1905 78.30 3.71 71.04 85.56 0.05 2.07 54061.
## 146 1906 70.94 3.88 63.33 78.54 0.05 1.24 35804.
## 147 1971 96.47 1.15 94.21 98.73 0.01 1.08 58196.
## 148 2101 85.89 3.39 79.24 92.54 0.04 2.76 61155.
## 149 2102 84.57 3.15 78.40 90.75 0.04 1.55 43029.
## 150 2103 91.53 2.34 86.94 96.12 0.03 0.72 21510.
## 151 2104 85.48 3.55 78.52 92.43 0.04 1.18 24472.
## 152 2105 24.65 4.61 15.61 33.68 0.19 0.60 11121.
## 153 2171 99.40 0.38 98.65 100.14 0.00 5.36 471419.
## 154 2172 91.96 1.93 88.18 95.74 0.02 1.38 57795.
## 155 3101 92.93 0.76 91.43 94.43 0.01 0.03 6739.
## 156 3171 100.00 0.00 100.00 100.00 0.00 4.23 609825.
## 157 3172 100.00 0.00 100.00 100.00 0.00 4.58 786125.
## 158 3173 100.00 0.00 100.00 100.00 0.00 4.50 253780.
## 159 3174 100.00 0.00 100.00 100.00 0.00 7.95 723310.
## 160 3175 99.85 0.15 99.57 100.14 0.00 3.46 500581.
## 161 3201 93.93 1.56 90.87 96.99 0.02 32.45 1595241.
## 162 3202 89.30 1.87 85.64 92.97 0.02 12.26 705558.
## 163 3203 91.62 1.63 88.41 94.82 0.02 10.78 653081.
## 164 3204 96.28 1.14 94.04 98.52 0.01 17.46 1005421.
## 165 3205 85.37 2.19 81.09 89.65 0.03 12.19 671674.
## 166 3206 86.78 2.34 82.20 91.37 0.03 11.75 518963.
## 167 3207 96.42 0.90 94.67 98.18 0.01 4.09 370607.
## 168 3208 97.92 0.98 96.00 99.85 0.01 6.64 295040.
## 169 3209 97.12 0.87 95.42 98.83 0.01 7.76 604921.
## 170 3210 94.38 1.67 91.12 97.65 0.02 9.15 367764.
## 171 3211 96.23 1.26 93.76 98.70 0.01 7.43 357724.
## 172 3212 97.39 0.92 95.59 99.20 0.01 8.42 531611.
## 173 3213 98.69 0.44 97.81 99.56 0.00 3.53 487121.
## 174 3214 91.02 2.26 86.59 95.44 0.02 7.59 256514.
## 175 3215 97.23 0.72 95.82 98.63 0.01 6.03 668717.
## 176 3216 96.76 1.09 94.63 98.89 0.01 20.85 1164693.
## 177 3217 89.25 2.34 84.67 93.83 0.03 12.43 459993.
## 178 3218 88.36 2.64 83.20 93.53 0.03 3.94 122856.
## 179 3271 99.11 0.35 98.43 99.79 0.00 1.85 286528.
## 180 3272 99.46 0.35 98.78 100.14 0.00 0.90 85597.
## 181 3273 97.62 1.04 95.58 99.67 0.01 15.84 709458.
## 182 3274 98.23 0.81 96.64 99.82 0.01 1.56 87235.
## 183 3275 98.89 0.84 97.24 100.55 0.01 25.84 838899.
## 184 3276 99.85 0.11 99.63 100.07 0.00 2.65 680217.
## 185 3277 99.60 0.30 99.02 100.19 0.00 1.77 165560.
## 186 3278 93.79 2.26 89.36 98.23 0.02 7.53 180489.
## 187 3279 96.97 1.18 94.67 99.28 0.01 1.19 52985.
## 188 3301 93.33 1.46 90.47 96.19 0.02 7.61 469454.
## 189 3302 94.33 1.53 91.34 97.33 0.02 9.73 469467.
## 190 3303 88.24 2.00 84.33 92.15 0.02 4.37 239946.
## 191 3304 91.75 1.99 87.85 95.65 0.02 6.05 243839.
## 192 3305 87.29 2.47 82.45 92.12 0.03 8.35 320660.
## 193 3306 90.04 2.54 85.06 95.02 0.03 7.11 207954.
## 194 3307 95.24 1.46 92.37 98.11 0.02 4.87 216925.
## 195 3308 93.27 1.85 89.65 96.90 0.02 9.29 358453.
## 196 3309 94.83 1.95 91.01 98.64 0.02 10.20 278539.
## 197 3310 99.48 0.35 98.79 100.17 0.00 3.90 340044.
## 198 3311 99.15 0.41 98.36 99.95 0.00 2.34 251446.
## 199 3312 96.96 1.19 94.62 99.30 0.01 6.06 263683.
## 200 3313 99.64 0.21 99.24 100.05 0.00 1.39 240515.
## 201 3314 98.45 0.60 97.27 99.62 0.01 2.88 257095.
## 202 3315 87.09 2.29 82.61 91.58 0.03 8.88 402401.
## 203 3316 97.02 1.15 94.78 99.27 0.01 5.36 247865.
## 204 3317 95.86 1.67 92.57 99.14 0.02 6.07 181188.
## 205 3318 98.32 1.06 96.23 100.40 0.01 11.99 369800.
## 206 3319 99.60 0.28 99.05 100.14 0.00 2.11 228174.
## 207 3320 89.46 2.69 84.18 94.74 0.03 12.79 349995.
## 208 3321 95.99 1.91 92.25 99.72 0.02 14.42 322424.
## 209 3322 99.60 0.20 99.22 99.99 0.00 1.38 297540.
## 210 3323 97.83 1.00 95.86 99.80 0.01 4.68 207027.
## 211 3324 97.65 1.36 94.99 100.32 0.01 10.29 269296.
## 212 3325 96.46 1.40 93.71 99.21 0.01 5.48 200547.
## 213 3326 96.91 1.39 94.17 99.64 0.01 6.59 214299.
## 214 3327 94.38 2.08 90.32 98.45 0.02 12.41 321360.
## 215 3328 98.28 0.92 96.48 100.07 0.01 8.70 369560.
## 216 3329 94.66 1.41 91.88 97.43 0.01 8.93 475705.
## 217 3371 99.93 0.07 99.78 100.07 0.00 0.11 32684.
## 218 3372 99.72 0.28 99.18 100.27 0.00 1.94 147404.
## 219 3373 98.98 0.67 97.67 100.28 0.01 1.17 56357.
## 220 3374 99.87 0.13 99.63 100.12 0.00 3.03 512063.
## 221 3375 99.27 0.37 98.55 99.99 0.00 0.70 80446.
## 222 3376 99.66 0.34 99.01 100.32 0.00 1.09 68223.
## 223 3401 92.93 1.68 89.63 96.22 0.02 2.61 127621.
## 224 3402 97.96 0.76 96.47 99.44 0.01 4.30 316489.
## 225 3403 98.47 0.47 97.54 99.40 0.00 1.59 226311.
## 226 3404 95.37 1.19 93.05 97.70 0.01 6.32 416432.
## 227 3471 99.94 0.06 99.83 100.06 0.00 0.43 151737.
## 228 3501 84.35 2.93 78.60 90.09 0.03 4.87 157678.
## 229 3502 95.57 1.58 92.48 98.65 0.02 6.90 248213.
## 230 3503 81.47 2.96 75.67 87.26 0.04 5.53 201105.
## 231 3504 98.62 0.77 97.11 100.13 0.01 6.14 297024.
## 232 3505 98.66 0.79 97.12 100.20 0.01 7.50 337500.
## 233 3506 96.72 1.45 93.88 99.57 0.02 13.69 432573.
## 234 3507 97.94 0.85 96.27 99.61 0.01 12.47 726036.
## 235 3508 97.31 0.87 95.59 99.02 0.01 4.05 292918.
## 236 3509 95.30 1.02 93.30 97.30 0.01 7.88 711964.
## 237 3510 97.06 0.95 95.20 98.93 0.01 7.37 487496.
## 238 3511 96.54 1.38 93.83 99.26 0.01 7.05 258914.
## 239 3512 95.17 2.14 90.97 99.36 0.02 10.64 225207.
## 240 3513 96.90 1.24 94.47 99.33 0.01 8.22 337782.
## 241 3514 98.62 0.60 97.44 99.80 0.01 5.76 458036.
## 242 3515 98.07 0.85 96.41 99.73 0.01 11.68 646783.
## 243 3516 97.46 1.45 94.62 100.31 0.01 12.54 310094.
## 244 3517 98.76 0.64 97.52 100.01 0.01 5.36 341556.
## 245 3518 99.71 0.23 99.25 100.16 0.00 2.53 296064.
## 246 3519 97.55 1.10 95.40 99.70 0.01 4.87 204687.
## 247 3520 98.56 0.66 97.27 99.85 0.01 2.58 178252.
## 248 3521 97.07 1.13 94.86 99.28 0.01 5.56 262231.
## 249 3522 96.13 1.68 92.84 99.42 0.02 12.50 347002.
## 250 3523 96.96 1.39 94.24 99.68 0.01 10.04 322850.
## 251 3524 84.60 3.45 77.84 91.37 0.04 13.49 310643.
## 252 3525 97.65 0.92 95.85 99.45 0.01 6.14 352286.
## 253 3526 98.00 0.90 96.23 99.76 0.01 4.92 250992.
## 254 3527 97.66 1.04 95.62 99.70 0.01 5.83 259827.
## 255 3528 96.98 0.92 95.19 98.77 0.01 3.28 241458.
## 256 3529 99.41 0.23 98.96 99.85 0.00 1.41 339395.
## 257 3571 99.73 0.20 99.34 100.12 0.00 0.55 77712.
## 258 3572 99.60 0.34 98.94 100.27 0.00 0.54 39104.
## 259 3573 99.48 0.34 98.82 100.14 0.00 2.49 240447.
## 260 3574 100.00 0.00 100.00 100.00 0.00 0.32 63050.
## 261 3575 99.94 0.06 99.83 100.06 0.00 0.14 52258.
## 262 3576 99.77 0.15 99.47 100.07 0.00 0.17 35393.
## 263 3577 100.00 0.00 100.00 100.00 0.00 0.43 51457.
## 264 3578 100.00 0.00 100.00 100.00 0.00 6.52 806086.
## 265 3579 98.27 1.06 96.19 100.36 0.01 1.78 56468.
## 266 3601 78.51 2.78 73.06 83.97 0.04 6.77 310686.
## 267 3602 77.92 3.31 71.43 84.41 0.04 10.39 343922.
## 268 3603 97.79 0.84 96.15 99.42 0.01 15.93 1042236.
## 269 3604 94.75 1.35 92.10 97.39 0.01 6.24 358024.
## 270 3671 98.04 0.59 96.89 99.20 0.01 5.45 635848.
## 271 3672 99.61 0.28 99.07 100.16 0.00 1.07 111203.
## 272 3673 99.27 0.31 98.65 99.88 0.00 1.07 166097.
## 273 3674 98.37 0.69 97.01 99.72 0.01 6.83 482304.
## 274 5101 95.14 1.55 92.11 98.17 0.02 1.95 79401.
## 275 5102 98.38 0.70 97.01 99.74 0.01 1.66 114950.
## 276 5103 100.00 0.00 100.00 100.00 0.00 3.95 203236.
## 277 5104 99.55 0.37 98.82 100.27 0.00 1.67 116176.
## 278 5105 99.02 0.94 97.19 100.86 0.01 2.01 46814.
## 279 5106 96.07 1.75 92.64 99.49 0.02 2.28 59599.
## 280 5107 99.50 0.35 98.82 100.18 0.00 1.27 110983.
## 281 5108 94.36 1.64 91.14 97.59 0.02 4.54 188280.
## 282 5171 99.96 0.04 99.89 100.04 0.00 0.52 293251.
## 283 5201 94.92 1.65 91.69 98.15 0.02 5.74 215024.
## 284 5202 94.66 1.66 91.40 97.91 0.02 7.84 302914.
## 285 5203 96.17 1.20 93.82 98.52 0.01 7.00 377938.
## 286 5204 95.75 1.34 93.12 98.38 0.01 2.71 129368.
## 287 5205 98.85 0.65 97.57 100.12 0.01 1.21 68892.
## 288 5206 97.13 1.17 94.84 99.42 0.01 3.10 133240.
## 289 5207 94.73 1.53 91.73 97.74 0.02 1.00 44579.
## 290 5208 95.02 2.03 91.04 99.00 0.02 2.75 66569.
## 291 5271 98.77 0.69 97.43 100.12 0.01 2.76 149571.
## 292 5272 99.57 0.30 98.98 100.15 0.00 0.48 48914.
## 293 5301 84.00 3.71 76.73 91.27 0.04 1.33 27488.
## 294 5302 73.31 4.08 65.31 81.32 0.06 2.26 55902.
## 295 5303 90.16 2.20 85.85 94.47 0.02 2.57 99433.
## 296 5304 75.48 3.38 68.85 82.11 0.04 3.51 119650.
## 297 5305 92.51 2.89 86.85 98.17 0.03 3.51 61470.
## 298 5306 86.36 3.07 80.34 92.38 0.04 2.00 52585.
## 299 5307 87.28 2.71 81.97 92.58 0.03 1.46 46514.
## 300 5308 96.72 1.99 92.82 100.62 0.02 2.22 37388.
## 301 5309 99.80 0.14 99.52 100.07 0.00 0.28 59704.
## 302 5310 93.95 2.41 89.23 98.68 0.03 3.40 70054.
## 303 5311 97.64 0.82 96.04 99.24 0.01 0.81 58886.
## 304 5312 97.20 1.39 94.48 99.93 0.01 1.22 36023.
## 305 5313 95.54 2.05 91.52 99.56 0.02 3.55 75802.
## 306 5314 82.48 3.69 75.26 89.71 0.04 2.20 49333.
## 307 5315 89.72 3.55 82.75 96.69 0.04 4.31 66248.
## 308 5316 59.59 4.15 51.47 67.72 0.07 0.51 15169.
## 309 5317 86.12 3.81 78.66 93.58 0.04 3.90 67847.
## 310 5318 93.46 2.23 89.09 97.84 0.02 1.15 29662.
## 311 5319 88.18 3.17 81.98 94.39 0.04 2.93 64136.
## 312 5320 72.95 4.04 65.04 80.86 0.06 1.05 26847.
## 313 5321 90.98 2.12 86.83 95.13 0.02 1.23 47562.
## 314 5371 93.68 3.56 86.70 100.66 0.04 11.71 115249.
## 315 6101 95.38 1.93 91.59 99.17 0.02 5.25 130377.
## 316 6102 86.41 4.43 77.72 95.09 0.05 5.09 64131.
## 317 6103 71.64 4.79 62.25 81.02 0.07 4.81 89779.
## 318 6104 90.19 2.38 85.52 94.86 0.03 1.93 63606.
## 319 6105 82.66 3.52 75.75 89.57 0.04 5.00 121602.
## 320 6106 69.49 4.63 60.42 78.56 0.07 6.68 139442.
## 321 6107 73.21 3.97 65.44 80.99 0.05 4.32 113453.
## 322 6108 67.25 4.65 58.13 76.37 0.07 3.31 70937.
## 323 6109 60.12 6.08 48.21 72.04 0.10 3.93 53747.
## 324 6110 51.24 5.05 41.35 61.14 0.10 2.63 54376.
## 325 6111 98.26 0.73 96.83 99.69 0.01 0.41 27639.
## 326 6112 92.84 1.90 89.12 96.56 0.02 3.62 141016.
## 327 6171 98.54 0.59 97.38 99.69 0.01 1.78 157007.
## 328 6172 91.96 1.89 88.26 95.67 0.02 1.22 53185.
## 329 6201 88.90 2.65 83.71 94.09 0.03 3.00 88951.
## 330 6202 83.59 3.68 76.38 90.80 0.04 6.02 128734.
## 331 6203 67.14 4.43 58.46 75.81 0.07 4.15 98499.
## 332 6204 59.88 5.74 48.64 71.13 0.10 2.40 36870.
## 333 6205 75.34 5.01 65.51 85.17 0.07 2.29 35630.
## 334 6206 89.48 2.37 84.84 94.13 0.03 0.55 19267.
## 335 6207 70.34 4.61 61.31 79.38 0.07 1.15 23748.
## 336 6208 70.99 5.37 60.47 81.52 0.08 4.29 64543.
## 337 6209 83.07 4.26 74.73 91.41 0.05 2.98 48813.
## 338 6210 85.05 3.39 78.40 91.70 0.04 1.51 35052.
## 339 6211 64.49 5.01 54.67 74.32 0.08 1.57 30053.
## 340 6212 72.46 3.69 65.22 79.70 0.05 1.17 35960.
## 341 6213 60.70 4.87 51.16 70.24 0.08 1.43 30367.
## 342 6271 99.46 0.47 98.54 100.38 0.00 1.57 80777.
## 343 6301 72.54 4.15 64.40 80.67 0.06 4.06 99083.
## 344 6302 71.50 3.96 63.75 79.26 0.06 3.51 96407.
## 345 6303 68.89 4.17 60.72 77.07 0.06 6.25 162115.
## 346 6304 57.63 5.17 47.49 67.78 0.09 4.59 88222.
## 347 6305 83.79 3.66 76.62 90.96 0.04 2.62 56029.
## 348 6306 69.89 5.16 59.77 80.00 0.07 4.10 68350.
## 349 6307 86.12 3.01 80.23 92.01 0.03 2.81 78323.
## 350 6308 85.99 3.31 79.50 92.47 0.04 2.75 63877.
## 351 6309 77.71 4.18 69.51 85.91 0.05 3.38 70554.
## 352 6310 64.37 3.97 56.59 72.15 0.06 3.64 111863.
## 353 6311 84.75 3.31 78.26 91.24 0.04 1.54 38204.
## 354 6371 99.53 0.25 99.04 100.02 0.00 1.25 193971.
## 355 6372 76.06 4.72 66.80 85.32 0.06 4.67 80353.
## 356 6401 72.31 4.82 62.86 81.77 0.07 4.21 76418.
## 357 6402 85.92 2.87 80.30 91.53 0.03 1.22 37839.
## 358 6403 86.76 3.26 80.36 93.16 0.04 8.99 204215.
## 359 6404 73.58 4.70 64.38 82.79 0.06 5.87 109053.
## 360 6405 89.17 2.87 83.55 94.80 0.03 2.38 58737.
## 361 6409 84.25 3.07 78.24 90.26 0.04 1.40 41521.
## 362 6411 52.45 4.57 43.49 61.41 0.09 0.28 6957.
## 363 6471 99.50 0.25 99.01 99.99 0.00 0.96 162378.
## 364 6472 97.50 0.88 95.78 99.22 0.01 3.29 218926.
## 365 6474 99.73 0.20 99.33 100.13 0.00 0.35 48111.
## 366 6501 72.60 4.70 63.38 81.81 0.06 1.08 20406.
## 367 6502 73.67 4.10 65.64 81.70 0.06 1.34 32625.
## 368 6503 94.40 1.40 91.65 97.15 0.01 0.14 8053.
## 369 6504 93.13 1.90 89.40 96.86 0.02 1.40 52076.
## 370 6571 99.23 0.63 98.00 100.47 0.01 1.57 63493.
## 371 7101 90.16 1.88 86.47 93.85 0.02 1.21 63632.
## 372 7102 91.50 1.79 87.99 95.00 0.02 1.90 97195.
## 373 7103 91.96 2.15 87.75 96.18 0.02 1.07 35884.
## 374 7104 92.61 1.78 89.11 96.11 0.02 0.52 23682.
## 375 7105 96.16 1.67 92.88 99.44 0.02 2.10 58600.
## 376 7106 97.29 1.27 94.80 99.78 0.01 1.62 55894.
## 377 7107 91.14 1.83 87.56 94.72 0.02 0.41 20965.
## 378 7108 97.34 0.98 95.42 99.26 0.01 0.31 17933.
## 379 7109 93.41 1.47 90.53 96.28 0.02 0.46 27802.
## 380 7110 89.09 2.71 83.78 94.40 0.03 0.59 16475.
## 381 7111 94.10 1.56 91.04 97.16 0.02 0.38 18294.
## 382 7171 97.28 1.24 94.86 99.70 0.01 3.12 113911.
## 383 7172 95.31 1.53 92.31 98.31 0.02 1.41 56867.
## 384 7173 98.21 0.69 96.86 99.56 0.01 0.38 29740.
## 385 7174 97.37 1.38 94.67 100.08 0.01 1.19 33766.
## 386 7201 95.78 1.63 92.59 98.97 0.02 0.99 31661.
## 387 7202 90.98 3.00 85.10 96.85 0.03 5.26 101415.
## 388 7203 88.92 3.32 82.40 95.43 0.04 1.81 33969.
## 389 7204 95.18 1.99 91.28 99.09 0.02 2.81 68388.
## 390 7205 89.41 2.71 84.10 94.72 0.03 2.68 72798.
## 391 7206 73.05 4.26 64.71 81.39 0.06 2.54 58082.
## 392 7207 80.61 3.47 73.81 87.41 0.04 1.45 39582.
## 393 7208 81.69 3.43 74.96 88.42 0.04 4.51 120573.
## 394 7209 95.08 1.64 91.86 98.30 0.02 1.01 36923.
## 395 7210 78.34 4.19 70.12 86.56 0.05 2.94 59791.
## 396 7211 85.95 3.90 78.31 93.59 0.05 1.26 21134.
## 397 7212 84.58 2.91 78.87 90.29 0.03 1.10 35683.
## 398 7271 99.25 0.62 98.03 100.47 0.01 2.44 99355.
## 399 7301 95.23 2.42 90.50 99.97 0.03 2.14 34952.
## 400 7302 86.77 2.73 81.41 92.12 0.03 3.32 107824.
## 401 7303 97.02 1.08 94.90 99.14 0.01 0.93 48460.
## 402 7304 97.43 1.16 95.16 99.69 0.01 2.23 87937.
## 403 7305 98.34 0.70 96.97 99.72 0.01 1.02 71277.
## 404 7306 94.88 1.92 91.11 98.65 0.02 6.93 192077.
## 405 7307 90.69 2.04 86.70 94.68 0.02 1.32 56510.
## 406 7308 91.34 2.75 85.95 96.74 0.03 3.76 82711.
## 407 7309 95.07 1.40 92.33 97.81 0.01 1.56 78955.
## 408 7310 98.26 0.71 96.88 99.65 0.01 0.60 43366.
## 409 7311 91.80 2.25 87.40 96.20 0.02 5.50 173044.
## 410 7312 95.09 3.31 88.60 101.58 0.03 6.70 60051.
## 411 7313 91.39 3.25 85.01 97.77 0.04 6.43 100571.
## 412 7314 96.78 1.41 94.03 99.54 0.01 2.25 74512.
## 413 7315 90.47 2.85 84.89 96.05 0.03 4.10 91953.
## 414 7316 83.34 3.83 75.84 90.85 0.05 2.50 49837.
## 415 7317 82.05 3.13 75.91 88.18 0.04 2.64 83475.
## 416 7318 66.50 4.41 57.86 75.15 0.07 2.40 57981.
## 417 7322 88.12 1.97 84.27 91.98 0.02 1.37 78026.
## 418 7325 88.27 3.21 81.98 94.56 0.04 3.51 74305.
## 419 7326 70.34 4.02 62.46 78.22 0.06 1.97 53529.
## 420 7371 98.57 0.53 97.53 99.61 0.01 3.53 370665.
## 421 7372 99.53 0.26 99.02 100.04 0.00 0.24 34296.
## 422 7373 99.26 0.38 98.52 100.00 0.00 0.40 43546.
## 423 7401 97.21 1.08 95.09 99.33 0.01 0.45 21908.
## 424 7402 97.66 1.06 95.58 99.75 0.01 1.24 52651.
## 425 7403 86.53 3.42 79.83 93.22 0.04 2.82 59247.
## 426 7404 93.67 1.99 89.76 97.57 0.02 2.07 65096.
## 427 7405 91.82 2.33 87.24 96.39 0.03 2.73 79472.
## 428 7406 96.88 1.67 93.61 100.15 0.02 2.04 46645.
## 429 7407 99.19 0.52 98.17 100.20 0.01 0.40 25556.
## 430 7408 94.08 1.80 90.56 97.60 0.02 1.05 38284.
## 431 7409 91.92 1.93 88.13 95.70 0.02 0.34 14403.
## 432 7410 85.87 2.56 80.86 90.88 0.03 0.38 14748.
## 433 7411 96.35 1.73 92.96 99.75 0.02 1.45 35815.
## 434 7412 96.42 0.95 94.56 98.27 0.01 0.10 8356.
## 435 7413 95.52 1.66 92.27 98.77 0.02 0.60 19768.
## 436 7414 97.67 1.27 95.17 100.17 0.01 0.69 20524.
## 437 7415 99.46 0.35 98.77 100.15 0.00 0.18 16600.
## 438 7471 99.70 0.27 99.17 100.23 0.00 1.18 102677.
## 439 7472 100.00 0.00 99.99 100.00 0.00 0.01 38924.
## 440 7501 95.37 1.15 93.12 97.63 0.01 0.62 43481.
## 441 7502 96.20 1.61 93.05 99.35 0.02 2.94 87718.
## 442 7503 96.20 1.13 93.99 98.41 0.01 0.66 40028.
## 443 7504 96.05 1.92 92.28 99.83 0.02 1.75 37780.
## 444 7505 90.71 2.58 85.65 95.76 0.03 1.07 28590.
## 445 7571 99.37 0.42 98.55 100.18 0.00 0.73 55939.
## 446 7601 98.50 0.79 96.94 100.05 0.01 0.75 37299.
## 447 7602 83.60 3.36 77.01 90.18 0.04 3.94 100920.
## 448 7603 36.70 5.48 25.95 47.45 0.15 2.41 39241.
## 449 7604 84.36 2.99 78.50 90.22 0.04 2.30 71613.
## 450 7605 87.28 2.37 82.64 91.92 0.03 1.08 45230.
## 451 7606 81.48 4.10 73.44 89.53 0.05 1.80 33902.
## 452 8101 92.89 2.46 88.08 97.71 0.03 1.01 23266.
## 453 8102 93.48 2.49 88.60 98.36 0.03 0.98 20302.
## 454 8103 93.40 2.54 88.42 98.39 0.03 4.14 83177.
## 455 8104 91.55 2.42 86.81 96.28 0.03 1.24 34571.
## 456 8105 83.85 3.28 77.42 90.28 0.04 0.72 19012.
## 457 8106 93.58 1.92 89.82 97.35 0.02 1.03 35433.
## 458 8107 80.48 3.85 72.92 88.04 0.05 1.03 23024.
## 459 8108 96.14 1.84 92.54 99.74 0.02 0.68 15730.
## 460 8109 84.81 3.45 78.04 91.57 0.04 0.55 12636.
## 461 8171 98.22 0.82 96.61 99.83 0.01 2.07 113107.
## 462 8172 95.56 1.28 93.04 98.07 0.01 0.29 15796.
## 463 8201 93.86 2.03 89.89 97.84 0.02 0.93 27549.
## 464 8202 89.34 2.97 83.51 95.16 0.03 0.56 12774.
## 465 8203 85.20 3.32 78.70 91.70 0.04 0.96 23096.
## 466 8204 84.34 3.56 77.36 91.32 0.04 2.25 49458.
## 467 8205 89.50 2.94 83.74 95.27 0.03 1.87 42860.
## 468 8206 71.17 4.43 62.49 79.85 0.06 1.06 23309.
## 469 8207 87.62 2.48 82.77 92.48 0.03 0.41 15094.
## 470 8208 85.77 3.19 79.51 92.03 0.04 0.52 13165.
## 471 8271 97.91 0.93 96.10 99.72 0.01 1.14 57496.
## 472 8272 93.42 2.02 89.45 97.39 0.02 0.72 22799.
## 473 9101 95.43 1.54 92.41 98.44 0.02 0.49 18964.
## 474 9102 85.15 2.81 79.64 90.66 0.03 0.46 15544.
## 475 9103 30.49 4.19 22.27 38.71 0.14 0.31 7888.
## 476 9104 89.90 2.53 84.93 94.87 0.03 0.54 16061.
## 477 9105 84.40 2.97 78.58 90.23 0.04 1.36 42721.
## 478 9111 83.12 2.83 77.58 88.66 0.03 0.18 6482.
## 479 9112 72.50 4.06 64.55 80.46 0.06 0.37 9417.
## 480 9201 88.48 2.04 84.49 92.47 0.02 0.22 11498.
## 481 9202 88.89 2.60 83.80 93.98 0.03 0.75 23122.
## 482 9203 57.76 5.17 47.62 67.90 0.09 0.54 10430.
## 483 9204 45.05 4.80 35.66 54.45 0.11 0.48 10894.
## 484 9205 57.96 7.91 42.45 73.47 0.14 0.45 3649.
## 485 9271 88.99 2.67 83.76 94.22 0.03 2.31 67055.
## 486 9403 81.58 3.62 74.49 88.67 0.04 1.32 31926.
## 487 9408 76.77 4.42 68.11 85.43 0.06 1.22 23568.
## 488 9409 87.17 3.21 80.88 93.47 0.04 1.53 35050.
## 489 9419 90.28 2.95 84.51 96.06 0.03 0.52 11103.
## 490 9420 60.69 3.47 53.88 67.50 0.06 0.39 16108.
## 491 9426 93.47 1.81 89.92 97.03 0.02 0.21 8066.
## 492 9427 68.84 2.78 63.40 74.29 0.04 0.08 4460.
## 493 9428 62.23 7.22 48.09 76.38 0.12 0.70 6664.
## 494 9471 96.86 1.07 94.75 98.97 0.01 1.36 75353.
## 495 9501 75.52 3.63 68.40 82.63 0.05 1.99 58882.
## 496 9502 67.26 3.89 59.64 74.89 0.06 0.61 18723.
## 497 9503 38.55 5.13 28.51 48.60 0.13 1.22 23121.
## 498 9504 98.53 0.63 97.29 99.77 0.01 0.32 24214.
## 499 9601 95.88 1.63 92.67 99.08 0.02 1.81 56362.
## 500 9602 70.66 4.64 61.56 79.76 0.07 1.32 26660.
## 501 9603 89.60 3.10 83.52 95.68 0.03 1.07 21853.
## 502 9604 87.49 3.14 81.33 93.64 0.04 1.58 36953.
## 503 9605 90.45 4.00 82.62 98.28 0.04 4.06 46233.
## 504 9606 90.15 2.80 84.67 95.63 0.03 0.58 13987.
## 505 9607 30.09 6.78 16.80 43.37 0.23 3.55 34309.
## 506 9608 88.90 3.69 81.66 96.14 0.04 2.81 42855.
## 507 9701 16.73 3.99 8.91 24.55 0.24 1.48 27341.
## 508 9702 47.45 5.86 35.97 58.94 0.12 3.89 59493.
## 509 9703 1.08 0.53 0.04 2.13 0.49 0.74 59071.
## 510 9704 18.37 4.54 9.48 27.25 0.25 2.54 39084.
## 511 9705 1.83 0.64 0.58 3.09 0.35 0.16 14354.
## 512 9706 44.58 4.33 36.10 53.05 0.10 0.64 17794.
## 513 9707 46.74 6.49 34.03 59.46 0.14 4.17 52006.
## 514 9708 78.31 3.99 70.49 86.14 0.05 0.87 19604.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. (2024). Konsep dan Definisi Susenas Maret 2024. 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
[3] Sekretariat SDGs Bappenas. 2024. Metadata Indikator Tujuan Pembangunan Berkelanjutan (TPB)/ Sustainable Development Goals (SDGs) Indonesia. Jakarta: Kementerian Perencanaan Pembangunan Nasional/Badan Perencanaan Pembangunan Nasional.
Direktorat Statistik Kesejahteraan Rakyat, BPS, saptahas@bps.go.id