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

Ejemplo de uso para obtener un muestreo en dos fases

Cargamos la librería sampling

if (!require(sampling)){
   install.packages("sampling",dep=TRUE)
   library(sampling)
}
## Loading required package: sampling

La 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)

Muestreo estratificado en la primera etapa y por conglomerados en la segunda

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

Ordenamos 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   0

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

Código

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

  • srswor : simple random sampling without replacement
  • srswr : simple random sampling with replacement
  • poisson : poisson sampling
  • systematic : systematic samplin

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 40

La 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       1
head(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.6666667

Para 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.5000000

Verificamos 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   0

Ejemplo con función strata

data(swissmunicipalities)
datos <- swissmunicipalities

datos <- datos[order(datos$REG),]
table(datos$REG)
## 
##   1   2   3   4   5   6   7 
## 589 913 321 171 471 186 245

Si 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