sms library tutorial: Preparing microdata
Dr. Dimitris Kavroudakis
Introduction
Spatial microsimulation is a methodology aiming to simulate entities such as households, individuals or business in the finer possible scale. This process requires the use of micro-datasets with individual based information. The sms library presented in this work, facilitates the production of small area population micro-data by combining various other datasets. This is a parallelized implementation of random selection with thoptimization to select a micro-dataset with characteristics that mach a macro description. The sms library contains functions for preparing micro-data from census and longitudinal datasets. This process is called fiting and is the primary step in spatial microsimulation methodology. It actualy fits small area datasets to census area descriptions. There is a need for a library in Statistcal Programming Language R for preparing small area microdata. The functions of the livbrary connect various data-sources and produce small area population dataset. The functions use multi-core approaches of modern personal desktop computers in order to simulate relatively big areas in less computational time. The sms library, uses the parallel interface of R platform to divide the main simulation process into smaller simulations which then run in parallel. The R library presented in this work, use parallel processing approaches for the efficient production of small area population micro-data, which can be later used for geographical analysis.
This is an example of using sms library. It uses the default data (census, survey) and prepares fitted microdata for the areas described in the census dataset.
Load data
In order to illustrate the use of sms library, it is useful to use a simple case study. The following case study will produce micro-data of the population of Lesvos island in Greece. The data-sources used are from the census of population of Greece and from an artificial panel dataset constructed for this case study. Similar panel datasets can be found from Eurostat.
library(sms)## Loading required package: doParallel
## Loading required package: foreach
## Loading required package: iterators
## Loading required package: paralleldata(survey)
head(survey)## pid he female agemature car_owner house_owner working annualIncome
## 1 6001 0 1 1 1 0 1 16567
## 2 6002 1 1 1 1 1 1 2458
## 3 6003 0 1 0 0 1 0 9437
## 4 6004 1 0 0 1 0 0 22130
## 5 6005 1 0 0 0 1 1 3936
## 6 6006 1 0 0 1 0 1 16695data(census)
census## areaid population he females
## 1 8301 56 46 42
## 2 8302 73 42 15
## 3 8303 58 12 10
## 4 8304 78 43 21
## 5 8305 73 17 60
## 6 8306 77 15 11
## 7 8307 66 37 20
## 8 8308 78 41 42
## 9 8309 77 56 10
## 10 8310 78 55 26
## 11 8311 58 19 60
## 12 8312 68 20 60
## 13 8313 60 56 51Now both datasets have been loaded as data.frames. We need another data.frame which will hold the associations between the census and the survey data.frames. It will hold the names of the columns from both datasets. This is called “data lexicon” and is just holding column names (Strings)
in.lexicon=createLexicon()
in.lexicon=addDataAssociation(in.lexicon, c("he","he"))
in.lexicon=addDataAssociation(in.lexicon, c("females","female"))
in.lexicon## con_1 con_2
## census_row he females
## survey_row he femaleNote about data
The first column of your census data should have a column named: areaid which should be a unique identifier of each area. ALso notice that the survey data sould use binary (0 or 1) values for all variables.
Fitting data
We need to construct a microsimulation object which will hold all necessery objects of our simulation. Thos objects are the imported data and after the simulation, the results will be stored inside the microsimulation object. This is the main object of every simulation process and it holds all necessary information and results of the fitting process. The following command creates a microsimulation object with the name insms and attach the census and survey data to it. It also defines the number of iterations (40) before ending the fitting process.
mysms=new("microsimulation", iterations=40,census=census, panel=survey,lexicon=in.lexicon)The object mysms now has 5 slots which contain various objects and information. We can explore the slot names with the following command:
slotNames(mysms)## [1] "census" "panel" "lexicon" "results" "iterations"The census slot, contains the census data and the lpane slot the panel data of the simulation. The lexicon slot contains the data association information and the results slot contains the results of the simulation process which is currently empty. Finally the iterations slot contains the number of iterations before the selection of the best combination of individuals for each area.
Now that we have a microsimulation object, we use it to start the fitting process.
try01 = run_parallel_HC( mysms )The object try01 is a new microsimulation object which holds the results of the fitting process as well as the initial data.
Explore the results
The results slot, inside the try01 object, holds the results of the fitting process.
try01## An object of class "microsimulation"
## Slot "census":
## areaid population he females
## 1 8301 56 46 42
## 2 8302 73 42 15
## 3 8303 58 12 10
## 4 8304 78 43 21
## 5 8305 73 17 60
## 6 8306 77 15 11
## 7 8307 66 37 20
## 8 8308 78 41 42
## 9 8309 77 56 10
## 10 8310 78 55 26
## 11 8311 58 19 60
## 12 8312 68 20 60
## 13 8313 60 56 51
##
## Slot "panel":
## pid he female agemature car_owner house_owner working annualIncome
## 1 6001 0 1 1 1 0 1 16567
## 2 6002 1 1 1 1 1 1 2458
## 3 6003 0 1 0 0 1 0 9437
## 4 6004 1 0 0 1 0 0 22130
## 5 6005 1 0 0 0 1 1 3936
## 6 6006 1 0 0 1 0 1 16695
## 7 6007 0 0 1 1 1 0 9535
## 8 6008 0 0 1 1 0 1 1136
## 9 6009 1 1 1 1 0 0 24170
## 10 6010 0 1 1 0 1 0 10119
## 11 6011 1 1 1 1 0 1 14821
## 12 6012 1 0 1 0 1 1 7016
## 13 6013 1 0 1 1 0 1 20227
## 14 6014 1 1 0 1 1 1 6014
## 15 6015 0 1 1 0 0 0 19546
## 16 6016 1 1 0 0 0 1 21714
## 17 6017 1 0 0 1 1 1 18431
## 18 6018 0 0 0 1 1 1 23950
## 19 6019 0 1 1 1 0 0 4546
## 20 6020 1 0 0 1 0 1 1010
## 21 6021 0 1 1 1 1 0 16588
## 22 6022 0 1 0 1 0 1 2871
## 23 6023 1 1 0 0 1 1 18244
## 24 6024 1 1 0 0 0 0 24865
## 25 6025 1 0 1 0 0 0 14752
## 26 6026 0 0 0 1 0 0 22364
## 27 6027 1 0 1 1 0 1 1347
## 28 6028 1 0 0 0 0 0 12338
## 29 6029 0 1 1 1 0 0 22979
## 30 6030 0 0 0 1 1 0 2245
## 31 6031 0 0 1 0 1 1 10229
## 32 6032 0 0 1 1 1 0 9894
## 33 6033 0 1 1 0 0 0 15822
## 34 6034 1 0 0 1 1 0 8897
## 35 6035 1 1 1 0 1 0 23377
## 36 6036 0 0 1 1 1 1 23880
## 37 6037 0 0 1 1 0 0 7093
## 38 6038 1 0 1 1 1 0 4508
## 39 6039 0 0 0 0 1 1 3470
## 40 6040 1 1 0 0 0 0 13271
## 41 6041 1 1 1 1 1 1 997
## 42 6042 0 0 1 1 0 1 19157
## 43 6043 0 0 0 1 0 1 3170
## 44 6044 1 0 1 0 0 0 1844
## 45 6045 0 0 0 0 1 1 13377
## 46 6046 1 0 0 1 1 0 24515
## 47 6047 1 1 0 0 0 0 14501
## 48 6048 0 1 1 0 0 1 4032
## 49 6049 0 1 1 1 0 1 7873
## 50 6050 0 1 1 1 1 0 12243
## 51 6051 0 0 1 1 0 1 6669
## 52 6052 0 1 0 1 1 1 21240
## 53 6053 1 0 1 0 0 0 8081
## 54 6054 1 0 0 1 1 1 20903
## 55 6055 0 0 1 1 0 0 9651
## 56 6056 0 1 0 1 1 0 19291
## 57 6057 1 0 1 1 1 1 24504
## 58 6058 1 0 0 0 0 0 22892
## 59 6059 1 0 1 0 0 1 18409
## 60 6060 0 0 0 1 0 0 7352
## 61 6061 0 0 0 1 1 0 1202
## 62 6062 0 0 0 0 0 0 5375
## 63 6063 0 0 1 1 0 0 22695
## 64 6064 1 1 0 0 1 0 19559
## 65 6065 1 0 0 1 0 0 11579
## 66 6066 0 0 1 1 1 1 9266
## 67 6067 1 1 1 1 0 1 19112
## 68 6068 0 1 0 1 1 0 20739
## 69 6069 1 1 1 1 1 1 15149
## 70 6070 1 0 0 1 1 1 16285
## 71 6071 0 0 0 1 0 0 23451
## 72 6072 0 1 0 1 0 1 861
## 73 6073 0 0 1 0 0 0 23682
## 74 6074 1 0 0 0 1 1 16511
## 75 6075 0 0 0 0 0 0 5855
## 76 6076 1 1 1 0 0 0 14320
## 77 6077 1 1 1 0 1 1 3660
## 78 6078 0 1 0 0 0 1 12626
## 79 6079 1 1 0 0 1 1 16117
## 80 6080 1 0 0 0 1 0 14511
## 81 6081 0 1 0 0 0 0 19338
## 82 6082 0 0 1 1 1 0 6233
## 83 6083 1 1 0 0 1 1 15286
## 84 6084 0 1 0 1 0 0 24604
## 85 6085 0 1 1 1 1 1 14183
## 86 6086 0 1 1 1 0 1 3552
## 87 6087 1 1 1 1 1 0 4265
## 88 6088 0 0 0 0 1 0 23350
## 89 6089 1 1 1 0 1 1 18804
## 90 6090 0 1 0 1 1 1 20528
## 91 6091 1 1 0 0 0 0 5695
## 92 6092 0 1 0 0 1 0 22133
## 93 6093 0 1 0 1 1 0 20926
## 94 6094 0 0 1 1 1 1 18017
## 95 6095 1 1 0 0 0 1 4143
## 96 6096 1 0 0 0 1 0 15223
## 97 6097 1 1 1 0 1 1 21521
## 98 6098 0 1 1 0 0 0 23423
## 99 6099 1 0 1 0 1 0 15471
## 100 6100 0 1 1 1 1 0 19540
## 101 6101 0 1 1 1 0 1 21034
## 102 6102 0 0 0 1 1 0 12107
## 103 6103 1 0 1 0 0 0 13352
## 104 6104 0 0 1 1 1 1 14046
## 105 6105 1 1 1 0 1 0 6742
## 106 6106 1 1 1 1 0 0 11994
## 107 6107 1 1 0 0 0 0 2070
## 108 6108 0 1 1 0 1 0 19083
## 109 6109 0 1 1 0 0 0 1120
## 110 6110 1 1 0 1 1 0 16503
## 111 6111 0 0 1 0 1 1 16075
## 112 6112 0 0 1 0 0 0 23382
## 113 6113 0 0 0 1 1 1 18582
## 114 6114 1 0 1 0 1 1 21039
## 115 6115 1 0 1 1 1 0 17012
## 116 6116 0 1 0 0 1 0 4414
## 117 6117 1 1 0 1 0 0 24753
## 118 6118 1 1 1 0 0 1 5610
## 119 6119 1 1 1 0 0 0 24288
## 120 6120 0 0 0 1 0 0 1025
## 121 6121 1 1 1 1 0 0 2352
## 122 6122 0 1 0 1 0 1 5706
## 123 6123 1 1 1 1 0 0 9458
## 124 6124 0 0 0 1 0 1 15521
## 125 6125 0 1 0 1 0 1 16596
## 126 6126 1 1 1 0 1 1 12386
## 127 6127 0 1 0 1 1 0 12893
## 128 6128 0 1 0 0 1 1 1066
## 129 6129 0 1 1 1 0 0 5539
## 130 6130 1 1 1 1 1 1 7340
## 131 6131 1 0 1 1 1 0 4416
## 132 6132 1 1 0 1 1 1 12931
## 133 6133 1 0 0 1 1 1 16596
## 134 6134 0 0 0 1 1 0 13271
## 135 6135 1 0 1 0 0 1 14497
## 136 6136 0 1 1 0 0 0 1716
## 137 6137 1 1 0 0 0 0 11019
## 138 6138 0 0 0 0 1 0 728
## 139 6139 1 1 0 1 0 0 24333
## 140 6140 0 1 0 1 0 1 17720
## 141 6141 0 0 0 0 0 1 8217
## 142 6142 0 1 1 1 0 0 24860
## 143 6143 0 0 1 0 0 0 4922
## 144 6144 1 0 0 0 0 0 10261
## 145 6145 1 1 0 1 0 1 9894
## 146 6146 0 0 0 0 1 0 17654
## 147 6147 1 0 1 0 1 0 8914
## 148 6148 1 0 1 1 0 1 4717
## 149 6149 0 0 0 0 0 1 21917
## 150 6150 0 1 0 1 1 0 23632
## 151 6151 0 0 1 0 0 0 18055
## 152 6152 1 1 0 0 1 0 11947
## 153 6153 0 0 0 1 0 1 13642
## 154 6154 0 0 0 1 0 0 24662
## 155 6155 1 0 0 1 0 0 13186
## 156 6156 0 1 1 0 1 0 12150
## 157 6157 0 0 1 1 1 1 18831
## 158 6158 1 0 0 1 0 0 755
## 159 6159 0 0 1 1 1 1 17199
## 160 6160 1 1 1 1 0 1 13746
## 161 6161 0 1 0 1 0 0 23118
## 162 6162 1 1 1 0 1 1 8955
## 163 6163 0 1 0 0 1 1 9332
## 164 6164 0 1 0 0 1 1 21865
## 165 6165 1 0 1 0 1 1 1125
## 166 6166 0 1 1 0 1 1 15312
## 167 6167 0 1 1 1 1 1 19045
## 168 6168 0 1 0 1 1 0 17495
## 169 6169 1 0 1 1 1 0 3719
## 170 6170 0 1 1 0 0 0 12948
## 171 6171 1 0 1 0 1 1 20752
## 172 6172 1 0 1 0 0 1 14346
## 173 6173 0 1 1 1 1 1 7025
## 174 6174 0 0 1 1 0 0 6798
## 175 6175 1 0 0 0 0 0 4067
## 176 6176 1 0 1 0 0 0 23390
## 177 6177 1 0 0 1 1 0 17102
## 178 6178 0 1 0 1 0 1 5271
## 179 6179 0 0 0 0 1 0 14058
## 180 6180 0 1 0 0 1 1 14699
## 181 6181 0 1 1 1 1 0 24356
## 182 6182 1 0 0 0 1 1 1369
## 183 6183 1 1 0 1 0 0 20091
## 184 6184 0 0 1 1 1 0 18460
## 185 6185 1 0 1 1 1 1 12862
## 186 6186 0 0 0 0 1 1 21279
## 187 6187 1 1 0 0 1 1 6766
## 188 6188 1 0 0 0 1 1 21100
## 189 6189 0 1 0 0 0 1 3754
## 190 6190 0 1 0 0 0 1 6535
## 191 6191 1 1 1 0 1 0 24296
## 192 6192 1 0 1 0 1 1 1649
## 193 6193 1 1 1 0 1 1 17777
## 194 6194 1 1 1 0 1 0 15198
## 195 6195 1 0 0 1 0 1 21223
## 196 6196 1 0 1 0 0 1 17683
## 197 6197 0 0 1 1 0 1 16944
## 198 6198 1 0 0 0 0 0 13894
## 199 6199 1 0 1 1 1 1 20173
## 200 6200 0 1 1 0 1 1 12623
##
## Slot "lexicon":
## con_1 con_2
## census_row he females
## survey_row he female
##
## Slot "results":
## [[1]]
## [[1]]$areaid
## [1] 8301
##
## [[1]]$selection
## pid he female agemature car_owner house_owner working annualIncome
## 42 6042 0 0 1 1 0 1 19157
## 118 6118 1 1 1 0 0 1 5610
## 148 6148 1 0 1 1 0 1 4717
## 2 6002 1 1 1 1 1 1 2458
## 108 6108 0 1 1 0 1 0 19083
## 14 6014 1 1 0 1 1 1 6014
## 123 6123 1 1 1 1 0 0 9458
## 125 6125 0 1 0 1 0 1 16596
## 129 6129 0 1 1 1 0 0 5539
## 72 6072 0 1 0 1 0 1 861
## 1 6001 0 1 1 1 0 1 16567
## 48 6048 0 1 1 0 0 1 4032
## 44 6044 1 0 1 0 0 0 1844
## 54 6054 1 0 0 1 1 1 20903
## 59 6059 1 0 1 0 0 1 18409
## 58 6058 1 0 0 0 0 0 22892
## 137 6137 1 1 0 0 0 0 11019
## 173 6173 0 1 1 1 1 1 7025
## 84 6084 0 1 0 1 0 0 24604
## 172 6172 1 0 1 0 0 1 14346
## 11 6011 1 1 1 1 0 1 14821
## 157 6157 0 0 1 1 1 1 18831
## 193 6193 1 1 1 0 1 1 17777
## 79 6079 1 1 0 0 1 1 16117
## 101 6101 0 1 1 1 0 1 21034
## 126 6126 1 1 1 0 1 1 12386
## 118.1 6118 1 1 1 0 0 1 5610
## 130 6130 1 1 1 1 1 1 7340
## 59.1 6059 1 0 1 0 0 1 18409
## 192 6192 1 0 1 0 1 1 1649
## 193.1 6193 1 1 1 0 1 1 17777
## 43 6043 0 0 0 1 0 1 3170
## 41 6041 1 1 1 1 1 1 997
## 16 6016 1 1 0 0 0 1 21714
## 137.1 6137 1 1 0 0 0 0 11019
## 21 6021 0 1 1 1 1 0 16588
## 28 6028 1 0 0 0 0 0 12338
## 154 6154 0 0 0 1 0 0 24662
## 44.1 6044 1 0 1 0 0 0 1844
## 44.2 6044 1 0 1 0 0 0 1844
## 84.1 6084 0 1 0 1 0 0 24604
## 5 6005 1 0 0 0 1 1 3936
## 160 6160 1 1 1 1 0 1 13746
## 190 6190 0 1 0 0 0 1 6535
## 82 6082 0 0 1 1 1 0 6233
## 153 6153 0 0 0 1 0 1 13642
## 119 6119 1 1 1 0 0 0 24288
## 121 6121 1 1 1 1 0 0 2352
## 75 6075 0 0 0 0 0 0 5855
## 122 6122 0 1 0 1 0 1 5706
## 14.1 6014 1 1 0 1 1 1 6014
## 103 6103 1 0 1 0 0 0 13352
## 198 6198 1 0 0 0 0 0 13894
## 54.1 6054 1 0 0 1 1 1 20903
## 132 6132 1 1 0 1 1 1 12931
## 26 6026 0 0 0 1 0 0 22364
##
## [[1]]$tae
## [1] 20
##
## [[1]]$tries
## [1] 35 29 25 26 32 33 36 38 33 42 32 28 41 41 36 29 33 25 20 27 32 32 33
## [24] 31 39 36 31 36 29 27 28 29 35 34 40 37 34 28 39 31
##
## [[1]]$error_states
## [1] 35 29 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 20 20 20 20 20
## [24] 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20
##
##
## [[2]]
## [[2]]$areaid
## [1] 8302
##
## [[2]]$selection
## pid he female agemature car_owner house_owner working annualIncome
## 137 6137 1 1 0 0 0 0 11019
## 176 6176 1 0 1 0 0 0 23390
## 152 6152 1 1 0 0 1 0 11947
## 146 6146 0 0 0 0 1 0 17654
## 191 6191 1 1 1 0 1 0 24296
## 63 6063 0 0 1 1 0 0 22695
## 199 6199 1 0 1 1 1 1 20173
## 2 6002 1 1 1 1 1 1 2458
## 94 6094 0 0 1 1 1 1 18017
## 148 6148 1 0 1 1 0 1 4717
## 137.1 6137 1 1 0 0 0 0 11019
## 4 6004 1 0 0 1 0 0 22130
## 13 6013 1 0 1 1 0 1 20227
## 57 6057 1 0 1 1 1 1 24504
## 23 6023 1 1 0 0 1 1 18244
## 27 6027 1 0 1 1 0 1 1347
## 17 6017 1 0 0 1 1 1 18431
## 92 6092 0 1 0 0 1 0 22133
## 51 6051 0 0 1 1 0 1 6669
## 163 6163 0 1 0 0 1 1 9332
## 66 6066 0 0 1 1 1 1 9266
## 107 6107 1 1 0 0 0 0 2070
## 41 6041 1 1 1 1 1 1 997
## 82 6082 0 0 1 1 1 0 6233
## 18 6018 0 0 0 1 1 1 23950
## 143 6143 0 0 1 0 0 0 4922
## 142 6142 0 1 1 1 0 0 24860
## 47 6047 1 1 0 0 0 0 14501
## 23.1 6023 1 1 0 0 1 1 18244
## 62 6062 0 0 0 0 0 0 5375
## 94.1 6094 0 0 1 1 1 1 18017
## 30 6030 0 0 0 1 1 0 2245
## 143.1 6143 0 0 1 0 0 0 4922
## 33 6033 0 1 1 0 0 0 15822
## 15 6015 0 1 1 0 0 0 19546
## 108 6108 0 1 1 0 1 0 19083
## 80 6080 1 0 0 0 1 0 14511
## 185 6185 1 0 1 1 1 1 12862
## 34 6034 1 0 0 1 1 0 8897
## 162 6162 1 1 1 0 1 1 8955
## 112 6112 0 0 1 0 0 0 23382
## 36 6036 0 0 1 1 1 1 23880
## 63.1 6063 0 0 1 1 0 0 22695
## 179 6179 0 0 0 0 1 0 14058
## 55 6055 0 0 1 1 0 0 9651
## 33.1 6033 0 1 1 0 0 0 15822
## 28 6028 1 0 0 0 0 0 12338
## 158 6158 1 0 0 1 0 0 755
## 74 6074 1 0 0 0 1 1 16511
## 171 6171 1 0 1 0 1 1 20752
## 59 6059 1 0 1 0 0 1 18409
## 192 6192 1 0 1 0 1 1 1649
## 62.1 6062 0 0 0 0 0 0 5375
## 20 6020 1 0 0 1 0 1 1010
## 67 6067 1 1 1 1 0 1 19112
## 45 6045 0 0 0 0 1 1 13377
## 105 6105 1 1 1 0 1 0 6742
## 186 6186 0 0 0 0 1 1 21279
## 161 6161 0 1 0 1 0 0 23118
## 77 6077 1 1 1 0 1 1 3660
## 61 6061 0 0 0 1 1 0 1202
## 143.2 6143 0 0 1 0 0 0 4922
## 173 6173 0 1 1 1 1 1 7025
## 188 6188 1 0 0 0 1 1 21100
## 106 6106 1 1 1 1 0 0 11994
## 177 6177 1 0 0 1 1 0 17102
## 121 6121 1 1 1 1 0 0 2352
## 111 6111 0 0 1 0 1 1 16075
## 73 6073 0 0 1 0 0 0 23682
## 126 6126 1 1 1 0 1 1 12386
## 140 6140 0 1 0 1 0 1 17720
## 156 6156 0 1 1 0 1 0 12150
## 88 6088 0 0 0 0 1 0 23350
##
## [[2]]$tae
## [1] 18
##
## [[2]]$tries
## [1] 22 22 23 27 39 29 25 27 31 23 28 29 19 31 21 40 26 26 18 32 26 27 31
## [24] 23 34 32 39 31 23 26 26 32 38 23 30 18 26 30 34 31
##
## [[2]]$error_states
## [1] 22 22 22 22 22 22 22 22 22 22 22 22 19 19 19 19 19 19 18 18 18 18 18
## [24] 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18 18
##
##
## [[3]]
## [[3]]$areaid
## [1] 8303
##
## [[3]]$selection
## pid he female agemature car_owner house_owner working annualIncome
## 158 6158 1 0 0 1 0 0 755
## 35 6035 1 1 1 0 1 0 23377
## 94 6094 0 0 1 1 1 1 18017
## 150 6150 0 1 0 1 1 0 23632
## 177 6177 1 0 0 1 1 0 17102
## 45 6045 0 0 0 0 1 1 13377
## 94.1 6094 0 0 1 1 1 1 18017
## 96 6096 1 0 0 0 1 0 15223
## 151 6151 0 0 1 0 0 0 18055
## 174 6174 0 0 1 1 0 0 6798
## 174.1 6174 0 0 1 1 0 0 6798
## 164 6164 0 1 0 0 1 1 21865
## 139 6139 1 1 0 1 0 0 24333
## 143 6143 0 0 1 0 0 0 4922
## 2 6002 1 1 1 1 1 1 2458
## 170 6170 0 1 1 0 0 0 12948
## 98 6098 0 1 1 0 0 0 23423
## 188 6188 1 0 0 0 1 1 21100
## 141 6141 0 0 0 0 0 1 8217
## 184 6184 0 0 1 1 1 0 18460
## 59 6059 1 0 1 0 0 1 18409
## 122 6122 0 1 0 1 0 1 5706
## 45.1 6045 0 0 0 0 1 1 13377
## 59.1 6059 1 0 1 0 0 1 18409
## 195 6195 1 0 0 1 0 1 21223
## 71 6071 0 0 0 1 0 0 23451
## 60 6060 0 0 0 1 0 0 7352
## 15 6015 0 1 1 0 0 0 19546
## 108 6108 0 1 1 0 1 0 19083
## 130 6130 1 1 1 1 1 1 7340
## 103 6103 1 0 1 0 0 0 13352
## 84 6084 0 1 0 1 0 0 24604
## 178 6178 0 1 0 1 0 1 5271
## 153 6153 0 0 0 1 0 1 13642
## 147 6147 1 0 1 0 1 0 8914
## 48 6048 0 1 1 0 0 1 4032
## 27 6027 1 0 1 1 0 1 1347
## 153.1 6153 0 0 0 1 0 1 13642
## 161 6161 0 1 0 1 0 0 23118
## 130.1 6130 1 1 1 1 1 1 7340
## 33 6033 0 1 1 0 0 0 15822
## 184.1 6184 0 0 1 1 1 0 18460
## 71.1 6071 0 0 0 1 0 0 23451
## 90 6090 0 1 0 1 1 1 20528
## 18 6018 0 0 0 1 1 1 23950
## 33.1 6033 0 1 1 0 0 0 15822
## 154 6154 0 0 0 1 0 0 24662
## 188.1 6188 1 0 0 0 1 1 21100
## 32 6032 0 0 1 1 1 0 9894
## 13 6013 1 0 1 1 0 1 20227
## 144 6144 1 0 0 0 0 0 10261
## 16 6016 1 1 0 0 0 1 21714
## 159 6159 0 0 1 1 1 1 17199
## 112 6112 0 0 1 0 0 0 23382
## 187 6187 1 1 0 0 1 1 6766
## 167 6167 0 1 1 1 1 1 19045
## 161.1 6161 0 1 0 1 0 0 23118
## 78 6078 0 1 0 0 0 1 12626
##
## [[3]]$tae
## [1] 22
##
## [[3]]$tries
## [1] 42 35 22 39 39 32 41 28 27 43 29 36 32 39 33 39 31 36 38 36 36 39 35
## [24] 39 37 41 39 27 27 33 41 24 31 33 25 30 28 34 35 41
##
## [[3]]$error_states
## [1] 42 35 22 22 22 22 22 22 22 22 22 22 22 22 22 22 22 22 22 22 22 22 22
## [24] 22 22 22 22 22 22 22 22 22 22 22 22 22 22 22 22 22
##
##
## [[4]]
## [[4]]$areaid
## [1] 8304
##
## [[4]]$selection
## pid he female agemature car_owner house_owner working annualIncome
## 39 6039 0 0 0 0 1 1 3470
## 113 6113 0 0 0 1 1 1 18582
## 34 6034 1 0 0 1 1 0 8897
## 187 6187 1 1 0 0 1 1 6766
## 64 6064 1 1 0 0 1 0 19559
## 151 6151 0 0 1 0 0 0 18055
## 48 6048 0 1 1 0 0 1 4032
## 158 6158 1 0 0 1 0 0 755
## 167 6167 0 1 1 1 1 1 19045
## 133 6133 1 0 0 1 1 1 16596
## 2 6002 1 1 1 1 1 1 2458
## 191 6191 1 1 1 0 1 0 24296
## 192 6192 1 0 1 0 1 1 1649
## 177 6177 1 0 0 1 1 0 17102
## 171 6171 1 0 1 0 1 1 20752
## 38 6038 1 0 1 1 1 0 4508
## 96 6096 1 0 0 0 1 0 15223
## 157 6157 0 0 1 1 1 1 18831
## 173 6173 0 1 1 1 1 1 7025
## 79 6079 1 1 0 0 1 1 16117
## 26 6026 0 0 0 1 0 0 22364
## 165 6165 1 0 1 0 1 1 1125
## 43 6043 0 0 0 1 0 1 3170
## 78 6078 0 1 0 0 0 1 12626
## 91 6091 1 1 0 0 0 0 5695
## 174 6174 0 0 1 1 0 0 6798
## 94 6094 0 0 1 1 1 1 18017
## 23 6023 1 1 0 0 1 1 18244
## 198 6198 1 0 0 0 0 0 13894
## 57 6057 1 0 1 1 1 1 24504
## 165.1 6165 1 0 1 0 1 1 1125
## 93 6093 0 1 0 1 1 0 20926
## 109 6109 0 1 1 0 0 0 1120
## 140 6140 0 1 0 1 0 1 17720
## 2.1 6002 1 1 1 1 1 1 2458
## 38.1 6038 1 0 1 1 1 0 4508
## 85 6085 0 1 1 1 1 1 14183
## 50 6050 0 1 1 1 1 0 12243
## 183 6183 1 1 0 1 0 0 20091
## 11 6011 1 1 1 1 0 1 14821
## 92 6092 0 1 0 0 1 0 22133
## 182 6182 1 0 0 0 1 1 1369
## 149 6149 0 0 0 0 0 1 21917
## 132 6132 1 1 0 1 1 1 12931
## 28 6028 1 0 0 0 0 0 12338
## 180 6180 0 1 0 0 1 1 14699
## 80 6080 1 0 0 0 1 0 14511
## 198.1 6198 1 0 0 0 0 0 13894
## 30 6030 0 0 0 1 1 0 2245
## 182.1 6182 1 0 0 0 1 1 1369
## 94.1 6094 0 0 1 1 1 1 18017
## 116 6116 0 1 0 0 1 0 4414
## 51 6051 0 0 1 1 0 1 6669
## 49 6049 0 1 1 1 0 1 7873
## 4 6004 1 0 0 1 0 0 22130
## 92.1 6092 0 1 0 0 1 0 22133
## 198.2 6198 1 0 0 0 0 0 13894
## 117 6117 1 1 0 1 0 0 24753
## 193 6193 1 1 1 0 1 1 17777
## 171.1 6171 1 0 1 0 1 1 20752
## 38.2 6038 1 0 1 1 1 0 4508
## 140.1 6140 0 1 0 1 0 1 17720
## 173.1 6173 0 1 1 1 1 1 7025
## 195 6195 1 0 0 1 0 1 21223
## 63 6063 0 0 1 1 0 0 22695
## 144 6144 1 0 0 0 0 0 10261
## 136 6136 0 1 1 0 0 0 1716
## 5 6005 1 0 0 0 1 1 3936
## 2.2 6002 1 1 1 1 1 1 2458
## 7 6007 0 0 1 1 1 0 9535
## 73 6073 0 0 1 0 0 0 23682
## 187.1 6187 1 1 0 0 1 1 6766
## 177.1 6177 1 0 0 1 1 0 17102
## 182.2 6182 1 0 0 0 1 1 1369
## 192.1 6192 1 0 1 0 1 1 1649
## 40 6040 1 1 0 0 0 0 13271
## 144.1 6144 1 0 0 0 0 0 10261
## 190 6190 0 1 0 0 0 1 6535
##
## [[4]]$tae
## [1] 15
##
## [[4]]$tries
## [1] 21 27 30 27 28 15 26 30 34 25 26 22 29 25 19 29 29 23 26 21 25 23 23
## [24] 26 39 20 26 30 28 29 27 25 28 28 25 23 30 24 27 24
##
## [[4]]$error_states
## [1] 21 21 21 21 21 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15
## [24] 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15
##
##
## [[5]]
## [[5]]$areaid
## [1] 8305
##
## [[5]]$selection
## pid he female agemature car_owner house_owner working annualIncome
## 158 6158 1 0 0 1 0 0 755
## 78 6078 0 1 0 0 0 1 12626
## 128 6128 0 1 0 0 1 1 1066
## 164 6164 0 1 0 0 1 1 21865
## 95 6095 1 1 0 0 0 1 4143
## 167 6167 0 1 1 1 1 1 19045
## 182 6182 1 0 0 0 1 1 1369
## 41 6041 1 1 1 1 1 1 997
## 65 6065 1 0 0 1 0 0 11579
## 140 6140 0 1 0 1 0 1 17720
## 46 6046 1 0 0 1 1 0 24515
## 176 6176 1 0 1 0 0 0 23390
## 186 6186 0 0 0 0 1 1 21279
## 190 6190 0 1 0 0 0 1 6535
## 51 6051 0 0 1 1 0 1 6669
## 112 6112 0 0 1 0 0 0 23382
## 85 6085 0 1 1 1 1 1 14183
## 90 6090 0 1 0 1 1 1 20528
## 143 6143 0 0 1 0 0 0 4922
## 127 6127 0 1 0 1 1 0 12893
## 101 6101 0 1 1 1 0 1 21034
## 63 6063 0 0 1 1 0 0 22695
## 125 6125 0 1 0 1 0 1 16596
## 48 6048 0 1 1 0 0 1 4032
## 70 6070 1 0 0 1 1 1 16285
## 43 6043 0 0 0 1 0 1 3170
## 29 6029 0 1 1 1 0 0 22979
## 70.1 6070 1 0 0 1 1 1 16285
## 50 6050 0 1 1 1 1 0 12243
## 178 6178 0 1 0 1 0 1 5271
## 81 6081 0 1 0 0 0 0 19338
## 21 6021 0 1 1 1 1 0 16588
## 13 6013 1 0 1 1 0 1 20227
## 91 6091 1 1 0 0 0 0 5695
## 80 6080 1 0 0 0 1 0 14511
## 129 6129 0 1 1 1 0 0 5539
## 69 6069 1 1 1 1 1 1 15149
## 173 6173 0 1 1 1 1 1 7025
## 15 6015 0 1 1 0 0 0 19546
## 162 6162 1 1 1 0 1 1 8955
## 11 6011 1 1 1 1 0 1 14821
## 124 6124 0 0 0 1 0 1 15521
## 155 6155 1 0 0 1 0 0 13186
## 7 6007 0 0 1 1 1 0 9535
## 74 6074 1 0 0 0 1 1 16511
## 72 6072 0 1 0 1 0 1 861
## 74.1 6074 1 0 0 0 1 1 16511
## 180 6180 0 1 0 0 1 1 14699
## 155.1 6155 1 0 0 1 0 0 13186
## 113 6113 0 0 0 1 1 1 18582
## 166 6166 0 1 1 0 1 1 15312
## 54 6054 1 0 0 1 1 1 20903
## 175 6175 1 0 0 0 0 0 4067
## 190.1 6190 0 1 0 0 0 1 6535
## 84 6084 0 1 0 1 0 0 24604
## 125.1 6125 0 1 0 1 0 1 16596
## 76 6076 1 1 1 0 0 0 14320
## 116 6116 0 1 0 0 1 0 4414
## 128.1 6128 0 1 0 0 1 1 1066
## 70.2 6070 1 0 0 1 1 1 16285
## 26 6026 0 0 0 1 0 0 22364
## 98 6098 0 1 1 0 0 0 23423
## 196 6196 1 0 1 0 0 1 17683
## 125.2 6125 0 1 0 1 0 1 16596
## 194 6194 1 1 1 0 1 0 15198
## 70.3 6070 1 0 0 1 1 1 16285
## 137 6137 1 1 0 0 0 0 11019
## 36 6036 0 0 1 1 1 1 23880
## 55 6055 0 0 1 1 0 0 9651
## 97 6097 1 1 1 0 1 1 21521
## 97.1 6097 1 1 1 0 1 1 21521
## 158.1 6158 1 0 0 1 0 0 755
## 177 6177 1 0 0 1 1 0 17102
##
## [[5]]$tae
## [1] 33
##
## [[5]]$tries
## [1] 39 36 55 40 45 45 50 43 41 43 33 40 41 38 45 37 47 38 42 50 45 54 46
## [24] 57 45 49 47 48 35 40 35 46 35 41 39 33 39 48 51 43
##
## [[5]]$error_states
## [1] 39 36 36 36 36 36 36 36 36 36 33 33 33 33 33 33 33 33 33 33 33 33 33
## [24] 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33
##
##
## [[6]]
## [[6]]$areaid
## [1] 8306
##
## [[6]]$selection
## pid he female agemature car_owner house_owner working annualIncome
## 43 6043 0 0 0 1 0 1 3170
## 44 6044 1 0 1 0 0 0 1844
## 154 6154 0 0 0 1 0 0 24662
## 144 6144 1 0 0 0 0 0 10261
## 7 6007 0 0 1 1 1 0 9535
## 46 6046 1 0 0 1 1 0 24515
## 157 6157 0 0 1 1 1 1 18831
## 171 6171 1 0 1 0 1 1 20752
## 94 6094 0 0 1 1 1 1 18017
## 6 6006 1 0 0 1 0 1 16695
## 71 6071 0 0 0 1 0 0 23451
## 200 6200 0 1 1 0 1 1 12623
## 177 6177 1 0 0 1 1 0 17102
## 83 6083 1 1 0 0 1 1 15286
## 90 6090 0 1 0 1 1 1 20528
## 148 6148 1 0 1 1 0 1 4717
## 192 6192 1 0 1 0 1 1 1649
## 46.1 6046 1 0 0 1 1 0 24515
## 32 6032 0 0 1 1 1 0 9894
## 54 6054 1 0 0 1 1 1 20903
## 185 6185 1 0 1 1 1 1 12862
## 107 6107 1 1 0 0 0 0 2070
## 97 6097 1 1 1 0 1 1 21521
## 105 6105 1 1 1 0 1 0 6742
## 44.1 6044 1 0 1 0 0 0 1844
## 31 6031 0 0 1 0 1 1 10229
## 18 6018 0 0 0 1 1 1 23950
## 167 6167 0 1 1 1 1 1 19045
## 109 6109 0 1 1 0 0 0 1120
## 168 6168 0 1 0 1 1 0 17495
## 42 6042 0 0 1 1 0 1 19157
## 91 6091 1 1 0 0 0 0 5695
## 37 6037 0 0 1 1 0 0 7093
## 28 6028 1 0 0 0 0 0 12338
## 77 6077 1 1 1 0 1 1 3660
## 45 6045 0 0 0 0 1 1 13377
## 124 6124 0 0 0 1 0 1 15521
## 8 6008 0 0 1 1 0 1 1136
## 174 6174 0 0 1 1 0 0 6798
## 129 6129 0 1 1 1 0 0 5539
## 70 6070 1 0 0 1 1 1 16285
## 73 6073 0 0 1 0 0 0 23682
## 136 6136 0 1 1 0 0 0 1716
## 45.1 6045 0 0 0 0 1 1 13377
## 180 6180 0 1 0 0 1 1 14699
## 150 6150 0 1 0 1 1 0 23632
## 186 6186 0 0 0 0 1 1 21279
## 43.1 6043 0 0 0 1 0 1 3170
## 156 6156 0 1 1 0 1 0 12150
## 54.1 6054 1 0 0 1 1 1 20903
## 75 6075 0 0 0 0 0 0 5855
## 185.1 6185 1 0 1 1 1 1 12862
## 152 6152 1 1 0 0 1 0 11947
## 171.1 6171 1 0 1 0 1 1 20752
## 152.1 6152 1 1 0 0 1 0 11947
## 55 6055 0 0 1 1 0 0 9651
## 57 6057 1 0 1 1 1 1 24504
## 16 6016 1 1 0 0 0 1 21714
## 40 6040 1 1 0 0 0 0 13271
## 72 6072 0 1 0 1 0 1 861
## 106 6106 1 1 1 1 0 0 11994
## 149 6149 0 0 0 0 0 1 21917
## 93 6093 0 1 0 1 1 0 20926
## 156.1 6156 0 1 1 0 1 0 12150
## 86 6086 0 1 1 1 0 1 3552
## 195 6195 1 0 0 1 0 1 21223
## 41 6041 1 1 1 1 1 1 997
## 163 6163 0 1 0 0 1 1 9332
## 18.1 6018 0 0 0 1 1 1 23950
## 124.1 6124 0 0 0 1 0 1 15521
## 165 6165 1 0 1 0 1 1 1125
## 99 6099 1 0 1 0 1 0 15471
## 20 6020 1 0 0 1 0 1 1010
## 177.1 6177 1 0 0 1 1 0 17102
## 148.1 6148 1 0 1 1 0 1 4717
## 163.1 6163 0 1 0 0 1 1 9332
## 174.1 6174 0 0 1 1 0 0 6798
##
## [[6]]$tae
## [1] 38
##
## [[6]]$tries
## [1] 51 49 56 57 52 54 53 49 57 51 50 54 47 51 45 46 38 52 46 50 57 54 45
## [24] 56 53 40 55 48 57 51 48 51 49 51 44 58 41 65 48 43
##
## [[6]]$error_states
## [1] 51 49 49 49 49 49 49 49 49 49 49 49 47 47 45 45 38 38 38 38 38 38 38
## [24] 38 38 38 38 38 38 38 38 38 38 38 38 38 38 38 38 38
##
##
## [[7]]
## [[7]]$areaid
## [1] 8307
##
## [[7]]$selection
## pid he female agemature car_owner house_owner working annualIncome
## 200 6200 0 1 1 0 1 1 12623
## 138 6138 0 0 0 0 1 0 728
## 35 6035 1 1 1 0 1 0 23377
## 115 6115 1 0 1 1 1 0 17012
## 51 6051 0 0 1 1 0 1 6669
## 39 6039 0 0 0 0 1 1 3470
## 188 6188 1 0 0 0 1 1 21100
## 169 6169 1 0 1 1 1 0 3719
## 188.1 6188 1 0 0 0 1 1 21100
## 153 6153 0 0 0 1 0 1 13642
## 16 6016 1 1 0 0 0 1 21714
## 71 6071 0 0 0 1 0 0 23451
## 6 6006 1 0 0 1 0 1 16695
## 169.1 6169 1 0 1 1 1 0 3719
## 5 6005 1 0 0 0 1 1 3936
## 114 6114 1 0 1 0 1 1 21039
## 130 6130 1 1 1 1 1 1 7340
## 120 6120 0 0 0 1 0 0 1025
## 134 6134 0 0 0 1 1 0 13271
## 7 6007 0 0 1 1 1 0 9535
## 102 6102 0 0 0 1 1 0 12107
## 158 6158 1 0 0 1 0 0 755
## 155 6155 1 0 0 1 0 0 13186
## 171 6171 1 0 1 0 1 1 20752
## 199 6199 1 0 1 1 1 1 20173
## 197 6197 0 0 1 1 0 1 16944
## 31 6031 0 0 1 0 1 1 10229
## 176 6176 1 0 1 0 0 0 23390
## 104 6104 0 0 1 1 1 1 14046
## 181 6181 0 1 1 1 1 0 24356
## 32 6032 0 0 1 1 1 0 9894
## 54 6054 1 0 0 1 1 1 20903
## 198 6198 1 0 0 0 0 0 13894
## 3 6003 0 1 0 0 1 0 9437
## 70 6070 1 0 0 1 1 1 16285
## 189 6189 0 1 0 0 0 1 3754
## 158.1 6158 1 0 0 1 0 0 755
## 182 6182 1 0 0 0 1 1 1369
## 30 6030 0 0 0 1 1 0 2245
## 66 6066 0 0 1 1 1 1 9266
## 92 6092 0 1 0 0 1 0 22133
## 61 6061 0 0 0 1 1 0 1202
## 195 6195 1 0 0 1 0 1 21223
## 48 6048 0 1 1 0 0 1 4032
## 132 6132 1 1 0 1 1 1 12931
## 159 6159 0 0 1 1 1 1 17199
## 153.1 6153 0 0 0 1 0 1 13642
## 181.1 6181 0 1 1 1 1 0 24356
## 147 6147 1 0 1 0 1 0 8914
## 123 6123 1 1 1 1 0 0 9458
## 85 6085 0 1 1 1 1 1 14183
## 98 6098 0 1 1 0 0 0 23423
## 46 6046 1 0 0 1 1 0 24515
## 59 6059 1 0 1 0 0 1 18409
## 83 6083 1 1 0 0 1 1 15286
## 95 6095 1 1 0 0 0 1 4143
## 147.1 6147 1 0 1 0 1 0 8914
## 89 6089 1 1 1 0 1 1 18804
## 41 6041 1 1 1 1 1 1 997
## 52 6052 0 1 0 1 1 1 21240
## 93 6093 0 1 0 1 1 0 20926
## 87 6087 1 1 1 1 1 0 4265
## 119 6119 1 1 1 0 0 0 24288
## 92.1 6092 0 1 0 0 1 0 22133
## 187 6187 1 1 0 0 1 1 6766
## 4 6004 1 0 0 1 0 0 22130
##
## [[7]]$tae
## [1] 5
##
## [[7]]$tries
## [1] 25 21 11 23 31 16 16 19 22 18 22 24 29 25 14 25 26 9 12 18 26 16 16
## [24] 20 21 8 17 14 9 17 16 25 23 26 15 15 19 14 5 25
##
## [[7]]$error_states
## [1] 25 21 11 11 11 11 11 11 11 11 11 11 11 11 11 11 11 9 9 9 9 9 9
## [24] 9 9 8 8 8 8 8 8 8 8 8 8 8 8 8 5 5
##
##
## [[8]]
## [[8]]$areaid
## [1] 8308
##
## [[8]]$selection
## pid he female agemature car_owner house_owner working annualIncome
## 191 6191 1 1 1 0 1 0 24296
## 12 6012 1 0 1 0 1 1 7016
## 151 6151 0 0 1 0 0 0 18055
## 3 6003 0 1 0 0 1 0 9437
## 30 6030 0 0 0 1 1 0 2245
## 106 6106 1 1 1 1 0 0 11994
## 81 6081 0 1 0 0 0 0 19338
## 41 6041 1 1 1 1 1 1 997
## 105 6105 1 1 1 0 1 0 6742
## 73 6073 0 0 1 0 0 0 23682
## 193 6193 1 1 1 0 1 1 17777
## 112 6112 0 0 1 0 0 0 23382
## 129 6129 0 1 1 1 0 0 5539
## 45 6045 0 0 0 0 1 1 13377
## 121 6121 1 1 1 1 0 0 2352
## 94 6094 0 0 1 1 1 1 18017
## 22 6022 0 1 0 1 0 1 2871
## 93 6093 0 1 0 1 1 0 20926
## 74 6074 1 0 0 0 1 1 16511
## 12.1 6012 1 0 1 0 1 1 7016
## 115 6115 1 0 1 1 1 0 17012
## 2 6002 1 1 1 1 1 1 2458
## 12.2 6012 1 0 1 0 1 1 7016
## 54 6054 1 0 0 1 1 1 20903
## 79 6079 1 1 0 0 1 1 16117
## 102 6102 0 0 0 1 1 0 12107
## 2.1 6002 1 1 1 1 1 1 2458
## 68 6068 0 1 0 1 1 0 20739
## 168 6168 0 1 0 1 1 0 17495
## 167 6167 0 1 1 1 1 1 19045
## 5 6005 1 0 0 0 1 1 3936
## 59 6059 1 0 1 0 0 1 18409
## 144 6144 1 0 0 0 0 0 10261
## 3.1 6003 0 1 0 0 1 0 9437
## 153 6153 0 0 0 1 0 1 13642
## 30.1 6030 0 0 0 1 1 0 2245
## 33 6033 0 1 1 0 0 0 15822
## 79.1 6079 1 1 0 0 1 1 16117
## 108 6108 0 1 1 0 1 0 19083
## 171 6171 1 0 1 0 1 1 20752
## 145 6145 1 1 0 1 0 1 9894
## 48 6048 0 1 1 0 0 1 4032
## 114 6114 1 0 1 0 1 1 21039
## 81.1 6081 0 1 0 0 0 0 19338
## 47 6047 1 1 0 0 0 0 14501
## 120 6120 0 0 0 1 0 0 1025
## 142 6142 0 1 1 1 0 0 24860
## 44 6044 1 0 1 0 0 0 1844
## 151.1 6151 0 0 1 0 0 0 18055
## 82 6082 0 0 1 1 1 0 6233
## 179 6179 0 0 0 0 1 0 14058
## 152 6152 1 1 0 0 1 0 11947
## 117 6117 1 1 0 1 0 0 24753
## 198 6198 1 0 0 0 0 0 13894
## 98 6098 0 1 1 0 0 0 23423
## 185 6185 1 0 1 1 1 1 12862
## 101 6101 0 1 1 1 0 1 21034
## 188 6188 1 0 0 0 1 1 21100
## 36 6036 0 0 1 1 1 1 23880
## 49 6049 0 1 1 1 0 1 7873
## 91 6091 1 1 0 0 0 0 5695
## 94.1 6094 0 0 1 1 1 1 18017
## 89 6089 1 1 1 0 1 1 18804
## 131 6131 1 0 1 1 1 0 4416
## 155 6155 1 0 0 1 0 0 13186
## 47.1 6047 1 1 0 0 0 0 14501
## 118 6118 1 1 1 0 0 1 5610
## 48.1 6048 0 1 1 0 0 1 4032
## 62 6062 0 0 0 0 0 0 5375
## 4 6004 1 0 0 1 0 0 22130
## 110 6110 1 1 0 1 1 0 16503
## 73.1 6073 0 0 1 0 0 0 23682
## 21 6021 0 1 1 1 1 0 16588
## 49.1 6049 0 1 1 1 0 1 7873
## 43 6043 0 0 0 1 0 1 3170
## 80 6080 1 0 0 0 1 0 14511
## 145.1 6145 1 1 0 1 0 1 9894
## 165 6165 1 0 1 0 1 1 1125
##
## [[8]]$tae
## [1] 3
##
## [[8]]$tries
## [1] 5 9 6 8 5 10 11 14 6 6 10 9 21 11 12 12 13 8 6 13 4 13 9
## [24] 7 9 16 7 9 8 15 13 12 6 16 4 3 14 12 4 8
##
## [[8]]$error_states
## [1] 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4
## [36] 3 3 3 3 3
##
##
## [[9]]
## [[9]]$areaid
## [1] 8309
##
## [[9]]$selection
## pid he female agemature car_owner house_owner working annualIncome
## 85 6085 0 1 1 1 1 1 14183
## 61 6061 0 0 0 1 1 0 1202
## 117 6117 1 1 0 1 0 0 24753
## 26 6026 0 0 0 1 0 0 22364
## 46 6046 1 0 0 1 1 0 24515
## 4 6004 1 0 0 1 0 0 22130
## 181 6181 0 1 1 1 1 0 24356
## 111 6111 0 0 1 0 1 1 16075
## 16 6016 1 1 0 0 0 1 21714
## 118 6118 1 1 1 0 0 1 5610
## 37 6037 0 0 1 1 0 0 7093
## 108 6108 0 1 1 0 1 0 19083
## 5 6005 1 0 0 0 1 1 3936
## 123 6123 1 1 1 1 0 0 9458
## 197 6197 0 0 1 1 0 1 16944
## 189 6189 0 1 0 0 0 1 3754
## 46.1 6046 1 0 0 1 1 0 24515
## 99 6099 1 0 1 0 1 0 15471
## 197.1 6197 0 0 1 1 0 1 16944
## 132 6132 1 1 0 1 1 1 12931
## 107 6107 1 1 0 0 0 0 2070
## 70 6070 1 0 0 1 1 1 16285
## 148 6148 1 0 1 1 0 1 4717
## 103 6103 1 0 1 0 0 0 13352
## 38 6038 1 0 1 1 1 0 4508
## 72 6072 0 1 0 1 0 1 861
## 108.1 6108 0 1 1 0 1 0 19083
## 135 6135 1 0 1 0 0 1 14497
## 172 6172 1 0 1 0 0 1 14346
## 106 6106 1 1 1 1 0 0 11994
## 49 6049 0 1 1 1 0 1 7873
## 128 6128 0 1 0 0 1 1 1066
## 106.1 6106 1 1 1 1 0 0 11994
## 105 6105 1 1 1 0 1 0 6742
## 105.1 6105 1 1 1 0 1 0 6742
## 160 6160 1 1 1 1 0 1 13746
## 13 6013 1 0 1 1 0 1 20227
## 62 6062 0 0 0 0 0 0 5375
## 197.2 6197 0 0 1 1 0 1 16944
## 167 6167 0 1 1 1 1 1 19045
## 133 6133 1 0 0 1 1 1 16596
## 3 6003 0 1 0 0 1 0 9437
## 38.1 6038 1 0 1 1 1 0 4508
## 46.2 6046 1 0 0 1 1 0 24515
## 76 6076 1 1 1 0 0 0 14320
## 63 6063 0 0 1 1 0 0 22695
## 186 6186 0 0 0 0 1 1 21279
## 104 6104 0 0 1 1 1 1 14046
## 55 6055 0 0 1 1 0 0 9651
## 169 6169 1 0 1 1 1 0 3719
## 74 6074 1 0 0 0 1 1 16511
## 135.1 6135 1 0 1 0 0 1 14497
## 65 6065 1 0 0 1 0 0 11579
## 138 6138 0 0 0 0 1 0 728
## 197.3 6197 0 0 1 1 0 1 16944
## 133.1 6133 1 0 0 1 1 1 16596
## 55.1 6055 0 0 1 1 0 0 9651
## 191 6191 1 1 1 0 1 0 24296
## 44 6044 1 0 1 0 0 0 1844
## 18 6018 0 0 0 1 1 1 23950
## 46.3 6046 1 0 0 1 1 0 24515
## 82 6082 0 0 1 1 1 0 6233
## 99.1 6099 1 0 1 0 1 0 15471
## 158 6158 1 0 0 1 0 0 755
## 57 6057 1 0 1 1 1 1 24504
## 157 6157 0 0 1 1 1 1 18831
## 96 6096 1 0 0 0 1 0 15223
## 130 6130 1 1 1 1 1 1 7340
## 172.1 6172 1 0 1 0 0 1 14346
## 162 6162 1 1 1 0 1 1 8955
## 1 6001 0 1 1 1 0 1 16567
## 144 6144 1 0 0 0 0 0 10261
## 121 6121 1 1 1 1 0 0 2352
## 54 6054 1 0 0 1 1 1 20903
## 40 6040 1 1 0 0 0 0 13271
## 114 6114 1 0 1 0 1 1 21039
## 111.1 6111 0 0 1 0 1 1 16075
##
## [[9]]$tae
## [1] 27
##
## [[9]]$tries
## [1] 55 49 53 42 27 43 44 48 54 47 68 47 50 51 44 67 49 43 48 40 50 36 47
## [24] 48 49 42 47 54 38 54 43 50 44 53 41 44 57 50 33 41
##
## [[9]]$error_states
## [1] 55 49 49 42 27 27 27 27 27 27 27 27 27 27 27 27 27 27 27 27 27 27 27
## [24] 27 27 27 27 27 27 27 27 27 27 27 27 27 27 27 27 27
##
##
## [[10]]
## [[10]]$areaid
## [1] 8310
##
## [[10]]$selection
## pid he female agemature car_owner house_owner working annualIncome
## 46 6046 1 0 0 1 1 0 24515
## 123 6123 1 1 1 1 0 0 9458
## 72 6072 0 1 0 1 0 1 861
## 193 6193 1 1 1 0 1 1 17777
## 141 6141 0 0 0 0 0 1 8217
## 198 6198 1 0 0 0 0 0 13894
## 83 6083 1 1 0 0 1 1 15286
## 38 6038 1 0 1 1 1 0 4508
## 148 6148 1 0 1 1 0 1 4717
## 4 6004 1 0 0 1 0 0 22130
## 19 6019 0 1 1 1 0 0 4546
## 58 6058 1 0 0 0 0 0 22892
## 145 6145 1 1 0 1 0 1 9894
## 193.1 6193 1 1 1 0 1 1 17777
## 199 6199 1 0 1 1 1 1 20173
## 80 6080 1 0 0 0 1 0 14511
## 136 6136 0 1 1 0 0 0 1716
## 175 6175 1 0 0 0 0 0 4067
## 155 6155 1 0 0 1 0 0 13186
## 180 6180 0 1 0 0 1 1 14699
## 140 6140 0 1 0 1 0 1 17720
## 122 6122 0 1 0 1 0 1 5706
## 115 6115 1 0 1 1 1 0 17012
## 184 6184 0 0 1 1 1 0 18460
## 113 6113 0 0 0 1 1 1 18582
## 69 6069 1 1 1 1 1 1 15149
## 44 6044 1 0 1 0 0 0 1844
## 161 6161 0 1 0 1 0 0 23118
## 113.1 6113 0 0 0 1 1 1 18582
## 28 6028 1 0 0 0 0 0 12338
## 127 6127 0 1 0 1 1 0 12893
## 184.1 6184 0 0 1 1 1 0 18460
## 139 6139 1 1 0 1 0 0 24333
## 57 6057 1 0 1 1 1 1 24504
## 100 6100 0 1 1 1 1 0 19540
## 53 6053 1 0 1 0 0 0 8081
## 142 6142 0 1 1 1 0 0 24860
## 100.1 6100 0 1 1 1 1 0 19540
## 65 6065 1 0 0 1 0 0 11579
## 77 6077 1 1 1 0 1 1 3660
## 175.1 6175 1 0 0 0 0 0 4067
## 175.2 6175 1 0 0 0 0 0 4067
## 39 6039 0 0 0 0 1 1 3470
## 44.1 6044 1 0 1 0 0 0 1844
## 34 6034 1 0 0 1 1 0 8897
## 192 6192 1 0 1 0 1 1 1649
## 4.1 6004 1 0 0 1 0 0 22130
## 124 6124 0 0 0 1 0 1 15521
## 48 6048 0 1 1 0 0 1 4032
## 157 6157 0 0 1 1 1 1 18831
## 117 6117 1 1 0 1 0 0 24753
## 52 6052 0 1 0 1 1 1 21240
## 125 6125 0 1 0 1 0 1 16596
## 38.1 6038 1 0 1 1 1 0 4508
## 44.2 6044 1 0 1 0 0 0 1844
## 131 6131 1 0 1 1 1 0 4416
## 199.1 6199 1 0 1 1 1 1 20173
## 191 6191 1 1 1 0 1 0 24296
## 81 6081 0 1 0 0 0 0 19338
## 182 6182 1 0 0 0 1 1 1369
## 170 6170 0 1 1 0 0 0 12948
## 200 6200 0 1 1 0 1 1 12623
## 145.1 6145 1 1 0 1 0 1 9894
## 194 6194 1 1 1 0 1 0 15198
## 71 6071 0 0 0 1 0 0 23451
## 143 6143 0 0 1 0 0 0 4922
## 174 6174 0 0 1 1 0 0 6798
## 80.1 6080 1 0 0 0 1 0 14511
## 191.1 6191 1 1 1 0 1 0 24296
## 34.1 6034 1 0 0 1 1 0 8897
## 20 6020 1 0 0 1 0 1 1010
## 7 6007 0 0 1 1 1 0 9535
## 103 6103 1 0 1 0 0 0 13352
## 29 6029 0 1 1 1 0 0 22979
## 105 6105 1 1 1 0 1 0 6742
## 21 6021 0 1 1 1 1 0 16588
## 171 6171 1 0 1 0 1 1 20752
## 28.1 6028 1 0 0 0 0 0 12338
##
## [[10]]$tae
## [1] 15
##
## [[10]]$tries
## [1] 34 31 31 28 30 22 25 35 22 34 34 39 30 28 15 25 28 40 31 25 28 24 30
## [24] 29 29 26 30 32 26 34 34 15 20 33 21 42 39 29 25 23
##
## [[10]]$error_states
## [1] 34 31 31 28 28 22 22 22 22 22 22 22 22 22 15 15 15 15 15 15 15 15 15
## [24] 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15
##
##
## [[11]]
## [[11]]$areaid
## [1] 8311
##
## [[11]]$selection
## pid he female agemature car_owner house_owner working annualIncome
## 88 6088 0 0 0 0 1 0 23350
## 200 6200 0 1 1 0 1 1 12623
## 173 6173 0 1 1 1 1 1 7025
## 51 6051 0 0 1 1 0 1 6669
## 79 6079 1 1 0 0 1 1 16117
## 166 6166 0 1 1 0 1 1 15312
## 136 6136 0 1 1 0 0 0 1716
## 139 6139 1 1 0 1 0 0 24333
## 65 6065 1 0 0 1 0 0 11579
## 39 6039 0 0 0 0 1 1 3470
## 115 6115 1 0 1 1 1 0 17012
## 3 6003 0 1 0 0 1 0 9437
## 49 6049 0 1 1 1 0 1 7873
## 93 6093 0 1 0 1 1 0 20926
## 111 6111 0 0 1 0 1 1 16075
## 8 6008 0 0 1 1 0 1 1136
## 75 6075 0 0 0 0 0 0 5855
## 67 6067 1 1 1 1 0 1 19112
## 173.1 6173 0 1 1 1 1 1 7025
## 127 6127 0 1 0 1 1 0 12893
## 96 6096 1 0 0 0 1 0 15223
## 92 6092 0 1 0 0 1 0 22133
## 18 6018 0 0 0 1 1 1 23950
## 191 6191 1 1 1 0 1 0 24296
## 128 6128 0 1 0 0 1 1 1066
## 168 6168 0 1 0 1 1 0 17495
## 184 6184 0 0 1 1 1 0 18460
## 116 6116 0 1 0 0 1 0 4414
## 75.1 6075 0 0 0 0 0 0 5855
## 95 6095 1 1 0 0 0 1 4143
## 96.1 6096 1 0 0 0 1 0 15223
## 96.2 6096 1 0 0 0 1 0 15223
## 42 6042 0 0 1 1 0 1 19157
## 48 6048 0 1 1 0 0 1 4032
## 180 6180 0 1 0 0 1 1 14699
## 35 6035 1 1 1 0 1 0 23377
## 124 6124 0 0 0 1 0 1 15521
## 47 6047 1 1 0 0 0 0 14501
## 157 6157 0 0 1 1 1 1 18831
## 164 6164 0 1 0 0 1 1 21865
## 151 6151 0 0 1 0 0 0 18055
## 140 6140 0 1 0 1 0 1 17720
## 113 6113 0 0 0 1 1 1 18582
## 181 6181 0 1 1 1 1 0 24356
## 126 6126 1 1 1 0 1 1 12386
## 6 6006 1 0 0 1 0 1 16695
## 68 6068 0 1 0 1 1 0 20739
## 107 6107 1 1 0 0 0 0 2070
## 69 6069 1 1 1 1 1 1 15149
## 135 6135 1 0 1 0 0 1 14497
## 99 6099 1 0 1 0 1 0 15471
## 145 6145 1 1 0 1 0 1 9894
## 19 6019 0 1 1 1 0 0 4546
## 198 6198 1 0 0 0 0 0 13894
## 59 6059 1 0 1 0 0 1 18409
## 69.1 6069 1 1 1 1 1 1 15149
## 129 6129 0 1 1 1 0 0 5539
## 47.1 6047 1 1 0 0 0 0 14501
##
## [[11]]$tae
## [1] 30
##
## [[11]]$tries
## [1] 35 39 35 35 46 36 44 41 44 33 45 37 45 38 42 40 36 37 45 48 51 44 46
## [24] 33 36 44 43 38 39 40 36 30 36 42 39 41 35 32 44 44
##
## [[11]]$error_states
## [1] 35 35 35 35 35 35 35 35 35 33 33 33 33 33 33 33 33 33 33 33 33 33 33
## [24] 33 33 33 33 33 33 33 33 30 30 30 30 30 30 30 30 30
##
##
## [[12]]
## [[12]]$areaid
## [1] 8312
##
## [[12]]$selection
## pid he female agemature car_owner house_owner working annualIncome
## 174 6174 0 0 1 1 0 0 6798
## 43 6043 0 0 0 1 0 1 3170
## 118 6118 1 1 1 0 0 1 5610
## 10 6010 0 1 1 0 1 0 10119
## 129 6129 0 1 1 1 0 0 5539
## 37 6037 0 0 1 1 0 0 7093
## 172 6172 1 0 1 0 0 1 14346
## 142 6142 0 1 1 1 0 0 24860
## 9 6009 1 1 1 1 0 0 24170
## 101 6101 0 1 1 1 0 1 21034
## 183 6183 1 1 0 1 0 0 20091
## 100 6100 0 1 1 1 1 0 19540
## 67 6067 1 1 1 1 0 1 19112
## 58 6058 1 0 0 0 0 0 22892
## 6 6006 1 0 0 1 0 1 16695
## 87 6087 1 1 1 1 1 0 4265
## 93 6093 0 1 0 1 1 0 20926
## 164 6164 0 1 0 0 1 1 21865
## 111 6111 0 0 1 0 1 1 16075
## 31 6031 0 0 1 0 1 1 10229
## 45 6045 0 0 0 0 1 1 13377
## 154 6154 0 0 0 1 0 0 24662
## 92 6092 0 1 0 0 1 0 22133
## 189 6189 0 1 0 0 0 1 3754
## 195 6195 1 0 0 1 0 1 21223
## 123 6123 1 1 1 1 0 0 9458
## 80 6080 1 0 0 0 1 0 14511
## 147 6147 1 0 1 0 1 0 8914
## 132 6132 1 1 0 1 1 1 12931
## 95 6095 1 1 0 0 0 1 4143
## 43.1 6043 0 0 0 1 0 1 3170
## 7 6007 0 0 1 1 1 0 9535
## 88 6088 0 0 0 0 1 0 23350
## 18 6018 0 0 0 1 1 1 23950
## 45.1 6045 0 0 0 0 1 1 13377
## 158 6158 1 0 0 1 0 0 755
## 38 6038 1 0 1 1 1 0 4508
## 149 6149 0 0 0 0 0 1 21917
## 141 6141 0 0 0 0 0 1 8217
## 19 6019 0 1 1 1 0 0 4546
## 41 6041 1 1 1 1 1 1 997
## 130 6130 1 1 1 1 1 1 7340
## 150 6150 0 1 0 1 1 0 23632
## 187 6187 1 1 0 0 1 1 6766
## 161 6161 0 1 0 1 0 0 23118
## 135 6135 1 0 1 0 0 1 14497
## 183.1 6183 1 1 0 1 0 0 20091
## 68 6068 0 1 0 1 1 0 20739
## 26 6026 0 0 0 1 0 0 22364
## 31.1 6031 0 0 1 0 1 1 10229
## 57 6057 1 0 1 1 1 1 24504
## 200 6200 0 1 1 0 1 1 12623
## 106 6106 1 1 1 1 0 0 11994
## 92.1 6092 0 1 0 0 1 0 22133
## 56 6056 0 1 0 1 1 0 19291
## 3 6003 0 1 0 0 1 0 9437
## 87.1 6087 1 1 1 1 1 0 4265
## 162 6162 1 1 1 0 1 1 8955
## 62 6062 0 0 0 0 0 0 5375
## 188 6188 1 0 0 0 1 1 21100
## 45.2 6045 0 0 0 0 1 1 13377
## 137 6137 1 1 0 0 0 0 11019
## 120 6120 0 0 0 1 0 0 1025
## 10.1 6010 0 1 1 0 1 0 10119
## 16 6016 1 1 0 0 0 1 21714
## 122 6122 0 1 0 1 0 1 5706
## 19.1 6019 0 1 1 1 0 0 4546
## 196 6196 1 0 1 0 0 1 17683
##
## [[12]]$tae
## [1] 32
##
## [[12]]$tries
## [1] 55 33 42 40 42 42 42 46 39 41 39 43 38 49 34 41 48 41 40 33 57 44 40
## [24] 32 45 34 46 35 35 35 36 49 38 33 34 40 40 40 44 42
##
## [[12]]$error_states
## [1] 55 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33 33
## [24] 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32
##
##
## [[13]]
## [[13]]$areaid
## [1] 8313
##
## [[13]]$selection
## pid he female agemature car_owner house_owner working annualIncome
## 99 6099 1 0 1 0 1 0 15471
## 164 6164 0 1 0 0 1 1 21865
## 197 6197 0 0 1 1 0 1 16944
## 81 6081 0 1 0 0 0 0 19338
## 41 6041 1 1 1 1 1 1 997
## 144 6144 1 0 0 0 0 0 10261
## 133 6133 1 0 0 1 1 1 16596
## 46 6046 1 0 0 1 1 0 24515
## 87 6087 1 1 1 1 1 0 4265
## 22 6022 0 1 0 1 0 1 2871
## 113 6113 0 0 0 1 1 1 18582
## 104 6104 0 0 1 1 1 1 14046
## 44 6044 1 0 1 0 0 0 1844
## 29 6029 0 1 1 1 0 0 22979
## 61 6061 0 0 0 1 1 0 1202
## 57 6057 1 0 1 1 1 1 24504
## 40 6040 1 1 0 0 0 0 13271
## 33 6033 0 1 1 0 0 0 15822
## 123 6123 1 1 1 1 0 0 9458
## 167 6167 0 1 1 1 1 1 19045
## 10 6010 0 1 1 0 1 0 10119
## 125 6125 0 1 0 1 0 1 16596
## 190 6190 0 1 0 0 0 1 6535
## 52 6052 0 1 0 1 1 1 21240
## 160 6160 1 1 1 1 0 1 13746
## 198 6198 1 0 0 0 0 0 13894
## 175 6175 1 0 0 0 0 0 4067
## 11 6011 1 1 1 1 0 1 14821
## 70 6070 1 0 0 1 1 1 16285
## 169 6169 1 0 1 1 1 0 3719
## 14 6014 1 1 0 1 1 1 6014
## 58 6058 1 0 0 0 0 0 22892
## 61.1 6061 0 0 0 1 1 0 1202
## 145 6145 1 1 0 1 0 1 9894
## 157 6157 0 0 1 1 1 1 18831
## 46.1 6046 1 0 0 1 1 0 24515
## 97 6097 1 1 1 0 1 1 21521
## 162 6162 1 1 1 0 1 1 8955
## 196 6196 1 0 1 0 0 1 17683
## 53 6053 1 0 1 0 0 0 8081
## 188 6188 1 0 0 0 1 1 21100
## 169.1 6169 1 0 1 1 1 0 3719
## 118 6118 1 1 1 0 0 1 5610
## 41.1 6041 1 1 1 1 1 1 997
## 134 6134 0 0 0 1 1 0 13271
## 178 6178 0 1 0 1 0 1 5271
## 140 6140 0 1 0 1 0 1 17720
## 162.1 6162 1 1 1 0 1 1 8955
## 100 6100 0 1 1 1 1 0 19540
## 79 6079 1 1 0 0 1 1 16117
## 59 6059 1 0 1 0 0 1 18409
## 178.1 6178 0 1 0 1 0 1 5271
## 163 6163 0 1 0 0 1 1 9332
## 123.1 6123 1 1 1 1 0 0 9458
## 151 6151 0 0 1 0 0 0 18055
## 185 6185 1 0 1 1 1 1 12862
## 11.1 6011 1 1 1 1 0 1 14821
## 4 6004 1 0 0 1 0 0 22130
## 41.2 6041 1 1 1 1 1 1 997
## 123.2 6123 1 1 1 1 0 0 9458
##
## [[13]]$tae
## [1] 37
##
## [[13]]$tries
## [1] 48 44 47 42 57 38 44 46 53 43 53 47 54 43 53 48 41 48 44 46 43 52 41
## [24] 47 51 59 43 38 48 55 55 52 51 51 58 37 43 39 40 49
##
## [[13]]$error_states
## [1] 48 44 44 42 42 38 38 38 38 38 38 38 38 38 38 38 38 38 38 38 38 38 38
## [24] 38 38 38 38 38 38 38 38 38 38 38 38 37 37 37 37 37
##
##
##
## Slot "iterations":
## [1] 40length( try01@results )## [1] 13The results object, holds the results for 4 areas. For each area, the microsimulation object holds a number of information. Let’s examine the third area of our microsimulation.
area03 = try01@results[[3]]
area03 $areaid
[1] 8303
$selection
pid he female agemature car_owner house_owner working annualIncome
158 6158 1 0 0 1 0 0 755
35 6035 1 1 1 0 1 0 23377
94 6094 0 0 1 1 1 1 18017
150 6150 0 1 0 1 1 0 23632
177 6177 1 0 0 1 1 0 17102
45 6045 0 0 0 0 1 1 13377
94.1 6094 0 0 1 1 1 1 18017
96 6096 1 0 0 0 1 0 15223
151 6151 0 0 1 0 0 0 18055
174 6174 0 0 1 1 0 0 6798
174.1 6174 0 0 1 1 0 0 6798
164 6164 0 1 0 0 1 1 21865
139 6139 1 1 0 1 0 0 24333
143 6143 0 0 1 0 0 0 4922
2 6002 1 1 1 1 1 1 2458
170 6170 0 1 1 0 0 0 12948
98 6098 0 1 1 0 0 0 23423
188 6188 1 0 0 0 1 1 21100
141 6141 0 0 0 0 0 1 8217
184 6184 0 0 1 1 1 0 18460
59 6059 1 0 1 0 0 1 18409
122 6122 0 1 0 1 0 1 5706
45.1 6045 0 0 0 0 1 1 13377
59.1 6059 1 0 1 0 0 1 18409
195 6195 1 0 0 1 0 1 21223
71 6071 0 0 0 1 0 0 23451
60 6060 0 0 0 1 0 0 7352
15 6015 0 1 1 0 0 0 19546
108 6108 0 1 1 0 1 0 19083
130 6130 1 1 1 1 1 1 7340
103 6103 1 0 1 0 0 0 13352
84 6084 0 1 0 1 0 0 24604
178 6178 0 1 0 1 0 1 5271
153 6153 0 0 0 1 0 1 13642
147 6147 1 0 1 0 1 0 8914
48 6048 0 1 1 0 0 1 4032
27 6027 1 0 1 1 0 1 1347
153.1 6153 0 0 0 1 0 1 13642
161 6161 0 1 0 1 0 0 23118
130.1 6130 1 1 1 1 1 1 7340
33 6033 0 1 1 0 0 0 15822
184.1 6184 0 0 1 1 1 0 18460
71.1 6071 0 0 0 1 0 0 23451
90 6090 0 1 0 1 1 1 20528
18 6018 0 0 0 1 1 1 23950
33.1 6033 0 1 1 0 0 0 15822
154 6154 0 0 0 1 0 0 24662
188.1 6188 1 0 0 0 1 1 21100
32 6032 0 0 1 1 1 0 9894
13 6013 1 0 1 1 0 1 20227
144 6144 1 0 0 0 0 0 10261
16 6016 1 1 0 0 0 1 21714
159 6159 0 0 1 1 1 1 17199
112 6112 0 0 1 0 0 0 23382
187 6187 1 1 0 0 1 1 6766
167 6167 0 1 1 1 1 1 19045
161.1 6161 0 1 0 1 0 0 23118
78 6078 0 1 0 0 0 1 12626
$tae
[1] 22
$tries
[1] 42 35 22 39 39 32 41 28 27 43 29 36 32 39 33 39 31 36 38 36 36 39 35
[24] 39 37 41 39 27 27 33 41 24 31 33 25 30 28 34 35 41
$error_states
[1] 42 35 22 22 22 22 22 22 22 22 22 22 22 22 22 22 22 22 22 22 22 22 22
[24] 22 22 22 22 22 22 22 22 22 22 22 22 22 22 22 22 22As can be see, the area03\(areaid_ is the name of this area. The _area03\)selection object is the generated microdata of the area is the following:
area03$selection pid he female agemature car_owner house_owner working annualIncome
158 6158 1 0 0 1 0 0 755
35 6035 1 1 1 0 1 0 23377
94 6094 0 0 1 1 1 1 18017
150 6150 0 1 0 1 1 0 23632
177 6177 1 0 0 1 1 0 17102
45 6045 0 0 0 0 1 1 13377
94.1 6094 0 0 1 1 1 1 18017
96 6096 1 0 0 0 1 0 15223
151 6151 0 0 1 0 0 0 18055
174 6174 0 0 1 1 0 0 6798
174.1 6174 0 0 1 1 0 0 6798
164 6164 0 1 0 0 1 1 21865
139 6139 1 1 0 1 0 0 24333
143 6143 0 0 1 0 0 0 4922
2 6002 1 1 1 1 1 1 2458
170 6170 0 1 1 0 0 0 12948
98 6098 0 1 1 0 0 0 23423
188 6188 1 0 0 0 1 1 21100
141 6141 0 0 0 0 0 1 8217
184 6184 0 0 1 1 1 0 18460
59 6059 1 0 1 0 0 1 18409
122 6122 0 1 0 1 0 1 5706
45.1 6045 0 0 0 0 1 1 13377
59.1 6059 1 0 1 0 0 1 18409
195 6195 1 0 0 1 0 1 21223
71 6071 0 0 0 1 0 0 23451
60 6060 0 0 0 1 0 0 7352
15 6015 0 1 1 0 0 0 19546
108 6108 0 1 1 0 1 0 19083
130 6130 1 1 1 1 1 1 7340
103 6103 1 0 1 0 0 0 13352
84 6084 0 1 0 1 0 0 24604
178 6178 0 1 0 1 0 1 5271
153 6153 0 0 0 1 0 1 13642
147 6147 1 0 1 0 1 0 8914
48 6048 0 1 1 0 0 1 4032
27 6027 1 0 1 1 0 1 1347
153.1 6153 0 0 0 1 0 1 13642
161 6161 0 1 0 1 0 0 23118
130.1 6130 1 1 1 1 1 1 7340
33 6033 0 1 1 0 0 0 15822
184.1 6184 0 0 1 1 1 0 18460
71.1 6071 0 0 0 1 0 0 23451
90 6090 0 1 0 1 1 1 20528
18 6018 0 0 0 1 1 1 23950
33.1 6033 0 1 1 0 0 0 15822
154 6154 0 0 0 1 0 0 24662
188.1 6188 1 0 0 0 1 1 21100
32 6032 0 0 1 1 1 0 9894
13 6013 1 0 1 1 0 1 20227
144 6144 1 0 0 0 0 0 10261
16 6016 1 1 0 0 0 1 21714
159 6159 0 0 1 1 1 1 17199
112 6112 0 0 1 0 0 0 23382
187 6187 1 1 0 0 1 1 6766
167 6167 0 1 1 1 1 1 19045
161.1 6161 0 1 0 1 0 0 23118
78 6078 0 1 0 0 0 1 12626which can be exported to a csv file with the following command:
write.table(area03$selection, "area03Microdata.csv", sep=",")The tae oject is the Total Absolute Error of the fitting process for this area. The tries object holds the TAE of the random selections during the fitting process. This is the error of each combination, even if the combination isn’t selected. Finally, the area03$error_states holds the progress of the TAE during the iterations. The length of this array, is the same as the number of iterations and it is used mainly for tracking and then plotting the progress of the TAE. The last two arrays, can be plotted in order to examine the progress of the fitting process:
plotTries(insms=try01, 3)
The above plot, shows the progress of TAE during the fitting process. The horizontal axis represents the iterations and the vertical the TAE. Each red circle, represents the TAE of a randomly selected combination of individuals, in each iteration. The dark filled dots, represent, the TAE of the currently selected best combination of individuals. This graph is very useful when the fitting algorithm is Simulated Annealing as it depicts the tolerance of the algorithm during the first iterations.
This was a small tutorial on preparing microdata with the use of R and the sms library. Please fill free to contact me regarding the sms library.