La librería sampling se puede utilizar para extraer muestras especificando el diseño. Si lo que se quiere es analizar una muestra con diseños complejos se utiliza la librería survey
Cargamos la librería sampling
if (!require(sampling)){
install.packages("sampling",dep=TRUE)
library(sampling)
}## Loading required package: samplingLa función es mstage (Multistage sampling), pero tiene un bug en la versión 2.6. Los autores me han enviado por email la función corregida mstage.R
Bajamos y cargamos la función
temp <- tempfile()
# uso método wget porque curl no me funciona
download.file(url = "https://www.dropbox.com/s/lzdootdhinjumkz/mstage.r?dl=0",
destfile = temp, method = "wget")
source(temp)Creamos unos datos artificiales
data=rbind(matrix(rep('n',165),165,1,byrow=TRUE),matrix(rep('s',70),70,1,byrow=TRUE))
data=cbind.data.frame(data,c(rep('A',115),rep('D',10),rep('E',40),rep('B',30),rep('C',40)),
100*runif(235))
names(data)=c("state","region","income")
head(data)## state region income
## 1 n A 99.521363
## 2 n A 1.801964
## 3 n A 35.176919
## 4 n A 93.349948
## 5 n A 56.204249
## 6 n A 19.615237Ordenamos de forma ascendente por las variables que vamos a tomar como estrato y como cluster (según pone en la ayuda de la función mstage)
data=data[order(data$state,data$region),]
table(data$state,data$region)##
## A B C D E
## n 115 0 0 10 40
## s 0 30 40 0 0La variable state tiene dos categorías y la variable region 5. La variable income ha sido generada por una uniforme (0,1)
Consideramos la variable state como variable de estratificación y la variable region como variable de conglomerados.
En el estrato n hay 165 unidades y en el s 70. Supongamos que queremos hacer muestreo estraficado y elegimos 2 conglomerados en el estrato n y 1 en el estrato s.
Hay que decirle en cada etapa (stage) el tipo de muestreo y la variable que se utiliza. En la opción size se especifica una lista con los tamaños, en este caso se pone la población total en cada estrato 165 y 70 y el número de conglomerados que seleccionamos en cada estrato (2 y 1). En method le indicamos como seleccionar las unidades
Aunque en este ejemplo no se usa también se pueden especificar probabilidades de selección desiguales de las unidades de cada etapa (argumento pik). Este argumento sólo se utiliza cuando en el anterior se utiliza muestreo de poisson o sistemático. Ver ejemplos en la función mstage
m <- mstage1(data, stage = list("stratified", "cluster"), varnames = list("state","region"),
size=list(c(165,70),c(2,1)),method="srswor", description=TRUE) ## STAGE 1
## STAGE 2
## Number of selected clusters: 2
## Population total and number of selected units 165 155
## Number of selected clusters: 1
## Population total and number of selected units 70 40La función mstage1 devuelve como valores una lista con tantos elementos como etapas. Da información sobre la unidad seleccionada, la probabilidad de inclusión en la etapa i y la probabilidad de inclusión en la última etapa
head(m[[1]])## state ID_unit Prob_ 1 _stage Stratum
## 1 n 1 1 1
## 2 n 2 1 1
## 3 n 3 1 1
## 4 n 4 1 1
## 5 n 5 1 1
## 6 n 6 1 1head(m[[2]])## region ID_unit Prob_ 2 _stage Prob
## 1 A 7 0.6666667 0.6666667
## 2 A 8 0.6666667 0.6666667
## 3 A 9 0.6666667 0.6666667
## 4 A 6 0.6666667 0.6666667
## 5 A 15 0.6666667 0.6666667
## 6 A 16 0.6666667 0.6666667Para obtener los datos originales según este muestreo podemos utilizar la función getdata, pasándole como argumento los datos originales y el objeto que nos devuelve mstage. getdatadevuelve una lista con las unidades seleccionadas en cada etapa, las unidades que se han seleccionado al final del proceso están en el último elemento de la lista, en este caso en el segundo.
xx=getdata(data,m)[[2]]
xx## state income region ID_unit Prob_ 2 _stage Prob
## 7 n 88.4715100 A 7 0.6666667 0.6666667
## 8 n 36.9091778 A 8 0.6666667 0.6666667
## 9 n 59.5332618 A 9 0.6666667 0.6666667
## 6 n 19.6152370 A 6 0.6666667 0.6666667
## 15 n 2.5127308 A 15 0.6666667 0.6666667
## 16 n 18.8063125 A 16 0.6666667 0.6666667
## 17 n 62.8124933 A 17 0.6666667 0.6666667
## 18 n 47.4494883 A 18 0.6666667 0.6666667
## 19 n 57.5449846 A 19 0.6666667 0.6666667
## 20 n 36.4627149 A 20 0.6666667 0.6666667
## 21 n 31.9632278 A 21 0.6666667 0.6666667
## 22 n 5.5820809 A 22 0.6666667 0.6666667
## 10 n 81.2333829 A 10 0.6666667 0.6666667
## 11 n 74.5630338 A 11 0.6666667 0.6666667
## 12 n 0.9358263 A 12 0.6666667 0.6666667
## 13 n 32.7977540 A 13 0.6666667 0.6666667
## 14 n 74.2639278 A 14 0.6666667 0.6666667
## 28 n 66.4795040 A 28 0.6666667 0.6666667
## 29 n 42.0094777 A 29 0.6666667 0.6666667
## 30 n 97.7673650 A 30 0.6666667 0.6666667
## 31 n 71.9226850 A 31 0.6666667 0.6666667
## 32 n 3.6643676 A 32 0.6666667 0.6666667
## 33 n 42.3470949 A 33 0.6666667 0.6666667
## 34 n 38.0503924 A 34 0.6666667 0.6666667
## 35 n 70.5297326 A 35 0.6666667 0.6666667
## 23 n 30.9508602 A 23 0.6666667 0.6666667
## 24 n 86.4545714 A 24 0.6666667 0.6666667
## 25 n 92.6925888 A 25 0.6666667 0.6666667
## 26 n 44.7693517 A 26 0.6666667 0.6666667
## 27 n 27.5694360 A 27 0.6666667 0.6666667
## 1 n 99.5213634 A 1 0.6666667 0.6666667
## 2 n 1.8019645 A 2 0.6666667 0.6666667
## 3 n 35.1769191 A 3 0.6666667 0.6666667
## 4 n 93.3499481 A 4 0.6666667 0.6666667
## 5 n 56.2042489 A 5 0.6666667 0.6666667
## 46 n 7.9515596 A 46 0.6666667 0.6666667
## 47 n 81.0193568 A 47 0.6666667 0.6666667
## 48 n 60.0569664 A 48 0.6666667 0.6666667
## 49 n 81.6952142 A 49 0.6666667 0.6666667
## 50 n 51.7822959 A 50 0.6666667 0.6666667
## 51 n 71.6329366 A 51 0.6666667 0.6666667
## 52 n 80.6264753 A 52 0.6666667 0.6666667
## 53 n 0.8064215 A 53 0.6666667 0.6666667
## 54 n 73.0552500 A 54 0.6666667 0.6666667
## 55 n 38.3441581 A 55 0.6666667 0.6666667
## 56 n 75.2735012 A 56 0.6666667 0.6666667
## 57 n 19.0246811 A 57 0.6666667 0.6666667
## 58 n 67.5058935 A 58 0.6666667 0.6666667
## 59 n 62.2341210 A 59 0.6666667 0.6666667
## 60 n 99.1942390 A 60 0.6666667 0.6666667
## 61 n 85.1055545 A 61 0.6666667 0.6666667
## 62 n 35.6849157 A 62 0.6666667 0.6666667
## 63 n 88.2010156 A 63 0.6666667 0.6666667
## 64 n 6.4637275 A 64 0.6666667 0.6666667
## 65 n 34.8643030 A 65 0.6666667 0.6666667
## 66 n 46.6413054 A 66 0.6666667 0.6666667
## 67 n 80.3911124 A 67 0.6666667 0.6666667
## 68 n 20.7018284 A 68 0.6666667 0.6666667
## 69 n 71.1877503 A 69 0.6666667 0.6666667
## 70 n 54.5215113 A 70 0.6666667 0.6666667
## 71 n 34.1398335 A 71 0.6666667 0.6666667
## 72 n 10.5526085 A 72 0.6666667 0.6666667
## 73 n 99.3685464 A 73 0.6666667 0.6666667
## 74 n 12.4523544 A 74 0.6666667 0.6666667
## 75 n 66.1307842 A 75 0.6666667 0.6666667
## 36 n 0.7626952 A 36 0.6666667 0.6666667
## 37 n 95.0506879 A 37 0.6666667 0.6666667
## 38 n 36.4105911 A 38 0.6666667 0.6666667
## 39 n 87.3448173 A 39 0.6666667 0.6666667
## 40 n 97.3343901 A 40 0.6666667 0.6666667
## 41 n 50.4162251 A 41 0.6666667 0.6666667
## 42 n 29.2536982 A 42 0.6666667 0.6666667
## 43 n 52.6080823 A 43 0.6666667 0.6666667
## 44 n 21.9438461 A 44 0.6666667 0.6666667
## 45 n 8.1710815 A 45 0.6666667 0.6666667
## 86 n 25.8601717 A 86 0.6666667 0.6666667
## 87 n 53.4591818 A 87 0.6666667 0.6666667
## 88 n 74.1454887 A 88 0.6666667 0.6666667
## 89 n 98.1671982 A 89 0.6666667 0.6666667
## 90 n 61.3992383 A 90 0.6666667 0.6666667
## 91 n 58.2479154 A 91 0.6666667 0.6666667
## 92 n 7.7950890 A 92 0.6666667 0.6666667
## 93 n 5.3062318 A 93 0.6666667 0.6666667
## 94 n 57.0284058 A 94 0.6666667 0.6666667
## 95 n 50.2401471 A 95 0.6666667 0.6666667
## 96 n 78.9598142 A 96 0.6666667 0.6666667
## 97 n 15.4643004 A 97 0.6666667 0.6666667
## 98 n 84.6005336 A 98 0.6666667 0.6666667
## 99 n 3.1739886 A 99 0.6666667 0.6666667
## 100 n 43.9629626 A 100 0.6666667 0.6666667
## 101 n 21.6785324 A 101 0.6666667 0.6666667
## 102 n 88.6651628 A 102 0.6666667 0.6666667
## 103 n 72.6373881 A 103 0.6666667 0.6666667
## 104 n 80.0055564 A 104 0.6666667 0.6666667
## 105 n 60.1750817 A 105 0.6666667 0.6666667
## 106 n 47.8393794 A 106 0.6666667 0.6666667
## 107 n 0.2397272 A 107 0.6666667 0.6666667
## 108 n 76.2090890 A 108 0.6666667 0.6666667
## 109 n 53.5498631 A 109 0.6666667 0.6666667
## 110 n 5.3196671 A 110 0.6666667 0.6666667
## 111 n 48.7121697 A 111 0.6666667 0.6666667
## 112 n 2.9144365 A 112 0.6666667 0.6666667
## 113 n 8.1327649 A 113 0.6666667 0.6666667
## 114 n 52.5297917 A 114 0.6666667 0.6666667
## 115 n 7.4287212 A 115 0.6666667 0.6666667
## 76 n 29.8141962 A 76 0.6666667 0.6666667
## 77 n 35.5496594 A 77 0.6666667 0.6666667
## 78 n 50.7121051 A 78 0.6666667 0.6666667
## 79 n 14.3348420 A 79 0.6666667 0.6666667
## 80 n 58.9701313 A 80 0.6666667 0.6666667
## 81 n 37.1326326 A 81 0.6666667 0.6666667
## 82 n 90.8194170 A 82 0.6666667 0.6666667
## 83 n 91.7900737 A 83 0.6666667 0.6666667
## 84 n 13.2420932 A 84 0.6666667 0.6666667
## 85 n 76.0856323 A 85 0.6666667 0.6666667
## 126 n 44.8293304 E 126 0.6666667 0.6666667
## 127 n 92.1129314 E 127 0.6666667 0.6666667
## 128 n 90.4778450 E 128 0.6666667 0.6666667
## 129 n 88.4983810 E 129 0.6666667 0.6666667
## 130 n 16.4212147 E 130 0.6666667 0.6666667
## 131 n 84.4840840 E 131 0.6666667 0.6666667
## 132 n 48.4419205 E 132 0.6666667 0.6666667
## 133 n 48.4527229 E 133 0.6666667 0.6666667
## 134 n 3.7589894 E 134 0.6666667 0.6666667
## 135 n 66.3089133 E 135 0.6666667 0.6666667
## 136 n 62.2732940 E 136 0.6666667 0.6666667
## 137 n 81.8737664 E 137 0.6666667 0.6666667
## 138 n 75.0281448 E 138 0.6666667 0.6666667
## 139 n 4.9235689 E 139 0.6666667 0.6666667
## 140 n 92.0513147 E 140 0.6666667 0.6666667
## 141 n 28.9218198 E 141 0.6666667 0.6666667
## 142 n 1.2979860 E 142 0.6666667 0.6666667
## 143 n 27.2616479 E 143 0.6666667 0.6666667
## 144 n 76.6855516 E 144 0.6666667 0.6666667
## 145 n 46.3674579 E 145 0.6666667 0.6666667
## 146 n 11.2265446 E 146 0.6666667 0.6666667
## 147 n 97.5223317 E 147 0.6666667 0.6666667
## 148 n 11.0521420 E 148 0.6666667 0.6666667
## 149 n 15.0415194 E 149 0.6666667 0.6666667
## 150 n 29.3264657 E 150 0.6666667 0.6666667
## 151 n 71.6697691 E 151 0.6666667 0.6666667
## 152 n 59.9324764 E 152 0.6666667 0.6666667
## 153 n 3.0451548 E 153 0.6666667 0.6666667
## 154 n 96.4667542 E 154 0.6666667 0.6666667
## 155 n 76.6638573 E 155 0.6666667 0.6666667
## 156 n 84.4972424 E 156 0.6666667 0.6666667
## 157 n 33.1655099 E 157 0.6666667 0.6666667
## 158 n 64.1573716 E 158 0.6666667 0.6666667
## 159 n 25.4832905 E 159 0.6666667 0.6666667
## 160 n 99.7938742 E 160 0.6666667 0.6666667
## 161 n 15.8569932 E 161 0.6666667 0.6666667
## 162 n 28.6158246 E 162 0.6666667 0.6666667
## 163 n 22.3719261 E 163 0.6666667 0.6666667
## 164 n 52.2040672 E 164 0.6666667 0.6666667
## 165 n 71.4342273 E 165 0.6666667 0.6666667
## 196 s 1.3032632 C 196 0.5000000 0.5000000
## 197 s 2.7865850 C 197 0.5000000 0.5000000
## 198 s 7.1243291 C 198 0.5000000 0.5000000
## 199 s 70.4815151 C 199 0.5000000 0.5000000
## 200 s 94.8395711 C 200 0.5000000 0.5000000
## 201 s 85.6252347 C 201 0.5000000 0.5000000
## 202 s 50.1201008 C 202 0.5000000 0.5000000
## 203 s 19.3657070 C 203 0.5000000 0.5000000
## 204 s 81.6607402 C 204 0.5000000 0.5000000
## 205 s 76.5675493 C 205 0.5000000 0.5000000
## 206 s 37.8898578 C 206 0.5000000 0.5000000
## 207 s 9.6171677 C 207 0.5000000 0.5000000
## 208 s 50.1739333 C 208 0.5000000 0.5000000
## 209 s 17.4049054 C 209 0.5000000 0.5000000
## 210 s 66.7061097 C 210 0.5000000 0.5000000
## 211 s 20.2554295 C 211 0.5000000 0.5000000
## 212 s 37.7490880 C 212 0.5000000 0.5000000
## 213 s 24.7867593 C 213 0.5000000 0.5000000
## 214 s 28.3039456 C 214 0.5000000 0.5000000
## 215 s 41.6149533 C 215 0.5000000 0.5000000
## 216 s 70.0464154 C 216 0.5000000 0.5000000
## 217 s 15.3257101 C 217 0.5000000 0.5000000
## 218 s 21.3272467 C 218 0.5000000 0.5000000
## 219 s 45.4484878 C 219 0.5000000 0.5000000
## 220 s 56.8538516 C 220 0.5000000 0.5000000
## 221 s 94.6382512 C 221 0.5000000 0.5000000
## 222 s 3.0915002 C 222 0.5000000 0.5000000
## 223 s 56.6517841 C 223 0.5000000 0.5000000
## 224 s 59.5928363 C 224 0.5000000 0.5000000
## 225 s 45.1554216 C 225 0.5000000 0.5000000
## 226 s 13.4272819 C 226 0.5000000 0.5000000
## 227 s 82.0646290 C 227 0.5000000 0.5000000
## 228 s 34.6941799 C 228 0.5000000 0.5000000
## 229 s 33.4355303 C 229 0.5000000 0.5000000
## 230 s 57.4564283 C 230 0.5000000 0.5000000
## 231 s 40.5542760 C 231 0.5000000 0.5000000
## 232 s 39.7777576 C 232 0.5000000 0.5000000
## 233 s 61.6708402 C 233 0.5000000 0.5000000
## 234 s 77.5421731 C 234 0.5000000 0.5000000
## 235 s 13.6150021 C 235 0.5000000 0.5000000Verificamos el resultado viendo una tabla
table(xx$state,xx$region)##
## A B C D E
## n 115 0 0 0 40
## s 0 0 40 0 0data(swissmunicipalities)
datos <- swissmunicipalities
datos <- datos[order(datos$REG),]
table(datos$REG)##
## 1 2 3 4 5 6 7
## 589 913 321 171 471 186 245Si queremos extraer una muestra de tamaños 30,20,45,15,20,11,44 en cada estrato por muestreo sin reemplazamiento, utilizamos la función strata
st <- strata(datos,stratanames = c("REG"), size = c(30,20,45,15,20,11,44),
method = "srswor")
# extracts the observed data
getdata(datos, st)## CT COM Nom HApoly Surfacesbois Surfacescult Alp
## 31 25 6630 Meyrin 996 66 392 0
## 187 22 5822 Payerne 2419 353 1613 0
## 211 23 6297 Visp 1317 788 224 0
## 276 23 6265 Saviese 7093 1078 663 909
## 502 22 5764 Vallorbe 2318 1568 404 114
## 554 23 6159 Vouvry 3354 1205 532 666
## 623 23 6232 Chalais 2446 1660 268 188
## 867 25 6606 Bellevue 440 85 216 0
## 952 23 6219 Vernayaz 561 310 73 0
## 1225 23 6293 Stalden (VS) 1043 753 96 9
## 1359 22 5415 Yvorne 1219 491 588 0
## 1733 23 6212 Dorenaz 1257 758 132 173
## 1810 22 5489 Mex (VD) 290 88 157 0
## 1967 23 6251 Vissoie 150 105 10 0
## 1993 22 5761 Romainmotier-Envy 701 482 192 0
## 2048 22 5654 Vullierens 682 60 583 0
## 2117 22 5815 Combremont-le-Petit 574 128 409 0
## 2119 23 6197 Kippel 1166 354 77 239
## 2234 23 6214 Finhaut 2288 600 33 41
## 2258 22 5751 Corcelles-sur-Chavornay 549 162 356 0
## 2294 22 5673 Hermenches 476 136 309 0
## 2346 23 6102 Albinen 1548 571 68 657
## 2347 23 6066 Oberwald 9692 755 152 1820
## 2523 22 5830 Villarzel 406 110 280 0
## 2563 22 5460 Mur (VD) 178 16 140 0
## 2646 22 5672 Forel-sur-Lucens 284 54 214 0
## 2679 22 5475 Chavannes-le-Veyron 258 70 176 0
## 2706 22 5683 Prevonloup 184 24 146 0
## 2798 22 5797 La Rogivue 106 20 69 0
## 2889 23 6172 Bister 584 233 64 43
## 596 2 886 Wattenwil 1454 702 611 0
## 633 11 2532 Subingen 625 190 325 0
## 711 10 2197 Givisiez 346 80 137 0
## 740 11 2494 Niedererlinsbach 585 226 249 0
## 922 2 731 Aegerten 217 62 85 0
## 940 11 2403 Kestenholz 858 327 447 0
## 1012 2 405 Ersigen 872 278 512 0
## 1083 2 971 Attiswil 764 275 386 18
## 1307 2 403 Bariswil 273 116 128 0
## 1478 2 357 Oberbalm 1242 262 922 0
## 1513 2 842 Lauenen 5872 1439 532 1736
## 1667 10 2171 Arconciel 617 177 358 0
## 1881 10 2176 Bonnefontaine 499 272 162 34
## 2080 26 6718 Movelier 807 401 271 98
## 2250 2 926 Forst 184 22 149 0
## 2335 11 2491 Hauenstein-Ifenthal 535 225 269 0
## 2342 10 2185 Corserey 344 66 262 0
## 2408 10 2062 Berlens 300 46 239 0
## 2512 2 325 Busswil bei Melchnau 287 121 148 0
## 2527 10 2110 Villarsiviriaux 416 186 185 35
## 30 12 2703 Riehen 1086 276 283 0
## 36 13 2762 Allschwil 892 242 273 0
## 78 19 4082 Wohlen (AG) 1251 349 494 0
## 149 13 2763 Arlesheim 694 372 76 0
## 188 19 4141 Reinach (AG) 948 257 442 0
## 200 19 4282 Rothrist 1184 387 445 0
## 265 19 4139 Menziken 639 152 316 0
## 464 19 4092 Birr 502 166 212 0
## 467 19 4083 Zufikon 480 130 239 0
## 511 19 4013 Unterentfelden 288 88 93 0
## 519 19 4283 Safenwil 596 277 181 0
## 539 19 4026 Ennetbaden 211 90 44 0
## 620 13 2764 Biel-Benken 412 83 241 0
## 713 19 4312 Lengnau (AG) 1264 482 680 0
## 720 19 4023 Bergdietikon 593 164 328 0
## 722 19 4140 Oberkulm 943 348 496 0
## 732 19 4210 Staufen 350 113 140 0
## 738 19 4313 Leuggern 1376 519 651 0
## 809 19 4068 Hagglingen 774 255 441 0
## 859 19 4104 Lupfig 515 125 265 0
## 988 19 4261 Wallbach 452 136 204 0
## 1073 13 2834 Ziefen 780 375 344 0
## 1107 19 4061 Arni (AG) 336 94 199 0
## 1142 19 4177 Sisseln 252 46 93 0
## 1154 19 4147 Zetzwil 580 223 307 0
## 1204 19 4007 Hirschthal 354 182 121 0
## 1387 13 2888 Langenbruck 1569 721 505 286
## 1388 19 4164 Gansingen 876 328 483 0
## 1522 19 4182 Wolflinswil 951 304 578 0
## 1590 19 4004 Densburen 1251 701 496 0
## 1644 13 2842 Bockten 229 99 100 0
## 1689 13 2832 Ramlinsburg 224 110 93 0
## 1732 19 4191 Ammerswil 318 185 108 0
## 1787 13 2885 Eptingen 1119 539 452 72
## 1856 19 4225 Benzenschwil 245 69 149 0
## 1897 13 2881 Arboldswil 348 149 174 0
## 1996 19 4308 Kaiserstuhl 32 4 12 0
## 2170 19 4322 Wislikofen 374 135 200 0
## 2223 13 2859 Rumlingen 228 94 113 0
## 2255 19 4174 Oberhofen (AG) 314 144 147 0
## 2286 13 2848 Hemmiken 339 83 237 0
## 2355 19 4314 Mellikon 271 143 74 0
## 2400 19 4097 Elfingen 423 207 202 0
## 2493 13 2854 Nusshof 172 70 98 0
## 2537 19 4302 Bobikon 262 99 145 0
## 41 1 62 Kloten 1922 510 658 0
## 57 1 131 Adliswil 778 248 182 0
## 156 1 155 Mannedorf 477 118 163 0
## 171 1 90 Niederhasli 1129 239 659 0
## 181 1 192 Egg 1449 334 842 0
## 248 1 245 Oberengstringen 216 56 50 0
## 286 1 159 Uetikon am See 345 59 165 0
## 349 1 64 Nurensdorf 1007 325 490 0
## 375 1 10 Obfelden 754 169 410 0
## 463 1 58 Glattfelden 1235 397 582 0
## 659 1 224 Pfungen 498 210 160 0
## 730 1 91 Niederweningen 687 260 334 0
## 1368 1 241 Aesch bei Birmensdorf 524 174 303 0
## 1760 1 21 Adlikon 661 151 438 0
## 1914 1 95 Regensberg 238 136 82 0
## 399 20 4781 Wangi 1641 269 1161 0
## 533 20 4501 Kradolf-Schonenberg 1095 259 675 0
## 606 18 3945 Trimmis 2840 1181 394 584
## 649 18 3732 Flims 5048 1457 454 1248
## 756 20 4711 Affeltrangen 1444 226 1063 0
## 773 20 4921 Bussnang 1889 361 1336 0
## 804 14 2971 Hallau 1532 618 774 0
## 906 15 3007 Waldstatt 674 174 417 0
## 1020 17 3375 Oberhelfenschwil 1266 427 784 1
## 1058 16 3102 Gonten 2471 791 1230 295
## 1126 20 4776 Tobel-Tagerschen 706 205 416 0
## 1181 18 3961 Grusch 1002 514 290 105
## 1185 18 3822 Mesocco 16477 3306 370 2936
## 1365 20 4811 Herdern 1367 383 922 0
## 1439 18 3640 Sils im Domleschg 932 639 80 59
## 1815 8 1608 Filzbach 1388 570 160 333
## 1835 14 2973 Trasadingen 415 52 313 0
## 2182 14 2961 Buch (SH) 379 75 286 0
## 2291 18 3592 Cumbel 446 203 174 33
## 2453 18 3805 Castaneda 396 302 15 0
## 238 5 1323 Wollerau 632 88 361 0
## 347 9 1704 Menzingen 2753 778 1719 0
## 577 3 1136 Menznau 3038 1056 1777 23
## 914 3 1088 Hildisrieden 705 83 545 0
## 1362 5 1366 Morschach 2081 752 367 721
## 1593 5 1368 Oberiberg 3294 1284 490 1158
## 1627 3 1096 Pfeffikon 248 126 96 0
## 1818 4 1211 Isenthal 6099 1598 350 1352
## 1977 3 1028 Hamikon 466 126 310 0
## 2103 3 1124 Buchs (LU) 467 98 327 0
## 2582 5 1343 Innerthal 5016 1709 173 1576
## 21 21 5192 Lugano 1170 673 22 11
## 47 21 5002 Bellinzona 1915 1196 173 29
## 167 21 5250 Chiasso 533 265 65 0
## 215 21 5118 Minusio 585 311 33 30
## 245 21 5115 Losone 953 619 83 7
## 258 21 5196 Massagno 73 1 0 0
## 441 21 5210 Paradiso 89 32 3 0
## 474 21 5242 Balerna 257 16 36 0
## 771 21 5201 Montagnola 314 166 29 0
## 857 21 5167 Canobbio 130 49 9 0
## 949 21 5224 Sonvico 1106 886 83 22
## 972 21 5147 Barbengo 266 139 25 0
## 976 21 5072 Faido 372 236 48 0
## 1013 21 5212 Ponte Capriasca 620 397 12 77
## 1085 21 5185 Gentilino 122 42 16 0
## 1192 21 5148 Bedano 187 100 29 2
## 1268 21 5112 Lavertezzo 5811 2934 42 403
## 1319 21 5241 Arzo 279 201 35 0
## 1396 21 5111 Intragna 2405 2078 60 111
## 1473 21 5043 Olivone 7609 2641 123 2107
## 1484 21 5252 Genestrerio 149 32 78 0
## 1571 21 5205 Muzzano 157 37 30 0
## 1587 21 5077 Pollegio 589 395 78 6
## 1800 21 5229 Valcolla 1134 589 30 350
## 1853 21 5181 Curio 287 235 26 0
## 1890 21 5190 Lopagno 537 322 11 120
## 1922 21 5015 Preonzo 1644 993 105 42
## 1963 21 5160 Brusino Arsizio 404 362 9 1
## 2128 21 5230 Vernate 151 125 7 0
## 2174 21 5246 Caneggio 390 346 15 2
## 2274 21 5150 Bidogno 349 183 12 77
## 2367 21 5322 Someo 3279 2197 22 180
## 2368 21 5114 Loco 907 663 19 82
## 2479 21 5259 Muggio 839 647 59 106
## 2676 21 5218 Roveredo (TI) 69 50 6 7
## 2712 21 5179 Cureggia 69 60 2 0
## 2732 21 5012 Moleno 749 594 46 5
## 2774 21 5068 Cavagnago 668 497 37 60
## 2796 21 5318 Menzonio 1064 584 17 31
## 2797 21 5063 Bedretto 7523 875 133 2169
## 2817 21 5173 Certara 273 222 16 14
## 2828 21 5248 Casima 99 89 5 0
## 2863 21 5067 Campello 396 114 35 121
## 2896 21 5102 Corippo 773 495 6 55
## Airbat Airind P00BMTOT P00BWTOT Pop020 Pop2040 Pop4065 Pop65P H00PTOT
## 31 232 77 9656 9892 4424 5800 7121 2203 7938
## 187 145 42 3475 3819 1724 2089 2261 1220 3271
## 211 70 90 3200 3350 1532 1891 2143 984 2536
## 276 183 1 2604 2737 1219 1519 1812 791 1878
## 502 78 14 1609 1638 787 841 991 628 1385
## 554 75 24 1441 1519 797 780 938 445 1121
## 623 107 8 1320 1337 639 728 868 422 1045
## 867 67 1 881 920 518 439 705 139 561
## 952 38 11 792 805 418 449 481 249 599
## 1225 23 2 541 608 307 288 353 201 432
## 1359 25 6 483 479 212 274 330 146 402
## 1733 21 8 299 308 176 164 179 88 181
## 1810 20 5 274 273 115 153 216 63 240
## 1967 12 1 227 224 127 123 145 56 170
## 1993 9 1 219 216 99 99 149 88 187
## 2048 22 0 197 205 108 111 117 66 163
## 2117 13 2 181 187 107 104 101 56 133
## 2119 5 2 173 195 82 96 116 74 134
## 2234 17 0 164 154 74 80 103 61 130
## 2258 13 0 143 162 77 67 113 48 110
## 2294 15 0 147 137 66 83 89 46 112
## 2346 13 0 121 140 61 64 71 65 107
## 2347 15 1 140 120 59 64 86 51 112
## 2523 9 0 96 94 59 42 57 32 72
## 2563 9 0 86 87 52 41 59 21 67
## 2646 6 0 67 70 44 26 42 25 46
## 2679 7 0 62 63 40 32 44 9 44
## 2706 7 0 53 61 21 21 46 26 51
## 2798 2 0 33 39 18 21 23 10 30
## 2889 2 0 15 18 3 3 16 11 14
## 596 59 5 1365 1405 727 743 811 489 1059
## 633 51 13 1305 1320 699 822 799 305 1046
## 711 39 28 1130 1161 535 905 622 229 828
## 740 50 3 1045 1146 524 575 671 421 890
## 922 37 2 821 842 365 395 670 233 704
## 940 41 7 807 810 421 496 497 203 647
## 1012 50 3 752 731 378 367 489 249 586
## 1083 38 2 665 669 301 328 464 241 567
## 1307 17 1 505 522 268 256 395 108 388
## 1478 34 0 444 391 196 240 267 132 336
## 1513 40 2 390 404 206 214 231 143 286
## 1667 24 0 325 331 206 185 201 64 227
## 1881 18 1 256 244 128 154 150 68 185
## 2080 21 0 195 195 98 94 138 60 157
## 2250 11 0 154 155 96 76 85 52 119
## 2335 12 0 133 132 80 68 85 32 95
## 2342 10 0 136 126 76 75 87 24 87
## 2408 12 0 126 111 84 67 67 19 75
## 2512 5 0 93 102 41 52 63 39 77
## 2527 8 0 103 86 58 48 58 25 72
## 30 302 2 9467 10903 4039 4191 7097 5043 9201
## 36 189 42 8714 9417 3284 4875 6506 3466 8751
## 78 224 48 6592 6737 3303 3862 4284 1880 5479
## 149 145 36 4004 4624 1772 2143 3055 1658 3830
## 188 146 32 3551 3707 1846 2016 2279 1117 2964
## 200 160 49 3430 3439 1618 1839 2381 1031 2748
## 265 92 19 2690 2821 1429 1601 1647 834 2111
## 464 51 30 1850 1679 973 1086 1170 300 1323
## 467 63 9 1757 1731 866 932 1342 348 1393
## 511 61 14 1589 1606 705 888 1041 561 1378
## 519 69 28 1577 1551 751 878 1051 448 1275
## 539 52 0 1516 1503 529 972 1068 450 1397
## 620 60 3 1350 1329 616 649 1047 367 1054
## 713 53 8 1124 1163 641 668 711 267 857
## 720 71 5 1112 1144 436 610 888 322 961
## 722 60 6 1109 1143 615 574 718 345 826
## 732 51 3 1104 1109 485 567 807 354 932
## 738 60 1 1089 1103 591 583 767 251 775
## 809 51 2 971 1019 531 547 626 286 766
## 859 33 28 895 924 390 562 646 221 749
## 988 37 6 758 768 376 403 558 189 623
## 1073 31 4 672 677 397 348 403 201 491
## 1107 30 0 642 660 358 312 536 96 496
## 1142 19 27 659 600 325 399 432 103 521
## 1154 35 1 639 606 315 351 419 160 479
## 1204 27 2 564 613 283 331 404 159 476
## 1387 29 0 455 474 244 248 325 112 362
## 1388 32 0 475 452 262 224 321 120 326
## 1522 32 1 411 373 200 239 255 90 294
## 1590 34 2 362 359 198 169 246 108 266
## 1644 16 2 340 329 158 182 240 89 270
## 1689 15 0 324 317 161 179 236 65 245
## 1732 18 1 307 302 161 197 197 54 234
## 1787 24 5 294 272 145 148 195 78 225
## 1856 13 0 277 241 158 160 149 51 185
## 1897 12 0 252 241 133 126 160 74 192
## 1996 7 0 224 210 99 112 164 59 185
## 2170 20 0 176 168 92 87 127 38 113
## 2223 8 1 162 159 68 97 98 58 130
## 2255 12 1 171 136 85 83 97 42 104
## 2286 6 0 150 140 63 83 95 49 123
## 2355 6 8 131 126 77 47 94 39 96
## 2400 7 0 116 124 61 53 89 37 92
## 2493 4 0 99 102 47 54 76 24 77
## 2537 7 0 97 86 58 45 62 18 61
## 41 239 30 8767 8423 3297 5984 5755 2154 8160
## 57 195 16 7838 7984 3193 4735 5464 2430 7436
## 156 115 6 3948 4400 1658 2366 2840 1484 3715
## 171 98 22 3847 3742 1965 2485 2589 550 3032
## 181 159 7 3592 3863 1775 1958 2806 916 3024
## 248 71 3 2902 2977 1154 1728 1973 1024 2851
## 286 83 8 2438 2772 1155 1426 1751 878 2071
## 349 79 5 2277 2198 1006 1221 1743 505 1866
## 375 78 9 2114 2068 1083 1197 1497 405 1625
## 463 64 11 1751 1793 925 1016 1213 390 1399
## 659 48 13 1278 1208 630 763 807 286 1040
## 730 49 7 1087 1133 600 666 774 180 863
## 1368 30 0 481 469 199 279 365 107 395
## 1760 26 1 301 286 173 159 189 66 204
## 1914 8 0 243 244 116 149 172 50 209
## 399 93 12 2013 1956 1117 1095 1248 509 1448
## 533 66 15 1524 1538 903 858 907 394 1165
## 606 45 24 1350 1374 739 858 902 225 1077
## 649 117 1 1262 1287 541 814 829 365 1139
## 756 68 9 1114 1021 639 582 664 250 789
## 773 85 5 1060 1025 625 545 600 315 708
## 804 63 6 964 1044 529 454 665 360 784
## 906 43 4 858 852 468 440 531 271 636
## 1020 29 2 734 722 492 355 397 212 461
## 1058 47 6 717 662 434 334 370 241 447
## 1126 38 3 627 653 357 389 380 154 476
## 1181 32 6 606 604 353 303 384 170 459
## 1185 72 8 606 595 215 323 387 276 503
## 1365 29 1 496 459 279 234 321 121 305
## 1439 21 11 451 429 218 240 275 147 367
## 1815 25 1 274 268 144 153 182 63 213
## 1835 22 2 262 268 135 118 183 94 226
## 2182 12 0 200 141 81 130 89 41 107
## 2291 8 0 143 145 71 57 61 99 93
## 2453 10 0 103 118 28 58 80 55 90
## 238 94 8 3127 2947 1348 1767 2286 673 2525
## 347 100 3 2085 2410 1108 1363 1270 754 1273
## 577 70 19 1449 1401 872 899 746 333 929
## 914 41 1 845 830 480 426 565 204 594
## 1362 30 0 521 437 274 331 247 106 293
## 1593 33 0 371 347 150 189 260 119 281
## 1627 17 2 335 351 168 163 201 154 289
## 1818 20 2 291 248 185 137 133 84 168
## 1977 15 0 233 211 147 122 141 34 148
## 2103 8 0 185 189 123 110 87 54 123
## 2582 7 0 81 86 46 44 56 21 63
## 21 304 3 11953 14607 4337 7846 8554 5823 13420
## 47 248 30 7701 8762 3170 4709 5432 3152 7294
## 167 98 9 3584 4136 1257 2186 2538 1739 3774
## 215 137 1 2940 3488 1126 1623 2278 1401 3080
## 245 112 19 2900 3007 1194 1828 2068 817 2545
## 258 56 0 2552 3006 993 1682 1834 1049 2685
## 441 29 0 1750 1944 625 1202 1235 632 1762
## 474 65 39 1647 1768 651 956 1104 704 1519
## 771 82 16 981 1111 594 478 696 324 810
## 857 29 9 885 940 336 601 627 261 806
## 949 51 0 774 826 325 444 550 281 646
## 972 44 25 744 815 404 495 497 163 613
## 976 39 4 725 823 316 365 501 366 614
## 1013 39 0 714 764 344 436 534 164 594
## 1085 47 0 644 684 254 345 512 217 576
## 1192 26 12 564 632 245 357 401 193 441
## 1268 28 14 545 553 263 361 326 148 450
## 1319 24 0 478 532 236 280 340 154 409
## 1396 49 0 416 499 173 214 309 219 395
## 1473 70 6 416 429 163 231 256 195 364
## 1484 13 10 400 427 146 283 278 120 360
## 1571 23 13 354 382 154 176 301 105 322
## 1587 18 5 393 330 149 207 249 118 284
## 1800 25 1 248 310 94 140 181 143 245
## 1853 15 0 263 258 128 137 167 89 211
## 1890 20 0 254 242 95 140 193 68 182
## 1922 17 7 237 247 117 137 173 57 195
## 1963 20 0 219 235 97 123 151 83 205
## 2128 16 0 169 194 76 104 128 55 172
## 2174 11 0 170 173 75 88 126 54 147
## 2274 12 0 142 154 42 86 120 48 138
## 2367 8 1 122 132 38 54 89 73 111
## 2368 5 0 116 138 32 55 91 76 107
## 2479 10 0 102 104 27 48 78 53 101
## 2676 2 0 55 71 30 32 43 21 49
## 2712 4 0 56 56 27 27 44 14 45
## 2732 7 0 50 55 23 30 31 21 41
## 2774 10 0 44 39 9 10 32 32 44
## 2796 6 0 31 42 13 22 27 11 34
## 2797 17 0 35 37 4 16 16 36 37
## 2817 3 0 31 34 9 17 20 19 31
## 2828 3 0 28 33 13 16 24 8 29
## 2863 12 0 25 20 8 7 18 12 22
## 2896 1 0 11 11 1 5 8 8 11
## H00P01 H00P02 H00P03 H00P04 POPTOT REG ID_unit Prob Stratum
## 31 2662 2391 1304 1581 19548 1 8 0.05093379 1
## 187 1221 973 472 605 7294 1 35 0.05093379 1
## 211 721 830 377 608 6550 1 39 0.05093379 1
## 276 458 624 334 462 5341 1 50 0.05093379 1
## 502 515 448 165 257 3247 1 70 0.05093379 1
## 554 327 346 172 276 2960 1 80 0.05093379 1
## 623 308 341 164 232 2657 1 89 0.05093379 1
## 867 129 144 99 189 1801 1 128 0.05093379 1
## 952 153 194 75 177 1597 1 141 0.05093379 1
## 1225 109 125 67 131 1149 1 178 0.05093379 1
## 1359 126 132 60 84 962 1 204 0.05093379 1
## 1733 48 40 31 62 607 1 275 0.05093379 1
## 1810 77 85 36 42 547 1 298 0.05093379 1
## 1967 41 47 28 54 451 1 326 0.05093379 1
## 1993 83 55 17 32 435 1 335 0.05093379 1
## 2048 47 53 23 40 402 1 350 0.05093379 1
## 2117 24 51 18 40 368 1 369 0.05093379 1
## 2119 36 39 20 39 368 1 371 0.05093379 1
## 2234 46 37 16 31 318 1 409 0.05093379 1
## 2258 24 38 15 33 305 1 418 0.05093379 1
## 2294 34 35 16 27 284 1 428 0.05093379 1
## 2346 35 37 13 22 261 1 440 0.05093379 1
## 2347 44 33 11 24 260 1 441 0.05093379 1
## 2523 16 28 9 19 190 1 486 0.05093379 1
## 2563 20 19 12 16 173 1 493 0.05093379 1
## 2646 6 20 4 16 137 1 521 0.05093379 1
## 2679 7 16 7 14 125 1 531 0.05093379 1
## 2706 11 28 6 6 114 1 540 0.05093379 1
## 2798 7 13 2 8 72 1 563 0.05093379 1
## 2889 7 3 1 3 33 1 588 0.05093379 1
## 596 354 321 137 247 2770 2 723 0.02190581 2
## 633 272 332 178 264 2625 2 734 0.02190581 2
## 711 241 238 151 198 2291 2 757 0.02190581 2
## 740 264 307 113 206 2191 2 765 0.02190581 2
## 922 191 268 103 142 1663 2 813 0.02190581 2
## 940 177 213 100 157 1617 2 822 0.02190581 2
## 1012 140 224 85 137 1483 2 846 0.02190581 2
## 1083 173 189 90 115 1334 2 871 0.02190581 2
## 1307 67 148 75 98 1027 2 935 0.02190581 2
## 1478 102 112 43 79 835 2 993 0.02190581 2
## 1513 77 83 46 80 794 2 1004 0.02190581 2
## 1667 39 76 32 80 656 2 1052 0.02190581 2
## 1881 41 65 28 51 500 2 1132 0.02190581 2
## 2080 46 54 24 33 390 2 1217 0.02190581 2
## 2250 36 38 14 31 309 2 1276 0.02190581 2
## 2335 27 31 8 29 265 2 1310 0.02190581 2
## 2342 13 23 18 33 262 2 1314 0.02190581 2
## 2408 9 23 12 31 237 2 1347 0.02190581 2
## 2512 17 32 10 18 195 2 1383 0.02190581 2
## 2527 21 22 12 17 189 2 1389 0.02190581 2
## 30 3248 3438 1070 1445 20370 3 1504 0.14018692 3
## 36 3440 3169 998 1144 18131 3 1506 0.14018692 3
## 78 1772 1672 825 1210 13329 3 1513 0.14018692 3
## 149 1354 1370 466 640 8628 3 1524 0.14018692 3
## 188 936 934 424 670 7258 3 1531 0.14018692 3
## 200 778 926 411 633 6869 3 1532 0.14018692 3
## 265 585 687 331 508 5511 3 1541 0.14018692 3
## 464 328 421 204 370 3529 3 1570 0.14018692 3
## 467 366 468 209 350 3488 3 1571 0.14018692 3
## 511 425 508 179 266 3195 3 1579 0.14018692 3
## 519 344 471 175 285 3128 3 1580 0.14018692 3
## 539 531 521 140 205 3019 3 1583 0.14018692 3
## 620 221 424 175 234 2679 3 1589 0.14018692 3
## 713 225 272 113 247 2287 3 1604 0.14018692 3
## 720 240 402 143 176 2256 3 1605 0.14018692 3
## 722 197 278 111 240 2252 3 1606 0.14018692 3
## 732 298 326 117 191 2213 3 1607 0.14018692 3
## 738 157 256 117 245 2192 3 1608 0.14018692 3
## 809 209 247 107 203 1990 3 1620 0.14018692 3
## 859 220 274 105 150 1819 3 1623 0.14018692 3
## 988 158 231 95 139 1526 3 1642 0.14018692 3
## 1073 102 165 74 150 1349 3 1653 0.14018692 3
## 1107 92 195 67 142 1302 3 1659 0.14018692 3
## 1142 145 181 80 115 1259 3 1663 0.14018692 3
## 1154 135 163 57 124 1245 3 1667 0.14018692 3
## 1204 122 187 54 113 1177 3 1681 0.14018692 3
## 1387 122 108 51 81 929 3 1709 0.14018692 3
## 1388 66 106 45 109 927 3 1710 0.14018692 3
## 1522 73 101 33 87 784 3 1728 0.14018692 3
## 1590 67 84 40 75 721 3 1736 0.14018692 3
## 1644 68 97 44 61 669 3 1748 0.14018692 3
## 1689 40 109 36 60 641 3 1757 0.14018692 3
## 1732 50 90 32 62 609 3 1759 0.14018692 3
## 1787 70 64 35 56 566 3 1765 0.14018692 3
## 1856 43 58 22 62 518 3 1770 0.14018692 3
## 1897 44 66 30 52 493 3 1775 0.14018692 3
## 1996 63 69 17 36 434 3 1781 0.14018692 3
## 2170 21 37 16 39 344 3 1795 0.14018692 3
## 2223 38 47 16 29 321 3 1797 0.14018692 3
## 2255 17 27 27 33 307 3 1799 0.14018692 3
## 2286 36 43 19 25 290 3 1803 0.14018692 3
## 2355 25 29 16 26 257 3 1807 0.14018692 3
## 2400 25 30 12 25 240 3 1810 0.14018692 3
## 2493 18 26 14 19 201 3 1816 0.14018692 3
## 2537 9 18 12 22 183 3 1817 0.14018692 3
## 41 3335 2733 931 1161 17190 4 1832 0.08771930 4
## 57 2997 2453 889 1097 15822 4 1833 0.08771930 4
## 156 1309 1366 446 594 8348 4 1855 0.08771930 4
## 171 801 1042 439 750 7589 4 1857 0.08771930 4
## 181 895 1022 404 703 7455 4 1859 0.08771930 4
## 248 1156 975 320 400 5879 4 1868 0.08771930 4
## 286 628 736 281 426 5210 4 1874 0.08771930 4
## 349 479 733 238 416 4475 4 1885 0.08771930 4
## 375 419 531 232 443 4182 4 1890 0.08771930 4
## 463 405 445 198 351 3544 4 1902 0.08771930 4
## 659 341 333 142 224 2486 4 1920 0.08771930 4
## 730 211 293 134 225 2220 4 1924 0.08771930 4
## 1368 96 156 59 84 950 4 1964 0.08771930 4
## 1760 39 68 31 66 587 4 1981 0.08771930 4
## 1914 78 79 27 25 487 4 1988 0.08771930 4
## 399 341 471 190 446 3969 5 2056 0.04246285 5
## 533 349 342 157 317 3062 5 2083 0.04246285 5
## 606 278 356 167 276 2724 5 2100 0.04246285 5
## 649 502 358 112 167 2549 5 2109 0.04246285 5
## 756 200 240 108 241 2135 5 2129 0.04246285 5
## 773 164 226 87 231 2085 5 2133 0.04246285 5
## 804 234 262 102 186 2008 5 2137 0.04246285 5
## 906 177 195 81 183 1710 5 2158 0.04246285 5
## 1020 100 153 46 162 1456 5 2172 0.04246285 5
## 1058 124 115 52 156 1379 5 2180 0.04246285 5
## 1126 129 154 64 129 1280 5 2193 0.04246285 5
## 1181 129 133 67 130 1210 5 2200 0.04246285 5
## 1185 164 158 88 93 1201 5 2203 0.04246285 5
## 1365 71 90 36 108 955 5 2234 0.04246285 5
## 1439 113 115 62 77 880 5 2244 0.04246285 5
## 1815 60 71 33 49 542 5 2287 0.04246285 5
## 1835 72 82 24 48 530 5 2292 0.04246285 5
## 2182 33 38 13 23 341 5 2336 0.04246285 5
## 2291 25 33 8 27 288 5 2350 0.04246285 5
## 2453 30 27 17 16 221 5 2370 0.04246285 5
## 238 732 871 385 537 6074 6 2495 0.05913978 6
## 347 349 384 196 344 4495 6 2508 0.05913978 6
## 577 215 243 126 345 2850 6 2535 0.05913978 6
## 914 120 195 94 185 1675 6 2573 0.05913978 6
## 1362 75 97 41 80 958 6 2598 0.05913978 6
## 1593 89 86 39 67 718 6 2614 0.05913978 6
## 1627 95 97 38 59 686 6 2616 0.05913978 6
## 1818 27 47 24 70 539 6 2625 0.05913978 6
## 1977 26 48 25 49 444 6 2630 0.05913978 6
## 2103 26 35 14 48 374 6 2634 0.05913978 6
## 2582 13 22 13 15 167 6 2648 0.05913978 6
## 21 6567 3668 1710 1475 26560 7 2652 0.17959184 7
## 47 2661 2171 1228 1234 16463 7 2653 0.17959184 7
## 167 1634 1077 576 487 7720 7 2655 0.17959184 7
## 215 1244 983 445 408 6428 7 2658 0.17959184 7
## 245 779 804 481 481 5907 7 2661 0.17959184 7
## 258 1127 790 399 369 5558 7 2663 0.17959184 7
## 441 816 498 226 222 3694 7 2669 0.17959184 7
## 474 511 488 301 219 3415 7 2673 0.17959184 7
## 771 269 265 135 141 2092 7 2685 0.17959184 7
## 857 290 240 129 147 1825 7 2689 0.17959184 7
## 949 188 205 120 133 1600 7 2693 0.17959184 7
## 972 161 176 127 149 1559 7 2698 0.17959184 7
## 976 214 170 106 124 1548 7 2700 0.17959184 7
## 1013 140 208 107 139 1478 7 2704 0.17959184 7
## 1085 174 193 105 104 1328 7 2712 0.17959184 7
## 1192 102 135 96 108 1196 7 2716 0.17959184 7
## 1268 138 130 75 107 1098 7 2724 0.17959184 7
## 1319 123 119 69 98 1010 7 2731 0.17959184 7
## 1396 148 114 65 68 915 7 2735 0.17959184 7
## 1473 133 115 53 63 845 7 2741 0.17959184 7
## 1484 121 90 82 67 827 7 2742 0.17959184 7
## 1571 110 106 47 59 736 7 2748 0.17959184 7
## 1587 84 90 41 69 723 7 2749 0.17959184 7
## 1800 86 85 40 34 558 7 2770 0.17959184 7
## 1853 68 56 44 43 521 7 2773 0.17959184 7
## 1890 62 43 28 49 496 7 2780 0.17959184 7
## 1922 52 51 55 37 484 7 2787 0.17959184 7
## 1963 73 71 31 30 454 7 2791 0.17959184 7
## 2128 70 51 21 30 363 7 2802 0.17959184 7
## 2174 47 40 29 31 343 7 2808 0.17959184 7
## 2274 50 42 26 20 296 7 2815 0.17959184 7
## 2367 45 32 18 16 254 7 2825 0.17959184 7
## 2368 47 32 13 15 254 7 2826 0.17959184 7
## 2479 39 34 16 12 206 7 2834 0.17959184 7
## 2676 12 14 9 14 126 7 2849 0.17959184 7
## 2712 11 17 8 9 112 7 2851 0.17959184 7
## 2732 9 15 7 10 105 7 2854 0.17959184 7
## 2774 18 18 6 2 83 7 2867 0.17959184 7
## 2796 14 7 8 5 73 7 2871 0.17959184 7
## 2797 12 17 6 2 72 7 2872 0.17959184 7
## 2817 13 10 4 4 65 7 2877 0.17959184 7
## 2828 10 11 4 4 61 7 2879 0.17959184 7
## 2863 10 4 6 2 45 7 2886 0.17959184 7
## 2896 4 4 3 0 22 7 2896 0.17959184 7# see the result using a contingency table
table(st$REG)##
## 1 2 3 4 5 6 7
## 30 20 45 15 20 11 44