Luttinen HW 6

#load libraries
library(haven)
library(gtsummary)
library(MASS)

## Warning: package 'MASS' was built under R version 4.1.3

## 
## Attaching package: 'MASS'

## The following object is masked from 'package:gtsummary':
## 
##     select

library(survey)

## Loading required package: grid

## Loading required package: Matrix

## Loading required package: survival

## 
## Attaching package: 'survey'

## The following object is masked from 'package:graphics':
## 
##     dotchart

#read in data
uganda16 <- read_dta("C:/Users/rlutt/Downloads/UGIR7BFL.DTA")
uganda16<-zap_labels(uganda16)

#recodes
#bank account
uganda16$bank<-car::Recode(uganda16$v170, recodes= "0= 'No'; 1= 'Yes'")
#age groups
uganda16$agegroup <- car::Recode(uganda16$v013, recodes = "1='15-19'; 2='20-24';3='25-29'; 4 = '30-34';5='35-39';6='40-44';7='45+' ", as.factor=T)
#fertility preferences
uganda16$wantanotherchild<-ifelse(uganda16$v602!=9&uganda16$v602==1,1,0)

#education level
uganda16$educationlevel <- car::recode(uganda16$v106, 
                            recodes = "0 = 'none'; 1 = 'primary'; 2:3='secondary and above' ",
                            as.factor=T)
#internet
uganda16$internet<-as.factor(uganda16$v171a)
uganda16$internet<-car::Recode(uganda16$v171a, recodes= "0='never'; 1= 'in the past year'; 2='over a year ago'; 3= 'yes, but unsure when'", as.factor=T)

#contraception
uganda16$contraception<-as.factor(uganda16$v313)
uganda16$contraception<-car::Recode(uganda16$v313, recodes= "0='none'; 1:2='folkloric method'; 3= 'modern method'", as.factor=T)

#usingmoderncontraception
uganda16$moderncon<-car::Recode(uganda16$v313, recodes= "0='1'; 1='2'; 2='2'; 3='2'", as.factor=F)

as.numeric(uganda16$moderncon)

##     [1] 1 1 1 1 1 1 1 1 2 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 2 1 2 1 1 1
##    [37] 2 1 2 1 2 1 2 2 2 1 1 1 1 2 1 1 1 1 1 1 1 1 1 2 1 2 1 1 1 1 2 1 2 1 1 1
##    [73] 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 2 2 2 2 2 2 1 1 1 1 2 1 1 1 1 1 1 1 1 1 2
##   [109] 1 1 2 2 2 1 2 1 2 2 1 1 1 1 1 1 2 2 1 1 1 1 1 1 2 1 1 1 2 2 1 1 1 1 1 2
##   [145] 1 1 2 2 2 2 1 2 1 1 2 1 1 1 1 1 1 1 2 1 1 1 1 1 2 1 1 1 2 1 1 2 1 2 1 2
##   [181] 1 1 1 1 1 2 1 2 1 2 1 2 1 1 1 1 1 1 1 1 1 2 1 1 1 1 2 1 1 1 1 2 1 1 1 2
##   [217] 2 2 1 1 2 1 1 1 1 1 1 2 1 1 2 1 1 1 1 1 1 2 1 1 1 1 1 1 1 2 2 1 1 1 1 1
##   [253] 2 1 1 1 1 1 1 2 2 1 1 1 1 2 1 1 1 2 1 1 2 1 1 1 1 1 1 2 1 2 2 2 1 2 1 1
##   [289] 1 2 1 1 1 2 1 1 1 1 1 2 1 1 2 1 1 1 1 2 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1 2
##   [325] 1 1 1 2 1 1 1 1 1 1 2 1 1 2 1 1 1 1 1 1 1 2 1 2 1 1 1 1 1 2 1 1 1 1 2 2
##   [361] 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 1 1 1 2 2 1 2 1 1 1 2 1 2 1
##   [397] 1 1 1 2 1 1 1 1 1 1 2 1 2 1 2 1 1 1 1 2 1 2 2 1 1 1 1 1 1 2 1 2 1 1 1 2
##   [433] 1 1 1 2 1 1 2 1 1 1 1 1 1 1 2 2 2 1 2 1 1 1 2 1 1 1 1 1 2 2 2 2 1 2 1 1
##   [469] 1 2 2 2 2 1 2 2 1 1 2 1 2 1 1 1 1 1 1 1 1 2 2 1 1 1 2 2 2 1 2 2 1 2 1 1
##   [505] 1 1 1 2 2 1 1 1 2 1 1 1 1 2 1 1 1 1 1 1 1 1 1 2 2 2 1 2 1 2 2 2 1 1 1 2
##   [541] 2 1 2 1 1 1 1 1 1 2 1 1 1 1 2 2 1 2 2 1 2 1 1 2 2 2 2 2 1 1 2 1 2 1 1 1
##   [577] 1 1 1 1 1 1 1 2 1 1 1 1 1 1 2 2 1 1 1 1 1 2 2 1 1 1 2 2 1 1 2 1 1 2 2 2
##   [613] 2 1 2 2 1 1 1 1 2 1 2 2 1 2 2 2 1 2 1 1 1 2 1 2 2 1 2 1 2 2 1 1 1 1 1 1
##   [649] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 2 2 2 1 1 1 1 1 1 1
##   [685] 1 1 1 1 1 1 1 2 1 1 1 1 1 2 1 1 1 1 1 1 1 2 2 2 1 2 2 2 2 1 2 1 2 1 1 1
##   [721] 2 1 1 2 2 1 1 2 1 2 2 1 1 1 2 1 2 1 1 1 2 2 1 2 1 2 2 1 1 2 1 2 1 1 1 1
##   [757] 1 1 2 2 1 2 2 1 2 1 1 1 1 1 1 2 1 1 1 1 1 1 2 1 2 1 1 1 2 2 1 2 2 1 1 1
##   [793] 1 2 1 1 1 1 1 2 1 1 1 2 2 1 1 2 1 1 2 1 2 1 1 1 1 2 1 1 2 2 2 1 1 1 1 1
##   [829] 1 2 2 1 1 2 1 1 1 1 1 1 1 2 1 2 2 1 2 1 2 1 1 1 1 2 1 2 1 2 1 1 1 2 1 2
##   [865] 1 2 1 1 1 2 2 2 2 1 2 1 1 2 1 1 1 2 2 1 2 1 2 1 1 1 1 2 1 1 1 2 1 1 2 1
##   [901] 1 1 2 2 1 1 1 1 1 2 2 2 1 2 1 1 1 2 1 1 1 1 1 1 1 1 2 2 1 2 2 2 1 1 1 2
##   [937] 1 1 1 1 1 1 1 2 2 1 1 2 1 1 1 1 1 1 1 1 1 2 1 2 2 2 1 1 1 1 1 1 1 1 1 2
##   [973] 1 2 1 1 1 1 2 1 1 1 1 2 1 1 1 1 1 1 2 1 1 1 1 1 2 1 1 1 2 2 2 1 1 2 2 1
##  [1009] 1 1 2 1 1 1 2 2 2 1 1 2 1 2 2 2 2 2 1 1 1 1 1 1 1 1 1 1 2 1 1 1 2 1 1 2
##  [1045] 1 2 1 2 2 1 1 1 1 2 2 2 1 1 2 1 1 2 2 2 2 1 1 1 2 2 1 1 1 2 1 1 1 2 1 2
##  [1081] 2 2 1 1 1 1 2 1 2 2 1 1 2 1 2 1 1 2 1 2 1 1 1 1 2 2 1 1 1 1 1 1 2 2 1 2
##  [1117] 1 2 2 1 2 2 1 2 2 1 1 1 2 1 1 2 1 1 2 1 2 1 1 1 1 2 1 2 1 2 2 1 1 1 1 1
##  [1153] 2 1 2 1 1 1 1 1 2 1 2 1 1 2 1 2 2 1 1 1 2 2 2 2 1 2 1 1 2 1 1 1 2 2 1 1
##  [1189] 2 1 1 1 1 2 2 2 1 2 1 2 1 1 1 2 1 1 1 1 1 1 1 1 1 1 2 1 2 2 2 1 1 1 2 2
##  [1225] 1 1 1 1 1 2 1 2 1 2 1 1 2 1 1 1 1 1 1 2 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 2
##  [1261] 2 2 1 2 1 1 1 1 1 2 2 1 2 2 2 2 1 1 1 2 1 2 2 1 1 2 1 1 1 2 2 2 1 2 1 2
##  [1297] 1 2 1 1 2 2 2 2 2 1 2 1 2 2 1 1 2 1 2 2 1 2 1 2 1 1 1 2 1 1 2 1 1 2 2 2
##  [1333] 2 2 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 2 1 2 2 1 1 1 1 1 1 2 1 2 2 2 1 2 1 1
##  [1369] 1 1 1 2 1 2 2 1 1 2 2 1 1 2 1 1 1 1 2 1 1 1 1 2 2 1 2 1 1 1 2 1 1 2 1 2
##  [1405] 2 1 1 2 1 1 1 2 1 2 1 1 1 2 2 1 1 1 2 1 1 1 2 1 1 1 1 1 1 1 1 2 1 2 1 1
##  [1441] 1 1 1 2 1 2 1 1 2 1 2 2 2 1 1 1 1 2 1 1 2 1 1 2 1 1 2 2 1 1 2 1 2 1 1 2
##  [1477] 1 1 2 2 1 1 2 1 1 1 1 2 1 1 1 1 1 1 1 2 1 1 1 2 1 2 1 1 1 1 1 1 1 2 1 2
##  [1513] 1 2 2 1 1 2 1 1 2 1 1 1 1 2 1 2 2 1 1 1 1 1 1 2 1 2 1 1 1 1 1 1 2 1 2 1
##  [1549] 2 1 1 1 1 2 2 2 1 2 1 1 2 1 1 1 1 1 1 2 1 2 2 1 1 1 1 1 2 2 1 1 1 1 2 1
##  [1585] 1 1 2 1 2 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 2 2 1 2 2 1 1 1 2 2 2 2 1 2 1
##  [1621] 1 1 1 1 1 2 2 2 1 1 2 1 2 1 1 2 1 1 2 2 2 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1
##  [1657] 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 2 2 1 1 1 1 1 1
##  [1693] 2 2 2 1 2 2 2 2 1 1 1 2 1 1 1 1 1 2 2 2 1 1 1 2 1 1 1 1 1 2 2 1 1 1 1 1
##  [1729] 1 1 1 2 1 1 1 2 2 2 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 1
##  [1765] 1 2 1 2 1 1 1 1 1 2 2 1 1 2 1 2 1 2 1 1 1 1 2 1 1 1 2 2 2 1 1 1 1 1 2 2
##  [1801] 1 2 2 2 2 1 1 1 1 1 1 1 1 2 1 1 1 2 1 2 2 2 2 2 2 1 2 1 2 1 2 1 2 1 1 1
##  [1837] 1 2 2 2 1 2 1 1 1 1 1 1 1 2 1 1 2 1 1 1 1 2 1 2 2 2 2 1 2 2 2 2 2 2 1 1
##  [1873] 1 1 1 1 1 1 1 1 2 1 1 1 1 1 2 1 2 1 1 2 1 1 2 2 2 2 1 1 1 2 1 1 1 1 1 1
##  [1909] 1 1 1 2 1 1 1 2 2 1 1 2 1 2 1 1 1 1 1 2 1 2 2 1 1 1 2 2 1 2 1 1 2 2 1 1
##  [1945] 1 1 1 1 2 2 1 1 1 1 1 1 1 1 1 1 1 2 2 1 1 2 1 1 2 1 1 2 1 2 1 1 2 1 1 2
##  [1981] 2 1 2 1 1 1 2 1 1 2 2 2 1 1 2 2 1 1 1 1 1 1 1 2 2 1 1 1 2 2 1 1 1 1 1 1
##  [2017] 1 1 1 1 2 1 2 1 1 1 1 1 1 2 2 2 2 1 2 2 2 2 1 2 1 2 2 2 2 2 2 1 2 2 1 1
##  [2053] 1 2 1 1 1 1 2 1 2 2 2 1 2 1 1 1 1 2 1 1 1 2 1 1 1 1 1 1 1 1 1 1 2 1 2 1
##  [2089] 2 1 1 1 1 1 2 1 2 2 1 1 1 1 1 1 2 1 2 2 1 1 1 1 2 1 2 1 1 1 2 1 1 1 2 2
##  [2125] 2 1 1 2 1 1 2 2 2 1 2 2 2 1 1 2 1 2 1 2 1 1 2 1 1 1 2 2 2 1 2 1 2 1 1 1
##  [2161] 1 1 1 2 1 1 1 2 2 2 1 2 1 2 1 1 1 1 1 1 1 1 1 1 2 2 2 1 1 1 2 1 1 1 1 1
##  [2197] 1 1 2 1 1 2 2 1 1 2 1 2 1 2 1 1 1 1 1 1 1 1 2 2 2 1 1 1 1 1 1 1 2 2 1 1
##  [2233] 1 2 1 1 1 1 2 1 1 2 1 1 2 2 2 1 2 2 1 1 1 1 2 1 1 1 1 2 2 1 1 1 1 2 1 1
##  [2269] 1 1 1 1 2 1 2 1 1 2 1 1 2 2 1 1 1 1 1 1 2 1 2 2 2 1 2 2 2 2 1 2 1 1 2 1
##  [2305] 1 1 1 1 1 2 2 1 2 1 1 2 2 2 1 2 1 2 2 2 1 1 2 1 1 2 1 2 1 1 2 2 1 1 1 1
##  [2341] 1 2 2 1 1 1 2 1 1 1 2 1 2 2 2 1 2 2 1 2 1 2 1 1 1 1 1 2 1 1 2 2 1 1 2 1
##  [2377] 1 1 2 2 2 2 1 2 1 1 1 1 1 1 1 1 1 1 2 1 1 2 1 2 1 1 1 1 2 2 2 2 1 2 1 1
##  [2413] 2 2 2 1 1 1 2 2 1 2 2 2 2 2 2 1 2 1 1 1 1 1 2 1 1 1 2 2 1 2 2 1 2 1 2 1
##  [2449] 1 1 1 1 2 2 2 1 1 1 1 1 1 1 2 2 1 1 2 1 1 1 1 1 1 2 2 1 1 2 1 1 2 2 2 2
##  [2485] 2 1 1 1 2 2 1 1 1 1 2 1 1 2 1 2 2 2 2 2 1 2 1 2 1 1 1 1 1 1 1 2 1 1 2 1
##  [2521] 1 1 1 1 2 2 1 1 1 1 1 1 1 2 2 2 1 1 1 1 2 1 1 1 1 1 2 2 2 1 2 2 1 1 2 1
##  [2557] 2 2 1 1 2 1 1 1 2 1 2 1 1 2 1 2 2 1 1 1 1 1 1 1 2 1 2 1 1 1 2 1 1 1 1 1
##  [2593] 2 1 1 2 2 2 2 1 1 1 1 2 1 2 2 1 1 1 1 1 2 1 1 1 1 1 1 2 1 2 1 1 1 1 1 2
##  [2629] 1 1 2 1 1 1 2 1 2 2 1 1 2 2 1 2 1 1 2 1 1 1 2 2 2 1 2 2 1 1 1 2 1 2 2 1
##  [2665] 2 1 2 2 1 1 1 1 2 2 1 1 1 1 1 2 1 2 1 1 1 1 1 2 1 1 1 2 1 1 1 1 1 1 2 1
##  [2701] 2 1 2 1 1 1 2 1 1 1 1 2 1 1 1 1 2 1 1 2 2 1 1 1 1 1 1 1 1 2 1 2 1 2 2 1
##  [2737] 1 1 2 2 1 2 1 1 1 2 1 2 1 2 1 2 1 2 1 1 1 1 1 1 1 1 1 2 1 1 1 2 1 1 1 2
##  [2773] 1 1 1 2 1 2 1 2 1 1 1 1 1 1 2 1 2 2 1 2 1 1 2 1 1 1 2 1 2 1 1 2 2 2 1 2
##  [2809] 2 1 2 2 1 2 1 1 2 1 1 1 1 2 1 2 2 1 1 1 1 2 2 1 1 1 1 2 2 1 1 1 1 2 2 1
##  [2845] 2 1 1 1 2 1 1 1 1 2 1 1 1 1 2 2 1 1 1 2 2 1 1 1 1 2 2 2 2 1 1 1 2 1 2 2
##  [2881] 1 2 1 2 2 2 2 1 1 1 1 1 1 1 1 2 1 2 1 1 1 1 2 1 2 1 1 2 1 1 1 1 1 2 1 1
##  [2917] 2 1 2 1 1 1 1 2 2 2 2 1 2 1 1 1 1 2 1 1 2 1 1 1 1 2 2 1 2 2 1 2 1 2 1 1
##  [2953] 1 2 2 2 1 1 2 2 1 2 2 2 1 2 1 1 1 1 1 1 1 1 2 1 1 1 1 1 2 1 2 2 1 1 2 2
##  [2989] 1 2 1 1 1 1 1 1 2 1 2 2 1 1 2 1 1 1 1 2 1 1 1 1 1 2 1 1 1 1 2 1 1 1 1 2
##  [3025] 1 1 2 1 1 1 2 1 1 1 2 1 2 1 2 2 1 1 2 1 2 1 1 1 2 1 2 1 1 1 1 1 2 1 2 1
##  [3061] 1 1 1 1 2 2 1 1 2 2 2 2 1 1 1 2 2 1 1 1 1 2 2 2 2 1 1 1 2 2 1 2 1 1 1 2
##  [3097] 2 2 2 1 1 1 1 1 2 2 2 2 1 2 2 2 1 1 1 2 1 1 2 2 2 1 1 1 2 2 1 2 2 2 2 1
##  [3133] 1 1 1 2 1 2 1 1 1 2 2 1 1 1 1 1 2 1 1 2 2 1 2 2 1 1 2 1 1 1 2 1 1 2 1 2
##  [3169] 2 2 2 1 2 2 2 1 1 1 1 2 1 2 1 1 1 2 2 1 1 1 2 1 1 1 1 1 1 1 1 2 1 1 1 1
##  [3205] 2 1 1 2 1 1 1 2 1 1 1 1 2 1 2 1 1 1 1 1 1 1 2 1 1 1 1 2 1 1 1 1 1 1 2 2
##  [3241] 1 1 2 2 2 2 1 1 1 1 2 1 2 1 1 1 1 2 2 2 2 1 1 1 1 2 1 2 1 2 1 1 2 1 1 1
##  [3277] 1 1 1 1 2 1 1 1 1 1 1 2 1 1 1 1 1 1 1 2 1 2 1 2 1 1 1 2 1 1 1 1 2 1 1 1
##  [3313] 1 1 1 1 2 2 1 1 2 2 1 2 2 1 2 1 1 1 2 1 1 1 1 2 1 1 1 1 1 1 1 1 1 2 1 1
##  [3349] 1 1 1 1 1 1 1 1 2 1 2 2 1 1 2 2 1 2 1 2 1 2 2 1 2 2 1 1 1 2 1 2 1 1 2 1
##  [3385] 2 1 1 2 1 1 2 1 1 2 2 2 1 1 1 1 2 2 2 1 1 2 1 2 1 2 1 1 1 1 2 1 2 1 1 1
##  [3421] 1 1 1 2 1 1 2 1 1 1 2 2 1 2 1 1 2 2 1 2 1 1 2 1 1 2 2 1 1 1 2 2 1 1 2 2
##  [3457] 1 1 1 1 1 1 1 1 1 1 1 2 1 1 2 2 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 1 2 1 1 1
##  [3493] 1 1 1 1 1 2 1 2 1 1 2 1 2 1 1 1 2 1 1 1 2 1 1 1 1 2 2 2 2 1 2 1 1 2 2 1
##  [3529] 2 1 1 1 2 2 1 1 1 2 2 2 2 2 1 2 1 2 1 2 2 2 1 1 2 1 1 1 2 2 1 1 1 2 1 1
##  [3565] 2 2 2 2 1 1 2 1 1 1 1 1 1 1 1 2 1 2 1 1 1 2 1 1 1 1 2 1 2 2 1 2 2 2 1 1
##  [3601] 2 1 2 1 2 1 1 1 1 2 2 2 2 1 2 1 2 1 2 1 1 2 1 2 1 1 2 1 2 2 1 2 2 1 2 1
##  [3637] 1 1 1 1 1 1 2 2 1 1 1 1 1 1 2 1 2 1 1 2 1 1 2 1 1 2 1 2 2 1 1 1 2 2 1 2
##  [3673] 2 1 1 1 1 1 1 2 1 1 1 1 1 2 1 2 1 2 2 2 1 1 2 1 1 2 1 1 1 2 2 2 1 1 1 2
##  [3709] 1 1 1 1 1 2 1 1 2 1 1 2 1 1 2 2 2 2 1 1 2 2 1 1 1 1 1 1 1 2 2 1 1 1 1 2
##  [3745] 1 2 1 1 1 1 1 1 1 2 1 2 1 1 2 1 2 2 1 2 1 1 1 1 2 2 1 2 2 1 2 2 2 1 1 1
##  [3781] 2 2 2 2 1 2 1 2 1 2 2 2 1 1 1 1 2 1 2 2 2 2 1 1 1 1 1 1 1 1 1 1 1 2 2 1
##  [3817] 1 2 2 2 2 1 2 1 1 2 1 1 1 1 2 1 2 1 1 1 1 1 1 2 1 1 1 1 1 1 1 1 2 2 2 1
##  [3853] 2 1 1 1 1 1 1 1 2 1 1 1 2 2 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
##  [3889] 1 1 1 1 1 1 1 1 1 2 2 1 1 1 2 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 1 1 1
##  [3925] 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 2 1 2 2 2 1 1 1 1 1 1 1 1 1 2 1
##  [3961] 1 1 1 2 1 1 2 1 2 1 1 2 1 2 2 1 1 1 1 1 1 1 2 1 1 1 1 2 1 2 2 1 1 1 1 2
##  [3997] 2 2 1 1 1 1 1 1 2 1 1 1 1 2 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 2
##  [4033] 2 1 1 1 1 1 1 1 1 1 1 1 2 2 1 1 2 1 2 1 1 1 1 1 1 1 2 2 1 2 1 1 1 2 1 2
##  [4069] 1 1 1 1 1 2 2 1 2 2 1 2 2 2 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 1 2 2 1 2 1 2
##  [4105] 1 1 2 1 2 1 1 2 2 1 2 1 1 2 1 1 1 1 1 2 2 1 1 1 1 1 1 1 1 1 2 1 1 1 1 2
##  [4141] 1 1 1 2 2 2 2 2 1 1 2 1 1 1 1 2 1 1 2 1 1 2 1 1 1 1 1 1 1 1 1 1 1 2 2 1
##  [4177] 1 1 1 2 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 2 1 1 1 2 1 1 1 1 1 2 1 1 1 2 1
##  [4213] 1 2 1 1 2 2 2 1 1 1 1 1 2 1 1 1 1 1 1 1 2 1 1 1 1 1 2 2 1 1 2 1 2 1 1 2
##  [4249] 1 1 1 1 1 1 1 1 1 2 1 1 1 1 2 1 2 1 1 1 1 2 1 2 1 1 1 2 2 1 1 1 1 1 1 1
##  [4285] 2 2 1 2 1 2 2 2 2 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 1 1 2 1 2 1 1 2 2 1 1
##  [4321] 2 2 1 1 2 1 2 2 1 1 1 1 1 1 1 1 2 1 1 1 2 2 2 1 1 1 2 1 1 1 2 1 2 1 2 2
##  [4357] 1 1 2 2 1 2 1 1 1 1 1 2 1 1 1 1 1 1 2 1 1 1 1 2 1 1 1 2 1 2 1 1 1 1 1 1
##  [4393] 1 2 1 2 1 2 1 1 2 2 1 1 1 2 2 1 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 2 1 2 2 2
##  [4429] 1 2 2 2 1 2 1 1 2 1 2 1 2 1 1 1 2 2 1 1 2 1 2 2 2 1 1 2 1 1 2 1 2 1 2 1
##  [4465] 2 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 2 2 1 1 1 1 1 1 2 2 1 2 1 2 1 2 1 2 1
##  [4501] 1 1 1 1 1 1 1 1 1 2 1 2 1 1 1 1 2 1 1 2 1 1 2 2 1 1 1 2 1 1 1 1 1 1 2 1
##  [4537] 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 2 1 1 1 1 1 1 1 2 1 1 1 1
##  [4573] 2 1 1 1 1 1 1 1 2 1 1 2 2 1 1 1 1 2 1 1 1 1 1 2 1 1 2 1 1 2 1 1 1 1 1 1
##  [4609] 1 1 2 1 1 1 1 1 1 1 1 2 1 1 1 2 1 2 1 1 2 1 1 1 1 1 1 2 1 1 2 2 2 2 2 1
##  [4645] 1 1 1 1 1 1 1 2 1 1 1 1 1 2 1 1 2 1 1 2 2 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1
##  [4681] 2 1 1 1 1 1 1 1 2 1 1 1 2 2 1 1 1 1 1 2 2 1 2 1 2 1 1 1 1 2 2 2 1 2 2 1
##  [4717] 1 1 2 2 1 1 2 1 2 1 2 1 1 2 1 1 2 2 1 1 1 1 1 2 1 1 1 1 2 1 1 1 1 2 1 1
##  [4753] 1 2 2 1 1 2 1 1 1 1 1 1 2 2 1 2 2 1 1 2 1 2 1 1 1 1 1 1 1 1 1 1 2 1 1 1
##  [4789] 1 1 1 2 1 1 1 1 1 1 2 2 2 1 1 1 2 1 1 1 1 1 1 1 1 2 1 2 1 1 1 2 1 1 2 1
##  [4825] 1 1 1 1 1 2 2 1 2 1 1 1 1 1 2 2 2 2 1 1 2 1 1 1 1 1 1 2 1 1 1 2 1 2 2 1
##  [4861] 2 2 2 2 2 2 2 1 2 2 2 1 2 1 1 2 2 1 1 1 2 1 1 1 2 2 1 1 1 2 2 1 1 2 1 1
##  [4897] 1 2 1 1 2 1 1 2 1 1 1 2 2 1 2 2 1 2 2 2 1 1 1 2 1 1 1 1 1 2 1 1 1 1 1 1
##  [4933] 2 1 1 1 2 1 1 1 2 2 2 1 1 1 1 1 2 1 1 2 1 1 1 1 2 1 1 1 2 2 1 1 2 2 1 1
##  [4969] 1 2 1 1 1 1 1 2 1 2 1 2 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 1 1 1 1 1 2
##  [5005] 1 1 2 2 1 1 2 1 1 2 1 1 2 2 1 1 2 1 2 1 1 1 2 2 1 2 1 1 1 2 2 2 1 1 1 1
##  [5041] 2 2 2 2 1 1 1 2 2 1 2 1 1 2 1 1 1 2 1 1 1 1 1 1 1 2 1 1 1 2 1 1 1 1 1 2
##  [5077] 1 1 1 1 1 2 1 1 2 1 1 1 2 1 2 1 1 1 2 1 2 2 1 2 2 2 1 1 1 2 2 1 1 1 2 1
##  [5113] 2 2 2 2 1 2 2 2 1 1 2 1 1 1 2 1 1 1 1 1 2 2 2 1 2 2 2 1 2 2 1 2 2 1 1 1
##  [5149] 1 1 1 1 2 2 1 2 2 1 1 1 1 1 2 1 2 1 1 1 1 1 2 1 2 1 1 1 1 2 1 1 2 1 1 1
##  [5185] 2 1 1 1 1 1 1 1 1 2 1 2 1 2 2 1 1 2 2 1 2 1 2 1 1 2 2 2 1 1 1 2 1 1 2 1
##  [5221] 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 1 2 1 1 2 1 2 1 1 2 2 1 1 1 1 1 1 1 2 1
##  [5257] 1 1 1 1 1 1 2 1 1 1 1 1 1 2 1 2 1 2 1 1 2 1 2 1 1 1 1 1 1 1 2 1 1 1 1 1
##  [5293] 1 1 1 2 2 1 1 2 1 2 1 1 1 2 2 1 2 1 2 1 2 1 1 1 1 1 1 1 1 1 1 2 1 2 2 2
##  [5329] 1 1 1 1 2 2 2 1 2 1 1 1 1 1 1 1 1 2 1 1 1 2 2 2 2 2 1 1 1 1 1 2 2 1 1 2
##  [5365] 1 1 1 1 2 2 1 1 1 1 2 1 1 1 1 1 1 2 2 1 2 1 1 1 1 2 1 2 1 2 2 1 2 1 2 2
##  [5401] 1 1 1 1 1 1 1 1 2 1 2 1 1 2 1 2 2 1 1 1 1 1 1 2 2 1 1 1 1 1 1 1 1 1 1 1
##  [5437] 1 1 2 2 1 1 1 1 1 1 1 1 2 2 1 1 1 1 1 1 2 2 1 1 1 2 2 1 1 1 2 1 1 1 1 2
##  [5473] 2 2 1 2 1 1 1 1 2 1 2 1 2 1 1 1 2 1 1 1 2 2 1 2 1 1 1 2 1 2 1 1 1 1 1 1
##  [5509] 1 2 2 1 1 2 1 2 2 1 1 1 1 1 2 1 2 1 1 1 1 1 1 2 1 1 1 2 1 1 1 1 1 1 1 1
##  [5545] 1 1 1 1 1 2 2 1 1 2 1 2 1 2 1 1 2 1 2 1 2 2 2 2 1 1 1 2 2 2 2 1 1 1 2 2
##  [5581] 2 1 1 2 2 1 2 1 1 1 2 2 1 2 2 2 2 1 1 1 1 1 1 2 2 1 1 1 1 2 1 1 1 1 2 1
##  [5617] 1 1 2 2 1 1 2 1 1 1 1 2 2 2 1 1 1 1 1 2 1 2 1 1 2 1 1 1 1 1 1 1 2 2 1 2
##  [5653] 2 1 1 2 2 1 2 2 1 1 1 1 1 2 2 2 1 1 1 2 1 2 1 1 2 1 1 1 2 1 1 1 1 2 1 1
##  [5689] 1 1 2 1 1 2 2 2 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1
##  [5725] 1 1 1 2 1 1 1 1 2 1 1 1 1 1 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2
##  [5761] 1 1 1 1 1 1 1 2 2 1 1 1 2 1 2 1 1 1 1 2 1 2 1 1 1 1 1 1 2 1 1 2 1 1 1 1
##  [5797] 2 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 2 1 2 1 2 2 2 1 1 2 1 1 1 2 2
##  [5833] 1 1 1 1 1 1 1 2 1 1 2 2 2 1 1 2 2 1 1 1 2 1 1 2 2 2 1 1 2 2 2 2 1 1 1 1
##  [5869] 2 2 1 2 1 1 2 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1 2 1 1 2 2 1 2 1 1 1 1 2 2 1
##  [5905] 2 2 1 1 1 2 2 1 1 1 1 1 1 1 2 2 1 2 2 1 1 2 1 2 2 2 2 1 1 1 1 1 1 2 1 1
##  [5941] 1 1 2 1 1 1 1 1 2 2 1 1 2 1 2 1 1 1 2 1 2 1 1 2 1 1 1 1 1 2 2 1 2 1 1 1
##  [5977] 2 2 1 2 1 1 2 1 2 1 1 1 1 2 2 2 1 1 2 1 1 1 1 1 2 1 1 2 2 1 2 2 1 1 1 1
##  [6013] 1 1 1 2 2 2 1 1 1 1 1 1 1 1 1 1 2 1 2 1 1 1 1 1 1 1 2 1 1 1 1 2 1 1 2 1
##  [6049] 2 2 2 2 2 2 2 1 1 1 2 1 1 1 2 1 2 1 2 1 1 2 2 2 2 2 2 1 1 2 2 1 2 2 1 2
##  [6085] 2 2 1 1 1 1 1 1 1 2 2 1 2 1 1 1 1 2 2 2 1 1 1 2 1 2 1 1 1 2 1 2 2 1 1 2
##  [6121] 2 1 1 1 2 2 2 1 1 2 1 2 1 2 2 1 2 2 1 1 1 2 1 1 1 1 1 2 2 2 2 1 1 1 1 2
##  [6157] 1 1 1 1 1 2 1 1 2 1 1 2 1 2 2 2 1 2 1 2 1 1 2 1 2 2 1 1 1 2 2 1 1 1 1 1
##  [6193] 2 1 1 1 2 2 1 2 2 1 1 1 2 1 1 1 2 1 2 2 2 1 1 1 1 2 2 2 1 1 1 2 2 1 1 1
##  [6229] 2 1 1 2 1 1 1 1 2 2 1 2 1 1 1 2 1 1 1 1 1 1 1 1 1 2 2 1 1 1 1 1 2 2 1 1
##  [6265] 1 2 1 1 1 1 1 1 1 1 2 1 1 1 2 1 1 2 1 2 1 1 1 2 1 1 2 1 1 1 1 1 2 2 1 1
##  [6301] 1 1 1 1 2 2 1 1 1 1 1 2 1 1 1 2 1 1 2 2 1 1 1 1 1 1 1 2 1 2 2 1 1 1 1 2
##  [6337] 1 2 2 2 2 2 2 1 2 1 1 2 2 2 1 2 1 2 2 1 1 1 2 1 1 2 2 2 2 1 1 2 2 2 1 2
##  [6373] 1 1 1 2 1 1 1 1 2 2 1 1 1 1 1 1 1 2 1 1 1 2 1 2 2 1 1 2 1 1 1 1 1 2 1 1
##  [6409] 1 1 2 1 1 1 2 1 1 2 1 1 2 2 1 1 2 1 1 1 1 2 2 1 2 1 1 1 2 1 1 1 1 1 1 1
##  [6445] 1 2 1 1 1 1 1 2 2 2 1 1 1 2 1 2 1 1 1 1 1 1 1 2 2 2 2 2 2 2 1 2 2 1 1 2
##  [6481] 2 2 1 2 1 1 1 1 1 1 1 1 1 1 2 2 2 1 2 2 2 2 2 1 2 1 1 1 2 1 1 1 1 2 2 1
##  [6517] 2 2 1 1 1 1 1 1 1 1 1 1 2 2 2 2 1 1 1 1 2 1 1 1 1 2 1 2 1 2 2 1 2 2 2 1
##  [6553] 1 1 2 1 1 1 1 1 2 1 1 1 2 2 1 2 1 2 2 1 1 1 2 1 1 1 2 1 2 1 1 1 1 2 2 2
##  [6589] 1 1 1 2 2 2 2 1 1 1 1 2 1 1 2 1 1 2 1 1 1 1 1 1 1 2 1 2 2 1 1 1 2 2 1 1
##  [6625] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 2 1 2 2 1 2 1 1 1 1 1 2
##  [6661] 1 1 1 1 1 1 1 1 1 2 2 1 1 1 2 1 2 1 2 2 1 2 1 1 2 1 2 1 1 2 2 2 1 1 1 1
##  [6697] 2 1 1 1 1 1 1 2 2 2 2 1 1 2 2 2 1 1 2 2 1 2 2 1 1 1 1 1 1 1 1 1 1 2 1 1
##  [6733] 2 1 1 1 1 1 2 1 2 1 2 1 1 2 1 1 1 2 2 1 1 1 2 2 2 2 2 1 1 1 1 2 1 1 2 1
##  [6769] 2 2 2 1 1 2 1 2 1 2 1 2 2 1 1 2 1 1 1 2 1 1 1 2 1 1 1 2 2 2 2 2 1 1 1 2
##  [6805] 1 1 2 1 1 1 1 1 1 1 2 2 1 1 2 2 1 1 1 2 2 2 1 1 1 1 2 2 1 2 1 1 1 1 1 2
##  [6841] 1 1 2 1 1 1 2 1 1 1 1 1 1 1 1 1 2 1 1 1 1 2 2 2 2 2 1 2 2 2 1 1 2 1 2 2
##  [6877] 1 1 1 2 1 2 2 1 1 1 1 1 2 1 2 1 1 1 1 1 1 2 2 1 2 2 1 1 1 1 1 1 2 2 1 1
##  [6913] 1 1 1 1 2 1 1 1 2 2 1 1 1 1 1 1 1 2 2 2 1 1 1 1 1 1 2 1 2 2 2 1 1 1 1 2
##  [6949] 2 1 2 1 2 2 1 2 1 1 2 2 1 1 1 1 2 1 1 2 1 1 2 1 1 1 2 2 2 2 1 1 1 1 1 1
##  [6985] 1 1 2 1 1 2 1 2 1 1 2 1 1 1 1 1 1 2 1 2 2 1 1 1 1 2 1 1 1 1 1 2 1 2 1 1
##  [7021] 1 2 1 1 2 1 1 2 2 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 1 2 2 1 2 2
##  [7057] 1 2 2 2 2 2 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 2 1 1 1 2 2 1 1 1 1
##  [7093] 1 1 1 1 1 1 2 1 1 1 1 2 1 1 1 1 2 1 1 1 1 2 1 1 1 1 1 2 2 2 1 2 1 1 1 1
##  [7129] 1 1 2 1 1 1 1 1 2 2 1 2 1 1 1 1 1 1 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2
##  [7165] 1 2 1 1 2 1 1 1 1 1 2 1 1 1 1 1 2 1 1 1 1 1 2 2 1 1 1 1 1 1 1 1 1 1 1 1
##  [7201] 1 1 2 1 1 2 1 1 1 2 1 1 1 1 1 1 1 1 2 1 1 1 1 2 1 1 1 1 1 1 2 1 2 1 1 1
##  [7237] 1 1 1 1 2 1 1 1 2 2 2 1 1 1 1 2 1 1 1 1 2 1 1 1 2 1 1 1 1 1 1 1 1 2 1 2
##  [7273] 2 2 2 1 1 1 1 1 2 1 2 2 1 1 2 1 2 2 2 2 1 1 1 1 2 1 1 2 1 1 1 1 1 1 1 2
##  [7309] 1 1 1 1 1 1 1 1 2 2 1 1 1 1 1 2 1 2 1 1 2 1 1 2 1 1 2 1 1 2 2 2 2 1 1 1
##  [7345] 1 1 2 1 1 1 1 2 2 1 2 1 1 1 1 2 2 1 1 1 2 1 2 1 1 1 1 1 1 2 1 1 1 1 2 1
##  [7381] 1 2 1 1 1 2 1 1 2 1 2 2 1 1 1 1 2 1 1 1 1 1 2 2 1 1 1 1 1 2 2 1 2 1 1 2
##  [7417] 1 1 2 1 2 1 2 2 2 1 1 1 1 1 1 1 1 2 1 1 2 1 2 1 2 1 1 1 1 1 1 1 1 1 1 1
##  [7453] 1 1 1 1 2 1 1 1 1 1 2 1 1 2 1 1 1 1 2 1 1 1 1 1 1 1 1 2 1 1 2 1 1 1 1 1
##  [7489] 1 2 1 1 1 1 1 1 1 1 1 1 1 2 2 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
##  [7525] 1 2 1 1 1 1 2 1 2 2 1 1 1 1 1 1 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1
##  [7561] 1 1 1 1 1 1 1 1 1 2 1 1 1 1 2 1 1 1 1 1 1 1 1 1 2 1 2 1 1 1 2 1 1 1 1 1
##  [7597] 1 1 1 1 1 1 1 1 2 2 2 1 2 1 1 1 1 1 1 2 2 1 2 1 1 1 2 1 1 1 1 1 1 1 1 1
##  [7633] 2 1 1 1 1 2 2 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 2 1
##  [7669] 1 1 1 1 1 1 2 1 1 2 2 1 1 1 1 1 1 1 2 2 1 1 1 1 2 1 2 1 2 1 2 2 1 1 1 2
##  [7705] 1 1 2 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 2 2 2 1 2 1 1 1 1 1
##  [7741] 1 2 1 1 2 1 1 1 1 1 1 2 2 1 2 2 1 1 1 1 1 1 2 2 1 1 1 2 1 1 1 1 1 1 1 1
##  [7777] 2 1 1 1 1 1 1 1 1 2 1 1 2 1 1 1 2 2 1 2 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 2
##  [7813] 2 1 1 1 2 2 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 1 1 1 2
##  [7849] 2 1 1 1 1 1 1 2 1 1 1 1 2 1 1 1 1 1 1 2 1 1 1 1 2 2 1 1 1 1 1 1 1 1 2 2
##  [7885] 2 1 2 1 1 1 1 1 2 1 1 2 1 1 1 2 2 1 1 1 1 2 1 1 1 1 1 1 1 1 2 1 1 1 1 1
##  [7921] 1 1 2 1 2 2 1 1 1 1 1 1 2 1 1 1 2 1 2 2 1 1 1 1 1 1 1 1 2 2 1 2 1 1 1 2
##  [7957] 1 2 1 1 1 1 2 2 1 1 1 1 1 2 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1
##  [7993] 1 1 1 1 1 2 2 1 2 1 1 1 2 2 1 2 1 1 1 2 1 2 2 2 1 2 2 1 1 2 1 1 1 1 1 1
##  [8029] 1 1 1 1 1 1 1 1 2 1 1 2 1 1 1 1 1 1 1 2 1 1 1 1 1 1 2 1 1 2 1 1 1 1 1 2
##  [8065] 2 1 1 2 2 2 2 1 2 2 2 2 1 1 1 2 1 1 1 1 2 1 1 2 1 1 1 1 1 1 1 1 1 1 2 1
##  [8101] 2 2 2 2 1 1 2 2 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 2 1 1 2 1 2 1 1 1 1 1 2 1
##  [8137] 2 1 2 1 1 1 1 1 2 1 2 2 1 1 1 1 2 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 2 1 1 2
##  [8173] 1 1 1 1 1 1 1 2 2 2 1 2 2 1 1 2 1 2 1 2 2 1 1 1 2 1 1 2 1 2 1 1 1 2 1 1
##  [8209] 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 2 2 1 1 2 1 1 2 1 1 1 1 1 1 1 1 1 2 1
##  [8245] 1 2 1 1 2 1 2 1 1 1 1 1 2 1 1 1 1 1 1 1 2 1 1 2 1 1 2 1 2 1 1 1 2 1 1 2
##  [8281] 2 2 1 2 1 1 1 2 1 1 1 2 1 2 2 2 2 1 1 1 2 2 2 1 1 1 2 1 2 1 1 1 2 2 1 2
##  [8317] 2 2 1 2 2 1 1 1 2 1 1 1 2 2 2 2 2 1 2 1 2 2 1 1 1 1 2 1 1 1 1 2 1 1 2 1
##  [8353] 1 1 1 1 1 2 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 1 1 1 1
##  [8389] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 2 2 1 2 1 2 1 1 2 1 1 1 2 2
##  [8425] 2 2 2 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1
##  [8461] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
##  [8497] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 2 1
##  [8533] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
##  [8569] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
##  [8605] 1 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
##  [8641] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
##  [8677] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 2 1 2 1 2 1 1 1 1
##  [8713] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
##  [8749] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
##  [8785] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
##  [8821] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
##  [8857] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
##  [8893] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 2 1 1 1 1 1 1 2
##  [8929] 1 1 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 1 1 1 1 1 2 1 1 1 2 2 2 1
##  [8965] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 2 2 1 1 2
##  [9001] 2 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 1 1 2 1 1 1 1 1
##  [9037] 1 1 1 2 1 1 2 1 2 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
##  [9073] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
##  [9109] 1 1 2 2 2 1 1 2 1 1 2 1 1 2 1 1 1 1 1 2 1 1 2 1 1 1 1 1 2 1 1 1 1 1 1 1
##  [9145] 1 1 1 1 1 1 1 1 1 1 2 1 1 1 2 2 2 1 1 1 1 1 1 1 1 2 1 1 1 2 1 2 1 1 2 1
##  [9181] 1 2 2 2 1 2 2 2 1 1 2 1 1 1 2 1 1 2 1 1 2 2 1 1 1 1 1 1 1 2 1 1 2 2 2 2
##  [9217] 2 1 2 2 1 1 1 1 1 2 1 2 2 2 2 1 1 2 1 1 1 1 2 2 1 1 1 1 1 1 1 2 1 1 1 1
##  [9253] 1 1 1 1 1 1 2 1 1 1 1 2 2 1 1 1 2 1 2 1 1 1 1 1 1 1 2 1 2 1 1 2 1 1 1 2
##  [9289] 2 1 2 2 1 1 1 2 1 2 2 2 2 2 1 1 1 2 2 1 2 2 1 1 1 1 2 2 1 2 1 1 1 1 1 2
##  [9325] 1 1 2 2 1 1 1 2 1 1 2 1 2 1 2 1 2 1 1 1 2 2 1 2 2 1 1 1 2 1 1 2 2 2 1 1
##  [9361] 1 1 1 2 1 1 2 1 1 1 1 1 1 1 2 2 1 1 1 1 1 2 1 2 1 1 2 1 1 1 1 2 2 1 2 1
##  [9397] 1 1 1 1 1 1 2 1 2 1 1 2 2 2 2 2 1 1 1 1 1 2 1 2 2 1 1 1 1 1 1 1 1 2 1 2
##  [9433] 2 2 1 1 2 1 1 1 1 1 2 1 1 2 1 1 1 2 1 1 2 1 1 1 2 1 2 2 2 2 2 1 2 2 1 1
##  [9469] 1 2 1 1 1 1 2 1 1 2 1 1 1 1 2 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1
##  [9505] 1 2 1 1 1 2 1 1 1 1 1 2 2 1 1 1 1 1 1 1 1 1 1 2 1 1 2 1 1 1 1 2 1 2 1 1
##  [9541] 1 1 1 2 2 1 2 2 1 1 2 1 1 2 2 2 2 2 1 1 2 1 2 1 1 1 1 2 1 2 2 2 2 1 1 1
##  [9577] 2 1 1 1 1 1 1 2 1 1 2 1 1 1 2 2 1 1 2 2 1 2 2 2 1 1 2 2 1 1 1 1 2 2 2 1
##  [9613] 1 1 2 1 1 1 2 1 1 1 1 1 1 1 1 2 2 2 2 2 1 2 2 2 1 2 1 1 1 2 1 1 2 1 1 1
##  [9649] 1 1 1 1 1 1 2 1 1 2 1 1 2 2 1 1 1 1 1 1 2 2 2 1 2 1 1 1 2 2 2 1 2 1 2 2
##  [9685] 1 1 1 2 1 1 1 1 1 2 2 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 1 1 2 1 2 1 1 1 2 2
##  [9721] 2 2 1 1 1 2 2 2 2 1 1 2 2 2 1 1 2 1 1 1 2 1 2 1 2 1 1 2 1 2 1 2 1 1 1 2
##  [9757] 2 1 2 1 1 1 1 1 1 1 1 2 1 2 1 1 1 1 1 2 1 1 1 2 1 2 2 1 1 1 1 1 1 1 1 1
##  [9793] 1 2 2 2 1 2 2 2 1 1 1 1 1 1 1 1 1 1 2 1 2 1 2 2 1 2 1 2 1 1 1 2 2 1 1 2
##  [9829] 1 1 1 2 1 1 1 2 1 1 1 2 2 1 1 1 1 2 2 2 2 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1
##  [9865] 1 1 1 1 1 1 1 1 1 2 2 1 1 1 2 1 1 1 1 1 1 1 1 1 2 1 2 1 1 1 2 2 1 1 1 2
##  [9901] 1 2 2 1 2 1 1 2 1 1 2 1 1 1 1 2 1 1 1 1 2 1 1 1 1 1 2 1 1 2 1 1 1 1 1 2
##  [9937] 1 2 1 1 1 1 2 1 2 1 1 1 1 1 1 2 1 2 1 1 1 1 2 2 2 1 1 1 2 2 2 1 2 2 2 2
##  [9973] 2 1 1 1 1 2 1 1 1 1 1 1 1 2 2 1 2 1 2 1 1 1 1 1 1 1 1 2 1 2 1 2 1 1 1 2
## [10009] 2 1 1 1 1 1 2 2 2 1 2 1 1 1 1 1 1 2 1 2 1 2 1 2 2 1 1 1 2 1 1 1 1 1 1 1
## [10045] 2 2 1 1 1 1 1 1 1 2 1 2 2 2 1 1 2 1 2 1 1 2 1 1 1 2 1 1 1 1 2 1 1 2 2 2
## [10081] 1 1 1 1 1 1 1 2 2 1 1 1 1 1 1 2 1 1 1 1 2 1 2 1 2 2 1 1 1 1 1 1 1 2 1 2
## [10117] 2 1 1 2 1 1 1 1 1 1 1 2 2 1 1 1 1 1 2 2 1 1 2 1 1 1 1 1 1 1 1 1 1 1 2 2
## [10153] 2 1 1 1 2 1 2 1 1 1 1 2 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 1 2 1 1 1 2 1 1 1
## [10189] 2 2 1 1 1 1 1 2 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 2 1 2 1 2 1 1 1 1 1 1 1
## [10225] 2 2 1 1 2 1 1 1 2 1 2 1 2 2 2 1 1 2 1 1 1 2 1 2 2 2 2 1 1 1 1 2 1 1 1 2
## [10261] 2 2 2 2 1 1 2 2 1 2 1 2 2 2 1 1 1 1 2 1 1 2 2 2 2 2 1 2 1 1 1 1 1 1 1 2
## [10297] 2 1 2 2 2 2 1 1 2 1 1 1 2 1 1 1 2 1 2 2 1 1 2 1 1 2 1 1 2 1 2 1 1 2 2 1
## [10333] 1 1 1 1 1 1 1 1 2 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
## [10369] 1 1 1 1 1 1 1 1 1 2 1 1 1 1 2 1 1 1 1 2 1 1 2 2 1 1 2 1 1 2 1 1 1 1 1 2
## [10405] 1 1 1 2 2 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 1
## [10441] 2 1 1 1 1 1 1 2 2 1 1 1 1 1 1 1 1 1 2 1 1 2 1 1 1 1 2 1 1 1 1 1 1 1 1 1
## [10477] 1 1 1 1 1 1 1 2 1 1 1 1 1 1 2 1 1 1 1 2 1 2 1 1 2 1 1 1 1 1 1 1 2 1 1 1
## [10513] 1 1 1 1 1 1 2 1 2 2 1 1 2 1 1 1 1 1 1 2 1 1 2 2 1 1 1 1 1 2 1 1 2 2 2 1
## [10549] 1 2 1 1 1 1 1 2 2 1 1 1 2 1 1 1 1 1 1 1 1 1 2 1 2 1 2 1 1 2 1 1 1 1 1 1
## [10585] 2 1 1 1 1 2 1 1 1 1 1 1 1 2 1 1 1 1 2 1 2 1 1 2 1 2 2 2 2 2 2 1 1 1 1 2
## [10621] 1 1 1 1 1 2 1 1 1 1 1 1 1 2 1 2 1 1 2 1 2 1 1 1 2 1 1 1 1 1 1 1 1 1 1 2
## [10657] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 2 1 2 1 1 2 1 1 1 1 2 2 1 2 2
## [10693] 1 1 2 1 1 1 1 1 1 1 1 2 2 1 1 1 1 2 1 2 1 1 1 2 2 1 1 1 2 1 1 1 1 1 1 1
## [10729] 1 1 1 2 1 1 1 2 1 1 1 1 2 1 1 1 2 2 1 1 1 1 2 2 2 1 1 2 1 2 2 2 2 1 2 2
## [10765] 2 1 1 2 1 1 1 1 1 2 1 1 1 2 1 1 1 1 2 2 1 2 1 1 1 1 2 1 1 2 2 2 1 1 1 1
## [10801] 2 1 1 2 1 2 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2
## [10837] 1 1 1 2 1 1 1 1 1 1 1 1 1 2 2 1 2 1 1 1 1 1 1 1 1 1 2 2 1 1 1 1 1 1 2 1
## [10873] 2 1 2 1 1 1 2 1 1 1 1 1 1 2 1 1 2 1 1 2 1 1 1 1 2 1 1 1 1 1 1 1 1 1 2 2
## [10909] 1 1 2 1 1 1 1 2 1 1 1 2 2 1 1 1 1 2 1 2 1 1 1 1 1 2 2 1 1 1 1 1 1 1 2 1
## [10945] 1 1 1 1 1 2 1 1 1 2 1 2 1 1 1 2 1 1 2 1 1 2 1 1 2 1 1 1 2 2 1 1 2 1 1 1
## [10981] 2 1 1 2 1 1 1 1 1 1 1 1 1 2 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1
## [11017] 1 1 2 2 2 2 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 1 1 1 1 1 1 1 2 2 1 1 1 1 1 1
## [11053] 1 2 1 1 1 1 1 1 1 1 2 2 1 2 2 2 1 1 2 1 1 1 1 2 1 1 2 2 1 2 2 2 2 2 1 1
## [11089] 1 1 2 1 1 1 2 2 1 1 1 1 1 2 1 1 2 1 1 2 2 2 2 1 1 1 1 1 2 1 1 1 1 1 1 1
## [11125] 1 1 2 2 2 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 2 1 1 2 1 1 1 1 1 2 1 1 1 1 1
## [11161] 1 1 1 1 1 1 2 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1 2 2 2 1 2 1 2 1 1 1 1 2 1 1
## [11197] 1 1 2 1 1 1 1 1 2 2 1 1 1 1 2 2 1 2 1 1 1 1 1 2 1 1 2 2 1 2 2 1 1 1 1 2
## [11233] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 2 1 1 2 1 1 1 1 2 1
## [11269] 1 2 2 1 2 1 1 2 2 1 2 1 1 1 1 1 2 2 1 2 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1
## [11305] 1 1 1 1 1 1 1 2 1 1 2 1 1 1 2 1 2 1 2 1 1 1 1 2 1 1 1 1 1 1 1 2 1 1 1 1
## [11341] 1 1 1 1 1 1 1 1 2 1 2 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 2 1 2 1 1 1 2 1 2 1
## [11377] 1 2 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 2 1 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2
## [11413] 2 2 2 1 1 1 1 1 2 1 2 1 1 1 2 2 1 2 2 1 1 2 1 1 2 1 2 1 2 1 2 1 2 2 1 1
## [11449] 1 2 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 1 1 1 1
## [11485] 2 1 2 1 1 1 1 1 1 1 2 2 2 1 1 1 2 2 2 1 1 1 1 1 1 2 2 1 1 1 1 1 1 1 1 1
## [11521] 1 2 2 2 1 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 1 1 2 1
## [11557] 2 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 2 1 2 1 1 2 1
## [11593] 2 1 1 1 2 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
## [11629] 2 2 1 1 2 2 1 1 1 2 1 1 1 2 1 1 2 1 1 1 1 2 2 1 1 2 2 1 2 1 1 1 1 1 1 1
## [11665] 1 1 1 1 1 1 1 1 1 1 2 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 2 1 1
## [11701] 2 1 1 1 1 2 1 1 1 1 1 1 1 1 1 2 1 1 1 2 1 1 1 1 1 1 2 1 1 1 2 2 2 2 1 1
## [11737] 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1
## [11773] 1 1 2 1 2 1 1 1 1 1 2 1 1 2 1 1 2 2 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1
## [11809] 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 1 1 1 1 1 1
## [11845] 1 1 1 1 1 1 1 1 1 1 2 2 2 1 2 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 2 1 1 1 1 1
## [11881] 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 2 1 1 1 1 1 1 1 1 2
## [11917] 1 1 1 2 1 1 1 1 1 1 1 2 1 2 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 1 2 1 1 2 1 1
## [11953] 1 1 2 1 1 1 2 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 2 1 1 1 1 1 2 1
## [11989] 2 2 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 2 1 1 1
## [12025] 1 1 1 1 1 1 1 2 2 1 2 1 1 1 2 2 1 1 1 1 1 1 1 1 1 2 1 1 1 2 1 1 2 1 1 1
## [12061] 1 1 1 1 1 2 1 2 2 2 1 1 1 1 1 1 1 2 2 1 2 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1
## [12097] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
## [12133] 2 1 2 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1
## [12169] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 2 1 2 2 1 1 1 2 1 1 2 1 1 1
## [12205] 1 1 2 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 1 2 2 1 1 2 1 1 2 1 1 1 2 1 2
## [12241] 1 1 1 2 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 2 1 1 1 1 1 1 1 1 1 1 2
## [12277] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
## [12313] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 2 1 1 1 2 1 1 1 1 1 2 1 1 1
## [12349] 1 1 1 1 2 1 1 2 1 1 2 1 1 2 1 1 1 1 2 2 1 2 2 1 2 1 1 1 1 1 2 2 2 1 2 1
## [12385] 1 1 1 1 1 1 1 1 1 1 1 1 2 1 2 1 1 2 2 1 1 1 1 1 1 1 1 2 1 1 1 2 1 1 1 2
## [12421] 1 1 1 1 1 1 1 1 1 2 1 2 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 2 1 1 1 1 1 1
## [12457] 1 2 2 2 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
## [12493] 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 1 2 1 2 1 2 1
## [12529] 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
## [12565] 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 2 1 1 1 2
## [12601] 1 1 1 2 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 2 1 1 1 1 1 2 1 1 1 1 1 1 1 2 1 1
## [12637] 1 1 1 1 1 1 1 2 1 2 2 2 1 1 2 1 1 2 1 1 1 1 2 1 1 2 1 1 1 1 1 1 1 1 2 1
## [12673] 1 1 1 1 1 2 2 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1 2 1 1 2 1 1 1 1 1 1 1 1
## [12709] 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1
## [12745] 1 1 2 1 1 1 1 1 2 1 1 1 1 1 2 1 1 1 1 1 1 1 1 2 1 1 1 1 2 1 1 2 2 1 2 1
## [12781] 1 1 2 1 1 1 1 1 2 1 1 1 1 1 2 2 2 1 2 1 1 2 2 1 1 1 1 2 2 1 1 1 2 1 1 1
## [12817] 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
## [12853] 2 2 1 1 1 2 1 1 2 2 1 1 2 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 2 1 2 1 1 1 1 1
## [12889] 1 1 1 1 2 1 1 1 1 2 1 1 1 1 1 2 2 1 2 2 1 1 1 2 1 1 2 2 1 2 1 1 2 2 2 1
## [12925] 1 2 2 1 1 2 2 2 1 1 1 1 2 1 1 2 1 2 2 1 2 1 2 2 1 1 2 1 1 1 1 2 1 1 1 1
## [12961] 2 1 1 2 2 1 1 2 2 2 1 1 1 1 1 1 2 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
## [12997] 1 1 1 1 1 1 1 2 1 1 1 1 2 1 1 1 1 1 1 1 1 1 2 2 1 2 2 2 1 1 1 2 2 1 2 2
## [13033] 1 2 2 1 1 1 1 1 2 1 1 1 1 1 2 1 1 1 1 1 1 1 2 2 1 1 1 1 1 1 1 1 1 1 1 1
## [13069] 1 1 2 1 1 1 1 1 2 1 1 1 1 2 2 2 2 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1 2 2 2 1
## [13105] 1 2 2 1 2 1 1 1 2 1 1 1 2 1 1 1 2 2 2 2 2 2 2 1 1 2 1 2 1 1 2 1 2 1 2 2
## [13141] 1 1 2 2 2 2 1 2 1 1 1 1 2 1 1 1 1 1 2 1 2 1 2 1 1 1 2 1 1 2 2 1 1 1 2 1
## [13177] 1 1 2 1 1 1 2 2 1 2 1 1 1 1 1 1 1 1 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
## [13213] 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1 2
## [13249] 1 2 2 1 2 1 1 1 1 1 1 1 2 1 1 1 2 1 1 1 2 2 2 1 1 2 1 1 2 1 2 2 1 1 2 1
## [13285] 2 1 1 1 1 1 1 1 1 2 1 1 1 2 1 1 2 2 1 2 1 1 1 2 1 1 2 1 1 2 1 2 1 1 2 1
## [13321] 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
## [13357] 1 1 2 1 1 1 2 1 2 1 1 2 1 2 1 2 1 2 1 2 2 1 1 1 2 1 1 1 2 1 1 1 1 1 1 1
## [13393] 1 1 1 2 1 1 2 2 1 1 1 1 1 1 1 1 1 1 2 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1
## [13429] 1 1 1 1 1 2 1 2 1 2 1 2 1 1 1 1 2 1 2 1 1 2 1 1 1 2 2 2 1 2 2 2 1 1 1 1
## [13465] 1 1 1 1 1 1 2 1 2 1 2 2 1 1 1 1 1 1 1 1 1 1 1 2 1 2 1 2 1 1 1 2 2 1 1 1
## [13501] 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 1 2 2 1 2 1 1 1 1 2 1 1 1
## [13537] 2 1 1 1 1 2 1 1 1 2 1 1 1 1 2 2 2 2 1 2 2 1 1 2 1 1 1 1 2 1 1 2 1 1 2 1
## [13573] 2 1 1 1 1 2 1 1 1 1 1 1 2 2 2 1 1 1 1 2 1 2 1 1 1 2 2 1 1 1 1 1 1 1 1 1
## [13609] 2 1 1 1 1 1 1 2 2 2 1 2 1 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1
## [13645] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 2 1 2 1 1 1 1 1 1
## [13681] 2 2 2 1 1 1 2 1 2 2 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 2 2 1 2 1 2 1 2 1 1
## [13717] 1 1 1 1 2 2 1 1 1 1 1 1 2 1 1 1 1 1 1 2 1 2 1 1 1 2 1 1 2 1 1 1 1 1 1 1
## [13753] 1 1 2 1 1 1 1 1 2 2 1 1 1 1 1 1 2 2 1 1 2 2 1 1 2 2 1 2 2 2 2 1 2 1 2 1
## [13789] 2 2 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 2 1 2 1 2 1 2 2 2 1 2 2 1 1 2 1 1 1
## [13825] 2 1 2 1 1 1 1 1 2 1 1 2 2 1 2 1 1 2 1 1 2 2 1 1 1 1 1 2 1 1 2 1 1 1 1 2
## [13861] 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 2 1 1 2 2 2 2 1 1 1 1 1 1 1 2 1 2
## [13897] 1 2 1 1 2 1 1 1 1 2 2 1 2 2 2 1 2 1 2 1 1 1 2 1 2 2 2 2 1 1 2 1 2 1 2 1
## [13933] 2 1 1 2 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 2 1 2 2 2 1 1 1 1 2 1 2 1 1 2 1 2
## [13969] 1 2 1 2 1 1 1 2 2 1 1 1 1 1 2 2 2 1 1 2 1 1 2 1 1 2 2 2 2 1 1 2 1 1 1 2
## [14005] 1 1 2 1 1 1 1 2 1 2 2 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 2 2 1 1
## [14041] 1 1 1 2 1 1 1 1 1 1 1 1 2 2 1 2 1 1 1 1 1 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1
## [14077] 1 1 1 1 1 2 1 1 1 1 1 2 1 2 1 1 1 2 1 1 1 1 1 1 1 2 1 2 1 2 2 2 2 1 1 2
## [14113] 1 1 1 1 1 1 2 2 2 1 2 1 1 1 2 1 2 2 1 1 1 2 1 1 1 2 1 1 1 1 1 1 1 1 2 2
## [14149] 2 2 1 1 1 2 2 1 2 1 2 2 2 1 1 1 1 2 2 2 1 1 2 2 1 2 2 1 1 1 2 1 2 2 1 1
## [14185] 2 1 1 2 2 2 1 2 1 1 2 1 2 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 2 2 2 1 1 1 1
## [14221] 1 1 1 2 1 1 2 1 2 1 1 1 1 1 2 1 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 2
## [14257] 1 1 2 1 2 1 1 1 1 1 1 1 1 1 1 1 2 2 1 2 1 2 2 2 2 2 1 1 2 2 2 2 1 1 2 1
## [14293] 1 1 1 2 2 1 1 2 2 1 2 1 1 2 1 1 2 2 1 1 2 2 1 1 2 2 1 2 2 2 2 1 2 2 1 1
## [14329] 2 2 1 2 1 1 1 1 1 2 2 2 1 1 2 1 1 1 1 2 1 1 1 1 1 2 2 1 1 2 1 2 2 2 1 2
## [14365] 1 1 2 1 2 2 1 2 2 1 1 2 1 1 2 1 1 2 1 1 1 2 2 2 1 2 1 1 1 2 1 1 2 2 1 1
## [14401] 1 1 1 1 1 2 1 1 2 2 2 1 2 2 2 2 2 2 2 1 2 1 1 2 1 2 1 1 2 1 2 1 1 1 1 1
## [14437] 1 2 2 1 1 1 1 2 2 1 1 2 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 1
## [14473] 2 2 2 2 2 2 1 1 1 2 1 2 2 2 2 2 1 1 1 2 2 1 2 1 2 2 1 2 1 1 1 1 1 1 2 2
## [14509] 1 2 2 2 2 1 2 2 1 2 1 1 1 2 2 2 2 1 1 1 1 1 1 1 2 2 2 1 2 1 2 1 2 1 2 2
## [14545] 2 1 1 1 2 1 1 2 1 1 2 1 1 2 2 1 1 1 2 2 1 2 1 2 1 1 2 2 1 1 1 1 1 1 1 1
## [14581] 2 1 1 1 1 1 2 2 1 1 2 2 1 2 1 2 1 1 1 1 1 2 2 2 2 1 1 1 2 1 1 2 1 1 1 1
## [14617] 2 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 1 2 1 1 1 1 1 1
## [14653] 1 2 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 2 1 1 2 1 1 1 2 2 2 1 2 1 1 1
## [14689] 1 1 1 1 1 1 2 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 2 1 1 2 1 1 1 2 1
## [14725] 1 1 1 2 1 2 1 1 2 2 1 1 1 1 1 2 1 2 1 2 1 1 1 2 1 2 1 1 1 1 1 1 1 1 1 1
## [14761] 2 2 2 1 1 1 1 2 2 1 2 1 1 1 2 2 2 2 1 1 1 2 2 2 2 1 1 1 1 1 1 1 1 1 1 2
## [14797] 2 2 1 1 1 1 1 2 1 1 2 1 2 1 1 1 2 1 1 2 1 2 1 1 2 2 1 1 1 1 1 1 1 1 1 1
## [14833] 1 2 1 1 2 1 1 2 2 1 1 2 1 2 2 1 1 1 1 2 1 1 2 1 1 1 2 2 2 1 2 1 1 1 1 1
## [14869] 1 1 1 2 1 1 1 1 1 1 2 1 1 1 2 1 1 2 1 2 1 2 2 2 1 1 2 1 1 2 2 1 2 2 2 2
## [14905] 1 1 1 1 2 1 2 1 1 1 1 1 1 1 1 1 1 1 2 1 1 2 1 1 1 1 1 1 1 1 1 2 1 1 1 1
## [14941] 2 2 2 2 2 1 2 1 1 1 2 1 2 2 1 2 2 1 1 1 2 1 1 1 2 1 2 1 1 1 1 1 2 2 2 1
## [14977] 2 2 1 2 2 2 2 1 1 1 1 2 2 1 1 2 1 2 1 2 1 2 1 2 2 2 2 1 2 2 1 1 1 1 1 1
## [15013] 1 1 1 1 1 1 1 1 1 1 1 2 1 1 2 2 1 1 2 1 1 2 2 2 1 1 1 1 1 2 1 1 1 2 1 1
## [15049] 1 1 1 2 1 2 1 1 1 2 1 2 1 1 2 2 1 2 2 1 2 1 2 2 1 2 2 1 1 2 1 2 2 1 2 1
## [15085] 1 2 1 1 2 1 2 1 2 1 1 1 1 1 2 1 1 2 2 1 1 2 1 2 2 2 2 1 2 1 2 2 1 2 2 1
## [15121] 2 1 2 2 2 1 2 2 2 1 1 2 2 2 2 2 1 1 1 1 2 2 1 2 2 1 2 2 1 1 1 2 1 2 1 1
## [15157] 1 2 1 2 1 2 2 1 2 1 1 2 1 2 2 1 1 1 2 1 1 1 1 1 1 1 1 1 1 2 1 1 2 1 1 1
## [15193] 2 2 1 2 2 2 1 1 1 1 2 1 1 1 1 1 2 1 1 1 2 1 2 2 2 1 1 1 1 2 1 2 1 1 1 2
## [15229] 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 2 1 2 1 1 2 2 1 2 1 1 1 1 2 2
## [15265] 1 1 2 1 1 1 1 1 1 1 2 1 1 2 1 1 1 1 1 1 1 2 1 2 1 1 1 1 1 1 2 2 1 1 1 1
## [15301] 2 1 1 1 2 1 2 2 1 1 2 1 1 2 2 1 1 2 1 1 2 1 1 1 1 2 1 1 1 1 1 1 1 2 1 2
## [15337] 2 1 2 1 2 1 2 1 2 2 1 2 2 2 1 2 1 1 2 1 1 2 1 2 2 1 2 1 1 2 1 1 1 1 1 1
## [15373] 2 1 1 1 1 2 2 1 1 1 2 1 1 1 1 2 1 2 2 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1
## [15409] 1 1 2 1 2 1 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 2 1 1 1 1 2 2 1 2 1
## [15445] 2 2 1 2 2 2 2 2 1 1 1 1 1 2 1 2 2 2 2 1 1 1 1 1 2 2 2 1 1 1 2 1 1 2 2 2
## [15481] 1 2 1 1 1 1 2 1 1 1 2 1 2 1 2 1 2 2 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1
## [15517] 1 1 1 1 1 1 1 1 2 1 1 1 1 2 1 1 1 1 1 1 1 1 2 1 1 2 2 1 1 1 1 1 1 1 1 1
## [15553] 1 1 1 1 1 2 2 2 1 1 1 1 1 1 1 1 1 2 1 2 1 1 2 2 1 2 1 1 1 2 1 2 1 2 2 1
## [15589] 1 1 1 2 1 2 2 2 2 1 2 2 1 1 2 1 1 1 2 1 2 1 2 2 2 1 1 1 2 1 1 1 1 1 2 2
## [15625] 1 1 1 1 1 1 1 2 1 1 1 1 2 1 2 2 1 2 2 2 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 2
## [15661] 1 2 1 1 2 2 2 1 1 2 2 1 1 1 2 1 1 2 1 1 2 2 1 1 2 2 1 1 1 2 2 2 1 1 1 1
## [15697] 1 2 2 2 1 1 2 1 1 1 2 1 1 1 1 1 1 1 1 2 1 1 2 2 1 2 1 1 2 1 1 1 1 2 2 1
## [15733] 2 1 1 2 1 2 1 1 1 1 2 1 1 1 2 2 1 1 2 1 1 1 1 2 2 2 1 1 1 2 2 1 1 1 1 1
## [15769] 2 1 1 2 2 1 2 1 1 1 1 1 1 1 1 2 1 2 1 1 2 1 2 1 2 1 1 1 1 1 1 1 1 1 1 1
## [15805] 1 1 1 2 1 2 1 1 1 1 1 1 2 2 1 2 2 1 1 2 1 2 1 1 1 1 2 1 1 2 1 1 2 2 1 2
## [15841] 2 1 2 1 2 2 1 2 1 1 2 1 1 1 2 1 1 2 2 1 1 1 2 2 1 1 1 2 1 2 1 1 1 2 1 1
## [15877] 1 1 1 2 1 1 2 1 1 1 1 1 2 1 2 1 1 1 1 1 1 1 2 2 1 1 2 1 2 1 2 1 2 1 2 2
## [15913] 2 1 2 2 2 2 1 1 1 1 1 1 1 2 1 2 2 1 2 1 1 2 2 1 2 1 2 2 1 1 1 1 2 1 2 2
## [15949] 1 1 2 1 1 2 2 1 1 1 1 1 2 1 1 1 1 1 2 2 2 1 2 2 2 1 1 2 1 1 1 1 1 2 1 1
## [15985] 1 1 1 2 1 2 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 1 1 1 1 1 1 2 2 1
## [16021] 1 1 1 2 1 1 1 1 2 2 1 2 2 2 1 2 1 1 1 2 1 2 2 1 1 1 2 1 2 2 2 1 1 1 2 1
## [16057] 1 1 1 1 1 1 1 1 1 1 2 1 1 1 2 1 2 1 2 1 1 1 2 1 1 2 2 1 1 1 1 1 2 1 1 1
## [16093] 1 1 2 1 1 2 1 1 1 1 2 1 2 1 1 1 1 1 1 1 1 2 1 1 2 1 1 1 1 1 1 1 1 1 1 2
## [16129] 1 2 1 1 1 2 2 1 2 1 2 1 2 1 2 1 1 2 1 1 2 1 1 2 1 2 1 1 1 2 1 1 1 1 2 1
## [16165] 1 1 1 1 1 1 1 1 1 1 1 2 1 2 2 1 2 1 2 2 2 1 2 1 1 2 1 2 1 1 1 1 1 2 2 2
## [16201] 1 1 1 2 1 1 2 2 1 1 1 1 2 2 1 1 1 2 2 2 1 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1
## [16237] 2 2 1 1 2 1 2 2 2 1 2 2 1 2 1 1 2 2 2 1 2 1 2 2 2 1 1 2 2 1 1 2 1 2 2 1
## [16273] 1 1 1 1 1 1 2 2 1 1 2 2 1 1 2 1 1 1 1 2 1 1 1 2 2 2 1 1 2 1 2 2 1 1 1 1
## [16309] 2 2 1 1 1 1 1 1 1 1 2 2 1 1 1 1 1 1 2 1 1 1 1 1 1 2 1 1 2 2 1 2 2 1 1 2
## [16345] 1 2 1 1 1 2 1 1 2 1 2 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1 2 1 1 1 2 1
## [16381] 1 1 2 1 1 2 2 2 1 2 1 1 2 1 1 1 2 1 1 1 2 1 1 1 2 2 2 2 1 1 2 1 1 1 2 1
## [16417] 1 2 1 2 2 1 1 1 1 1 1 1 2 1 1 1 1 1 1 2 1 1 1 1 1 1 2 1 1 2 1 2 1 1 1 1
## [16453] 2 1 1 2 2 1 2 2 1 1 1 2 1 2 1 1 1 1 2 2 1 1 1 1 1 2 1 1 1 1 1 1 2 2 1 2
## [16489] 2 1 1 1 1 1 2 2 2 1 2 1 2 1 1 1 2 1 1 2 1 2 1 1 2 1 1 1 1 1 1 1 1 2 2 1
## [16525] 1 2 1 1 2 1 1 1 2 2 1 2 1 2 1 1 2 1 2 1 2 1 2 1 1 1 1 1 1 2 2 1 1 2 1 1
## [16561] 1 1 2 2 2 1 1 1 2 1 1 2 1 2 1 2 1 2 1 1 1 1 2 2 2 1 1 2 1 1 2 1 1 2 1 1
## [16597] 1 1 1 2 1 1 2 1 2 2 1 1 1 1 1 2 1 2 1 2 2 1 2 1 1 1 2 1 1 1 2 2 1 2 1 1
## [16633] 2 1 1 1 1 2 2 1 2 1 1 1 1 1 2 1 2 1 1 1 1 2 2 1 1 2 2 1 1 1 1 2 2 1 2 2
## [16669] 1 1 1 2 1 1 1 1 1 1 1 1 2 2 1 1 1 1 1 1 1 2 2 1 1 1 2 2 2 1 1 1 1 1 2 2
## [16705] 1 1 1 2 2 1 1 1 1 1 1 2 2 2 2 1 1 2 2 2 1 1 2 1 1 2 1 1 1 1 2 2 2 1 2 2
## [16741] 2 2 2 1 1 1 2 1 1 1 1 1 1 1 1 1 1 2 1 2 2 1 1 1 1 2 1 1 1 2 1 1 1 2 1 1
## [16777] 1 2 1 2 1 1 2 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 1 2 2 1 1 1 1 2 1 1 1
## [16813] 1 2 1 2 1 1 1 1 1 1 2 1 2 1 1 2 1 2 1 1 2 1 1 1 1 2 2 2 1 2 1 2 2 1 1 1
## [16849] 1 2 2 1 1 1 2 2 1 1 2 1 2 1 2 2 1 1 1 2 1 2 2 2 1 1 2 2 2 1 1 2 2 2 1 1
## [16885] 1 2 1 1 1 1 1 1 2 2 1 2 1 1 1 2 1 1 1 1 1 1 1 1 1 2 1 1 2 1 1 1 1 1 1 1
## [16921] 1 2 2 1 1 1 1 2 2 1 2 1 2 2 2 1 2 1 1 2 2 2 1 1 2 1 1 2 1 1 1 2 1 1 1 1
## [16957] 2 2 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 2 2 1 2 2 1 1 1 1 1 1 1 1 1 1
## [16993] 2 1 1 2 2 1 1 1 1 2 1 1 2 1 2 1 1 2 1 2 1 1 1 2 2 2 1 1 1 1 2 1 1 2 2 1
## [17029] 1 1 2 1 2 1 2 2 1 1 2 1 1 1 2 1 1 2 2 1 2 1 1 1 1 1 2 1 1 1 1 2 2 2 1 1
## [17065] 1 1 2 1 2 1 2 1 1 1 1 1 2 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 1 2 1 2 1 2 2 2
## [17101] 1 2 1 1 2 1 2 1 2 2 1 1 1 1 1 1 2 1 1 2 1 2 1 1 2 1 1 2 2 1 1 2 1 1 1 1
## [17137] 2 1 2 1 1 1 2 2 2 1 1 2 1 1 2 1 1 1 2 1 1 1 1 1 1 2 2 2 1 1 1 1 1 1 2 2
## [17173] 1 2 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 2 1 1 1 2 2 1 1 1 1 1 1 1 1 1 1 1
## [17209] 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1
## [17245] 1 1 1 1 1 1 1 2 2 1 1 1 2 1 1 1 1 1 1 2 1 1 1 1 1 2 2 1 2 1 1 1 1 1 1 2
## [17281] 1 2 1 1 2 2 2 2 2 1 1 2 1 1 1 1 1 2 1 2 2 1 2 1 1 1 2 1 1 1 1 1 1 2 1 1
## [17317] 1 2 1 1 2 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 2 1 1 1 2 1 1 1 1 1 1
## [17353] 1 1 1 1 1 2 1 2 1 2 1 2 2 2 1 1 1 1 2 1 1 2 1 1 1 1 2 1 1 2 2 1 1 1 1 2
## [17389] 1 1 2 1 1 2 2 1 1 2 1 1 2 1 1 1 1 1 1 1 1 1 2 2 1 2 1 2 2 2 1 2 2 2 1 1
## [17425] 2 1 1 2 2 2 2 2 1 2 1 2 2 2 1 1 1 2 1 2 1 2 1 2 2 2 2 1 2 2 2 2 2 1 1 1
## [17461] 2 2 1 1 2 2 1 2 1 2 1 1 1 1 1 1 2 2 1 1 1 1 1 2 1 2 2 1 1 2 1 1 1 1 2 2
## [17497] 1 1 1 1 2 2 1 1 2 2 2 1 1 2 2 1 2 1 1 1 1 1 1 1 2 1 1 1 2 2 1 1 2 2 1 1
## [17533] 1 1 2 1 1 2 2 1 1 2 1 1 1 1 1 1 1 1 1 1 2 2 1 2 2 1 1 2 1 1 1 1 2 2 2 1
## [17569] 1 1 2 2 2 1 2 2 1 2 1 2 1 2 2 1 1 1 1 2 1 2 1 2 2 2 2 1 2 2 2 1 2 1 2 2
## [17605] 1 1 1 2 1 1 1 1 2 1 2 1 2 2 1 2 1 1 1 1 1 1 2 1 1 2 2 2 2 1 1 2 1 1 1 1
## [17641] 1 1 1 1 2 2 2 1 2 2 1 2 1 2 2 2 1 1 1 2 2 1 1 2 1 1 2 1 1 1 2 1 1 1 1 1
## [17677] 1 1 1 2 1 1 2 1 2 2 1 1 2 2 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 2 1 2 1
## [17713] 1 2 1 2 1 2 2 1 2 1 1 1 1 2 1 1 1 2 1 1 1 1 1 1 1 2 1 1 1 1 2 1 2 1 2 1
## [17749] 1 1 1 1 1 2 2 1 2 2 1 1 1 2 1 1 1 1 1 2 2 1 1 2 1 2 1 1 2 1 1 1 1 1 2 2
## [17785] 1 1 2 2 1 1 1 2 1 2 2 1 2 2 2 1 1 1 1 1 2 1 1 2 1 2 2 1 1 1 1 2 1 1 1 2
## [17821] 2 1 1 1 2 1 1 2 1 2 1 1 1 1 1 1 1 1 1 1 1 2 1 1 2 1 1 1 1 1 1 1 1 1 1 1
## [17857] 2 1 1 1 1 1 1 1 2 1 2 1 2 1 1 2 1 2 1 2 2 1 2 2 1 1 1 2 1 1 2 2 2 1 1 1
## [17893] 2 2 1 2 1 1 2 1 1 2 1 1 1 1 2 2 2 2 1 1 1 2 1 2 2 2 2 1 1 1 1 1 1 1 1 2
## [17929] 1 1 2 2 1 1 1 1 2 1 1 1 1 2 1 2 2 1 1 1 2 1 2 2 1 1 2 1 1 1 1 1 2 1 2 1
## [17965] 2 1 1 1 2 2 2 1 2 2 1 2 1 1 1 1 2 1 1 1 1 2 1 1 2 2 1 1 1 1 1 1 1 1 1 1
## [18001] 1 1 1 1 1 1 2 2 1 1 1 1 1 2 2 1 1 1 2 1 1 1 1 2 1 2 1 1 1 2 1 1 1 1 1 1
## [18037] 1 2 2 2 2 2 1 1 2 2 1 1 1 1 1 1 1 1 2 1 1 2 1 1 2 2 2 1 1 2 1 2 2 1 1 1
## [18073] 2 1 1 1 1 2 1 2 2 2 2 1 1 1 1 1 2 1 2 1 2 1 1 1 1 1 2 1 2 1 1 2 1 1 1 1
## [18109] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 1 2 2 1 1 1 1 1 2
## [18145] 1 1 1 1 1 1 1 1 1 1 1 2 2 2 1 2 1 1 1 1 1 2 2 1 1 1 1 1 1 1 1 1 2 2 1 2
## [18181] 1 1 1 2 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 2 2 1 1 1 1 1 1 1 1 1 1 2 1 1 1
## [18217] 1 1 2 2 2 1 1 2 1 2 1 1 1 1 1 2 1 1 1 1 1 1 1 1 2 2 1 1 1 1 1 1 1 1 1 1
## [18253] 1 2 1 1 1 1 2 2 2 1 1 1 1 1 1 1 1 1 1 1 1 2 2 1 1 1 1 1 1 1 1 2 1 1 1 1
## [18289] 2 2 2 1 1 1 2 1 2 1 2 1 1 1 1 1 1 2 1 2 1 1 1 1 2 1 1 1 1 1 1 1 1 2 1 1
## [18325] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 2 1 1 1 1 1 1 2 2 1 1 2 1 2 1 1 1 1 1 2
## [18361] 2 1 1 1 2 1 2 2 1 1 1 1 1 1 2 2 1 2 2 1 1 1 1 1 1 2 2 1 1 2 1 1 2 1 1 1
## [18397] 1 1 2 1 1 1 2 2 1 2 1 1 1 1 2 2 2 2 2 1 2 1 1 1 1 1 2 1 1 2 1 1 1 2 1 1
## [18433] 1 1 2 1 1 1 2 1 1 1 1 1 2 2 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 1
## [18469] 1 2 1 2 1 1 1 2 1 1 1 1 1 2 1 2 2 2 1 1 1 1 2 1 1 2 1 1 1 2 1 1 1 1 1 2
## [18505] 1 2

#has another kid
uganda16$currentchildren<-cut(uganda16$v201, breaks = c(0,2,4, 20))
#regions
uganda16$region<-as.factor(uganda16$v024)
uganda16$regions<-car::Recode(uganda16$v024, recodes= "0='Kampala'; 1= 'South Buganda'; 2='North Buganda';3='Busoga';4='Bukedi'; 5='Bugisu'; 6='Teso';7='Karamoja'; 8='Lango'; 9='Acholi'; 10='West Nile'; 11='Bunyoro';12='Tooro'; 13='Ankole';14='Kigezi'", as.factor=T)

# survey design variables
uganda16$psu <- uganda16$v021
uganda16$strata <- uganda16$v022
uganda16$pwt <- uganda16$v005/1000000
desi<-svydesign(ids = ~ psu, strata = ~ strata, weights =~ pwt, data=uganda16)

1) Define a count outcome for the dataset of your choosing, the Area Resource File used in class provides many options here

State a research question about your outcome

At what rate are women in Uganda actively using contraception? In other words: how many women per women in child-rearing years are actively using contraception per region?

Is an offset term necessary? why or why not? If an offset term is need, what is the appropriate offset?

Yes, the offset term that is necessary in this context is the population of women in Uganda in their child-rearing years (ages 15-49). Offset terms are used in the Poisson model, to convert a count outcome into a rate.

#offset term
library(survey)
library(srvyr)

## 
## Attaching package: 'srvyr'

## The following object is masked from 'package:MASS':
## 
##     select

## The following object is masked from 'package:stats':
## 
##     filter

test<-uganda16%>%
  select(pwt, psu, strata,regions, agegroup, currentchildren, contraception,v313, moderncon)%>%
  filter(complete.cases(pwt, psu, strata, currentchildren, agegroup, regions, contraception, v313, moderncon))%>%
  #na.omit()%>%
  as_survey_design( ids = psu, 
                    strata =strata,
                    weights =pwt, nest=T)%>%
  
  group_by(regions, agegroup, currentchildren, contraception, .drop = FALSE)%>%
  summarise(nwom = round(survey_total()) )%>%
  ungroup()%>%
  filter(nwom>0)

contra <- test%>%
  filter(contraception=="modern method")%>%
  select(!contraception, womenmoderncon=nwom)

nocontra <- test%>%
  filter(contraception!="modern method")%>%
  select(!contraception, womennomoderncon = nwom)

joined<- dplyr::left_join(contra, nocontra,by=c( "currentchildren" ="currentchildren", "agegroup"="agegroup", "regions"= "regions"))

joined$nwomen = joined$womenmoderncon+joined$womennomoderncon

Consider a Poisson regression model for the outcome
1. Evaluate the level of dispersion in the outcome
2. Is the Poisson model a good choice?

#filter na
library(tidyverse)

## -- Attaching packages --------------------------------------- tidyverse 1.3.1 --

## v ggplot2 3.3.5     v purrr   0.3.4
## v tibble  3.1.6     v dplyr   1.0.7
## v tidyr   1.1.4     v stringr 1.4.0
## v readr   2.1.1     v forcats 0.5.1

## -- Conflicts ------------------------------------------ tidyverse_conflicts() --
## x tidyr::expand() masks Matrix::expand()
## x dplyr::filter() masks srvyr::filter(), stats::filter()
## x dplyr::lag()    masks stats::lag()
## x tidyr::pack()   masks Matrix::pack()
## x dplyr::select() masks srvyr::select(), MASS::select(), gtsummary::select()
## x tidyr::unpack() masks Matrix::unpack()

test %>% 
  drop_na()

## # A tibble: 678 x 6
##    regions agegroup currentchildren contraception     nwom nwom_se
##    <fct>   <fct>    <fct>           <fct>            <dbl>   <dbl>
##  1 Acholi  15-19    (0,2]           modern method       13       4
##  2 Acholi  15-19    (0,2]           none                33       4
##  3 Acholi  20-24    (0,2]           modern method       40       7
##  4 Acholi  20-24    (0,2]           none                72       7
##  5 Acholi  20-24    (2,4]           modern method        6       2
##  6 Acholi  20-24    (2,4]           none                23       4
##  7 Acholi  20-24    (4,20]          none                 2       1
##  8 Acholi  25-29    (0,2]           folkloric method     2       1
##  9 Acholi  25-29    (0,2]           modern method       10       2
## 10 Acholi  25-29    (0,2]           none                30       5
## # ... with 668 more rows

options(na.action = "na.omit")

#Here I use the original variables because the model will not run with factor varibles, which I how I recoded them


#poisson model
glmp_s<-glm(womenmoderncon ~offset(log(nwomen))+factor(regions)+currentchildren +factor(agegroup), data=joined, family=poisson)

summary(glmp_s)

## 
## Call:
## glm(formula = womenmoderncon ~ offset(log(nwomen)) + factor(regions) + 
##     currentchildren + factor(agegroup), family = poisson, data = joined)
## 
## Deviance Residuals: 
##     Min       1Q   Median       3Q      Max  
## -4.4227  -1.5117  -0.4502   1.8493   5.6648  
## 
## Coefficients:
##                              Estimate Std. Error z value Pr(>|z|)    
## (Intercept)                  -1.20101    0.07803 -15.392  < 2e-16 ***
## factor(regions)Ankole         0.27621    0.06873   4.019 5.85e-05 ***
## factor(regions)Bugisu         0.30183    0.07588   3.978 6.96e-05 ***
## factor(regions)Bukedi         0.29788    0.07170   4.154 3.26e-05 ***
## factor(regions)Bunyoro        0.05347    0.07870   0.679 0.496877    
## factor(regions)Busoga         0.08315    0.07087   1.173 0.240646    
## factor(regions)Kampala        0.30618    0.07579   4.040 5.34e-05 ***
## factor(regions)Karamoja      -1.21685    0.20137  -6.043 1.51e-09 ***
## factor(regions)Kigezi         0.36691    0.07849   4.675 2.94e-06 ***
## factor(regions)Lango          0.35625    0.07228   4.929 8.26e-07 ***
## factor(regions)North Buganda  0.36734    0.06602   5.564 2.64e-08 ***
## factor(regions)South Buganda  0.38240    0.06472   5.908 3.46e-09 ***
## factor(regions)Teso           0.13262    0.07578   1.750 0.080112 .  
## factor(regions)Tooro          0.30153    0.06937   4.347 1.38e-05 ***
## factor(regions)West Nile     -0.26410    0.08120  -3.252 0.001145 ** 
## currentchildren(2,4]          0.10376    0.03170   3.273 0.001065 ** 
## currentchildren(4,20]         0.15658    0.03794   4.127 3.67e-05 ***
## factor(agegroup)20-24         0.21772    0.05830   3.734 0.000188 ***
## factor(agegroup)25-29         0.25179    0.06054   4.159 3.19e-05 ***
## factor(agegroup)30-34         0.23185    0.06417   3.613 0.000303 ***
## factor(agegroup)35-39         0.17022    0.06787   2.508 0.012138 *  
## factor(agegroup)40-44         0.15445    0.07008   2.204 0.027524 *  
## factor(agegroup)45+          -0.36244    0.08123  -4.462 8.12e-06 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## (Dispersion parameter for poisson family taken to be 1)
## 
##     Null deviance: 2236.8  on 389  degrees of freedom
## Residual deviance: 1754.5  on 367  degrees of freedom
##   (5 observations deleted due to missingness)
## AIC: 3529.4
## 
## Number of Fisher Scoring iterations: 5

#check dispersion

1-pchisq(glmp_s$deviance, df=glmp_s$df.residual)

## [1] 0

#The p value is 0 which means this model does not fit the data.

#3) Consider a Negative binomial model

#negative binomial

init.coef = c(1, 14)
glmb<-glm.nb(womenmoderncon ~offset(log(nwomen))+ factor(regions) + factor(agegroup) +currentchildren, data=joined)
glmb%>%
  tbl_regression(exp=T)

Characteristic	IRR¹	95% CI¹	p-value
factor(regions)
Acholi	—	—
Ankole	1.29	0.95, 1.74	0.092
Bugisu	1.29	0.93, 1.78	0.12
Bukedi	1.35	0.99, 1.84	0.055
Bunyoro	1.08	0.78, 1.50	0.6
Busoga	1.13	0.82, 1.54	0.4
Kampala	1.28	0.93, 1.74	0.12
Karamoja	0.32	0.18, 0.53	<0.001
Kigezi	1.39	1.01, 1.91	0.042
Lango	1.34	0.98, 1.84	0.063
North Buganda	1.42	1.05, 1.90	0.020
South Buganda	1.44	1.07, 1.92	0.014
Teso	1.18	0.86, 1.62	0.3
Tooro	1.36	1.00, 1.85	0.045
West Nile	0.89	0.64, 1.24	0.5
factor(agegroup)
15-19	—	—
20-24	1.21	0.93, 1.56	0.14
25-29	1.13	0.87, 1.46	0.3
30-34	1.08	0.83, 1.41	0.5
35-39	0.95	0.72, 1.26	0.7
40-44	0.94	0.70, 1.24	0.6
45+	0.63	0.46, 0.86	0.003
currentchildren
(0,2]	—	—
(2,4]	1.20	1.04, 1.39	0.010
(4,20]	1.33	1.14, 1.56	<0.001
¹ IRR = Incidence Rate Ratio, CI = Confidence Interval

#4) Compare the model fits of the alternative models using AIC

#calculate AICs

AIC(glmp_s, glmb)

##        df      AIC
## glmp_s 23 3529.379
## glmb   24 2723.333

#In summary, the negative binomial is a better choice because of a lower AIC.

Luttinen HW 6

Rebecca Luttinen

3/10/2022

1) Define a count outcome for the dataset of your choosing, the Area Resource File used in class provides many options here