SEM immunisation data in Australia, 2015-16 and 2016-17
Alex Brooks
10/15/2018
First explain the model using 2016-17 immunisation data at postcode level, then build 2015-16 immunisation data at postcode level
We are trying to answer the question: Do geographic location, socio-economics, political affiliation or government funding have a relationship with immunisation coverage?
We will do this by using SEM on two years of postcode-level immunisation data to either confirm or deny whether a relationship exists.
Latent variables % immunisation - this is confusing because this is a manifest variable when used as pc_immun_class , but also a constructed latent variable as it’s the ‘target’ of interest. x =immunisation coverage (as measured by manifest variables of PC immun code
y = geo (Geographic areas or clusters, as measured by state, postcode, sa3, electorate) z = socio-economics (as measured by SEIFA score and mean tax) p = political (as measured by political score) f = funding (as measured by immunisation per person in that postcode or sa3)
#read in the data for 2016 only
Data <-read.csv('../cleaned_data/2016_data_REDUCED.csv')
str(Data)## 'data.frame': 7474 obs. of 18 variables:
## $ postcode : int 800 800 800 810 810 810 812 812 812 820 ...
## $ state : Factor w/ 12 levels "ACT","NSW","NSW/ACT",..: 6 6 6 6 6 6 6 6 6 6 ...
## $ year : int 2016 2016 2016 2016 2016 2016 2016 2016 2016 2016 ...
## $ age : int 1 2 5 5 1 2 1 5 2 2 ...
## $ pc_immun : Factor w/ 9 levels "<70.0","70.0-74.9",..: 6 3 6 6 7 6 6 7 5 4 ...
## $ caution : int 1 1 1 0 0 0 0 0 0 0 ...
## $ pc_immun_class : int 6 3 6 6 7 6 6 7 5 4 ...
## $ PHN_code : Factor w/ 31 levels "PHN101","PHN102",..: 30 30 30 30 30 30 30 30 30 30 ...
## $ PHN_number : int 701 701 701 701 701 701 701 701 701 701 ...
## $ IEO_SCORE : int 1089 1089 1089 1045 1045 1045 997 997 997 1085 ...
## $ IER_SCORE : int 946 946 946 1014 1014 1014 1013 1013 1013 1011 ...
## $ IRSAD_SCORE : int 1096 1096 1096 1052 1052 1052 1020 1020 1020 1094 ...
## $ IRSD_SCORE : int 1066 1066 1066 1037 1037 1037 1013 1013 1013 1073 ...
## $ SEIFA_URP : int 6464 6464 6464 33302 33302 33302 18873 18873 18873 19447 ...
## $ electorate : Factor w/ 150 levels "Adelaide","Aston",..: 136 136 136 136 136 136 136 136 136 136 ...
## $ mean_tax_000s : num 71.3 71.3 71.3 64.7 64.7 ...
## $ political_score : num 4 4 4 4 4 4 4 4 4 4 ...
## $ govt_funds_per_person: Factor w/ 2229 levels "#VALUE!","100",..: 837 837 837 249 249 249 1913 1913 1913 1947 ...View(Data)#calculate government funding per postcode
#Budget papers say the Feds spent $6.4m in 2016 on improving childhood immunisation rates
#Total population of Australia in 2016 was 24 598 900 according to the ABS
#We have the SEIFA URP which has the usual resident population per postcode
#First get per person amount
total_population <- (24598900)
total_immun_spend <- (6400000)
per_capita <- total_population/total_immun_spend
per_capita## [1] 3.843578#3.843578 spend per person# Remove NA values
Data <- Data[complete.cases(Data), ]How many postcodes are there in the 2016 set
2484 postcodes in 2016 150 electorates in 2016 31 PHNs in 2016 12 states in 2016
#count postcodes
length(unique(Data$postcode))## [1] 2484#count eletorates
length(unique(Data$electorate))## [1] 150#count PHNs
length(unique(Data$PHN_number))## [1] 31#count states - because there are doubles from the weird AIHW counting
length(unique(Data$state))## [1] 12#Need to make the factors ordered factors to work in Lavaan
Data$postcode = postcode_factor_ordered <- factor(Data$postcode, ordered = TRUE)
Data$state = state_factor_ordered <- factor(Data$state, ordered = TRUE)
Data$pc_immun = pc_immun_factor_ordered <- factor(Data$pc_immun, ordered = TRUE)
Data$PHN_code = PHN_code_factor_ordered <- factor(Data$PHN_code, ordered = TRUE)
Data$electorate = electorate_factor_ordered <- factor(Data$electorate, ordered = TRUE)
Data$govt_funds_per_person = funds_ordered <- as.numeric(Data$govt_funds_per_person)
glimpse(Data)## Observations: 7,286
## Variables: 18
## $ postcode <ord> 800, 800, 800, 810, 810, 810, 812, 812, ...
## $ state <ord> NT, NT, NT, NT, NT, NT, NT, NT, NT, NT, ...
## $ year <int> 2016, 2016, 2016, 2016, 2016, 2016, 2016...
## $ age <int> 1, 2, 5, 5, 1, 2, 1, 5, 2, 2, 5, 1, 1, 5...
## $ pc_immun <ord> 90.0-92.4, 75.0-79.9, 90.0-92.4, 90.0-92...
## $ caution <int> 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0...
## $ pc_immun_class <int> 6, 3, 6, 6, 7, 6, 6, 7, 5, 4, 7, 6, 7, 7...
## $ PHN_code <ord> PHN701, PHN701, PHN701, PHN701, PHN701, ...
## $ PHN_number <int> 701, 701, 701, 701, 701, 701, 701, 701, ...
## $ IEO_SCORE <int> 1089, 1089, 1089, 1045, 1045, 1045, 997,...
## $ IER_SCORE <int> 946, 946, 946, 1014, 1014, 1014, 1013, 1...
## $ IRSAD_SCORE <int> 1096, 1096, 1096, 1052, 1052, 1052, 1020...
## $ IRSD_SCORE <int> 1066, 1066, 1066, 1037, 1037, 1037, 1013...
## $ SEIFA_URP <int> 6464, 6464, 6464, 33302, 33302, 33302, 1...
## $ electorate <ord> Solomon, Solomon, Solomon, Solomon, Solo...
## $ mean_tax_000s <dbl> 71.263, 71.263, 71.263, 64.719, 64.719, ...
## $ political_score <dbl> 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4...
## $ govt_funds_per_person <dbl> 837, 837, 837, 249, 249, 249, 1913, 1913...#filter just by 5 year olds - as the per capita info in spreadsheet is duplicating by age group
Five_Data = Data%>%
filter(age == 5)
str(Five_Data)## 'data.frame': 2434 obs. of 18 variables:
## $ postcode : Ord.factor w/ 2484 levels "800"<"810"<"812"<..: 1 2 3 4 5 6 7 8 9 10 ...
## $ state : Ord.factor w/ 12 levels "ACT"<"NSW"<"NSW/ACT"<..: 6 6 6 6 6 6 6 6 6 6 ...
## $ year : int 2016 2016 2016 2016 2016 2016 2016 2016 2016 2016 ...
## $ age : int 5 5 5 5 5 5 5 5 5 5 ...
## $ pc_immun : Ord.factor w/ 9 levels "<70.0"<"70.0-74.9"<..: 6 6 7 7 7 9 7 6 6 8 ...
## $ caution : int 1 0 0 0 0 0 0 0 1 0 ...
## $ pc_immun_class : int 6 6 7 7 7 0 7 6 6 8 ...
## $ PHN_code : Ord.factor w/ 31 levels "PHN101"<"PHN102"<..: 30 30 30 30 30 30 30 30 30 30 ...
## $ PHN_number : int 701 701 701 701 701 701 701 701 701 701 ...
## $ IEO_SCORE : int 1089 1045 997 1085 880 993 967 1006 1000 969 ...
## $ IER_SCORE : int 946 1014 1013 1011 658 942 990 1060 1120 1114 ...
## $ IRSAD_SCORE : int 1096 1052 1020 1094 754 979 997 1075 1070 1055 ...
## $ IRSD_SCORE : int 1066 1037 1013 1073 650 970 989 1071 1072 1049 ...
## $ SEIFA_URP : int 6464 33302 18873 19447 25560 1201 18207 15570 5132 7756 ...
## $ electorate : Ord.factor w/ 150 levels "Adelaide"<"Aston"<..: 136 136 136 136 90 136 136 136 90 90 ...
## $ mean_tax_000s : num 71.3 64.7 63.3 79.6 45.9 ...
## $ political_score : num 4 4 4 4 4 4 4 4 4 4 ...
## $ govt_funds_per_person: num 837 249 1913 1947 2214 ...First 2016 model - assessing a real value against the other real values (IE, not latent yet)
# single factor model imm =~ pc_immun_class + IRSAD_SCORE + political_score + govt_funds_per_person + mean_tax_000s
#let's build 2016pc.model the standard way
pc2016.model <- 'imm =~ pc_immun_class + IRSAD_SCORE + political_score + govt_funds_per_person + mean_tax_000s'
pc2016.fit <- cfa(model=pc2016.model,
data = Five_Data, std.lv = TRUE)## Warning in lav_data_full(data = data, group = group, cluster = cluster, :
## lavaan WARNING: some observed variances are (at least) a factor 1000 times
## larger than others; use varTable(fit) to investigate## Warning in lav_object_post_check(object): lavaan WARNING: some estimated ov
## variances are negativesummary(pc2016.fit, standardized = TRUE, fit.measures = TRUE)## lavaan 0.6-3 ended normally after 39 iterations
##
## Optimization method NLMINB
## Number of free parameters 10
##
## Number of observations 2434
##
## Estimator ML
## Model Fit Test Statistic 105.776
## Degrees of freedom 5
## P-value (Chi-square) 0.000
##
## Model test baseline model:
##
## Minimum Function Test Statistic 1692.805
## Degrees of freedom 10
## P-value 0.000
##
## User model versus baseline model:
##
## Comparative Fit Index (CFI) 0.940
## Tucker-Lewis Index (TLI) 0.880
##
## Loglikelihood and Information Criteria:
##
## Loglikelihood user model (H0) -53841.964
## Loglikelihood unrestricted model (H1) -53789.076
##
## Number of free parameters 10
## Akaike (AIC) 107703.928
## Bayesian (BIC) 107761.901
## Sample-size adjusted Bayesian (BIC) 107730.128
##
## Root Mean Square Error of Approximation:
##
## RMSEA 0.091
## 90 Percent Confidence Interval 0.076 0.106
## P-value RMSEA <= 0.05 0.000
##
## Standardized Root Mean Square Residual:
##
## SRMR 0.046
##
## Parameter Estimates:
##
## Information Expected
## Information saturated (h1) model Structured
## Standard Errors Standard
##
## Latent Variables:
## Estimate Std.Err z-value P(>|z|) Std.lv
## imm =~
## pc_immun_class 0.775 0.077 10.013 0.000 0.775
## IRSAD_SCORE 50.348 2.628 19.158 0.000 50.348
## political_scor -0.077 0.033 -2.297 0.022 -0.077
## gvt_fnds_pr_pr 60.321 13.239 4.556 0.000 60.321
## mean_tax_000s 16.809 0.788 21.318 0.000 16.809
## Std.all
##
## 0.224
## 0.664
## -0.047
## 0.094
## 1.005
##
## Variances:
## Estimate Std.Err z-value P(>|z|) Std.lv
## .pc_immun_class 11.380 0.331 34.425 0.000 11.380
## .IRSAD_SCORE 3212.642 244.649 13.132 0.000 3212.642
## .political_scor 2.702 0.077 34.885 0.000 2.702
## .gvt_fnds_pr_pr 410281.632 11764.479 34.875 0.000 410281.632
## .mean_tax_000s -2.577 25.263 -0.102 0.919 -2.577
## imm 1.000 1.000
## Std.all
## 0.950
## 0.559
## 0.998
## 0.991
## -0.009
## 1.000#let's look at the variance of this dataset
varTable(pc2016.fit)## name idx nobs type exo user mean var nlev
## 1 pc_immun_class 7 2434 numeric 0 0 4.173 11.985 0
## 2 IRSAD_SCORE 12 2434 numeric 0 0 986.465 5749.967 0
## 3 political_score 17 2434 numeric 0 0 5.095 2.709 0
## 4 govt_funds_per_person 18 2434 numeric 0 0 1103.384 414090.352 0
## 5 mean_tax_000s 16 2434 numeric 0 0 53.675 280.083 0
## lnam
## 1
## 2
## 3
## 4
## 5lavInspect(pc2016.fit, what="estimates")## $lambda
## imm
## pc_immun_class 0.775
## IRSAD_SCORE 50.348
## political_score -0.077
## govt_funds_per_person 60.321
## mean_tax_000s 16.809
##
## $theta
## pc_mm_ IRSAD_ pltcl_ gvt___
## pc_immun_class 11.380
## IRSAD_SCORE 0.000 3212.642
## political_score 0.000 0.000 2.702
## govt_funds_per_person 0.000 0.000 0.000 410281.632
## mean_tax_000s 0.000 0.000 0.000 0.000
## m__000
## pc_immun_class
## IRSAD_SCORE
## political_score
## govt_funds_per_person
## mean_tax_000s -2.577
##
## $psi
## imm
## imm 1semPaths(pc2016.fit, title=TRUE, exoVar = TRUE, exoCov = TRUE)
semPaths(pc2016.fit, what = "stand", rotation = 4)
#
semCors(pc2016.fit, titles = TRUE)
#let's build to account for missing pc immun data
pc2016.fit2 <- cfa(model=pc2016.model,
data = Five_Data, estimator = "MLR", missing = "ml")## Warning in lav_data_full(data = data, group = group, cluster = cluster, :
## lavaan WARNING: some observed variances are (at least) a factor 1000 times
## larger than others; use varTable(fit) to investigate## Warning in lav_object_post_check(object): lavaan WARNING: some estimated ov
## variances are negativesummary(pc2016.fit2, standardized = TRUE, fit.measures = TRUE)## lavaan 0.6-3 ended normally after 51 iterations
##
## Optimization method NLMINB
## Number of free parameters 15
##
## Number of observations 2434
## Number of missing patterns 1
##
## Estimator ML Robust
## Model Fit Test Statistic 105.776 101.520
## Degrees of freedom 5 5
## P-value (Chi-square) 0.000 0.000
## Scaling correction factor 1.042
## for the Yuan-Bentler correction (Mplus variant)
##
## Model test baseline model:
##
## Minimum Function Test Statistic 1692.805 1115.088
## Degrees of freedom 10 10
## P-value 0.000 0.000
##
## User model versus baseline model:
##
## Comparative Fit Index (CFI) 0.940 0.913
## Tucker-Lewis Index (TLI) 0.880 0.825
##
## Robust Comparative Fit Index (CFI) 0.940
## Robust Tucker-Lewis Index (TLI) 0.880
##
## Loglikelihood and Information Criteria:
##
## Loglikelihood user model (H0) -53841.964 -53841.964
## Scaling correction factor 1.734
## for the MLR correction
## Loglikelihood unrestricted model (H1) -53789.076 -53789.076
## Scaling correction factor 1.561
## for the MLR correction
##
## Number of free parameters 15 15
## Akaike (AIC) 107713.928 107713.928
## Bayesian (BIC) 107800.887 107800.887
## Sample-size adjusted Bayesian (BIC) 107753.229 107753.229
##
## Root Mean Square Error of Approximation:
##
## RMSEA 0.091 0.089
## 90 Percent Confidence Interval 0.076 0.106 0.075 0.104
## P-value RMSEA <= 0.05 0.000 0.000
##
## Robust RMSEA 0.091
## 90 Percent Confidence Interval 0.076 0.107
##
## Standardized Root Mean Square Residual:
##
## SRMR 0.040 0.040
##
## Parameter Estimates:
##
## Information Observed
## Observed information based on Hessian
## Standard Errors Robust.huber.white
##
## Latent Variables:
## Estimate Std.Err z-value P(>|z|) Std.lv
## imm =~
## pc_immun_class 1.000 0.775
## IRSAD_SCORE 64.955 5.147 12.621 0.000 50.348
## political_scor -0.099 0.061 -1.615 0.106 -0.077
## gvt_fnds_pr_pr 77.820 15.656 4.971 0.000 60.321
## mean_tax_000s 21.686 3.092 7.014 0.000 16.809
## Std.all
##
## 0.224
## 0.664
## -0.047
## 0.094
## 1.005
##
## Intercepts:
## Estimate Std.Err z-value P(>|z|) Std.lv
## .pc_immun_class 4.173 0.070 59.480 0.000 4.173
## .IRSAD_SCORE 986.465 1.537 641.946 0.000 986.465
## .political_scor 5.095 0.033 152.736 0.000 5.095
## .gvt_fnds_pr_pr 1103.384 13.041 84.611 0.000 1103.384
## .mean_tax_000s 53.675 0.339 158.262 0.000 53.675
## imm 0.000 0.000
## Std.all
## 1.206
## 13.012
## 3.096
## 1.715
## 3.208
## 0.000
##
## Variances:
## Estimate Std.Err z-value P(>|z|) Std.lv
## .pc_immun_class 11.380 0.158 71.935 0.000 11.380
## .IRSAD_SCORE 3212.642 347.557 9.243 0.000 3212.642
## .political_scor 2.702 0.116 23.378 0.000 2.702
## .gvt_fnds_pr_pr 410281.636 7628.110 53.785 0.000 410281.636
## .mean_tax_000s -2.577 37.371 -0.069 0.945 -2.577
## imm 0.601 0.115 5.207 0.000 1.000
## Std.all
## 0.950
## 0.559
## 0.998
## 0.991
## -0.009
## 1.000#THIS IS THE ONE TO PUBLISH AS MODEL FIT#Let's bootstrap the model
pc2016.fit3 <- cfa(model=pc2016.model,
data = Five_Data, se = "bootstrap", verbose = TRUE, bootstrap = 1000)## Warning in lav_data_full(data = data, group = group, cluster = cluster, :
## lavaan WARNING: some observed variances are (at least) a factor 1000 times
## larger than others; use varTable(fit) to investigate## Quasi-Newton steps using NLMINB:
## Objective function = 1.0166881031435544
## Objective function = 0.2867231627086966
## Objective function = 0.1531022526498234
## Objective function = 0.0776357860488197
## Objective function = 0.0312178394323652
## Objective function = 0.0281607306461993
## Objective function = 0.0273120835943246
## Objective function = 0.0269954144078639
## Objective function = 0.0267756853596808
## Objective function = 0.0259536809907708
## Objective function = 0.0246474581795404
## Objective function = 0.0239694236687100
## Objective function = 0.0234478974084453
## Objective function = 0.0233481242915055
## Objective function = 0.0232638422559717
## Objective function = 0.0231289163456303
## Objective function = 0.0229430889609077
## Objective function = 0.0227408280796837
## Objective function = 0.0225357089505600
## Objective function = 0.0223407311577404
## Objective function = 0.0221711646327165
## Objective function = 0.0221178137388520
## Objective function = 0.0220226188845274
## Objective function = 0.0219620891379222
## Objective function = 0.0218897960714095
## Objective function = 0.0218546166280547
## Objective function = 0.0218379200646588
## Objective function = 0.0218269488756420
## Objective function = 0.0218206000170582
## Objective function = 0.0218129429156022
## Objective function = 0.0218092931344280
## Objective function = 0.0218048735948670
## Objective function = 0.0218002499188508
## Objective function = 0.0217947282285031
## Objective function = 0.0217899912266191
## Objective function = 0.0217846001301467
## Objective function = 0.0217728166456830
## Objective function = 0.0217615706568779
## Objective function = 0.0217418819328650
## Objective function = 0.0217330868238896
## Objective function = 0.0217297350647208
## Objective function = 0.0217291528116554
## Objective function = 0.0217289353600094
## Objective function = 0.0217288264538453
## Objective function = 0.0217287795261107
## Objective function = 0.0217287664417221
## Objective function = 0.0217287614934278
## Objective function = 0.0217287585929213
## Objective function = 0.0217287571105373
## Objective function = 0.0217287567487148
## Objective function = 0.0217287566879492
## Objective function = 0.0217287566864073
## Objective function = 0.0217287566864073
## convergence status (0=ok): 0
## nlminb message says: relative convergence (4)
## number of iterations: 51
## number of function evaluations [objective, gradient]: 52 52
## Computing VCOV for se = bootstrap ...
## ... bootstrap draw number: 1 OK -- niter = 42 fx = 0.029982193
## ... bootstrap draw number: 2 OK -- niter = 42 fx = 0.029384332
## ... bootstrap draw number: 3 OK -- niter = 46 fx = 0.018054687
## ... bootstrap draw number: 4 OK -- niter = 43 fx = 0.024094823
## ... bootstrap draw number: 5 OK -- niter = 46 fx = 0.024348647
## ... bootstrap draw number: 6 OK -- niter = 46 fx = 0.023893611
## ... bootstrap draw number: 7 OK -- niter = 43 fx = 0.021292824
## ... bootstrap draw number: 8 OK -- niter = 46 fx = 0.023640633
## ... bootstrap draw number: 9 OK -- niter = 48 fx = 0.017592670
## ... bootstrap draw number: 10 OK -- niter = 43 fx = 0.027638660
## ... bootstrap draw number: 11 OK -- niter = 50 fx = 0.026402997
## ... bootstrap draw number: 12 OK -- niter = 49 fx = 0.027095849
## ... bootstrap draw number: 13 OK -- niter = 47 fx = 0.026910756
## ... bootstrap draw number: 14 OK -- niter = 53 fx = 0.024617091
## ... bootstrap draw number: 15 OK -- niter = 46 fx = 0.022211748
## ... bootstrap draw number: 16 OK -- niter = 41 fx = 0.031776607
## ... bootstrap draw number: 17 OK -- niter = 42 fx = 0.017727955
## ... bootstrap draw number: 18 OK -- niter = 44 fx = 0.018346384
## ... bootstrap draw number: 19 OK -- niter = 61 fx = 0.027105879
## ... bootstrap draw number: 20 OK -- niter = 52 fx = 0.023123475
## ... bootstrap draw number: 21 OK -- niter = 47 fx = 0.022047154
## ... bootstrap draw number: 22 OK -- niter = 59 fx = 0.022066303
## ... bootstrap draw number: 23 OK -- niter = 45 fx = 0.023022773
## ... bootstrap draw number: 24 OK -- niter = 47 fx = 0.017048090
## ... bootstrap draw number: 25 OK -- niter = 47 fx = 0.027010227
## ... bootstrap draw number: 26 OK -- niter = 43 fx = 0.025097586
## ... bootstrap draw number: 27 OK -- niter = 46 fx = 0.028399291
## ... bootstrap draw number: 28 OK -- niter = 54 fx = 0.022796773
## ... bootstrap draw number: 29 OK -- niter = 48 fx = 0.024925931
## ... bootstrap draw number: 30 OK -- niter = 48 fx = 0.018317072
## ... bootstrap draw number: 31 OK -- niter = 46 fx = 0.032989417
## ... bootstrap draw number: 32 OK -- niter = 44 fx = 0.023747637
## ... bootstrap draw number: 33 OK -- niter = 44 fx = 0.017307917
## ... bootstrap draw number: 34 OK -- niter = 45 fx = 0.021400171
## ... bootstrap draw number: 35 OK -- niter = 43 fx = 0.021162171
## ... bootstrap draw number: 36 OK -- niter = 48 fx = 0.014612214
## ... bootstrap draw number: 37 OK -- niter = 42 fx = 0.013774878
## ... bootstrap draw number: 38 OK -- niter = 53 fx = 0.019808431
## ... bootstrap draw number: 39 OK -- niter = 50 fx = 0.016288481
## ... bootstrap draw number: 40 OK -- niter = 44 fx = 0.023556173
## ... bootstrap draw number: 41 OK -- niter = 45 fx = 0.022998350
## ... bootstrap draw number: 42 OK -- niter = 43 fx = 0.021380551
## ... bootstrap draw number: 43 OK -- niter = 47 fx = 0.022925923
## ... bootstrap draw number: 44 OK -- niter = 45 fx = 0.022190647
## ... bootstrap draw number: 45 OK -- niter = 47 fx = 0.027332594
## ... bootstrap draw number: 46 OK -- niter = 47 fx = 0.022468579
## ... bootstrap draw number: 47 OK -- niter = 51 fx = 0.015335799
## ... bootstrap draw number: 48 OK -- niter = 47 fx = 0.014256009
## ... bootstrap draw number: 49 OK -- niter = 42 fx = 0.022597414
## ... bootstrap draw number: 50 OK -- niter = 45 fx = 0.019745830
## ... bootstrap draw number: 51 OK -- niter = 43 fx = 0.024567894
## ... bootstrap draw number: 52 OK -- niter = 188 fx = 0.036778471
## ... bootstrap draw number: 53 OK -- niter = 46 fx = 0.022398670
## ... bootstrap draw number: 54 OK -- niter = 43 fx = 0.032006588
## ... bootstrap draw number: 55 OK -- niter = 48 fx = 0.021156661
## ... bootstrap draw number: 56 OK -- niter = 42 fx = 0.016272245
## ... bootstrap draw number: 57 OK -- niter = 47 fx = 0.011176298
## ... bootstrap draw number: 58 OK -- niter = 44 fx = 0.025470513
## ... bootstrap draw number: 59 OK -- niter = 43 fx = 0.020411646
## ... bootstrap draw number: 60 OK -- niter = 47 fx = 0.022611200
## ... bootstrap draw number: 61 OK -- niter = 58 fx = 0.025683491
## ... bootstrap draw number: 62 OK -- niter = 47 fx = 0.019698802
## ... bootstrap draw number: 63 OK -- niter = 47 fx = 0.021417546
## ... bootstrap draw number: 64 OK -- niter = 44 fx = 0.024706439
## ... bootstrap draw number: 65 OK -- niter = 41 fx = 0.023703386
## ... bootstrap draw number: 66 OK -- niter = 45 fx = 0.022244373
## ... bootstrap draw number: 67 OK -- niter = 46 fx = 0.024493148
## ... bootstrap draw number: 68 OK -- niter = 49 fx = 0.021760268
## ... bootstrap draw number: 69 OK -- niter = 46 fx = 0.026972388
## ... bootstrap draw number: 70 OK -- niter = 48 fx = 0.014556179
## ... bootstrap draw number: 71 OK -- niter = 47 fx = 0.024517165
## ... bootstrap draw number: 72 OK -- niter = 52 fx = 0.025870059
## ... bootstrap draw number: 73 OK -- niter = 45 fx = 0.030357971
## ... bootstrap draw number: 74 OK -- niter = 48 fx = 0.019858045
## ... bootstrap draw number: 75 OK -- niter = 49 fx = 0.021380493
## ... bootstrap draw number: 76 OK -- niter = 42 fx = 0.028572569
## ... bootstrap draw number: 77 OK -- niter = 46 fx = 0.021112907
## ... bootstrap draw number: 78 OK -- niter = 45 fx = 0.026307720
## ... bootstrap draw number: 79 OK -- niter = 42 fx = 0.022857750
## ... bootstrap draw number: 80 OK -- niter = 40 fx = 0.026411133
## ... bootstrap draw number: 81 OK -- niter = 50 fx = 0.021519750
## ... bootstrap draw number: 82 OK -- niter = 46 fx = 0.021278257
## ... bootstrap draw number: 83 OK -- niter = 58 fx = 0.026191285
## ... bootstrap draw number: 84 OK -- niter = 51 fx = 0.017180759
## ... bootstrap draw number: 85 OK -- niter = 46 fx = 0.018447436
## ... bootstrap draw number: 86 OK -- niter = 51 fx = 0.027233075
## ... bootstrap draw number: 87 OK -- niter = 44 fx = 0.014956792
## ... bootstrap draw number: 88 OK -- niter = 49 fx = 0.019138559
## ... bootstrap draw number: 89 OK -- niter = 48 fx = 0.032661419
## ... bootstrap draw number: 90 OK -- niter = 44 fx = 0.021761808
## ... bootstrap draw number: 91 OK -- niter = 48 fx = 0.019266479
## ... bootstrap draw number: 92 OK -- niter = 56 fx = 0.025045980
## ... bootstrap draw number: 93 OK -- niter = 48 fx = 0.024892243
## ... bootstrap draw number: 94 OK -- niter = 45 fx = 0.022436821
## ... bootstrap draw number: 95 OK -- niter = 50 fx = 0.026855130
## ... bootstrap draw number: 96 OK -- niter = 49 fx = 0.019180605
## ... bootstrap draw number: 97 OK -- niter = 53 fx = 0.018299653
## ... bootstrap draw number: 98 OK -- niter = 44 fx = 0.024664375
## ... bootstrap draw number: 99 OK -- niter = 44 fx = 0.023783673
## ... bootstrap draw number: 100 OK -- niter = 53 fx = 0.016499351
## ... bootstrap draw number: 101 OK -- niter = 49 fx = 0.011473554
## ... bootstrap draw number: 102 OK -- niter = 48 fx = 0.022172165
## ... bootstrap draw number: 103 OK -- niter = 42 fx = 0.019082588
## ... bootstrap draw number: 104 OK -- niter = 46 fx = 0.027200536
## ... bootstrap draw number: 105 OK -- niter = 46 fx = 0.019838468
## ... bootstrap draw number: 106 OK -- niter = 48 fx = 0.018939141
## ... bootstrap draw number: 107 OK -- niter = 49 fx = 0.027115429
## ... bootstrap draw number: 108 OK -- niter = 45 fx = 0.020825009
## ... bootstrap draw number: 109 OK -- niter = 45 fx = 0.025581218
## ... bootstrap draw number: 110 OK -- niter = 44 fx = 0.020316512
## ... bootstrap draw number: 111 OK -- niter = 43 fx = 0.019989066
## ... bootstrap draw number: 112 OK -- niter = 46 fx = 0.016412617
## ... bootstrap draw number: 113 OK -- niter = 41 fx = 0.015867154
## ... bootstrap draw number: 114 OK -- niter = 50 fx = 0.021039412
## ... bootstrap draw number: 115 OK -- niter = 49 fx = 0.019035208
## ... bootstrap draw number: 116 OK -- niter = 47 fx = 0.017212385
## ... bootstrap draw number: 117 OK -- niter = 48 fx = 0.024253714
## ... bootstrap draw number: 118 OK -- niter = 48 fx = 0.018858927
## ... bootstrap draw number: 119 OK -- niter = 46 fx = 0.021142933
## ... bootstrap draw number: 120 OK -- niter = 48 fx = 0.026816033
## ... bootstrap draw number: 121 OK -- niter = 43 fx = 0.022557792
## ... bootstrap draw number: 122 OK -- niter = 43 fx = 0.026271338
## ... bootstrap draw number: 123 OK -- niter = 55 fx = 0.019493871
## ... bootstrap draw number: 124 OK -- niter = 47 fx = 0.022095892
## ... bootstrap draw number: 125 OK -- niter = 52 fx = 0.024842061
## ... bootstrap draw number: 126 OK -- niter = 43 fx = 0.020241280
## ... bootstrap draw number: 127 OK -- niter = 49 fx = 0.021839287
## ... bootstrap draw number: 128 OK -- niter = 45 fx = 0.021765424
## ... bootstrap draw number: 129 OK -- niter = 26 fx = 0.012828980
## ... bootstrap draw number: 130 OK -- niter = 46 fx = 0.013515362
## ... bootstrap draw number: 131 OK -- niter = 42 fx = 0.034741067
## ... bootstrap draw number: 132 OK -- niter = 46 fx = 0.019035237
## ... bootstrap draw number: 133 OK -- niter = 53 fx = 0.025206171
## ... bootstrap draw number: 134 OK -- niter = 43 fx = 0.020261790
## ... bootstrap draw number: 135 OK -- niter = 53 fx = 0.023875574
## ... bootstrap draw number: 136 OK -- niter = 50 fx = 0.025101183
## ... bootstrap draw number: 137 OK -- niter = 46 fx = 0.019531643
## ... bootstrap draw number: 138 OK -- niter = 52 fx = 0.019982620
## ... bootstrap draw number: 139 OK -- niter = 43 fx = 0.031400463
## ... bootstrap draw number: 140 OK -- niter = 47 fx = 0.029205199
## ... bootstrap draw number: 141 OK -- niter = 50 fx = 0.019377641
## ... bootstrap draw number: 142 OK -- niter = 46 fx = 0.023956728
## ... bootstrap draw number: 143 OK -- niter = 46 fx = 0.018931945
## ... bootstrap draw number: 144 OK -- niter = 58 fx = 0.022264552
## ... bootstrap draw number: 145 OK -- niter = 42 fx = 0.022710189
## ... bootstrap draw number: 146 OK -- niter = 44 fx = 0.022431688
## ... bootstrap draw number: 147 OK -- niter = 48 fx = 0.022581355
## ... bootstrap draw number: 148 OK -- niter = 51 fx = 0.017762831
## ... bootstrap draw number: 149 OK -- niter = 49 fx = 0.032647767
## ... bootstrap draw number: 150 OK -- niter = 53 fx = 0.019998221
## ... bootstrap draw number: 151 OK -- niter = 51 fx = 0.017612827
## ... bootstrap draw number: 152 OK -- niter = 45 fx = 0.013681569
## ... bootstrap draw number: 153 OK -- niter = 46 fx = 0.027851723
## ... bootstrap draw number: 154 OK -- niter = 45 fx = 0.026388994
## ... bootstrap draw number: 155 OK -- niter = 49 fx = 0.023068147
## ... bootstrap draw number: 156 OK -- niter = 42 fx = 0.022299088
## ... bootstrap draw number: 157 OK -- niter = 45 fx = 0.023127189
## ... bootstrap draw number: 158 OK -- niter = 23 fx = 0.025824096
## ... bootstrap draw number: 159 OK -- niter = 42 fx = 0.024808715
## ... bootstrap draw number: 160 OK -- niter = 45 fx = 0.019806886
## ... bootstrap draw number: 161 OK -- niter = 44 fx = 0.019731264
## ... bootstrap draw number: 162 OK -- niter = 52 fx = 0.021027886
## ... bootstrap draw number: 163 OK -- niter = 47 fx = 0.023212963
## ... bootstrap draw number: 164 OK -- niter = 45 fx = 0.023891533
## ... bootstrap draw number: 165 OK -- niter = 47 fx = 0.017494031
## ... bootstrap draw number: 166 OK -- niter = 44 fx = 0.026106735
## ... bootstrap draw number: 167 OK -- niter = 43 fx = 0.023595589
## ... bootstrap draw number: 168 OK -- niter = 45 fx = 0.020448029
## ... bootstrap draw number: 169 OK -- niter = 44 fx = 0.028124098
## ... bootstrap draw number: 170 OK -- niter = 52 fx = 0.028017131
## ... bootstrap draw number: 171 OK -- niter = 45 fx = 0.028104945
## ... bootstrap draw number: 172 OK -- niter = 52 fx = 0.024853135
## ... bootstrap draw number: 173 OK -- niter = 48 fx = 0.024149470
## ... bootstrap draw number: 174 OK -- niter = 44 fx = 0.027876517
## ... bootstrap draw number: 175 OK -- niter = 48 fx = 0.021557563
## ... bootstrap draw number: 176 OK -- niter = 47 fx = 0.020953052
## ... bootstrap draw number: 177 OK -- niter = 44 fx = 0.015949309
## ... bootstrap draw number: 178 OK -- niter = 52 fx = 0.032406138
## ... bootstrap draw number: 179 OK -- niter = 42 fx = 0.022423502
## ... bootstrap draw number: 180 OK -- niter = 43 fx = 0.022294598
## ... bootstrap draw number: 181 OK -- niter = 44 fx = 0.016344643
## ... bootstrap draw number: 182 OK -- niter = 53 fx = 0.022936600
## ... bootstrap draw number: 183 OK -- niter = 50 fx = 0.021956779
## ... bootstrap draw number: 184 OK -- niter = 46 fx = 0.023709670
## ... bootstrap draw number: 185 OK -- niter = 45 fx = 0.029392422
## ... bootstrap draw number: 186 OK -- niter = 47 fx = 0.021092099
## ... bootstrap draw number: 187 OK -- niter = 45 fx = 0.025498630
## ... bootstrap draw number: 188 OK -- niter = 47 fx = 0.020074246
## ... bootstrap draw number: 189 OK -- niter = 42 fx = 0.028109698
## ... bootstrap draw number: 190 OK -- niter = 47 fx = 0.021762687
## ... bootstrap draw number: 191 OK -- niter = 45 fx = 0.021314121
## ... bootstrap draw number: 192 OK -- niter = 49 fx = 0.025886428
## ... bootstrap draw number: 193 OK -- niter = 49 fx = 0.023913106
## ... bootstrap draw number: 194 OK -- niter = 54 fx = 0.026085207
## ... bootstrap draw number: 195 OK -- niter = 43 fx = 0.027523750
## ... bootstrap draw number: 196 OK -- niter = 47 fx = 0.024329394
## ... bootstrap draw number: 197 OK -- niter = 43 fx = 0.015940493
## ... bootstrap draw number: 198 OK -- niter = 49 fx = 0.024774647
## ... bootstrap draw number: 199 OK -- niter = 46 fx = 0.018052549
## ... bootstrap draw number: 200 OK -- niter = 46 fx = 0.024309469
## ... bootstrap draw number: 201 OK -- niter = 42 fx = 0.019757950
## ... bootstrap draw number: 202 OK -- niter = 45 fx = 0.020983579
## ... bootstrap draw number: 203 OK -- niter = 46 fx = 0.024453789
## ... bootstrap draw number: 204 OK -- niter = 44 fx = 0.024619201
## ... bootstrap draw number: 205 OK -- niter = 45 fx = 0.026794802
## ... bootstrap draw number: 206 OK -- niter = 46 fx = 0.017184942
## ... bootstrap draw number: 207 OK -- niter = 45 fx = 0.021354185
## ... bootstrap draw number: 208 OK -- niter = 28 fx = 0.025439349
## ... bootstrap draw number: 209 OK -- niter = 44 fx = 0.018487643
## ... bootstrap draw number: 210 OK -- niter = 43 fx = 0.027386957
## ... bootstrap draw number: 211 OK -- niter = 51 fx = 0.025718442
## ... bootstrap draw number: 212 OK -- niter = 45 fx = 0.022490433
## ... bootstrap draw number: 213 OK -- niter = 41 fx = 0.018075633
## ... bootstrap draw number: 214 OK -- niter = 43 fx = 0.032481973
## ... bootstrap draw number: 215 OK -- niter = 47 fx = 0.027094255
## ... bootstrap draw number: 216 OK -- niter = 45 fx = 0.025598597
## ... bootstrap draw number: 217 OK -- niter = 26 fx = 0.026739610
## ... bootstrap draw number: 218 OK -- niter = 45 fx = 0.020745360
## ... bootstrap draw number: 219 OK -- niter = 46 fx = 0.024510428
## ... bootstrap draw number: 220 OK -- niter = 45 fx = 0.020609321
## ... bootstrap draw number: 221 OK -- niter = 53 fx = 0.020661657
## ... bootstrap draw number: 222 OK -- niter = 48 fx = 0.020617437
## ... bootstrap draw number: 223 OK -- niter = 48 fx = 0.028909526
## ... bootstrap draw number: 224 OK -- niter = 48 fx = 0.020056859
## ... bootstrap draw number: 225 OK -- niter = 44 fx = 0.020180939
## ... bootstrap draw number: 226 OK -- niter = 44 fx = 0.016073913
## ... bootstrap draw number: 227 OK -- niter = 46 fx = 0.020479435
## ... bootstrap draw number: 228 OK -- niter = 51 fx = 0.020748504
## ... bootstrap draw number: 229 OK -- niter = 28 fx = 0.019443020
## ... bootstrap draw number: 230 OK -- niter = 41 fx = 0.029326588
## ... bootstrap draw number: 231 OK -- niter = 45 fx = 0.026181294
## ... bootstrap draw number: 232 OK -- niter = 51 fx = 0.012780612
## ... bootstrap draw number: 233 OK -- niter = 46 fx = 0.025434745
## ... bootstrap draw number: 234 OK -- niter = 46 fx = 0.020878297
## ... bootstrap draw number: 235 OK -- niter = 60 fx = 0.026974308
## ... bootstrap draw number: 236 OK -- niter = 47 fx = 0.024639604
## ... bootstrap draw number: 237 OK -- niter = 46 fx = 0.019060600
## ... bootstrap draw number: 238 OK -- niter = 47 fx = 0.024534943
## ... bootstrap draw number: 239 OK -- niter = 42 fx = 0.017342762
## ... bootstrap draw number: 240 OK -- niter = 50 fx = 0.024538282
## ... bootstrap draw number: 241 OK -- niter = 59 fx = 0.023557684
## ... bootstrap draw number: 242 OK -- niter = 51 fx = 0.025772589
## ... bootstrap draw number: 243 OK -- niter = 44 fx = 0.026406148
## ... bootstrap draw number: 244 OK -- niter = 49 fx = 0.027303660
## ... bootstrap draw number: 245 OK -- niter = 52 fx = 0.031539282
## ... bootstrap draw number: 246 OK -- niter = 44 fx = 0.022057794
## ... bootstrap draw number: 247 OK -- niter = 46 fx = 0.016996203
## ... bootstrap draw number: 248 OK -- niter = 54 fx = 0.016135393
## ... bootstrap draw number: 249 OK -- niter = 47 fx = 0.020432575
## ... bootstrap draw number: 250 OK -- niter = 44 fx = 0.029606427
## ... bootstrap draw number: 251 OK -- niter = 46 fx = 0.021684283
## ... bootstrap draw number: 252 OK -- niter = 47 fx = 0.020583596
## ... bootstrap draw number: 253 OK -- niter = 47 fx = 0.027908091
## ... bootstrap draw number: 254 OK -- niter = 57 fx = 0.033356116
## ... bootstrap draw number: 255 OK -- niter = 42 fx = 0.017034473
## ... bootstrap draw number: 256 OK -- niter = 43 fx = 0.021661078
## ... bootstrap draw number: 257 OK -- niter = 47 fx = 0.031725281
## ... bootstrap draw number: 258 OK -- niter = 51 fx = 0.019095862
## ... bootstrap draw number: 259 OK -- niter = 48 fx = 0.022343840
## ... bootstrap draw number: 260 OK -- niter = 42 fx = 0.018932971
## ... bootstrap draw number: 261 OK -- niter = 46 fx = 0.020383875
## ... bootstrap draw number: 262 OK -- niter = 47 fx = 0.022422956
## ... bootstrap draw number: 263 OK -- niter = 44 fx = 0.014093132
## ... bootstrap draw number: 264 OK -- niter = 42 fx = 0.020985630
## ... bootstrap draw number: 265 OK -- niter = 56 fx = 0.019000772
## ... bootstrap draw number: 266 OK -- niter = 45 fx = 0.025093333
## ... bootstrap draw number: 267 OK -- niter = 48 fx = 0.030563961
## ... bootstrap draw number: 268 OK -- niter = 53 fx = 0.019652139
## ... bootstrap draw number: 269 OK -- niter = 65 fx = 0.022870947
## ... bootstrap draw number: 270 OK -- niter = 42 fx = 0.031765057
## ... bootstrap draw number: 271 OK -- niter = 50 fx = 0.017455937
## ... bootstrap draw number: 272 OK -- niter = 43 fx = 0.018120158
## ... bootstrap draw number: 273 OK -- niter = 54 fx = 0.017019507
## ... bootstrap draw number: 274 OK -- niter = 56 fx = 0.020176608
## ... bootstrap draw number: 275 OK -- niter = 43 fx = 0.021323797
## ... bootstrap draw number: 276 OK -- niter = 44 fx = 0.019320505
## ... bootstrap draw number: 277 OK -- niter = 43 fx = 0.021033102
## ... bootstrap draw number: 278 OK -- niter = 51 fx = 0.015055382
## ... bootstrap draw number: 279 OK -- niter = 44 fx = 0.020099355
## ... bootstrap draw number: 280 OK -- niter = 45 fx = 0.022445479
## ... bootstrap draw number: 281 OK -- niter = 42 fx = 0.025229279
## ... bootstrap draw number: 282 OK -- niter = 46 fx = 0.016075439
## ... bootstrap draw number: 283 OK -- niter = 45 fx = 0.021394856
## ... bootstrap draw number: 284 OK -- niter = 43 fx = 0.028716242
## ... bootstrap draw number: 285 OK -- niter = 48 fx = 0.020754251
## ... bootstrap draw number: 286 OK -- niter = 46 fx = 0.023775074
## ... bootstrap draw number: 287 OK -- niter = 48 fx = 0.023641208
## ... bootstrap draw number: 288 OK -- niter = 43 fx = 0.022117844
## ... bootstrap draw number: 289 OK -- niter = 50 fx = 0.027553242
## ... bootstrap draw number: 290 OK -- niter = 51 fx = 0.024756997
## ... bootstrap draw number: 291 OK -- niter = 46 fx = 0.021926080
## ... bootstrap draw number: 292 OK -- niter = 50 fx = 0.018239796
## ... bootstrap draw number: 293 OK -- niter = 48 fx = 0.019545442
## ... bootstrap draw number: 294 OK -- niter = 45 fx = 0.014499653
## ... bootstrap draw number: 295 OK -- niter = 47 fx = 0.021217849
## ... bootstrap draw number: 296 OK -- niter = 43 fx = 0.022140307
## ... bootstrap draw number: 297 OK -- niter = 43 fx = 0.024112722
## ... bootstrap draw number: 298 OK -- niter = 42 fx = 0.018712728
## ... bootstrap draw number: 299 OK -- niter = 45 fx = 0.019938227
## ... bootstrap draw number: 300 OK -- niter = 42 fx = 0.016896216
## ... bootstrap draw number: 301 OK -- niter = 51 fx = 0.021904604
## ... bootstrap draw number: 302 OK -- niter = 51 fx = 0.021143666
## ... bootstrap draw number: 303 OK -- niter = 49 fx = 0.015638355
## ... bootstrap draw number: 304 OK -- niter = 48 fx = 0.024020897
## ... bootstrap draw number: 305 OK -- niter = 48 fx = 0.019242930
## ... bootstrap draw number: 306 OK -- niter = 44 fx = 0.022118980
## ... bootstrap draw number: 307 OK -- niter = 44 fx = 0.024154152
## ... bootstrap draw number: 308 OK -- niter = 45 fx = 0.025801044
## ... bootstrap draw number: 309 OK -- niter = 44 fx = 0.022085692
## ... bootstrap draw number: 310 OK -- niter = 50 fx = 0.024239708
## ... bootstrap draw number: 311 OK -- niter = 49 fx = 0.023660884
## ... bootstrap draw number: 312 OK -- niter = 49 fx = 0.027554763
## ... bootstrap draw number: 313 OK -- niter = 41 fx = 0.028365492
## ... bootstrap draw number: 314 OK -- niter = 43 fx = 0.015274548
## ... bootstrap draw number: 315 OK -- niter = 46 fx = 0.022912526
## ... bootstrap draw number: 316 OK -- niter = 52 fx = 0.017801920
## ... bootstrap draw number: 317 OK -- niter = 44 fx = 0.025179889
## ... bootstrap draw number: 318 OK -- niter = 49 fx = 0.024351225
## ... bootstrap draw number: 319 OK -- niter = 42 fx = 0.020857003
## ... bootstrap draw number: 320 OK -- niter = 53 fx = 0.025514983
## ... bootstrap draw number: 321 OK -- niter = 25 fx = 0.030319395
## ... bootstrap draw number: 322 OK -- niter = 47 fx = 0.029784522
## ... bootstrap draw number: 323 OK -- niter = 45 fx = 0.024166109
## ... bootstrap draw number: 324 OK -- niter = 47 fx = 0.020764105
## ... bootstrap draw number: 325 OK -- niter = 45 fx = 0.019945848
## ... bootstrap draw number: 326 OK -- niter = 45 fx = 0.023460953
## ... bootstrap draw number: 327 OK -- niter = 44 fx = 0.019079797
## ... bootstrap draw number: 328 OK -- niter = 47 fx = 0.021214486
## ... bootstrap draw number: 329 OK -- niter = 51 fx = 0.022122662
## ... bootstrap draw number: 330 OK -- niter = 43 fx = 0.019019288
## ... bootstrap draw number: 331 OK -- niter = 46 fx = 0.018034141
## ... bootstrap draw number: 332 OK -- niter = 45 fx = 0.017677684
## ... bootstrap draw number: 333 OK -- niter = 42 fx = 0.024812853
## ... bootstrap draw number: 334 OK -- niter = 43 fx = 0.027278964
## ... bootstrap draw number: 335 OK -- niter = 47 fx = 0.016383790
## ... bootstrap draw number: 336 OK -- niter = 48 fx = 0.024615733
## ... bootstrap draw number: 337 OK -- niter = 46 fx = 0.012750360
## ... bootstrap draw number: 338 OK -- niter = 42 fx = 0.029977596
## ... bootstrap draw number: 339 OK -- niter = 49 fx = 0.025082931
## ... bootstrap draw number: 340 OK -- niter = 47 fx = 0.018947866
## ... bootstrap draw number: 341 OK -- niter = 43 fx = 0.014763188
## ... bootstrap draw number: 342 OK -- niter = 41 fx = 0.017706887
## ... bootstrap draw number: 343 OK -- niter = 49 fx = 0.025633787
## ... bootstrap draw number: 344 OK -- niter = 45 fx = 0.018388407
## ... bootstrap draw number: 345 OK -- niter = 40 fx = 0.021972756
## ... bootstrap draw number: 346 OK -- niter = 44 fx = 0.020826615
## ... bootstrap draw number: 347 OK -- niter = 47 fx = 0.028764379
## ... bootstrap draw number: 348 OK -- niter = 59 fx = 0.030373596
## ... bootstrap draw number: 349 OK -- niter = 50 fx = 0.022876143
## ... bootstrap draw number: 350 OK -- niter = 46 fx = 0.027610841
## ... bootstrap draw number: 351 OK -- niter = 45 fx = 0.027179818
## ... bootstrap draw number: 352 OK -- niter = 43 fx = 0.035760762
## ... bootstrap draw number: 353 OK -- niter = 54 fx = 0.015732110
## ... bootstrap draw number: 354 OK -- niter = 51 fx = 0.022349496
## ... bootstrap draw number: 355 OK -- niter = 42 fx = 0.020329977
## ... bootstrap draw number: 356 OK -- niter = 43 fx = 0.015366705
## ... bootstrap draw number: 357 OK -- niter = 48 fx = 0.026834371
## ... bootstrap draw number: 358 OK -- niter = 45 fx = 0.034276915
## ... bootstrap draw number: 359 OK -- niter = 43 fx = 0.018360365
## ... bootstrap draw number: 360 OK -- niter = 45 fx = 0.022644518
## ... bootstrap draw number: 361 OK -- niter = 44 fx = 0.024782843
## ... bootstrap draw number: 362 OK -- niter = 45 fx = 0.024745266
## ... bootstrap draw number: 363 OK -- niter = 47 fx = 0.027950347
## ... bootstrap draw number: 364 OK -- niter = 45 fx = 0.028522586
## ... bootstrap draw number: 365 OK -- niter = 59 fx = 0.019679578
## ... bootstrap draw number: 366 OK -- niter = 48 fx = 0.018393228
## ... bootstrap draw number: 367 OK -- niter = 45 fx = 0.020554339
## ... bootstrap draw number: 368 OK -- niter = 51 fx = 0.023259083
## ... bootstrap draw number: 369 OK -- niter = 47 fx = 0.027113343
## ... bootstrap draw number: 370 OK -- niter = 24 fx = 0.023597672
## ... bootstrap draw number: 371 OK -- niter = 49 fx = 0.014750327
## ... bootstrap draw number: 372 OK -- niter = 44 fx = 0.024600048
## ... bootstrap draw number: 373 OK -- niter = 54 fx = 0.023555877
## ... bootstrap draw number: 374 OK -- niter = 45 fx = 0.010603490
## ... bootstrap draw number: 375 OK -- niter = 47 fx = 0.023709940
## ... bootstrap draw number: 376 OK -- niter = 45 fx = 0.018137021
## ... bootstrap draw number: 377 OK -- niter = 46 fx = 0.023878786
## ... bootstrap draw number: 378 OK -- niter = 47 fx = 0.026788148
## ... bootstrap draw number: 379 OK -- niter = 44 fx = 0.023536575
## ... bootstrap draw number: 380 OK -- niter = 47 fx = 0.025753685
## ... bootstrap draw number: 381 OK -- niter = 43 fx = 0.023571636
## ... bootstrap draw number: 382 OK -- niter = 44 fx = 0.029438156
## ... bootstrap draw number: 383 OK -- niter = 41 fx = 0.031900727
## ... bootstrap draw number: 384 OK -- niter = 49 fx = 0.018354953
## ... bootstrap draw number: 385 OK -- niter = 44 fx = 0.013686895
## ... bootstrap draw number: 386 OK -- niter = 49 fx = 0.017046074
## ... bootstrap draw number: 387 OK -- niter = 44 fx = 0.026466542
## ... bootstrap draw number: 388 OK -- niter = 41 fx = 0.020960250
## ... bootstrap draw number: 389 OK -- niter = 49 fx = 0.022127807
## ... bootstrap draw number: 390 OK -- niter = 45 fx = 0.024401761
## ... bootstrap draw number: 391 OK -- niter = 47 fx = 0.026017036
## ... bootstrap draw number: 392 OK -- niter = 51 fx = 0.025105538
## ... bootstrap draw number: 393 OK -- niter = 47 fx = 0.019356782
## ... bootstrap draw number: 394 OK -- niter = 47 fx = 0.038243940
## ... bootstrap draw number: 395 OK -- niter = 42 fx = 0.018856991
## ... bootstrap draw number: 396 OK -- niter = 46 fx = 0.022953421
## ... bootstrap draw number: 397 OK -- niter = 50 fx = 0.018991856
## ... bootstrap draw number: 398 OK -- niter = 42 fx = 0.022345206
## ... bootstrap draw number: 399 OK -- niter = 51 fx = 0.022837344
## ... bootstrap draw number: 400 OK -- niter = 50 fx = 0.022227486
## ... bootstrap draw number: 401 OK -- niter = 52 fx = 0.020776247
## ... bootstrap draw number: 402 OK -- niter = 47 fx = 0.013798111
## ... bootstrap draw number: 403 OK -- niter = 43 fx = 0.019576688
## ... bootstrap draw number: 404 OK -- niter = 45 fx = 0.030317964
## ... bootstrap draw number: 405 OK -- niter = 42 fx = 0.023086753
## ... bootstrap draw number: 406 OK -- niter = 42 fx = 0.028810797
## ... bootstrap draw number: 407 OK -- niter = 45 fx = 0.021260954
## ... bootstrap draw number: 408 OK -- niter = 42 fx = 0.026571459
## ... bootstrap draw number: 409 OK -- niter = 49 fx = 0.018530238
## ... bootstrap draw number: 410 OK -- niter = 41 fx = 0.021702562
## ... bootstrap draw number: 411 OK -- niter = 47 fx = 0.019375174
## ... bootstrap draw number: 412 OK -- niter = 44 fx = 0.022917424
## ... bootstrap draw number: 413 OK -- niter = 48 fx = 0.029297679
## ... bootstrap draw number: 414 OK -- niter = 50 fx = 0.025919253
## ... bootstrap draw number: 415 OK -- niter = 44 fx = 0.033353343
## ... bootstrap draw number: 416 OK -- niter = 50 fx = 0.023436504
## ... bootstrap draw number: 417 OK -- niter = 48 fx = 0.016090388
## ... bootstrap draw number: 418 OK -- niter = 42 fx = 0.028524964
## ... bootstrap draw number: 419 OK -- niter = 43 fx = 0.021444920
## ... bootstrap draw number: 420 OK -- niter = 44 fx = 0.028023662
## ... bootstrap draw number: 421 OK -- niter = 46 fx = 0.022768928
## ... bootstrap draw number: 422 OK -- niter = 48 fx = 0.021905520
## ... bootstrap draw number: 423 OK -- niter = 44 fx = 0.022221888
## ... bootstrap draw number: 424 OK -- niter = 45 fx = 0.024789259
## ... bootstrap draw number: 425 OK -- niter = 43 fx = 0.020255759
## ... bootstrap draw number: 426 OK -- niter = 47 fx = 0.026326920
## ... bootstrap draw number: 427 OK -- niter = 41 fx = 0.025736978
## ... bootstrap draw number: 428 OK -- niter = 45 fx = 0.018149901
## ... bootstrap draw number: 429 OK -- niter = 43 fx = 0.023290390
## ... bootstrap draw number: 430 OK -- niter = 42 fx = 0.020251239
## ... bootstrap draw number: 431 OK -- niter = 47 fx = 0.029128937
## ... bootstrap draw number: 432 OK -- niter = 49 fx = 0.017996595
## ... bootstrap draw number: 433 OK -- niter = 47 fx = 0.021669125
## ... bootstrap draw number: 434 OK -- niter = 43 fx = 0.028455968
## ... bootstrap draw number: 435 OK -- niter = 51 fx = 0.029206820
## ... bootstrap draw number: 436 OK -- niter = 47 fx = 0.018448585
## ... bootstrap draw number: 437 OK -- niter = 47 fx = 0.025227299
## ... bootstrap draw number: 438 OK -- niter = 45 fx = 0.027507822
## ... bootstrap draw number: 439 OK -- niter = 52 fx = 0.028802625
## ... bootstrap draw number: 440 OK -- niter = 44 fx = 0.018189833
## ... bootstrap draw number: 441 OK -- niter = 43 fx = 0.018967004
## ... bootstrap draw number: 442 OK -- niter = 51 fx = 0.029706838
## ... bootstrap draw number: 443 OK -- niter = 51 fx = 0.020880430
## ... bootstrap draw number: 444 OK -- niter = 49 fx = 0.026322057
## ... bootstrap draw number: 445 OK -- niter = 47 fx = 0.012659889
## ... bootstrap draw number: 446 OK -- niter = 43 fx = 0.022650879
## ... bootstrap draw number: 447 OK -- niter = 44 fx = 0.030617417
## ... bootstrap draw number: 448 OK -- niter = 48 fx = 0.023648316
## ... bootstrap draw number: 449 OK -- niter = 43 fx = 0.021083462
## ... bootstrap draw number: 450 OK -- niter = 42 fx = 0.022384318
## ... bootstrap draw number: 451 OK -- niter = 48 fx = 0.023402208
## ... bootstrap draw number: 452 OK -- niter = 47 fx = 0.019071355
## ... bootstrap draw number: 453 OK -- niter = 49 fx = 0.020339048
## ... bootstrap draw number: 454 OK -- niter = 43 fx = 0.018571777
## ... bootstrap draw number: 455 OK -- niter = 47 fx = 0.030775766
## ... bootstrap draw number: 456 OK -- niter = 53 fx = 0.019933246
## ... bootstrap draw number: 457 OK -- niter = 52 fx = 0.020830367
## ... bootstrap draw number: 458 OK -- niter = 43 fx = 0.027166922
## ... bootstrap draw number: 459 OK -- niter = 49 fx = 0.025022896
## ... bootstrap draw number: 460 OK -- niter = 44 fx = 0.018623727
## ... bootstrap draw number: 461 OK -- niter = 49 fx = 0.019628669
## ... bootstrap draw number: 462 OK -- niter = 44 fx = 0.017713838
## ... bootstrap draw number: 463 OK -- niter = 42 fx = 0.021403477
## ... bootstrap draw number: 464 OK -- niter = 43 fx = 0.022734433
## ... bootstrap draw number: 465 OK -- niter = 26 fx = 0.014026112
## ... bootstrap draw number: 466 OK -- niter = 51 fx = 0.022221247
## ... bootstrap draw number: 467 OK -- niter = 49 fx = 0.016412063
## ... bootstrap draw number: 468 OK -- niter = 43 fx = 0.025209490
## ... bootstrap draw number: 469 OK -- niter = 47 fx = 0.023693983
## ... bootstrap draw number: 470 OK -- niter = 43 fx = 0.017058374
## ... bootstrap draw number: 471 OK -- niter = 48 fx = 0.016620060
## ... bootstrap draw number: 472 OK -- niter = 49 fx = 0.019226318
## ... bootstrap draw number: 473 OK -- niter = 42 fx = 0.028490431
## ... bootstrap draw number: 474 OK -- niter = 43 fx = 0.028190122
## ... bootstrap draw number: 475 OK -- niter = 46 fx = 0.022721695
## ... bootstrap draw number: 476 OK -- niter = 46 fx = 0.026297550
## ... bootstrap draw number: 477 OK -- niter = 43 fx = 0.018947760
## ... bootstrap draw number: 478 OK -- niter = 41 fx = 0.020166949
## ... bootstrap draw number: 479 OK -- niter = 43 fx = 0.027828734
## ... bootstrap draw number: 480 OK -- niter = 43 fx = 0.023106708
## ... bootstrap draw number: 481 OK -- niter = 44 fx = 0.021402371
## ... bootstrap draw number: 482 OK -- niter = 46 fx = 0.019028147
## ... bootstrap draw number: 483 OK -- niter = 49 fx = 0.019378792
## ... bootstrap draw number: 484 OK -- niter = 45 fx = 0.029479426
## ... bootstrap draw number: 485 OK -- niter = 47 fx = 0.020730218
## ... bootstrap draw number: 486 OK -- niter = 44 fx = 0.021055110
## ... bootstrap draw number: 487 OK -- niter = 47 fx = 0.020114664
## ... bootstrap draw number: 488 OK -- niter = 51 fx = 0.018944188
## ... bootstrap draw number: 489 OK -- niter = 45 fx = 0.019485017
## ... bootstrap draw number: 490 OK -- niter = 43 fx = 0.018801675
## ... bootstrap draw number: 491 OK -- niter = 42 fx = 0.020015259
## ... bootstrap draw number: 492 OK -- niter = 42 fx = 0.019863783
## ... bootstrap draw number: 493 OK -- niter = 42 fx = 0.021241093
## ... bootstrap draw number: 494 OK -- niter = 44 fx = 0.028282592
## ... bootstrap draw number: 495 OK -- niter = 46 fx = 0.018380784
## ... bootstrap draw number: 496 OK -- niter = 47 fx = 0.019308840
## ... bootstrap draw number: 497 OK -- niter = 47 fx = 0.031565753
## ... bootstrap draw number: 498 OK -- niter = 43 fx = 0.023975558
## ... bootstrap draw number: 499 OK -- niter = 43 fx = 0.027340778
## ... bootstrap draw number: 500 OK -- niter = 40 fx = 0.024037676
## ... bootstrap draw number: 501 OK -- niter = 46 fx = 0.021673180
## ... bootstrap draw number: 502 OK -- niter = 47 fx = 0.017896417
## ... bootstrap draw number: 503 OK -- niter = 45 fx = 0.025313601
## ... bootstrap draw number: 504 OK -- niter = 45 fx = 0.020934162
## ... bootstrap draw number: 505 OK -- niter = 24 fx = 0.018720165
## ... bootstrap draw number: 506 OK -- niter = 46 fx = 0.025847772
## ... bootstrap draw number: 507 OK -- niter = 42 fx = 0.020439667
## ... bootstrap draw number: 508 OK -- niter = 45 fx = 0.026978531
## ... bootstrap draw number: 509 OK -- niter = 43 fx = 0.019452696
## ... bootstrap draw number: 510 OK -- niter = 55 fx = 0.023390826
## ... bootstrap draw number: 511 OK -- niter = 49 fx = 0.020636392
## ... bootstrap draw number: 512 OK -- niter = 41 fx = 0.019208848
## ... bootstrap draw number: 513 OK -- niter = 48 fx = 0.025304296
## ... bootstrap draw number: 514 OK -- niter = 49 fx = 0.026095891
## ... bootstrap draw number: 515 OK -- niter = 46 fx = 0.017166805
## ... bootstrap draw number: 516 OK -- niter = 41 fx = 0.023922227
## ... bootstrap draw number: 517 OK -- niter = 49 fx = 0.023208396
## ... bootstrap draw number: 518 OK -- niter = 46 fx = 0.026291700
## ... bootstrap draw number: 519 OK -- niter = 48 fx = 0.018888878
## ... bootstrap draw number: 520 OK -- niter = 44 fx = 0.019668516
## ... bootstrap draw number: 521 OK -- niter = 45 fx = 0.023162369
## ... bootstrap draw number: 522 OK -- niter = 48 fx = 0.016044658
## ... bootstrap draw number: 523 OK -- niter = 49 fx = 0.022388704
## ... bootstrap draw number: 524 OK -- niter = 47 fx = 0.026407923
## ... bootstrap draw number: 525 OK -- niter = 43 fx = 0.023104730
## ... bootstrap draw number: 526 OK -- niter = 48 fx = 0.022894009
## ... bootstrap draw number: 527 OK -- niter = 48 fx = 0.026974215
## ... bootstrap draw number: 528 OK -- niter = 48 fx = 0.021400196
## ... bootstrap draw number: 529 OK -- niter = 46 fx = 0.023049136
## ... bootstrap draw number: 530 OK -- niter = 44 fx = 0.013719780
## ... bootstrap draw number: 531 OK -- niter = 44 fx = 0.018961778
## ... bootstrap draw number: 532 OK -- niter = 47 fx = 0.021593237
## ... bootstrap draw number: 533 OK -- niter = 46 fx = 0.029452597
## ... bootstrap draw number: 534 OK -- niter = 45 fx = 0.035853318
## ... bootstrap draw number: 535 OK -- niter = 45 fx = 0.024226959
## ... bootstrap draw number: 536 OK -- niter = 48 fx = 0.027329603
## ... bootstrap draw number: 537 OK -- niter = 52 fx = 0.020142300
## ... bootstrap draw number: 538 OK -- niter = 53 fx = 0.020270980
## ... bootstrap draw number: 539 OK -- niter = 44 fx = 0.018344924
## ... bootstrap draw number: 540 OK -- niter = 50 fx = 0.025378265
## ... bootstrap draw number: 541 OK -- niter = 46 fx = 0.028900823
## ... bootstrap draw number: 542 OK -- niter = 47 fx = 0.020519619
## ... bootstrap draw number: 543 OK -- niter = 42 fx = 0.017681667
## ... bootstrap draw number: 544 OK -- niter = 45 fx = 0.026838978
## ... bootstrap draw number: 545 OK -- niter = 43 fx = 0.020624679
## ... bootstrap draw number: 546 OK -- niter = 54 fx = 0.016322920
## ... bootstrap draw number: 547 OK -- niter = 51 fx = 0.032983978
## ... bootstrap draw number: 548 OK -- niter = 42 fx = 0.032550984
## ... bootstrap draw number: 549 OK -- niter = 43 fx = 0.026139420
## ... bootstrap draw number: 550 OK -- niter = 54 fx = 0.026533877
## ... bootstrap draw number: 551 OK -- niter = 47 fx = 0.028832424
## ... bootstrap draw number: 552 OK -- niter = 43 fx = 0.029462805
## ... bootstrap draw number: 553 OK -- niter = 46 fx = 0.029105386
## ... bootstrap draw number: 554 OK -- niter = 46 fx = 0.023084588
## ... bootstrap draw number: 555 OK -- niter = 49 fx = 0.023028315
## ... bootstrap draw number: 556 OK -- niter = 50 fx = 0.021327260
## ... bootstrap draw number: 557 OK -- niter = 43 fx = 0.028815551
## ... bootstrap draw number: 558 OK -- niter = 47 fx = 0.031080187
## ... bootstrap draw number: 559 OK -- niter = 49 fx = 0.027695517
## ... bootstrap draw number: 560 OK -- niter = 49 fx = 0.020551968
## ... bootstrap draw number: 561 OK -- niter = 46 fx = 0.026464679
## ... bootstrap draw number: 562 OK -- niter = 46 fx = 0.036192667
## ... bootstrap draw number: 563 OK -- niter = 47 fx = 0.032586064
## ... bootstrap draw number: 564 OK -- niter = 42 fx = 0.019531499
## ... bootstrap draw number: 565 OK -- niter = 26 fx = 0.019374348
## ... bootstrap draw number: 566 OK -- niter = 47 fx = 0.019251744
## ... bootstrap draw number: 567 OK -- niter = 43 fx = 0.019472378
## ... bootstrap draw number: 568 OK -- niter = 46 fx = 0.015229737
## ... bootstrap draw number: 569 OK -- niter = 50 fx = 0.024668196
## ... bootstrap draw number: 570 OK -- niter = 55 fx = 0.019040948
## ... bootstrap draw number: 571 OK -- niter = 45 fx = 0.030598955
## ... bootstrap draw number: 572 OK -- niter = 42 fx = 0.027768681
## ... bootstrap draw number: 573 OK -- niter = 53 fx = 0.021600795
## ... bootstrap draw number: 574 OK -- niter = 48 fx = 0.036487747
## ... bootstrap draw number: 575 OK -- niter = 42 fx = 0.026762078
## ... bootstrap draw number: 576 OK -- niter = 60 fx = 0.029647917
## ... bootstrap draw number: 577 OK -- niter = 42 fx = 0.022451759
## ... bootstrap draw number: 578 OK -- niter = 43 fx = 0.026878974
## ... bootstrap draw number: 579 OK -- niter = 47 fx = 0.031715864
## ... bootstrap draw number: 580 OK -- niter = 44 fx = 0.016460189
## ... bootstrap draw number: 581 OK -- niter = 43 fx = 0.017695485
## ... bootstrap draw number: 582 OK -- niter = 43 fx = 0.018478166
## ... bootstrap draw number: 583 OK -- niter = 49 fx = 0.020177343
## ... bootstrap draw number: 584 OK -- niter = 49 fx = 0.023113627
## ... bootstrap draw number: 585 OK -- niter = 42 fx = 0.023077243
## ... bootstrap draw number: 586 OK -- niter = 49 fx = 0.024654945
## ... bootstrap draw number: 587 OK -- niter = 47 fx = 0.021374386
## ... bootstrap draw number: 588 OK -- niter = 43 fx = 0.020507361
## ... bootstrap draw number: 589 OK -- niter = 48 fx = 0.021102263
## ... bootstrap draw number: 590 OK -- niter = 25 fx = 0.016497242
## ... bootstrap draw number: 591 OK -- niter = 44 fx = 0.026803098
## ... bootstrap draw number: 592 OK -- niter = 48 fx = 0.024923639
## ... bootstrap draw number: 593 OK -- niter = 46 fx = 0.027432570
## ... bootstrap draw number: 594 OK -- niter = 45 fx = 0.020426533
## ... bootstrap draw number: 595 OK -- niter = 49 fx = 0.024114839
## ... bootstrap draw number: 596 OK -- niter = 54 fx = 0.038681119
## ... bootstrap draw number: 597 OK -- niter = 48 fx = 0.020725706
## ... bootstrap draw number: 598 OK -- niter = 49 fx = 0.026687846
## ... bootstrap draw number: 599 OK -- niter = 46 fx = 0.023034747
## ... bootstrap draw number: 600 OK -- niter = 43 fx = 0.021932770
## ... bootstrap draw number: 601 OK -- niter = 50 fx = 0.021530049
## ... bootstrap draw number: 602 OK -- niter = 45 fx = 0.020327939
## ... bootstrap draw number: 603 OK -- niter = 43 fx = 0.021524182
## ... bootstrap draw number: 604 OK -- niter = 50 fx = 0.019392020
## ... bootstrap draw number: 605 OK -- niter = 45 fx = 0.019321240
## ... bootstrap draw number: 606 OK -- niter = 49 fx = 0.029857882
## ... bootstrap draw number: 607 OK -- niter = 44 fx = 0.019477239
## ... bootstrap draw number: 608 OK -- niter = 47 fx = 0.023841065
## ... bootstrap draw number: 609 OK -- niter = 46 fx = 0.022534044
## ... bootstrap draw number: 610 OK -- niter = 47 fx = 0.029470490
## ... bootstrap draw number: 611 OK -- niter = 43 fx = 0.015606744
## ... bootstrap draw number: 612 OK -- niter = 46 fx = 0.009161857
## ... bootstrap draw number: 613 OK -- niter = 47 fx = 0.018893479
## ... bootstrap draw number: 614 OK -- niter = 42 fx = 0.024417720
## ... bootstrap draw number: 615 OK -- niter = 42 fx = 0.022367449
## ... bootstrap draw number: 616 OK -- niter = 47 fx = 0.020533226
## ... bootstrap draw number: 617 OK -- niter = 48 fx = 0.023740487
## ... bootstrap draw number: 618 OK -- niter = 44 fx = 0.024687543
## ... bootstrap draw number: 619 OK -- niter = 46 fx = 0.024566207
## ... bootstrap draw number: 620 OK -- niter = 45 fx = 0.023806696
## ... bootstrap draw number: 621 OK -- niter = 42 fx = 0.022047295
## ... bootstrap draw number: 622 OK -- niter = 50 fx = 0.023443089
## ... bootstrap draw number: 623 OK -- niter = 49 fx = 0.026603272
## ... bootstrap draw number: 624 OK -- niter = 47 fx = 0.033841299
## ... bootstrap draw number: 625 OK -- niter = 47 fx = 0.029778410
## ... bootstrap draw number: 626 OK -- niter = 45 fx = 0.024625020
## ... bootstrap draw number: 627 OK -- niter = 51 fx = 0.020721423
## ... bootstrap draw number: 628 OK -- niter = 46 fx = 0.024429662
## ... bootstrap draw number: 629 OK -- niter = 48 fx = 0.032851495
## ... bootstrap draw number: 630 OK -- niter = 45 fx = 0.034523129
## ... bootstrap draw number: 631 OK -- niter = 46 fx = 0.018043318
## ... bootstrap draw number: 632 OK -- niter = 46 fx = 0.025396401
## ... bootstrap draw number: 633 OK -- niter = 45 fx = 0.021299276
## ... bootstrap draw number: 634 OK -- niter = 44 fx = 0.026678754
## ... bootstrap draw number: 635 OK -- niter = 47 fx = 0.015474323
## ... bootstrap draw number: 636 OK -- niter = 47 fx = 0.026421281
## ... bootstrap draw number: 637 OK -- niter = 41 fx = 0.024504763
## ... bootstrap draw number: 638 OK -- niter = 47 fx = 0.033144640
## ... bootstrap draw number: 639 OK -- niter = 58 fx = 0.020490647
## ... bootstrap draw number: 640 OK -- niter = 43 fx = 0.022988378
## ... bootstrap draw number: 641 OK -- niter = 49 fx = 0.014496710
## ... bootstrap draw number: 642 OK -- niter = 48 fx = 0.020088387
## ... bootstrap draw number: 643 OK -- niter = 45 fx = 0.017043504
## ... bootstrap draw number: 644 OK -- niter = 45 fx = 0.023219343
## ... bootstrap draw number: 645 OK -- niter = 44 fx = 0.018327179
## ... bootstrap draw number: 646 OK -- niter = 46 fx = 0.026189321
## ... bootstrap draw number: 647 OK -- niter = 43 fx = 0.018623208
## ... bootstrap draw number: 648 OK -- niter = 43 fx = 0.026765894
## ... bootstrap draw number: 649 OK -- niter = 43 fx = 0.014595338
## ... bootstrap draw number: 650 OK -- niter = 47 fx = 0.026702377
## ... bootstrap draw number: 651 OK -- niter = 41 fx = 0.019867040
## ... bootstrap draw number: 652 OK -- niter = 47 fx = 0.031108242
## ... bootstrap draw number: 653 OK -- niter = 49 fx = 0.025517579
## ... bootstrap draw number: 654 OK -- niter = 45 fx = 0.024935353
## ... bootstrap draw number: 655 OK -- niter = 45 fx = 0.024127369
## ... bootstrap draw number: 656 OK -- niter = 52 fx = 0.028933293
## ... bootstrap draw number: 657 OK -- niter = 42 fx = 0.024776136
## ... bootstrap draw number: 658 OK -- niter = 47 fx = 0.018169977
## ... bootstrap draw number: 659 OK -- niter = 43 fx = 0.023641902
## ... bootstrap draw number: 660 OK -- niter = 57 fx = 0.034901543
## ... bootstrap draw number: 661 OK -- niter = 45 fx = 0.023239081
## ... bootstrap draw number: 662 OK -- niter = 49 fx = 0.021655104
## ... bootstrap draw number: 663 OK -- niter = 45 fx = 0.019189277
## ... bootstrap draw number: 664 OK -- niter = 49 fx = 0.016836400
## ... bootstrap draw number: 665 OK -- niter = 43 fx = 0.019286989
## ... bootstrap draw number: 666 OK -- niter = 54 fx = 0.031835907
## ... bootstrap draw number: 667 OK -- niter = 43 fx = 0.031570841
## ... bootstrap draw number: 668 OK -- niter = 48 fx = 0.014579678
## ... bootstrap draw number: 669 OK -- niter = 24 fx = 0.018949937
## ... bootstrap draw number: 670 OK -- niter = 48 fx = 0.026014973
## ... bootstrap draw number: 671 OK -- niter = 44 fx = 0.022366708
## ... bootstrap draw number: 672 OK -- niter = 51 fx = 0.027045115
## ... bootstrap draw number: 673 OK -- niter = 46 fx = 0.021109824
## ... bootstrap draw number: 674 OK -- niter = 49 fx = 0.015193904
## ... bootstrap draw number: 675 OK -- niter = 48 fx = 0.028622932
## ... bootstrap draw number: 676 OK -- niter = 26 fx = 0.028342197
## ... bootstrap draw number: 677 OK -- niter = 25 fx = 0.022999498
## ... bootstrap draw number: 678 OK -- niter = 49 fx = 0.021788864
## ... bootstrap draw number: 679 OK -- niter = 45 fx = 0.027730885
## ... bootstrap draw number: 680 OK -- niter = 45 fx = 0.025429320
## ... bootstrap draw number: 681 OK -- niter = 45 fx = 0.022516471
## ... bootstrap draw number: 682 OK -- niter = 43 fx = 0.027535737
## ... bootstrap draw number: 683 OK -- niter = 50 fx = 0.022002246
## ... bootstrap draw number: 684 OK -- niter = 46 fx = 0.021805769
## ... bootstrap draw number: 685 OK -- niter = 44 fx = 0.028691320
## ... bootstrap draw number: 686 OK -- niter = 45 fx = 0.024469234
## ... bootstrap draw number: 687 OK -- niter = 46 fx = 0.027766880
## ... bootstrap draw number: 688 OK -- niter = 46 fx = 0.021475500
## ... bootstrap draw number: 689 OK -- niter = 46 fx = 0.023427175
## ... bootstrap draw number: 690 OK -- niter = 51 fx = 0.018500747
## ... bootstrap draw number: 691 OK -- niter = 46 fx = 0.023715448
## ... bootstrap draw number: 692 OK -- niter = 44 fx = 0.021416302
## ... bootstrap draw number: 693 OK -- niter = 42 fx = 0.019280995
## ... bootstrap draw number: 694 OK -- niter = 42 fx = 0.030695446
## ... bootstrap draw number: 695 OK -- niter = 42 fx = 0.021960891
## ... bootstrap draw number: 696 OK -- niter = 49 fx = 0.017284015
## ... bootstrap draw number: 697 OK -- niter = 56 fx = 0.019603784
## ... bootstrap draw number: 698 OK -- niter = 46 fx = 0.023228897
## ... bootstrap draw number: 699 OK -- niter = 49 fx = 0.015947847
## ... bootstrap draw number: 700 OK -- niter = 47 fx = 0.025154669
## ... bootstrap draw number: 701 OK -- niter = 53 fx = 0.027609127
## ... bootstrap draw number: 702 OK -- niter = 52 fx = 0.025115796
## ... bootstrap draw number: 703 OK -- niter = 41 fx = 0.021978053
## ... bootstrap draw number: 704 OK -- niter = 48 fx = 0.026817498
## ... bootstrap draw number: 705 OK -- niter = 50 fx = 0.024808076
## ... bootstrap draw number: 706 OK -- niter = 43 fx = 0.017324939
## ... bootstrap draw number: 707 OK -- niter = 44 fx = 0.017492474
## ... bootstrap draw number: 708 OK -- niter = 42 fx = 0.022406567
## ... bootstrap draw number: 709 OK -- niter = 51 fx = 0.017331505
## ... bootstrap draw number: 710 OK -- niter = 43 fx = 0.020247440
## ... bootstrap draw number: 711 OK -- niter = 49 fx = 0.023272706
## ... bootstrap draw number: 712 OK -- niter = 45 fx = 0.035401251
## ... bootstrap draw number: 713 OK -- niter = 44 fx = 0.017626523
## ... bootstrap draw number: 714 OK -- niter = 43 fx = 0.017493254
## ... bootstrap draw number: 715 OK -- niter = 50 fx = 0.024691064
## ... bootstrap draw number: 716 OK -- niter = 45 fx = 0.024246129
## ... bootstrap draw number: 717 OK -- niter = 53 fx = 0.021056680
## ... bootstrap draw number: 718 OK -- niter = 42 fx = 0.014013086
## ... bootstrap draw number: 719 OK -- niter = 49 fx = 0.022465680
## ... bootstrap draw number: 720 OK -- niter = 47 fx = 0.024703437
## ... bootstrap draw number: 721 OK -- niter = 43 fx = 0.025383403
## ... bootstrap draw number: 722 OK -- niter = 58 fx = 0.017231563
## ... bootstrap draw number: 723 OK -- niter = 52 fx = 0.023194455
## ... bootstrap draw number: 724 OK -- niter = 46 fx = 0.020783321
## ... bootstrap draw number: 725 OK -- niter = 43 fx = 0.025993745
## ... bootstrap draw number: 726 OK -- niter = 44 fx = 0.021292070
## ... bootstrap draw number: 727 OK -- niter = 46 fx = 0.023885750
## ... bootstrap draw number: 728 OK -- niter = 55 fx = 0.022054084
## ... bootstrap draw number: 729 OK -- niter = 50 fx = 0.027463600
## ... bootstrap draw number: 730 OK -- niter = 48 fx = 0.023833457
## ... bootstrap draw number: 731 OK -- niter = 43 fx = 0.023459048
## ... bootstrap draw number: 732 OK -- niter = 51 fx = 0.023608346
## ... bootstrap draw number: 733 OK -- niter = 48 fx = 0.020452460
## ... bootstrap draw number: 734 OK -- niter = 47 fx = 0.020533034
## ... bootstrap draw number: 735 OK -- niter = 60 fx = 0.032582603
## ... bootstrap draw number: 736 OK -- niter = 47 fx = 0.021364609
## ... bootstrap draw number: 737 OK -- niter = 47 fx = 0.016463805
## ... bootstrap draw number: 738 OK -- niter = 46 fx = 0.031069227
## ... bootstrap draw number: 739 OK -- niter = 48 fx = 0.019496397
## ... bootstrap draw number: 740 OK -- niter = 42 fx = 0.020306223
## ... bootstrap draw number: 741 OK -- niter = 79 fx = 0.028999472
## ... bootstrap draw number: 742 OK -- niter = 48 fx = 0.022890219
## ... bootstrap draw number: 743 OK -- niter = 45 fx = 0.015805772
## ... bootstrap draw number: 744 OK -- niter = 48 fx = 0.019354747
## ... bootstrap draw number: 745 OK -- niter = 27 fx = 0.017312067
## ... bootstrap draw number: 746 OK -- niter = 44 fx = 0.017350708
## ... bootstrap draw number: 747 OK -- niter = 44 fx = 0.031867947
## ... bootstrap draw number: 748 OK -- niter = 48 fx = 0.024947325
## ... bootstrap draw number: 749 OK -- niter = 25 fx = 0.014905812
## ... bootstrap draw number: 750 OK -- niter = 47 fx = 0.017826283
## ... bootstrap draw number: 751 OK -- niter = 51 fx = 0.019474181
## ... bootstrap draw number: 752 OK -- niter = 51 fx = 0.022226718
## ... bootstrap draw number: 753 OK -- niter = 46 fx = 0.018050553
## ... bootstrap draw number: 754 OK -- niter = 48 fx = 0.031591681
## ... bootstrap draw number: 755 OK -- niter = 47 fx = 0.018633784
## ... bootstrap draw number: 756 OK -- niter = 45 fx = 0.017325971
## ... bootstrap draw number: 757 OK -- niter = 47 fx = 0.023650575
## ... bootstrap draw number: 758 OK -- niter = 45 fx = 0.019003829
## ... bootstrap draw number: 759 OK -- niter = 43 fx = 0.021923336
## ... bootstrap draw number: 760 OK -- niter = 45 fx = 0.022839694
## ... bootstrap draw number: 761 OK -- niter = 44 fx = 0.030276024
## ... bootstrap draw number: 762 OK -- niter = 53 fx = 0.021003295
## ... bootstrap draw number: 763 OK -- niter = 43 fx = 0.018673681
## ... bootstrap draw number: 764 OK -- niter = 52 fx = 0.017599180
## ... bootstrap draw number: 765 OK -- niter = 46 fx = 0.019940558
## ... bootstrap draw number: 766 OK -- niter = 52 fx = 0.026551261
## ... bootstrap draw number: 767 OK -- niter = 47 fx = 0.015645038
## ... bootstrap draw number: 768 OK -- niter = 39 fx = 0.012785124
## ... bootstrap draw number: 769 OK -- niter = 48 fx = 0.030594890
## ... bootstrap draw number: 770 OK -- niter = 47 fx = 0.033776248
## ... bootstrap draw number: 771 OK -- niter = 44 fx = 0.025093395
## ... bootstrap draw number: 772 OK -- niter = 43 fx = 0.017705945
## ... bootstrap draw number: 773 OK -- niter = 48 fx = 0.028757915
## ... bootstrap draw number: 774 OK -- niter = 42 fx = 0.021689454
## ... bootstrap draw number: 775 OK -- niter = 46 fx = 0.022919980
## ... bootstrap draw number: 776 OK -- niter = 24 fx = 0.020367214
## ... bootstrap draw number: 777 OK -- niter = 42 fx = 0.018410228
## ... bootstrap draw number: 778 OK -- niter = 50 fx = 0.025229904
## ... bootstrap draw number: 779 OK -- niter = 48 fx = 0.017798835
## ... bootstrap draw number: 780 OK -- niter = 44 fx = 0.027341077
## ... bootstrap draw number: 781 OK -- niter = 46 fx = 0.019629100
## ... bootstrap draw number: 782 OK -- niter = 45 fx = 0.014348748
## ... bootstrap draw number: 783 OK -- niter = 49 fx = 0.030550594
## ... bootstrap draw number: 784 OK -- niter = 43 fx = 0.022642955
## ... bootstrap draw number: 785 OK -- niter = 45 fx = 0.031822269
## ... bootstrap draw number: 786 OK -- niter = 44 fx = 0.020283612
## ... bootstrap draw number: 787 OK -- niter = 40 fx = 0.018469998
## ... bootstrap draw number: 788 OK -- niter = 52 fx = 0.019947465
## ... bootstrap draw number: 789 OK -- niter = 45 fx = 0.025918081
## ... bootstrap draw number: 790 OK -- niter = 49 fx = 0.026817604
## ... bootstrap draw number: 791 OK -- niter = 50 fx = 0.014209226
## ... bootstrap draw number: 792 OK -- niter = 46 fx = 0.023892914
## ... bootstrap draw number: 793 OK -- niter = 44 fx = 0.026619491
## ... bootstrap draw number: 794 OK -- niter = 44 fx = 0.021785947
## ... bootstrap draw number: 795 OK -- niter = 42 fx = 0.025340718
## ... bootstrap draw number: 796 OK -- niter = 46 fx = 0.024863850
## ... bootstrap draw number: 797 OK -- niter = 48 fx = 0.029994986
## ... bootstrap draw number: 798 OK -- niter = 48 fx = 0.032459049
## ... bootstrap draw number: 799 OK -- niter = 44 fx = 0.018549450
## ... bootstrap draw number: 800 OK -- niter = 49 fx = 0.021324089
## ... bootstrap draw number: 801 OK -- niter = 43 fx = 0.024187583
## ... bootstrap draw number: 802 OK -- niter = 48 fx = 0.019643225
## ... bootstrap draw number: 803 OK -- niter = 46 fx = 0.025168559
## ... bootstrap draw number: 804 OK -- niter = 49 fx = 0.020765404
## ... bootstrap draw number: 805 OK -- niter = 46 fx = 0.014786932
## ... bootstrap draw number: 806 OK -- niter = 46 fx = 0.010468514
## ... bootstrap draw number: 807 OK -- niter = 45 fx = 0.020641484
## ... bootstrap draw number: 808 OK -- niter = 48 fx = 0.023049346
## ... bootstrap draw number: 809 OK -- niter = 50 fx = 0.026458254
## ... bootstrap draw number: 810 OK -- niter = 51 fx = 0.025622765
## ... bootstrap draw number: 811 OK -- niter = 45 fx = 0.021556901
## ... bootstrap draw number: 812 OK -- niter = 54 fx = 0.023764893
## ... bootstrap draw number: 813 OK -- niter = 53 fx = 0.028313930
## ... bootstrap draw number: 814 OK -- niter = 48 fx = 0.023068685
## ... bootstrap draw number: 815 OK -- niter = 42 fx = 0.025719112
## ... bootstrap draw number: 816 OK -- niter = 51 fx = 0.023351299
## ... bootstrap draw number: 817 OK -- niter = 45 fx = 0.015979773
## ... bootstrap draw number: 818 OK -- niter = 49 fx = 0.028266021
## ... bootstrap draw number: 819 OK -- niter = 44 fx = 0.029248847
## ... bootstrap draw number: 820 OK -- niter = 42 fx = 0.029267005
## ... bootstrap draw number: 821 OK -- niter = 46 fx = 0.020842062
## ... bootstrap draw number: 822 OK -- niter = 44 fx = 0.016824363
## ... bootstrap draw number: 823 OK -- niter = 44 fx = 0.028335012
## ... bootstrap draw number: 824 OK -- niter = 43 fx = 0.026845700
## ... bootstrap draw number: 825 OK -- niter = 44 fx = 0.027891114
## ... bootstrap draw number: 826 OK -- niter = 49 fx = 0.025569251
## ... bootstrap draw number: 827 OK -- niter = 48 fx = 0.025731722
## ... bootstrap draw number: 828 OK -- niter = 48 fx = 0.022745301
## ... bootstrap draw number: 829 OK -- niter = 59 fx = 0.029920427
## ... bootstrap draw number: 830 OK -- niter = 72 fx = 0.020967712
## ... bootstrap draw number: 831 OK -- niter = 52 fx = 0.027014802
## ... bootstrap draw number: 832 OK -- niter = 50 fx = 0.022150353
## ... bootstrap draw number: 833 OK -- niter = 46 fx = 0.017326506
## ... bootstrap draw number: 834 OK -- niter = 52 fx = 0.024747110
## ... bootstrap draw number: 835 OK -- niter = 43 fx = 0.022781640
## ... bootstrap draw number: 836 OK -- niter = 47 fx = 0.018773210
## ... bootstrap draw number: 837 OK -- niter = 46 fx = 0.020621972
## ... bootstrap draw number: 838 OK -- niter = 49 fx = 0.017218724
## ... bootstrap draw number: 839 OK -- niter = 45 fx = 0.026835937
## ... bootstrap draw number: 840 OK -- niter = 47 fx = 0.014822838
## ... bootstrap draw number: 841 OK -- niter = 46 fx = 0.018693712
## ... bootstrap draw number: 842 OK -- niter = 43 fx = 0.018483526
## ... bootstrap draw number: 843 OK -- niter = 46 fx = 0.030389317
## ... bootstrap draw number: 844 OK -- niter = 42 fx = 0.022448585
## ... bootstrap draw number: 845 OK -- niter = 43 fx = 0.020384353
## ... bootstrap draw number: 846 OK -- niter = 46 fx = 0.021109983
## ... bootstrap draw number: 847 OK -- niter = 48 fx = 0.017327491
## ... bootstrap draw number: 848 OK -- niter = 48 fx = 0.019176974
## ... bootstrap draw number: 849 OK -- niter = 51 fx = 0.032376002
## ... bootstrap draw number: 850 OK -- niter = 51 fx = 0.016796147
## ... bootstrap draw number: 851 OK -- niter = 42 fx = 0.013643078
## ... bootstrap draw number: 852 OK -- niter = 56 fx = 0.025025533
## ... bootstrap draw number: 853 OK -- niter = 51 fx = 0.019104976
## ... bootstrap draw number: 854 OK -- niter = 43 fx = 0.024832757
## ... bootstrap draw number: 855 OK -- niter = 44 fx = 0.015304258
## ... bootstrap draw number: 856 OK -- niter = 43 fx = 0.027357611
## ... bootstrap draw number: 857 OK -- niter = 43 fx = 0.026310603
## ... bootstrap draw number: 858 OK -- niter = 55 fx = 0.019101042
## ... bootstrap draw number: 859 OK -- niter = 51 fx = 0.025398326
## ... bootstrap draw number: 860 OK -- niter = 50 fx = 0.018020815
## ... bootstrap draw number: 861 OK -- niter = 44 fx = 0.021089010
## ... bootstrap draw number: 862 OK -- niter = 45 fx = 0.024401618
## ... bootstrap draw number: 863 OK -- niter = 44 fx = 0.024594704
## ... bootstrap draw number: 864 OK -- niter = 48 fx = 0.020921016
## ... bootstrap draw number: 865 OK -- niter = 46 fx = 0.028370712
## ... bootstrap draw number: 866 OK -- niter = 46 fx = 0.025968666
## ... bootstrap draw number: 867 OK -- niter = 44 fx = 0.024381187
## ... bootstrap draw number: 868 OK -- niter = 46 fx = 0.013923800
## ... bootstrap draw number: 869 OK -- niter = 44 fx = 0.017832182
## ... bootstrap draw number: 870 OK -- niter = 43 fx = 0.012732329
## ... bootstrap draw number: 871 OK -- niter = 43 fx = 0.020661606
## ... bootstrap draw number: 872 OK -- niter = 46 fx = 0.020999931
## ... bootstrap draw number: 873 OK -- niter = 44 fx = 0.019286663
## ... bootstrap draw number: 874 OK -- niter = 43 fx = 0.026694606
## ... bootstrap draw number: 875 OK -- niter = 44 fx = 0.026721686
## ... bootstrap draw number: 876 OK -- niter = 47 fx = 0.025211295
## ... bootstrap draw number: 877 OK -- niter = 44 fx = 0.025373213
## ... bootstrap draw number: 878 OK -- niter = 44 fx = 0.039911022
## ... bootstrap draw number: 879 OK -- niter = 55 fx = 0.023322597
## ... bootstrap draw number: 880 OK -- niter = 44 fx = 0.019298174
## ... bootstrap draw number: 881 OK -- niter = 55 fx = 0.016829822
## ... bootstrap draw number: 882 OK -- niter = 45 fx = 0.024768558
## ... bootstrap draw number: 883 OK -- niter = 51 fx = 0.022958922
## ... bootstrap draw number: 884 OK -- niter = 47 fx = 0.021248976
## ... bootstrap draw number: 885 OK -- niter = 50 fx = 0.026565693
## ... bootstrap draw number: 886 OK -- niter = 55 fx = 0.027243404
## ... bootstrap draw number: 887 OK -- niter = 45 fx = 0.014565514
## ... bootstrap draw number: 888 OK -- niter = 47 fx = 0.029320049
## ... bootstrap draw number: 889 OK -- niter = 46 fx = 0.020588835
## ... bootstrap draw number: 890 OK -- niter = 47 fx = 0.022794794
## ... bootstrap draw number: 891 OK -- niter = 46 fx = 0.026452775
## ... bootstrap draw number: 892 OK -- niter = 55 fx = 0.024822249
## ... bootstrap draw number: 893 OK -- niter = 53 fx = 0.029665482
## ... bootstrap draw number: 894 OK -- niter = 47 fx = 0.026671061
## ... bootstrap draw number: 895 OK -- niter = 48 fx = 0.011195522
## ... bootstrap draw number: 896 OK -- niter = 44 fx = 0.026777680
## ... bootstrap draw number: 897 OK -- niter = 42 fx = 0.025439168
## ... bootstrap draw number: 898 OK -- niter = 45 fx = 0.021855233
## ... bootstrap draw number: 899 OK -- niter = 47 fx = 0.027499366
## ... bootstrap draw number: 900 OK -- niter = 42 fx = 0.030329078
## ... bootstrap draw number: 901 OK -- niter = 49 fx = 0.015795437
## ... bootstrap draw number: 902 OK -- niter = 53 fx = 0.025197776
## ... bootstrap draw number: 903 OK -- niter = 44 fx = 0.028251798
## ... bootstrap draw number: 904 OK -- niter = 50 fx = 0.019027738
## ... bootstrap draw number: 905 OK -- niter = 41 fx = 0.019589166
## ... bootstrap draw number: 906 OK -- niter = 50 fx = 0.027029705
## ... bootstrap draw number: 907 OK -- niter = 52 fx = 0.016943924
## ... bootstrap draw number: 908 OK -- niter = 42 fx = 0.019245185
## ... bootstrap draw number: 909 OK -- niter = 43 fx = 0.021250387
## ... bootstrap draw number: 910 OK -- niter = 41 fx = 0.017583888
## ... bootstrap draw number: 911 OK -- niter = 57 fx = 0.017317714
## ... bootstrap draw number: 912 OK -- niter = 47 fx = 0.028776050
## ... bootstrap draw number: 913 OK -- niter = 44 fx = 0.012236675
## ... bootstrap draw number: 914 OK -- niter = 44 fx = 0.018277825
## ... bootstrap draw number: 915 OK -- niter = 45 fx = 0.018414232
## ... bootstrap draw number: 916 OK -- niter = 67 fx = 0.025386236
## ... bootstrap draw number: 917 OK -- niter = 53 fx = 0.026526334
## ... bootstrap draw number: 918 OK -- niter = 41 fx = 0.030030550
## ... bootstrap draw number: 919 OK -- niter = 42 fx = 0.023154014
## ... bootstrap draw number: 920 OK -- niter = 43 fx = 0.020139166
## ... bootstrap draw number: 921 OK -- niter = 43 fx = 0.028325375
## ... bootstrap draw number: 922 OK -- niter = 54 fx = 0.024754559
## ... bootstrap draw number: 923 OK -- niter = 54 fx = 0.020194499
## ... bootstrap draw number: 924 OK -- niter = 45 fx = 0.021933642
## ... bootstrap draw number: 925 OK -- niter = 44 fx = 0.022552356
## ... bootstrap draw number: 926 OK -- niter = 43 fx = 0.028492935
## ... bootstrap draw number: 927 OK -- niter = 46 fx = 0.027884761
## ... bootstrap draw number: 928 OK -- niter = 46 fx = 0.032347635
## ... bootstrap draw number: 929 OK -- niter = 49 fx = 0.028936197
## ... bootstrap draw number: 930 OK -- niter = 44 fx = 0.016713929
## ... bootstrap draw number: 931 OK -- niter = 52 fx = 0.023251784
## ... bootstrap draw number: 932 OK -- niter = 50 fx = 0.025157629
## ... bootstrap draw number: 933 OK -- niter = 27 fx = 0.023517232
## ... bootstrap draw number: 934 OK -- niter = 43 fx = 0.021016921
## ... bootstrap draw number: 935 OK -- niter = 45 fx = 0.027317168
## ... bootstrap draw number: 936 OK -- niter = 44 fx = 0.016917742
## ... bootstrap draw number: 937 OK -- niter = 47 fx = 0.018600358
## ... bootstrap draw number: 938 OK -- niter = 45 fx = 0.027398933
## ... bootstrap draw number: 939 OK -- niter = 42 fx = 0.019882643
## ... bootstrap draw number: 940 OK -- niter = 45 fx = 0.022700338
## ... bootstrap draw number: 941 OK -- niter = 43 fx = 0.025995514
## ... bootstrap draw number: 942 OK -- niter = 44 fx = 0.024117838
## ... bootstrap draw number: 943 OK -- niter = 44 fx = 0.025650776
## ... bootstrap draw number: 944 OK -- niter = 44 fx = 0.022382917
## ... bootstrap draw number: 945 OK -- niter = 43 fx = 0.015421655
## ... bootstrap draw number: 946 OK -- niter = 47 fx = 0.021797376
## ... bootstrap draw number: 947 OK -- niter = 49 fx = 0.025921381
## ... bootstrap draw number: 948 OK -- niter = 44 fx = 0.020657052
## ... bootstrap draw number: 949 OK -- niter = 47 fx = 0.024346856
## ... bootstrap draw number: 950 OK -- niter = 42 fx = 0.024694735
## ... bootstrap draw number: 951 OK -- niter = 44 fx = 0.025322615
## ... bootstrap draw number: 952 OK -- niter = 44 fx = 0.016429603
## ... bootstrap draw number: 953 OK -- niter = 45 fx = 0.025944306
## ... bootstrap draw number: 954 OK -- niter = 47 fx = 0.026413586
## ... bootstrap draw number: 955 OK -- niter = 52 fx = 0.016273678
## ... bootstrap draw number: 956 OK -- niter = 49 fx = 0.020750222
## ... bootstrap draw number: 957 OK -- niter = 46 fx = 0.027741326
## ... bootstrap draw number: 958 OK -- niter = 53 fx = 0.017584960
## ... bootstrap draw number: 959 OK -- niter = 42 fx = 0.031840406
## ... bootstrap draw number: 960 OK -- niter = 48 fx = 0.018568451
## ... bootstrap draw number: 961 OK -- niter = 47 fx = 0.023273185
## ... bootstrap draw number: 962 OK -- niter = 45 fx = 0.025014357
## ... bootstrap draw number: 963 OK -- niter = 42 fx = 0.013502737
## ... bootstrap draw number: 964 OK -- niter = 54 fx = 0.018791448
## ... bootstrap draw number: 965 OK -- niter = 44 fx = 0.018208004
## ... bootstrap draw number: 966 OK -- niter = 47 fx = 0.018084117
## ... bootstrap draw number: 967 OK -- niter = 46 fx = 0.015603945
## ... bootstrap draw number: 968 OK -- niter = 46 fx = 0.028443308
## ... bootstrap draw number: 969 OK -- niter = 48 fx = 0.027766728
## ... bootstrap draw number: 970 OK -- niter = 44 fx = 0.017903094
## ... bootstrap draw number: 971 OK -- niter = 42 fx = 0.027404737
## ... bootstrap draw number: 972 OK -- niter = 44 fx = 0.022987483
## ... bootstrap draw number: 973 OK -- niter = 49 fx = 0.021096962
## ... bootstrap draw number: 974 OK -- niter = 45 fx = 0.017614530
## ... bootstrap draw number: 975 OK -- niter = 47 fx = 0.023142346
## ... bootstrap draw number: 976 OK -- niter = 51 fx = 0.030840103
## ... bootstrap draw number: 977 OK -- niter = 42 fx = 0.016417971
## ... bootstrap draw number: 978 OK -- niter = 49 fx = 0.023166072
## ... bootstrap draw number: 979 OK -- niter = 47 fx = 0.018849862
## ... bootstrap draw number: 980 OK -- niter = 47 fx = 0.020233867
## ... bootstrap draw number: 981 OK -- niter = 24 fx = 0.025956899
## ... bootstrap draw number: 982 OK -- niter = 53 fx = 0.027068788
## ... bootstrap draw number: 983 OK -- niter = 46 fx = 0.023211580
## ... bootstrap draw number: 984 OK -- niter = 43 fx = 0.019761567
## ... bootstrap draw number: 985 OK -- niter = 49 fx = 0.025568093
## ... bootstrap draw number: 986 OK -- niter = 47 fx = 0.019560174
## ... bootstrap draw number: 987 OK -- niter = 45 fx = 0.024789723
## ... bootstrap draw number: 988 OK -- niter = 41 fx = 0.035223256
## ... bootstrap draw number: 989 OK -- niter = 52 fx = 0.023340883
## ... bootstrap draw number: 990 OK -- niter = 41 fx = 0.019146866
## ... bootstrap draw number: 991 OK -- niter = 51 fx = 0.024284213
## ... bootstrap draw number: 992 OK -- niter = 44 fx = 0.029974395
## ... bootstrap draw number: 993 OK -- niter = 48 fx = 0.028438065
## ... bootstrap draw number: 994 OK -- niter = 48 fx = 0.021096509
## ... bootstrap draw number: 995 OK -- niter = 43 fx = 0.021194098
## ... bootstrap draw number: 996 OK -- niter = 47 fx = 0.025295991
## ... bootstrap draw number: 997 OK -- niter = 45 fx = 0.024140955
## ... bootstrap draw number: 998 OK -- niter = 75 fx = 0.037265406
## ... bootstrap draw number: 999 OK -- niter = 51 fx = 0.020652497
## ... bootstrap draw number: 1000 OK -- niter = 43 fx = 0.029454065
## Number of successful bootstrap draws: 1000
## done.
## Computing TEST for test = standard ... done.## Warning in lav_object_post_check(object): lavaan WARNING: some estimated ov
## variances are negative#fitmeasures
fitMeasures(pc2016.fit2)## npar fmin
## 15.000 0.022
## chisq df
## 105.776 5.000
## pvalue chisq.scaled
## 0.000 101.520
## df.scaled pvalue.scaled
## 5.000 0.000
## chisq.scaling.factor baseline.chisq
## 1.042 1692.805
## baseline.df baseline.pvalue
## 10.000 0.000
## baseline.chisq.scaled baseline.df.scaled
## 1115.088 10.000
## baseline.pvalue.scaled baseline.chisq.scaling.factor
## 0.000 1.518
## cfi tli
## 0.940 0.880
## nnfi rfi
## 0.880 0.875
## nfi pnfi
## 0.938 0.469
## ifi rni
## 0.940 0.940
## cfi.scaled tli.scaled
## 0.913 0.825
## cfi.robust tli.robust
## 0.940 0.880
## nnfi.scaled nnfi.robust
## 0.825 0.880
## rfi.scaled nfi.scaled
## 0.818 0.909
## ifi.scaled rni.scaled
## 0.913 0.913
## rni.robust logl
## 0.940 -53841.964
## unrestricted.logl aic
## -53789.076 107713.928
## bic ntotal
## 107800.887 2434.000
## bic2 scaling.factor.h1
## 107753.229 1.561
## scaling.factor.h0 rmsea
## 1.734 0.091
## rmsea.ci.lower rmsea.ci.upper
## 0.076 0.106
## rmsea.pvalue rmsea.scaled
## 0.000 0.089
## rmsea.ci.lower.scaled rmsea.ci.upper.scaled
## 0.075 0.104
## rmsea.pvalue.scaled rmsea.robust
## 0.000 0.091
## rmsea.ci.lower.robust rmsea.ci.upper.robust
## 0.076 0.107
## rmsea.pvalue.robust rmr
## NA 255.803
## rmr_nomean srmr
## 295.376 0.040
## srmr_bentler srmr_bentler_nomean
## 0.040 0.046
## crmr crmr_nomean
## 0.046 0.056
## srmr_mplus srmr_mplus_nomean
## 0.040 0.046
## cn_05 cn_01
## 255.743 348.150
## gfi agfi
## 1.000 0.999
## pgfi mfi
## 0.250 0.980
## ecvi
## 0.056#fit measures for first simple model)
fitMeasures(pc2016.fit)## npar fmin chisq
## 10.000 0.022 105.776
## df pvalue baseline.chisq
## 5.000 0.000 1692.805
## baseline.df baseline.pvalue cfi
## 10.000 0.000 0.940
## tli nnfi rfi
## 0.880 0.880 0.875
## nfi pnfi ifi
## 0.938 0.469 0.940
## rni logl unrestricted.logl
## 0.940 -53841.964 -53789.076
## aic bic ntotal
## 107703.928 107761.901 2434.000
## bic2 rmsea rmsea.ci.lower
## 107730.128 0.091 0.076
## rmsea.ci.upper rmsea.pvalue rmr
## 0.106 0.000 295.377
## rmr_nomean srmr srmr_bentler
## 295.377 0.046 0.046
## srmr_bentler_nomean crmr crmr_nomean
## 0.046 0.056 0.056
## srmr_mplus srmr_mplus_nomean cn_05
## 0.046 0.046 255.743
## cn_01 gfi agfi
## 348.150 0.983 0.950
## pgfi mfi ecvi
## 0.328 0.980 0.052#parameter table
parameterTable(pc2016.fit)## id lhs op rhs user block group free
## 1 1 imm =~ pc_immun_class 1 1 1 1
## 2 2 imm =~ IRSAD_SCORE 1 1 1 2
## 3 3 imm =~ political_score 1 1 1 3
## 4 4 imm =~ govt_funds_per_person 1 1 1 4
## 5 5 imm =~ mean_tax_000s 1 1 1 5
## 6 6 pc_immun_class ~~ pc_immun_class 0 1 1 6
## 7 7 IRSAD_SCORE ~~ IRSAD_SCORE 0 1 1 7
## 8 8 political_score ~~ political_score 0 1 1 8
## 9 9 govt_funds_per_person ~~ govt_funds_per_person 0 1 1 9
## 10 10 mean_tax_000s ~~ mean_tax_000s 0 1 1 10
## 11 11 imm ~~ imm 0 1 1 0
## ustart exo label plabel start est se
## 1 NA 0 .p1. 1.275 0.775 0.077
## 2 NA 0 .p2. 53.194 50.348 2.628
## 3 NA 0 .p3. -0.108 -0.077 0.033
## 4 NA 0 .p4. 78.262 60.321 13.239
## 5 NA 0 .p5. 14.385 16.809 0.788
## 6 NA 0 .p6. 5.990 11.380 0.331
## 7 NA 0 .p7. 2873.802 3212.642 244.649
## 8 NA 0 .p8. 1.354 2.702 0.077
## 9 NA 0 .p9. 206960.112 410281.632 11764.479
## 10 NA 0 .p10. 139.984 -2.577 25.263
## 11 1 0 .p11. 1.000 1.000 0.000 #read in the data for 2015 only
#try this to stop data confusion
char_data <- read.csv('../cleaned_data/2015_dataREDUCED.csv')
str(char_data)## 'data.frame': 7559 obs. of 18 variables:
## $ postcode : int 800 800 800 810 810 810 812 812 812 820 ...
## $ state : Factor w/ 13 levels "","ACT","NSW",..: 7 7 7 7 7 7 7 7 7 7 ...
## $ year : int 2015 2015 2015 2015 2015 2015 2015 2015 2015 2015 ...
## $ age : int 1 5 2 2 1 5 2 5 1 2 ...
## $ pc_immun : Factor w/ 10 levels "","<70.0","70.0-74.9",..: 6 6 4 7 7 7 7 7 8 7 ...
## $ caution : int 1 1 1 0 0 0 0 0 0 0 ...
## $ pc_immun_class : int 5 5 3 6 6 6 6 6 7 6 ...
## $ PHN_code : Factor w/ 32 levels "","PHN101","PHN102",..: 31 31 31 31 31 31 31 31 31 31 ...
## $ PHN_number : int 701 701 701 701 701 701 701 701 701 701 ...
## $ IEO_SCORE : int 1077 1077 1077 1051 1051 1051 1007 1007 1007 1081 ...
## $ IER_SCORE : int 952 952 952 1008 1008 1008 1015 1015 1015 1008 ...
## $ IRSAD_SCORE : int 1072 1072 1072 1037 1037 1037 1019 1019 1019 1076 ...
## $ IRSAD_URP : int 4564 4564 4564 29725 29725 29725 19334 19334 19334 19057 ...
## $ IRSD_SCORE : int 1060 1060 1060 1027 1027 1027 1013 1013 1013 1067 ...
## $ electorate : Factor w/ 151 levels "","Adelaide",..: 137 137 137 137 137 137 137 137 137 137 ...
## $ mean_tax_000s : num 71.3 71.3 71.3 64.7 64.7 ...
## $ political_score: num 6.5 6.5 6.5 6.5 6.5 6.5 6.5 6.5 6.5 6.5 ...
## $ per_capita : Factor w/ 2178 levels "","#VALUE!","0",..: 223 223 223 2018 2018 2018 1597 1597 1597 1582 ...#Need to make the factors ordered factors to work in Lavaan
char_data$postcode = postcode_factor_ordered <- factor(char_data$postcode, ordered = TRUE)
char_data$state = state_factor_ordered <- factor(char_data$state, ordered = TRUE)
char_data$pc_immun = pc_immun_factor_ordered <- factor(char_data$pc_immun, ordered = TRUE)
char_data$PHN_code = PHN_code_factor_ordered <- factor(char_data$PHN_code, ordered = TRUE)
char_data$electorate = electorate_factor_ordered <- factor(char_data$electorate, ordered = TRUE)
char_data$per_capita = funds_ordered <- as.numeric(char_data$per_capita)
glimpse(char_data)## Observations: 7,559
## Variables: 18
## $ postcode <ord> 800, 800, 800, 810, 810, 810, 812, 812, 812, 8...
## $ state <ord> NT, NT, NT, NT, NT, NT, NT, NT, NT, NT, NT, NT...
## $ year <int> 2015, 2015, 2015, 2015, 2015, 2015, 2015, 2015...
## $ age <int> 1, 5, 2, 2, 1, 5, 2, 5, 1, 2, 5, 1, 1, 5, 2, 1...
## $ pc_immun <ord> 85.0-89.9, 85.0-89.9, 75.0-79.9, 90.0-92.4, 90...
## $ caution <int> 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0...
## $ pc_immun_class <int> 5, 5, 3, 6, 6, 6, 6, 6, 7, 6, 5, 6, 7, 7, 7, 0...
## $ PHN_code <ord> PHN701, PHN701, PHN701, PHN701, PHN701, PHN701...
## $ PHN_number <int> 701, 701, 701, 701, 701, 701, 701, 701, 701, 7...
## $ IEO_SCORE <int> 1077, 1077, 1077, 1051, 1051, 1051, 1007, 1007...
## $ IER_SCORE <int> 952, 952, 952, 1008, 1008, 1008, 1015, 1015, 1...
## $ IRSAD_SCORE <int> 1072, 1072, 1072, 1037, 1037, 1037, 1019, 1019...
## $ IRSAD_URP <int> 4564, 4564, 4564, 29725, 29725, 29725, 19334, ...
## $ IRSD_SCORE <int> 1060, 1060, 1060, 1027, 1027, 1027, 1013, 1013...
## $ electorate <ord> Solomon, Solomon, Solomon, Solomon, Solomon, S...
## $ mean_tax_000s <dbl> 71.263, 71.263, 71.263, 64.719, 64.719, 64.719...
## $ political_score <dbl> 6.50, 6.50, 6.50, 6.50, 6.50, 6.50, 6.50, 6.50...
## $ per_capita <dbl> 223, 223, 223, 2018, 2018, 2018, 1597, 1597, 1...Calculate per capita spend Govt spent $8.4m, population total was 24 127 200 in 2015-16
calculate government funding per postcode
#We have the SEIFA URP which has the usual resident population per postcode #First get per person amount
#total_population15 <- (24127200)
#total_immun_spend15 <- (8400000)
#per_capita15 <- total_population15/total_immun_spend15
#per_capita15
#2.87228 spend per personFilter for 5 year olds
#filter just by 5 year olds - as the per capita info in spreadsheet is duplicating by age group
Five_Data15 = char_data%>%
filter(age == 5)
#make sure all are ordered factors, as data seems to have changed from previous str
Five_Data15$per_capita = funds_ordered <- as.numeric(Five_Data15$per_capita)
glimpse(Five_Data15)## Observations: 2,422
## Variables: 18
## $ postcode <ord> 800, 810, 812, 820, 822, 828, 830, 832, 835, 8...
## $ state <ord> NT, NT, NT, NT, NT, NT, NT, NT, NT, NT, NT, NT...
## $ year <int> 2015, 2015, 2015, 2015, 2015, 2015, 2015, 2015...
## $ age <int> 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5...
## $ pc_immun <ord> 85.0-89.9, 90.0-92.4, 90.0-92.4, 85.0-89.9, 92...
## $ caution <int> 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0...
## $ pc_immun_class <int> 5, 6, 6, 5, 7, 0, 7, 7, 8, 7, 0, 0, 0, 0, 0, 0...
## $ PHN_code <ord> PHN701, PHN701, PHN701, PHN701, PHN701, PHN701...
## $ PHN_number <int> 701, 701, 701, 701, 701, 701, 701, 701, 701, 7...
## $ IEO_SCORE <int> 1077, 1051, 1007, 1081, 882, 973, 966, 1013, 9...
## $ IER_SCORE <int> 952, 1008, 1015, 1008, 645, 958, 991, 1060, 10...
## $ IRSAD_SCORE <int> 1072, 1037, 1019, 1076, 702, 972, 996, 1076, 1...
## $ IRSAD_URP <int> 4564, 29725, 19334, 19057, 24120, 1699, 17990,...
## $ IRSD_SCORE <int> 1060, 1027, 1013, 1067, 630, 970, 991, 1069, 1...
## $ electorate <ord> Solomon, Solomon, Solomon, Solomon, Lingiari, ...
## $ mean_tax_000s <dbl> 71.263, 64.719, 63.273, 79.588, 45.882, 71.860...
## $ political_score <dbl> 6.50, 6.50, 6.50, 6.50, 5.25, 6.50, 6.50, 6.50...
## $ per_capita <dbl> 223, 2018, 1597, 1582, 1821, 1460, 1517, 877, ...str(Five_Data15)## 'data.frame': 2422 obs. of 18 variables:
## $ postcode : Ord.factor w/ 2462 levels "800"<"810"<"812"<..: 1 2 3 4 5 6 7 8 9 10 ...
## $ state : Ord.factor w/ 13 levels ""<"ACT"<"NSW"<..: 7 7 7 7 7 7 7 7 7 7 ...
## $ year : int 2015 2015 2015 2015 2015 2015 2015 2015 2015 2015 ...
## $ age : int 5 5 5 5 5 5 5 5 5 5 ...
## $ pc_immun : Ord.factor w/ 10 levels ""<"<70.0"<"70.0-74.9"<..: 6 7 7 6 8 10 8 8 9 8 ...
## $ caution : int 1 0 0 0 0 0 0 0 1 1 ...
## $ pc_immun_class : int 5 6 6 5 7 0 7 7 8 7 ...
## $ PHN_code : Ord.factor w/ 32 levels ""<"PHN101"<"PHN102"<..: 31 31 31 31 31 31 31 31 31 31 ...
## $ PHN_number : int 701 701 701 701 701 701 701 701 701 701 ...
## $ IEO_SCORE : int 1077 1051 1007 1081 882 973 966 1013 991 968 ...
## $ IER_SCORE : int 952 1008 1015 1008 645 958 991 1060 1087 1088 ...
## $ IRSAD_SCORE : int 1072 1037 1019 1076 702 972 996 1076 1054 1041 ...
## $ IRSAD_URP : int 4564 29725 19334 19057 24120 1699 17990 9671 5194 7092 ...
## $ IRSD_SCORE : int 1060 1027 1013 1067 630 970 991 1069 1054 1035 ...
## $ electorate : Ord.factor w/ 151 levels ""<"Adelaide"<..: 137 137 137 137 91 137 137 137 91 91 ...
## $ mean_tax_000s : num 71.3 64.7 63.3 79.6 45.9 ...
## $ political_score: num 6.5 6.5 6.5 6.5 5.25 6.5 6.5 6.5 4 4 ...
## $ per_capita : num 223 2018 1597 1582 1821 ...#Our geo factor is measured by PHN, IRSAD score and pc immun class
geo_2015.model <- 'imm =~ mean_tax_000s + political_score + IRSAD_SCORE + per_capita + pc_immun_class'
geo_2015.fit <- cfa(model = geo_2015.model,
data = Five_Data15)## Warning in lav_data_full(data = data, group = group, cluster = cluster, :
## lavaan WARNING: some observed variances are (at least) a factor 1000 times
## larger than others; use varTable(fit) to investigatesummary(geo_2015.fit, standardized = TRUE, fit.measures = TRUE)## lavaan 0.6-3 ended normally after 117 iterations
##
## Optimization method NLMINB
## Number of free parameters 10
##
## Used Total
## Number of observations 2388 2422
##
## Estimator ML
## Model Fit Test Statistic 141.382
## Degrees of freedom 5
## P-value (Chi-square) 0.000
##
## Model test baseline model:
##
## Minimum Function Test Statistic 1440.151
## Degrees of freedom 10
## P-value 0.000
##
## User model versus baseline model:
##
## Comparative Fit Index (CFI) 0.905
## Tucker-Lewis Index (TLI) 0.809
##
## Loglikelihood and Information Criteria:
##
## Loglikelihood user model (H0) -52778.945
## Loglikelihood unrestricted model (H1) -52708.254
##
## Number of free parameters 10
## Akaike (AIC) 105577.890
## Bayesian (BIC) 105635.672
## Sample-size adjusted Bayesian (BIC) 105603.900
##
## Root Mean Square Error of Approximation:
##
## RMSEA 0.107
## 90 Percent Confidence Interval 0.092 0.122
## P-value RMSEA <= 0.05 0.000
##
## Standardized Root Mean Square Residual:
##
## SRMR 0.052
##
## Parameter Estimates:
##
## Information Expected
## Information saturated (h1) model Structured
## Standard Errors Standard
##
## Latent Variables:
## Estimate Std.Err z-value P(>|z|) Std.lv
## imm =~
## mean_tax_000s 1.000 15.756
## political_scor -0.004 0.002 -1.954 0.051 -0.069
## IRSAD_SCORE 3.177 0.307 10.345 0.000 50.052
## per_capita 1.720 0.879 1.955 0.051 27.093
## pc_immun_class 0.052 0.006 8.280 0.000 0.826
## Std.all
##
## 0.939
## -0.043
## 0.660
## 0.043
## 0.246
##
## Variances:
## Estimate Std.Err z-value P(>|z|) Std.lv
## .mean_tax_000s 33.022 23.330 1.415 0.157 33.022
## .political_scor 2.566 0.074 34.546 0.000 2.566
## .IRSAD_SCORE 3248.726 253.341 12.824 0.000 3248.726
## .per_capita 398235.153 11527.582 34.546 0.000 398235.153
## .pc_immun_class 10.565 0.312 33.807 0.000 10.565
## imm 248.255 24.673 10.062 0.000 1.000
## Std.all
## 0.117
## 0.998
## 0.565
## 0.998
## 0.939
## 1.000lavInspect(geo_2015.fit, what="estimates")## $lambda
## imm
## mean_tax_000s 1.000
## political_score -0.004
## IRSAD_SCORE 3.177
## per_capita 1.720
## pc_immun_class 0.052
##
## $theta
## m__000 pltcl_ IRSAD_ pr_cpt pc_mm_
## mean_tax_000s 33.022
## political_score 0.000 2.566
## IRSAD_SCORE 0.000 0.000 3248.726
## per_capita 0.000 0.000 0.000 398235.153
## pc_immun_class 0.000 0.000 0.000 0.000 10.565
##
## $psi
## imm
## imm 248.255semPaths(geo_2015.fit, what = "stand", rotation = 4)
semCors(geo_2015.fit)
semCors(geo_2015.fit, titles = TRUE)
Now let’s build a more complete model using the two years to see if things can be confirmed
# we will expect the mean tax data to be correlated, as it was repeated across the years due to lack of availability in 2016