library(magrittr)
library(Zelig)
library(dplyr)
library(texreg)
library(radiant.data)
data(titanic)
head(titanic)

Use Zelig 5 syntax to recreate slides 5.2 -> 5.18

titanic <- mutate(titanic, survival = ifelse(survived == "No",0,1))
min(titanic$age)
[1] 0.1667
max(titanic$age)
[1] 80
min(titanic$fare)
[1] 0
max(titanic$fare)
[1] 512.3292

In this chunk I created a new variable called survival, which is the numerical version of the variable “survived”, where 0 is equivlenet to not surviving, and 1 is equivelent to surviving. Then, I obtained the maximum and minimum values for age and fare so that I could use them in my range for the following analyses.

Age effect Using Zelig 5:

zhw.tit <- zlogit$new()
zhw.tit$zelig(survival ~ age + sex*pclass + fare, data = titanic)
summary(zhw.tit)
Model: 

Call:
zhw.tit$zelig(formula = survival ~ age + sex * pclass + fare, 
    data = titanic)

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-3.0634  -0.6641  -0.4943   0.4336   2.4941  

Coefficients:
                    Estimate Std. Error z value            Pr(>|z|)
(Intercept)        4.8978215  0.6131092   7.988 0.00000000000000137
age               -0.0387245  0.0068044  -5.691 0.00000001262563680
sexmale           -3.8996177  0.5015659  -7.775 0.00000000000000755
pclass2nd         -1.5923247  0.6024844  -2.643             0.00822
pclass3rd         -4.1382715  0.5601819  -7.387 0.00000000000014976
fare              -0.0009058  0.0020509  -0.442             0.65874
sexmale:pclass2nd -0.0600076  0.6372949  -0.094             0.92498
sexmale:pclass3rd  2.5019110  0.5479901   4.566 0.00000498035051247

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 1409.99  on 1042  degrees of freedom
Residual deviance:  931.45  on 1035  degrees of freedom
AIC: 947.45

Number of Fisher Scoring iterations: 5

Next step: Use 'setx' method
zhw.tit$setrange(age = 0:80)
zhw.tit$sim()
summary(zhw.tit)

[1] 0


 sim range :
 -----
ev
          mean         sd      50%      2.5%     97.5%
[1,] 0.3408924 0.04764033 0.338573 0.2506044 0.4325741
pv
         0     1
[1,] 0.638 0.362


[1] 1


 sim range :
 -----
ev
          mean         sd       50%    2.5%     97.5%
[1,] 0.3322111 0.04606172 0.3299581 0.24498 0.4214041
pv
         0     1
[1,] 0.666 0.334


[1] 2


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.3236374 0.04450769 0.3209356 0.2393563 0.4096313
pv
         0     1
[1,] 0.673 0.327


[1] 3


 sim range :
 -----
ev
          mean         sd       50%      2.5%    97.5%
[1,] 0.3151762 0.04298203 0.3128584 0.2341375 0.398604
pv
         0     1
[1,] 0.687 0.313


[1] 4


 sim range :
 -----
ev
          mean        sd       50%      2.5%     97.5%
[1,] 0.3068322 0.0414884 0.3049756 0.2286534 0.3889693
pv
       0   1
[1,] 0.7 0.3


[1] 5


 sim range :
 -----
ev
          mean         sd      50%      2.5%     97.5%
[1,] 0.2986097 0.04003024 0.296884 0.2235053 0.3789114
pv
        0    1
[1,] 0.68 0.32


[1] 6


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.2905127 0.03861081 0.2883608 0.2185785 0.3673043
pv
         0     1
[1,] 0.719 0.281


[1] 7


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.2825449 0.03723314 0.2804607 0.2136908 0.3564031
pv
         0     1
[1,] 0.725 0.275


[1] 8


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.2747096 0.03590001 0.2726663 0.2087691 0.3485728
pv
         0     1
[1,] 0.717 0.283


[1] 9


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.2670098 0.03461397 0.2651861 0.2038288 0.3404221
pv
         0     1
[1,] 0.757 0.243


[1] 10


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.2594483 0.03337724 0.2580928 0.1982964 0.3297163
pv
        0    1
[1,] 0.74 0.26


[1] 11


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.2520274 0.03219178 0.2510753 0.1931566 0.3187075
pv
         0     1
[1,] 0.721 0.279


[1] 12


 sim range :
 -----
ev
          mean         sd       50%     2.5%    97.5%
[1,] 0.2447491 0.03105919 0.2440137 0.188478 0.308571
pv
         0     1
[1,] 0.769 0.231


[1] 13


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.2376151 0.02998078 0.2371124 0.1829962 0.2988492
pv
         0     1
[1,] 0.756 0.244


[1] 14


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.2306268 0.02895746 0.2295916 0.1778751 0.2906243
pv
         0     1
[1,] 0.781 0.219


[1] 15


 sim range :
 -----
ev
          mean        sd       50%      2.5%     97.5%
[1,] 0.2237855 0.0279898 0.2228629 0.1725218 0.2812586
pv
         0     1
[1,] 0.777 0.223


[1] 16


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.2170917 0.02707798 0.2161289 0.1672924 0.2717886
pv
         0     1
[1,] 0.762 0.238


[1] 17


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.2105462 0.02622181 0.2099987 0.1621451 0.2632305
pv
       0   1
[1,] 0.8 0.2


[1] 18


 sim range :
 -----
ev
          mean         sd       50%      2.5%    97.5%
[1,] 0.2041491 0.02542071 0.2034513 0.1566427 0.254988
pv
         0     1
[1,] 0.799 0.201


[1] 19


 sim range :
 -----
ev
          mean        sd       50%      2.5%     97.5%
[1,] 0.1979004 0.0246737 0.1974186 0.1513898 0.2475363
pv
       0   1
[1,] 0.8 0.2


[1] 20


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1917999 0.02397947 0.1910858 0.1472841 0.2391261
pv
         0     1
[1,] 0.823 0.177


[1] 21


 sim range :
 -----
ev
         mean         sd      50%     2.5%    97.5%
[1,] 0.185847 0.02333633 0.184897 0.142136 0.232736
pv
         0     1
[1,] 0.814 0.186


[1] 22


 sim range :
 -----
ev
          mean         sd       50%      2.5%    97.5%
[1,] 0.1800411 0.02274227 0.1792745 0.1376233 0.224976
pv
         0     1
[1,] 0.821 0.179


[1] 23


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1743813 0.02219502 0.1736403 0.1329147 0.2186248
pv
         0     1
[1,] 0.818 0.182


[1] 24


 sim range :
 -----
ev
          mean         sd       50%     2.5%     97.5%
[1,] 0.1688663 0.02169203 0.1683477 0.128156 0.2132391
pv
         0     1
[1,] 0.815 0.185


[1] 25


 sim range :
 -----
ev
         mean         sd       50%     2.5%     97.5%
[1,] 0.163495 0.02123055 0.1624726 0.123659 0.2073548
pv
         0     1
[1,] 0.853 0.147


[1] 26


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1582658 0.02080771 0.1571929 0.1194621 0.2018142
pv
         0     1
[1,] 0.852 0.148


[1] 27


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1531771 0.02042048 0.1524171 0.1151712 0.1964334
pv
         0     1
[1,] 0.848 0.152


[1] 28


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1482271 0.02006583 0.1473966 0.1108243 0.1907328
pv
         0     1
[1,] 0.855 0.145


[1] 29


 sim range :
 -----
ev
         mean         sd       50%      2.5%    97.5%
[1,] 0.143414 0.01974069 0.1425403 0.1067078 0.184952
pv
        0    1
[1,] 0.86 0.14


[1] 30


 sim range :
 -----
ev
          mean         sd      50%      2.5%     97.5%
[1,] 0.1387357 0.01944203 0.137839 0.1027338 0.1798054
pv
         0     1
[1,] 0.854 0.146


[1] 31


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1341901 0.01916692 0.1332829 0.09892269 0.1750962
pv
         0     1
[1,] 0.863 0.137


[1] 32


 sim range :
 -----
ev
         mean         sd       50%       2.5%    97.5%
[1,] 0.129775 0.01891252 0.1290234 0.09514884 0.170216
pv
         0     1
[1,] 0.855 0.145


[1] 33


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1254881 0.01867614 0.1245432 0.09098166 0.1663792
pv
         0     1
[1,] 0.893 0.107


[1] 34


 sim range :
 -----
ev
         mean         sd       50%       2.5%    97.5%
[1,] 0.121327 0.01845525 0.1202472 0.08701208 0.162725
pv
        0    1
[1,] 0.89 0.11


[1] 35


 sim range :
 -----
ev
          mean        sd       50%      2.5%     97.5%
[1,] 0.1172892 0.0182475 0.1162327 0.0839042 0.1589066
pv
         0     1
[1,] 0.884 0.116


[1] 36


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1133724 0.01805073 0.1121403 0.08099865 0.1553696
pv
         0     1
[1,] 0.866 0.134


[1] 37


 sim range :
 -----
ev
         mean         sd       50%       2.5%     97.5%
[1,] 0.109574 0.01786295 0.1083219 0.07798937 0.1514367
pv
         0     1
[1,] 0.904 0.096


[1] 38


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1058914 0.01768239 0.1048028 0.07468479 0.1475735
pv
        0    1
[1,] 0.89 0.11


[1] 39


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1023221 0.01750744 0.1012666 0.07145976 0.1438625
pv
         0     1
[1,] 0.901 0.099


[1] 40


 sim range :
 -----
ev
           mean         sd        50%       2.5%    97.5%
[1,] 0.09886341 0.01733671 0.09771301 0.06852918 0.140397
pv
         0     1
[1,] 0.892 0.108


[1] 41


 sim range :
 -----
ev
           mean         sd        50%       2.5%    97.5%
[1,] 0.09551277 0.01716895 0.09429206 0.06554115 0.137124
pv
         0     1
[1,] 0.905 0.095


[1] 42


 sim range :
 -----
ev
           mean         sd        50%       2.5%     97.5%
[1,] 0.09226756 0.01700309 0.09126092 0.06267465 0.1329059
pv
        0    1
[1,] 0.92 0.08


[1] 43


 sim range :
 -----
ev
           mean        sd        50%       2.5%     97.5%
[1,] 0.08912517 0.0168382 0.08798469 0.05971491 0.1292827
pv
         0     1
[1,] 0.922 0.078


[1] 44


 sim range :
 -----
ev
           mean         sd        50%       2.5%     97.5%
[1,] 0.08608299 0.01667349 0.08493439 0.05666966 0.1256791
pv
        0    1
[1,] 0.92 0.08


[1] 45


 sim range :
 -----
ev
           mean         sd       50%       2.5%     97.5%
[1,] 0.08313845 0.01650833 0.0819605 0.05405285 0.1223051
pv
         0     1
[1,] 0.923 0.077


[1] 46


 sim range :
 -----
ev
           mean         sd        50%       2.5%     97.5%
[1,] 0.08028896 0.01634217 0.07915197 0.05172897 0.1189409
pv
         0     1
[1,] 0.917 0.083


[1] 47


 sim range :
 -----
ev
           mean         sd        50%       2.5%     97.5%
[1,] 0.07753198 0.01617458 0.07648748 0.04960747 0.1160638
pv
         0     1
[1,] 0.907 0.093


[1] 48


 sim range :
 -----
ev
           mean         sd        50%       2.5%     97.5%
[1,] 0.07486498 0.01600523 0.07371709 0.04760905 0.1127564
pv
        0    1
[1,] 0.91 0.09


[1] 49


 sim range :
 -----
ev
           mean         sd        50%       2.5%     97.5%
[1,] 0.07228547 0.01583386 0.07116528 0.04564096 0.1095407
pv
         0     1
[1,] 0.945 0.055


[1] 50


 sim range :
 -----
ev
           mean         sd      50%       2.5%     97.5%
[1,] 0.06979099 0.01566029 0.068587 0.04364744 0.1066857
pv
         0     1
[1,] 0.931 0.069


[1] 51


 sim range :
 -----
ev
           mean         sd        50%       2.5%     97.5%
[1,] 0.06737911 0.01548441 0.06604288 0.04169408 0.1033801
pv
        0    1
[1,] 0.93 0.07


[1] 52


 sim range :
 -----
ev
           mean         sd        50%       2.5%     97.5%
[1,] 0.06504744 0.01530617 0.06360829 0.03982451 0.1003921
pv
         0     1
[1,] 0.933 0.067


[1] 53


 sim range :
 -----
ev
           mean         sd        50%       2.5%      97.5%
[1,] 0.06279364 0.01512555 0.06135914 0.03761477 0.09765145
pv
        0    1
[1,] 0.93 0.07


[1] 54


 sim range :
 -----
ev
          mean         sd        50%       2.5%      97.5%
[1,] 0.0606154 0.01494259 0.05913281 0.03557768 0.09520395
pv
         0     1
[1,] 0.931 0.069


[1] 55


 sim range :
 -----
ev
           mean         sd        50%       2.5%      97.5%
[1,] 0.05851046 0.01475736 0.05708026 0.03407753 0.09281149
pv
         0     1
[1,] 0.944 0.056


[1] 56


 sim range :
 -----
ev
           mean         sd        50%       2.5%      97.5%
[1,] 0.05647659 0.01456996 0.05496943 0.03264022 0.09047314
pv
         0     1
[1,] 0.945 0.055


[1] 57


 sim range :
 -----
ev
           mean        sd        50%       2.5%      97.5%
[1,] 0.05451163 0.0143805 0.05288393 0.03125827 0.08818798
pv
         0     1
[1,] 0.936 0.064


[1] 58


 sim range :
 -----
ev
           mean         sd        50%      2.5%      97.5%
[1,] 0.05261344 0.01418914 0.05103357 0.0299346 0.08602259
pv
         0     1
[1,] 0.941 0.059


[1] 59


 sim range :
 -----
ev
           mean         sd        50%       2.5%      97.5%
[1,] 0.05077995 0.01399602 0.04920966 0.02866532 0.08436905
pv
         0     1
[1,] 0.955 0.045


[1] 60


 sim range :
 -----
ev
           mean         sd        50%       2.5%      97.5%
[1,] 0.04900914 0.01380133 0.04744387 0.02735175 0.08239429
pv
         0     1
[1,] 0.949 0.051


[1] 61


 sim range :
 -----
ev
           mean         sd        50%       2.5%      97.5%
[1,] 0.04729901 0.01360524 0.04568749 0.02610953 0.08035663
pv
         0     1
[1,] 0.953 0.047


[1] 62


 sim range :
 -----
ev
           mean         sd        50%       2.5%      97.5%
[1,] 0.04564763 0.01340794 0.04399322 0.02485707 0.07812202
pv
         0     1
[1,] 0.963 0.037


[1] 63


 sim range :
 -----
ev
           mean         sd        50%       2.5%      97.5%
[1,] 0.04405312 0.01320961 0.04239603 0.02369232 0.07598909
pv
         0     1
[1,] 0.945 0.055


[1] 64


 sim range :
 -----
ev
           mean         sd        50%      2.5%      97.5%
[1,] 0.04251365 0.01301046 0.04081985 0.0225031 0.07386137
pv
         0     1
[1,] 0.952 0.048


[1] 65


 sim range :
 -----
ev
           mean         sd        50%       2.5%      97.5%
[1,] 0.04102741 0.01281068 0.03932936 0.02137849 0.07209718
pv
         0     1
[1,] 0.961 0.039


[1] 66


 sim range :
 -----
ev
           mean         sd        50%       2.5%      97.5%
[1,] 0.03959267 0.01261045 0.03784948 0.02028571 0.07036413
pv
         0     1
[1,] 0.965 0.035


[1] 67


 sim range :
 -----
ev
           mean         sd       50%       2.5%      97.5%
[1,] 0.03820774 0.01240996 0.0364661 0.01924816 0.06845363
pv
         0     1
[1,] 0.958 0.042


[1] 68


 sim range :
 -----
ev
           mean         sd        50%       2.5%      97.5%
[1,] 0.03687096 0.01220941 0.03515598 0.01831083 0.06659131
pv
         0     1
[1,] 0.968 0.032


[1] 69


 sim range :
 -----
ev
           mean         sd        50%       2.5%      97.5%
[1,] 0.03558075 0.01200896 0.03391262 0.01733884 0.06477614
pv
         0     1
[1,] 0.967 0.033


[1] 70


 sim range :
 -----
ev
           mean         sd        50%       2.5%      97.5%
[1,] 0.03433554 0.01180881 0.03270166 0.01645521 0.06300337
pv
         0     1
[1,] 0.969 0.031


[1] 71


 sim range :
 -----
ev
           mean        sd        50%      2.5%      97.5%
[1,] 0.03313382 0.0116091 0.03149383 0.0156159 0.06129065
pv
         0     1
[1,] 0.973 0.027


[1] 72


 sim range :
 -----
ev
           mean         sd        50%      2.5%      97.5%
[1,] 0.03197414 0.01141002 0.03029042 0.0148571 0.05968853
pv
         0     1
[1,] 0.974 0.026


[1] 73


 sim range :
 -----
ev
           mean         sd        50%       2.5%      97.5%
[1,] 0.03085506 0.01121171 0.02915377 0.01418001 0.05812936
pv
         0     1
[1,] 0.966 0.034


[1] 74


 sim range :
 -----
ev
           mean         sd        50%       2.5%      97.5%
[1,] 0.02977521 0.01101432 0.02806034 0.01353337 0.05661157
pv
         0     1
[1,] 0.964 0.036


[1] 75


 sim range :
 -----
ev
           mean         sd        50%       2.5%      97.5%
[1,] 0.02873326 0.01081799 0.02709099 0.01288839 0.05513708
pv
         0     1
[1,] 0.978 0.022


[1] 76


 sim range :
 -----
ev
          mean         sd        50%       2.5%      97.5%
[1,] 0.0277279 0.01062286 0.02611991 0.01224973 0.05369711
pv
         0     1
[1,] 0.978 0.022


[1] 77


 sim range :
 -----
ev
           mean         sd        50%      2.5%      97.5%
[1,] 0.02675789 0.01042905 0.02517552 0.0115985 0.05229188
pv
         0     1
[1,] 0.974 0.026


[1] 78


 sim range :
 -----
ev
           mean         sd        50%       2.5%      97.5%
[1,] 0.02582201 0.01023669 0.02425094 0.01097764 0.05095567
pv
         0     1
[1,] 0.962 0.038


[1] 79


 sim range :
 -----
ev
           mean         sd        50%       2.5%      97.5%
[1,] 0.02491908 0.01004588 0.02334606 0.01038967 0.04977603
pv
         0     1
[1,] 0.976 0.024


[1] 80


 sim range :
 -----
ev
           mean          sd       50%        2.5%      97.5%
[1,] 0.02404796 0.009856723 0.0224781 0.009834989 0.04862231
pv
         0     1
[1,] 0.986 0.014
ci.plot(zhw.tit)

In the above analysis, I estimated the age effect of surviving the titanic using zelig 5. First, I estimated the model, then I set used the command “set range” to set the maximum and minimum values for age, and ultimately run the simulation that estimates the probability of surviving based on age.

The results show that there is a difference in the probability of surviving based on ones age, such that the older one gets, the less likely that they are to survive the titanic. There is an approximate .316 (31.6%) difference in the mean survival probability in 0 year olds as compared to 80 year olds.

Fare effect using zelig 5:

zhw.tit <- zlogit$new()
zhw.tit$zelig(survival ~ age + sex*pclass + fare, data = titanic)
zhw.tit$setrange(fare = 0:512)
zhw.tit$sim()
summary(zhw.tit)

[1] 0


 sim range :
 -----
ev
          mean         sd       50%     2.5%     97.5%
[1,] 0.1442755 0.01906302 0.1433357 0.109647 0.1845318
pv
         0     1
[1,] 0.848 0.152


[1] 1


 sim range :
 -----
ev
          mean         sd      50%      2.5%     97.5%
[1,] 0.1441492 0.01901685 0.143177 0.1094387 0.1841484
pv
         0     1
[1,] 0.845 0.155


[1] 2


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1440232 0.01897395 0.1429603 0.1094303 0.1837656
pv
         0     1
[1,] 0.858 0.142


[1] 3


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1438977 0.01893433 0.1429317 0.1093683 0.1833821
pv
         0     1
[1,] 0.854 0.146


[1] 4


 sim range :
 -----
ev
          mean       sd       50%      2.5%     97.5%
[1,] 0.1437727 0.018898 0.1428719 0.1093475 0.1832709
pv
         0     1
[1,] 0.857 0.143


[1] 5


 sim range :
 -----
ev
          mean         sd      50%      2.5%     97.5%
[1,] 0.1436481 0.01886494 0.142819 0.1093267 0.1832169
pv
         0     1
[1,] 0.863 0.137


[1] 6


 sim range :
 -----
ev
          mean         sd       50%      2.5%    97.5%
[1,] 0.1435239 0.01883518 0.1427518 0.1089889 0.183142
pv
        0    1
[1,] 0.85 0.15


[1] 7


 sim range :
 -----
ev
          mean        sd       50%      2.5%     97.5%
[1,] 0.1434001 0.0188087 0.1425846 0.1088829 0.1828144
pv
         0     1
[1,] 0.844 0.156


[1] 8


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1432768 0.01878551 0.1424181 0.1089386 0.1828405
pv
        0    1
[1,] 0.86 0.14


[1] 9


 sim range :
 -----
ev
          mean        sd       50%      2.5%     97.5%
[1,] 0.1431539 0.0187656 0.1423067 0.1089013 0.1824621
pv
         0     1
[1,] 0.855 0.145


[1] 10


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1430315 0.01874897 0.1422091 0.1086962 0.1820812
pv
         0     1
[1,] 0.858 0.142


[1] 11


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1429095 0.01873561 0.1421165 0.1085023 0.1817011
pv
         0     1
[1,] 0.863 0.137


[1] 12


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1427879 0.01872551 0.1420775 0.1084628 0.1813399
pv
        0    1
[1,] 0.87 0.13


[1] 13


 sim range :
 -----
ev
          mean         sd       50%     2.5%     97.5%
[1,] 0.1426667 0.01871867 0.1419845 0.108489 0.1815834
pv
         0     1
[1,] 0.852 0.148


[1] 14


 sim range :
 -----
ev
         mean         sd       50%      2.5%     97.5%
[1,] 0.142546 0.01871506 0.1421199 0.1085197 0.1814228
pv
         0     1
[1,] 0.854 0.146


[1] 15


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1424258 0.01871469 0.1419527 0.1085283 0.1814929
pv
         0     1
[1,] 0.845 0.155


[1] 16


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1423059 0.01871754 0.1417979 0.1082246 0.1816744
pv
         0     1
[1,] 0.866 0.134


[1] 17


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1421865 0.01872358 0.1418154 0.1082301 0.1819223
pv
         0     1
[1,] 0.867 0.133


[1] 18


 sim range :
 -----
ev
          mean        sd       50%      2.5%     97.5%
[1,] 0.1420675 0.0187328 0.1415938 0.1082877 0.1820273
pv
         0     1
[1,] 0.884 0.116


[1] 19


 sim range :
 -----
ev
          mean         sd       50%      2.5%    97.5%
[1,] 0.1419489 0.01874519 0.1414432 0.1080748 0.182212
pv
         0     1
[1,] 0.843 0.157


[1] 20


 sim range :
 -----
ev
          mean         sd       50%      2.5%    97.5%
[1,] 0.1418308 0.01876072 0.1412979 0.1079171 0.182383
pv
         0     1
[1,] 0.864 0.136


[1] 21


 sim range :
 -----
ev
          mean         sd       50%      2.5%    97.5%
[1,] 0.1417131 0.01877937 0.1410688 0.1080168 0.182335
pv
        0    1
[1,] 0.88 0.12


[1] 22


 sim range :
 -----
ev
          mean         sd      50%      2.5%    97.5%
[1,] 0.1415959 0.01880112 0.140884 0.1081166 0.182287
pv
        0    1
[1,] 0.87 0.13


[1] 23


 sim range :
 -----
ev
         mean         sd       50%      2.5%    97.5%
[1,] 0.141479 0.01882594 0.1407761 0.1081109 0.182237
pv
         0     1
[1,] 0.865 0.135


[1] 24


 sim range :
 -----
ev
          mean         sd      50%      2.5%     97.5%
[1,] 0.1413626 0.01885381 0.140691 0.1079614 0.1821821
pv
         0     1
[1,] 0.851 0.149


[1] 25


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1412466 0.01888471 0.1405695 0.1078121 0.1822489
pv
        0    1
[1,] 0.86 0.14


[1] 26


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1411311 0.01891859 0.1404214 0.1076965 0.1824328
pv
         0     1
[1,] 0.833 0.167


[1] 27


 sim range :
 -----
ev
         mean         sd       50%      2.5%     97.5%
[1,] 0.141016 0.01895544 0.1404051 0.1074984 0.1821396
pv
         0     1
[1,] 0.872 0.128


[1] 28


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1409013 0.01899523 0.1404247 0.1073377 0.1822926
pv
         0     1
[1,] 0.859 0.141


[1] 29


 sim range :
 -----
ev
         mean         sd       50%      2.5%     97.5%
[1,] 0.140787 0.01903791 0.1402814 0.1071772 0.1819073
pv
         0     1
[1,] 0.844 0.156


[1] 30


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1406731 0.01908347 0.1400835 0.1069562 0.1818524
pv
         0     1
[1,] 0.864 0.136


[1] 31


 sim range :
 -----
ev
          mean         sd      50%    2.5%     97.5%
[1,] 0.1405597 0.01913187 0.140079 0.10666 0.1817975
pv
         0     1
[1,] 0.855 0.145


[1] 32


 sim range :
 -----
ev
          mean         sd      50%      2.5%     97.5%
[1,] 0.1404467 0.01918306 0.139728 0.1063644 0.1817427
pv
         0     1
[1,] 0.838 0.162


[1] 33


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1403341 0.01923703 0.1395777 0.1060702 0.1816878
pv
         0     1
[1,] 0.847 0.153


[1] 34


 sim range :
 -----
ev
         mean         sd       50%      2.5%     97.5%
[1,] 0.140222 0.01929374 0.1395325 0.1058045 0.1816329
pv
        0    1
[1,] 0.84 0.16


[1] 35


 sim range :
 -----
ev
          mean         sd       50%     2.5%    97.5%
[1,] 0.1401103 0.01935314 0.1395436 0.105472 0.181657
pv
         0     1
[1,] 0.868 0.132


[1] 36


 sim range :
 -----
ev
         mean        sd       50%      2.5%     97.5%
[1,] 0.139999 0.0194152 0.1394828 0.1051732 0.1815944
pv
         0     1
[1,] 0.856 0.144


[1] 37


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1398881 0.01947989 0.1395641 0.1048751 0.1815074
pv
         0     1
[1,] 0.864 0.136


[1] 38


 sim range :
 -----
ev
          mean         sd       50%      2.5%    97.5%
[1,] 0.1397776 0.01954717 0.1394029 0.1045713 0.181445
pv
         0     1
[1,] 0.859 0.141


[1] 39


 sim range :
 -----
ev
          mean       sd       50%      2.5%    97.5%
[1,] 0.1396676 0.019617 0.1394221 0.1042856 0.181383
pv
         0     1
[1,] 0.858 0.142


[1] 40


 sim range :
 -----
ev
          mean         sd     50%      2.5%     97.5%
[1,] 0.1395579 0.01968934 0.13923 0.1039897 0.1813201
pv
         0     1
[1,] 0.856 0.144


[1] 41


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1394487 0.01976415 0.1390724 0.1036946 0.1812573
pv
         0     1
[1,] 0.871 0.129


[1] 42


 sim range :
 -----
ev
        mean         sd      50%      2.5%    97.5%
[1,] 0.13934 0.01984139 0.138904 0.1034001 0.181004
pv
         0     1
[1,] 0.849 0.151


[1] 43


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1392316 0.01992103 0.1387651 0.1031086 0.1809076
pv
         0     1
[1,] 0.865 0.135


[1] 44


 sim range :
 -----
ev
          mean         sd       50%     2.5%     97.5%
[1,] 0.1391237 0.02000303 0.1385741 0.102819 0.1808113
pv
         0     1
[1,] 0.862 0.138


[1] 45


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1390161 0.02008735 0.1384046 0.1025303 0.1809317
pv
         0     1
[1,] 0.863 0.137


[1] 46


 sim range :
 -----
ev
         mean         sd       50%      2.5%     97.5%
[1,] 0.138909 0.02017394 0.1382545 0.1023053 0.1809896
pv
         0     1
[1,] 0.863 0.137


[1] 47


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1388023 0.02026277 0.1380303 0.1022621 0.1809432
pv
         0     1
[1,] 0.867 0.133


[1] 48


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1386961 0.02035381 0.1378295 0.1022188 0.1808947
pv
         0     1
[1,] 0.857 0.143


[1] 49


 sim range :
 -----
ev
          mean       sd       50%      2.5%     97.5%
[1,] 0.1385902 0.020447 0.1376953 0.1018195 0.1808371
pv
         0     1
[1,] 0.857 0.143


[1] 50


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1384848 0.02054231 0.1376597 0.1014941 0.1807796
pv
         0     1
[1,] 0.887 0.113


[1] 51


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1383797 0.02063971 0.1375414 0.1013617 0.1807222
pv
         0     1
[1,] 0.863 0.137


[1] 52


 sim range :
 -----
ev
          mean         sd       50%      2.5%    97.5%
[1,] 0.1382751 0.02073915 0.1373773 0.1012295 0.180702
pv
        0    1
[1,] 0.86 0.14


[1] 53


 sim range :
 -----
ev
          mean        sd       50%      2.5%     97.5%
[1,] 0.1381709 0.0208406 0.1371845 0.1009367 0.1807691
pv
         0     1
[1,] 0.846 0.154


[1] 54


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1380671 0.02094402 0.1371252 0.1006217 0.1810383
pv
         0     1
[1,] 0.882 0.118


[1] 55


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1379637 0.02104936 0.1369604 0.1003798 0.1813062
pv
         0     1
[1,] 0.864 0.136


[1] 56


 sim range :
 -----
ev
          mean        sd      50%      2.5%    97.5%
[1,] 0.1378608 0.0211566 0.136747 0.1001732 0.181565
pv
         0     1
[1,] 0.863 0.137


[1] 57


 sim range :
 -----
ev
          mean        sd       50%       2.5%     97.5%
[1,] 0.1377582 0.0212657 0.1366485 0.09992644 0.1818367
pv
         0     1
[1,] 0.854 0.146


[1] 58


 sim range :
 -----
ev
          mean         sd       50%       2.5%    97.5%
[1,] 0.1376561 0.02137661 0.1365102 0.09963561 0.182106
pv
         0     1
[1,] 0.868 0.132


[1] 59


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1375544 0.02148931 0.1363722 0.09915953 0.1823602
pv
         0     1
[1,] 0.855 0.145


[1] 60


 sim range :
 -----
ev
         mean         sd       50%       2.5%     97.5%
[1,] 0.137453 0.02160375 0.1362812 0.09873359 0.1821857
pv
         0     1
[1,] 0.867 0.133


[1] 61


 sim range :
 -----
ev
          mean        sd       50%       2.5%     97.5%
[1,] 0.1373521 0.0217199 0.1361913 0.09842106 0.1819767
pv
         0     1
[1,] 0.839 0.161


[1] 62


 sim range :
 -----
ev
          mean         sd       50%       2.5%    97.5%
[1,] 0.1372516 0.02183773 0.1361844 0.09813641 0.181935
pv
        0    1
[1,] 0.86 0.14


[1] 63


 sim range :
 -----
ev
          mean        sd       50%       2.5%     97.5%
[1,] 0.1371515 0.0219572 0.1360874 0.09789608 0.1816895
pv
         0     1
[1,] 0.877 0.123


[1] 64


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1370518 0.02207828 0.1359951 0.09754171 0.1814971
pv
         0     1
[1,] 0.853 0.147


[1] 65


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1369526 0.02220093 0.1358353 0.0971836 0.1815612
pv
         0     1
[1,] 0.863 0.137


[1] 66


 sim range :
 -----
ev
          mean         sd       50%     2.5%     97.5%
[1,] 0.1368537 0.02232512 0.1357026 0.096873 0.1816254
pv
         0     1
[1,] 0.871 0.129


[1] 67


 sim range :
 -----
ev
          mean         sd      50%       2.5%     97.5%
[1,] 0.1367552 0.02245082 0.135553 0.09651861 0.1816896
pv
        0    1
[1,] 0.87 0.13


[1] 68


 sim range :
 -----
ev
          mean       sd       50%       2.5%     97.5%
[1,] 0.1366572 0.022578 0.1354519 0.09615835 0.1817539
pv
         0     1
[1,] 0.854 0.146


[1] 69


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1365595 0.02270662 0.1353298 0.0957993 0.1819427
pv
         0     1
[1,] 0.872 0.128


[1] 70


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1364622 0.02283665 0.1352336 0.09549434 0.1822918
pv
         0     1
[1,] 0.855 0.145


[1] 71


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1363654 0.02296806 0.1351505 0.09518596 0.1826415
pv
         0     1
[1,] 0.857 0.143


[1] 72


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1362689 0.02310082 0.1350507 0.09475259 0.1829917
pv
         0     1
[1,] 0.868 0.132


[1] 73


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1361729 0.02323491 0.1349006 0.09436603 0.1835872
pv
         0     1
[1,] 0.859 0.141


[1] 74


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1360773 0.02337029 0.1347403 0.09400892 0.1841754
pv
         0     1
[1,] 0.846 0.154


[1] 75


 sim range :
 -----
ev
         mean         sd       50%       2.5%     97.5%
[1,] 0.135982 0.02350693 0.1345561 0.09365301 0.1843711
pv
         0     1
[1,] 0.864 0.136


[1] 76


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1358872 0.02364482 0.1345109 0.09325865 0.1845671
pv
         0     1
[1,] 0.857 0.143


[1] 77


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1357927 0.02378391 0.1344696 0.09284827 0.1847874
pv
         0     1
[1,] 0.882 0.118


[1] 78


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1356987 0.02392418 0.1343604 0.09259776 0.1851526
pv
         0     1
[1,] 0.866 0.134


[1] 79


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1356051 0.02406561 0.1342803 0.09221738 0.1855173
pv
         0     1
[1,] 0.887 0.113


[1] 80


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1355118 0.02420817 0.1341811 0.09186334 0.1858685
pv
         0     1
[1,] 0.863 0.137


[1] 81


 sim range :
 -----
ev
         mean         sd       50%       2.5%     97.5%
[1,] 0.135419 0.02435184 0.1341089 0.09147888 0.1862195
pv
         0     1
[1,] 0.866 0.134


[1] 82


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1353265 0.02449659 0.1340148 0.0910624 0.1865674
pv
         0     1
[1,] 0.873 0.127


[1] 83


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1352345 0.02464239 0.1339753 0.09070423 0.1869159
pv
         0     1
[1,] 0.854 0.146


[1] 84


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1351428 0.02478922 0.1339078 0.09049159 0.1872164
pv
         0     1
[1,] 0.876 0.124


[1] 85


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1350516 0.02493707 0.1338123 0.09017395 0.1871131
pv
         0     1
[1,] 0.867 0.133


[1] 86


 sim range :
 -----
ev
          mean        sd       50%      2.5%     97.5%
[1,] 0.1349607 0.0250859 0.1337168 0.0898311 0.1870071
pv
         0     1
[1,] 0.877 0.123


[1] 87


 sim range :
 -----
ev
          mean        sd       50%       2.5%     97.5%
[1,] 0.1348702 0.0252357 0.1336491 0.08949529 0.1868963
pv
         0     1
[1,] 0.865 0.135


[1] 88


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1347801 0.02538643 0.1335631 0.08916178 0.1867939
pv
         0     1
[1,] 0.867 0.133


[1] 89


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1346905 0.02553809 0.1335439 0.0889535 0.1869767
pv
         0     1
[1,] 0.845 0.155


[1] 90


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1346012 0.02569065 0.1335048 0.08857193 0.1876384
pv
         0     1
[1,] 0.867 0.133


[1] 91


 sim range :
 -----
ev
          mean         sd      50%       2.5%     97.5%
[1,] 0.1345123 0.02584409 0.133399 0.08818941 0.1881846
pv
         0     1
[1,] 0.869 0.131


[1] 92


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1344238 0.02599839 0.1332468 0.08781109 0.1882462
pv
         0     1
[1,] 0.868 0.132


[1] 93


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1343356 0.02615353 0.1331485 0.08743095 0.1882967
pv
         0     1
[1,] 0.855 0.145


[1] 94


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1342479 0.02630949 0.1329488 0.0870523 0.1883603
pv
         0     1
[1,] 0.864 0.136


[1] 95


 sim range :
 -----
ev
          mean         sd      50%       2.5%     97.5%
[1,] 0.1341606 0.02646626 0.132835 0.08667943 0.1884357
pv
         0     1
[1,] 0.871 0.129


[1] 96


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1340736 0.02662381 0.1328517 0.0864558 0.1885111
pv
         0     1
[1,] 0.875 0.125


[1] 97


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1339871 0.02678213 0.1327584 0.08610243 0.1885866
pv
         0     1
[1,] 0.868 0.132


[1] 98


 sim range :
 -----
ev
          mean         sd       50%       2.5%    97.5%
[1,] 0.1339009 0.02694119 0.1326529 0.08577682 0.188709
pv
         0     1
[1,] 0.881 0.119


[1] 99


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1338151 0.02710099 0.1325476 0.0855058 0.1891595
pv
         0     1
[1,] 0.875 0.125


[1] 100


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1337297 0.02726151 0.1324423 0.08517963 0.1896117
pv
         0     1
[1,] 0.866 0.134


[1] 101


 sim range :
 -----
ev
          mean         sd      50%      2.5%     97.5%
[1,] 0.1336447 0.02742273 0.132337 0.0848216 0.1900647
pv
         0     1
[1,] 0.858 0.142


[1] 102


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1335601 0.02758463 0.1321555 0.08442945 0.1905186
pv
         0     1
[1,] 0.848 0.152


[1] 103


 sim range :
 -----
ev
          mean        sd       50%       2.5%     97.5%
[1,] 0.1334758 0.0277472 0.1320124 0.08403894 0.1909733
pv
        0    1
[1,] 0.87 0.13


[1] 104


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1333919 0.02791043 0.1318438 0.08365008 0.1914289
pv
         0     1
[1,] 0.871 0.129


[1] 105


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1333085 0.02807429 0.1318277 0.08326436 0.1918853
pv
         0     1
[1,] 0.868 0.132


[1] 106


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1332254 0.02823879 0.1317592 0.08288055 0.1923426
pv
        0    1
[1,] 0.87 0.13


[1] 107


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1331427 0.02840389 0.1315663 0.08249834 0.1928006
pv
         0     1
[1,] 0.855 0.145


[1] 108


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1330603 0.02856959 0.1314583 0.08220561 0.1932596
pv
         0     1
[1,] 0.872 0.128


[1] 109


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1329784 0.02873588 0.1313067 0.0818954 0.1934985
pv
         0     1
[1,] 0.863 0.137


[1] 110


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1328968 0.02890274 0.1311946 0.08161345 0.1936146
pv
         0     1
[1,] 0.875 0.125


[1] 111


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1328156 0.02907015 0.1310923 0.08132996 0.1937308
pv
         0     1
[1,] 0.862 0.138


[1] 112


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1327348 0.02923812 0.1309948 0.08104598 0.1938471
pv
         0     1
[1,] 0.869 0.131


[1] 113


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1326544 0.02940662 0.1307863 0.08076291 0.1943215
pv
         0     1
[1,] 0.867 0.133


[1] 114


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1325743 0.02957564 0.1307373 0.08048075 0.1948981
pv
         0     1
[1,] 0.889 0.111


[1] 115


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1324946 0.02974517 0.1306103 0.08020143 0.1954759
pv
         0     1
[1,] 0.871 0.129


[1] 116


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1324153 0.02991521 0.1304487 0.07992527 0.1960551
pv
         0     1
[1,] 0.876 0.124


[1] 117


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1323364 0.03008574 0.1302719 0.07964997 0.1966309
pv
         0     1
[1,] 0.884 0.116


[1] 118


 sim range :
 -----
ev
          mean         sd      50%       2.5%     97.5%
[1,] 0.1322578 0.03025674 0.130207 0.07937555 0.1972024
pv
         0     1
[1,] 0.854 0.146


[1] 119


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1321797 0.03042821 0.1301536 0.07902131 0.1974594
pv
         0     1
[1,] 0.873 0.127


[1] 120


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1321018 0.03060015 0.1300831 0.0786386 0.1975382
pv
         0     1
[1,] 0.867 0.133


[1] 121


 sim range :
 -----
ev
          mean         sd       50%       2.5%    97.5%
[1,] 0.1320244 0.03077253 0.1298889 0.07840951 0.197617
pv
         0     1
[1,] 0.884 0.116


[1] 122


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1319474 0.03094535 0.1298028 0.07818462 0.1979748
pv
         0     1
[1,] 0.864 0.136


[1] 123


 sim range :
 -----
ev
          mean        sd       50%       2.5%     97.5%
[1,] 0.1318707 0.0311186 0.1296971 0.07796033 0.1984476
pv
         0     1
[1,] 0.879 0.121


[1] 124


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1317944 0.03129227 0.1295333 0.07770242 0.1989183
pv
         0     1
[1,] 0.881 0.119


[1] 125


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1317184 0.03146636 0.1294643 0.07740252 0.1993898
pv
         0     1
[1,] 0.859 0.141


[1] 126


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1316428 0.03164085 0.1294413 0.07710368 0.1998622
pv
         0     1
[1,] 0.864 0.136


[1] 127


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1315676 0.03181573 0.1292945 0.07678952 0.2003354
pv
         0     1
[1,] 0.868 0.132


[1] 128


 sim range :
 -----
ev
          mean       sd       50%       2.5%     97.5%
[1,] 0.1314928 0.031991 0.1291234 0.07642593 0.2008094
pv
         0     1
[1,] 0.869 0.131


[1] 129


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1314183 0.03216665 0.1290101 0.07593839 0.2012843
pv
         0     1
[1,] 0.882 0.118


[1] 130


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1313442 0.03234267 0.1288978 0.07537793 0.2017006
pv
         0     1
[1,] 0.865 0.135


[1] 131


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1312705 0.03251906 0.1287441 0.07483005 0.2019609
pv
         0     1
[1,] 0.868 0.132


[1] 132


 sim range :
 -----
ev
          mean        sd       50%       2.5%     97.5%
[1,] 0.1311971 0.0326958 0.1286109 0.07452659 0.2022182
pv
         0     1
[1,] 0.877 0.123


[1] 133


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1311241 0.03287289 0.1285454 0.07423136 0.2024742
pv
         0     1
[1,] 0.878 0.122


[1] 134


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1310515 0.03305032 0.1284513 0.07393721 0.2027304
pv
         0     1
[1,] 0.871 0.129


[1] 135


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1309792 0.03322808 0.1283532 0.07358035 0.2029868
pv
         0     1
[1,] 0.868 0.132


[1] 136


 sim range :
 -----
ev
          mean         sd      50%      2.5%     97.5%
[1,] 0.1309073 0.03340618 0.128184 0.0732161 0.2032435
pv
         0     1
[1,] 0.844 0.156


[1] 137


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1308357 0.03358459 0.1279954 0.07295158 0.2035004
pv
         0     1
[1,] 0.873 0.127


[1] 138


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1307645 0.03376332 0.1278961 0.07268796 0.2037575
pv
         0     1
[1,] 0.883 0.117


[1] 139


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1306937 0.03394236 0.1277771 0.07242522 0.2040149
pv
         0     1
[1,] 0.886 0.114


[1] 140


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1306232 0.03412171 0.1276711 0.07216336 0.2042726
pv
        0    1
[1,] 0.86 0.14


[1] 141


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1305531 0.03430135 0.1276155 0.07182321 0.2045305
pv
         0     1
[1,] 0.872 0.128


[1] 142


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1304834 0.03448128 0.1275276 0.07147863 0.2047886
pv
         0     1
[1,] 0.879 0.121


[1] 143


 sim range :
 -----
ev
         mean        sd       50%       2.5%    97.5%
[1,] 0.130414 0.0346615 0.1273873 0.07113558 0.205047
pv
        0    1
[1,] 0.87 0.13


[1] 144


 sim range :
 -----
ev
         mean         sd       50%       2.5%     97.5%
[1,] 0.130345 0.03484199 0.1272344 0.07073266 0.2053056
pv
        0    1
[1,] 0.89 0.11


[1] 145


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1302763 0.03502277 0.1271353 0.07030353 0.2055644
pv
         0     1
[1,] 0.861 0.139


[1] 146


 sim range :
 -----
ev
         mean         sd      50%       2.5%     97.5%
[1,] 0.130208 0.03520381 0.127086 0.06987749 0.2058214
pv
         0     1
[1,] 0.876 0.124


[1] 147


 sim range :
 -----
ev
        mean         sd       50%       2.5%     97.5%
[1,] 0.13014 0.03538512 0.1270367 0.06945385 0.2060784
pv
         0     1
[1,] 0.864 0.136


[1] 148


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1300724 0.03556668 0.1269219 0.06903259 0.2063357
pv
         0     1
[1,] 0.881 0.119


[1] 149


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1300052 0.03574851 0.1267579 0.06861369 0.2065932
pv
         0     1
[1,] 0.886 0.114


[1] 150


 sim range :
 -----
ev
          mean         sd      50%       2.5%     97.5%
[1,] 0.1299383 0.03593058 0.126594 0.06819716 0.2069249
pv
         0     1
[1,] 0.872 0.128


[1] 151


 sim range :
 -----
ev
          mean        sd       50%       2.5%     97.5%
[1,] 0.1298718 0.0361129 0.1264868 0.06788004 0.2075433
pv
         0     1
[1,] 0.879 0.121


[1] 152


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1298056 0.03629546 0.1264136 0.06774489 0.2077961
pv
         0     1
[1,] 0.885 0.115


[1] 153


 sim range :
 -----
ev
          mean         sd       50%    2.5%     97.5%
[1,] 0.1297397 0.03647825 0.1262488 0.06761 0.2080165
pv
        0    1
[1,] 0.89 0.11


[1] 154


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1296743 0.03666128 0.1261052 0.06747537 0.2084029
pv
         0     1
[1,] 0.865 0.135


[1] 155


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1296091 0.03684454 0.1260066 0.06715577 0.2087898
pv
        0    1
[1,] 0.86 0.14


[1] 156


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1295443 0.03702802 0.1259485 0.06683382 0.2091773
pv
         0     1
[1,] 0.863 0.137


[1] 157


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1294799 0.03721173 0.1258737 0.06645527 0.2095653
pv
         0     1
[1,] 0.874 0.126


[1] 158


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1294158 0.03739565 0.1257772 0.06605088 0.2099538
pv
         0     1
[1,] 0.876 0.124


[1] 159


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1293521 0.03757978 0.1256722 0.06564879 0.2103429
pv
         0     1
[1,] 0.847 0.153


[1] 160


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1292887 0.03776412 0.1255549 0.06524897 0.2107325
pv
         0     1
[1,] 0.857 0.143


[1] 161


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1292257 0.03794867 0.1254091 0.06493905 0.2110866
pv
         0     1
[1,] 0.876 0.124


[1] 162


 sim range :
 -----
ev
         mean         sd       50%       2.5%     97.5%
[1,] 0.129163 0.03813342 0.1252914 0.06464953 0.2113863
pv
         0     1
[1,] 0.865 0.135


[1] 163


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1291006 0.03831837 0.1251599 0.06436121 0.2116864
pv
         0     1
[1,] 0.857 0.143


[1] 164


 sim range :
 -----
ev
          mean         sd      50%       2.5%     97.5%
[1,] 0.1290386 0.03850352 0.125063 0.06407668 0.2119868
pv
         0     1
[1,] 0.878 0.122


[1] 165


 sim range :
 -----
ev
         mean         sd       50%       2.5%     97.5%
[1,] 0.128977 0.03868885 0.1249836 0.06379474 0.2122875
pv
         0     1
[1,] 0.863 0.137


[1] 166


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1289157 0.03887438 0.1248248 0.06351396 0.2125885
pv
         0     1
[1,] 0.884 0.116


[1] 167


 sim range :
 -----
ev
          mean        sd       50%       2.5%     97.5%
[1,] 0.1288547 0.0390601 0.1247922 0.06323433 0.2128898
pv
         0     1
[1,] 0.862 0.138


[1] 168


 sim range :
 -----
ev
          mean       sd       50%       2.5%     97.5%
[1,] 0.1287941 0.039246 0.1247201 0.06295585 0.2131915
pv
         0     1
[1,] 0.877 0.123


[1] 169


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1287338 0.03943207 0.1246002 0.06267851 0.2135826
pv
         0     1
[1,] 0.869 0.131


[1] 170


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1286739 0.03961833 0.1245007 0.06240232 0.2140847
pv
         0     1
[1,] 0.862 0.138


[1] 171


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1286143 0.03980476 0.1243897 0.06215503 0.2145876
pv
         0     1
[1,] 0.864 0.136


[1] 172


 sim range :
 -----
ev
         mean         sd       50%       2.5%     97.5%
[1,] 0.128555 0.03999137 0.1242752 0.06187401 0.2150914
pv
         0     1
[1,] 0.868 0.132


[1] 173


 sim range :
 -----
ev
          mean         sd       50%       2.5%    97.5%
[1,] 0.1284961 0.04017814 0.1241589 0.06157736 0.215596
pv
        0    1
[1,] 0.87 0.13


[1] 174


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1284375 0.04036508 0.1240243 0.06130388 0.2161015
pv
         0     1
[1,] 0.876 0.124


[1] 175


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1283792 0.04055218 0.1238985 0.06103154 0.2166079
pv
         0     1
[1,] 0.856 0.144


[1] 176


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1283213 0.04073945 0.1237376 0.06076033 0.2170762
pv
         0     1
[1,] 0.861 0.139


[1] 177


 sim range :
 -----
ev
          mean         sd       50%       2.5%   97.5%
[1,] 0.1282638 0.04092688 0.1235946 0.06049025 0.21748
pv
        0    1
[1,] 0.88 0.12


[1] 178


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1282065 0.04111446 0.1234333 0.06022129 0.2178844
pv
         0     1
[1,] 0.864 0.136


[1] 179


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1281496 0.04130221 0.1232668 0.05995345 0.2183564
pv
         0     1
[1,] 0.858 0.142


[1] 180


 sim range :
 -----
ev
          mean        sd      50%       2.5%     97.5%
[1,] 0.1280931 0.0414901 0.123159 0.05968673 0.2188806
pv
         0     1
[1,] 0.855 0.145


[1] 181


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1280368 0.04167814 0.1230391 0.05942112 0.2194058
pv
         0     1
[1,] 0.873 0.127


[1] 182


 sim range :
 -----
ev
          mean         sd     50%       2.5%     97.5%
[1,] 0.1279809 0.04186634 0.12288 0.05915662 0.2199319
pv
         0     1
[1,] 0.887 0.113


[1] 183


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1279253 0.04205468 0.1227401 0.05889322 0.2204589
pv
         0     1
[1,] 0.871 0.129


[1] 184


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1278701 0.04224316 0.1225419 0.05863092 0.2209868
pv
         0     1
[1,] 0.864 0.136


[1] 185


 sim range :
 -----
ev
          mean         sd      50%       2.5%     97.5%
[1,] 0.1278152 0.04243179 0.122475 0.05836971 0.2215109
pv
         0     1
[1,] 0.873 0.127


[1] 186


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1277606 0.04262056 0.1223813 0.0581096 0.2217536
pv
         0     1
[1,] 0.863 0.137


[1] 187


 sim range :
 -----
ev
          mean         sd      50%       2.5%     97.5%
[1,] 0.1277064 0.04280947 0.122169 0.05785058 0.2219523
pv
         0     1
[1,] 0.877 0.123


[1] 188


 sim range :
 -----
ev
          mean         sd      50%       2.5%     97.5%
[1,] 0.1276525 0.04299851 0.121944 0.05759264 0.2223987
pv
        0    1
[1,] 0.87 0.13


[1] 189


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1275989 0.04318769 0.1217861 0.05733578 0.2229499
pv
         0     1
[1,] 0.862 0.138


[1] 190


 sim range :
 -----
ev
          mean         sd       50%    2.5%     97.5%
[1,] 0.1275456 0.04337701 0.1216297 0.05708 0.2235021
pv
         0     1
[1,] 0.867 0.133


[1] 191


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1274927 0.04356646 0.1215155 0.0567797 0.2240463
pv
         0     1
[1,] 0.874 0.126


[1] 192


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1274401 0.04375603 0.1213989 0.05646485 0.2245957
pv
        0    1
[1,] 0.89 0.11


[1] 193


 sim range :
 -----
ev
          mean         sd      50%       2.5%     97.5%
[1,] 0.1273878 0.04394574 0.121237 0.05615165 0.2251545
pv
         0     1
[1,] 0.869 0.131


[1] 194


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1273359 0.04413557 0.1210437 0.05584008 0.2257143
pv
         0     1
[1,] 0.884 0.116


[1] 195


 sim range :
 -----
ev
          mean         sd      50%       2.5%     97.5%
[1,] 0.1272842 0.04432553 0.120953 0.05553014 0.2262751
pv
         0     1
[1,] 0.874 0.126


[1] 196


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1272329 0.04451562 0.1208983 0.05522202 0.2268356
pv
         0     1
[1,] 0.867 0.133


[1] 197


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1271819 0.04470582 0.1208058 0.05491614 0.2273945
pv
         0     1
[1,] 0.891 0.109


[1] 198


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1271313 0.04489615 0.1207272 0.05461185 0.2279545
pv
         0     1
[1,] 0.873 0.127


[1] 199


 sim range :
 -----
ev
          mean        sd       50%       2.5%     97.5%
[1,] 0.1270809 0.0450866 0.1206487 0.05436659 0.2285154
pv
         0     1
[1,] 0.867 0.133


[1] 200


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1270309 0.04527717 0.1205845 0.05413258 0.2290773
pv
         0     1
[1,] 0.861 0.139


[1] 201


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1269812 0.04546785 0.1205381 0.05390032 0.2296402
pv
         0     1
[1,] 0.878 0.122


[1] 202


 sim range :
 -----
ev
          mean         sd       50%       2.5%    97.5%
[1,] 0.1269318 0.04565865 0.1205009 0.05366916 0.230204
pv
         0     1
[1,] 0.858 0.142


[1] 203


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1268828 0.04584957 0.1204409 0.05343893 0.2307688
pv
         0     1
[1,] 0.886 0.114


[1] 204


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1268341 0.04604059 0.1203445 0.05320964 0.2313346
pv
         0     1
[1,] 0.867 0.133


[1] 205


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1267856 0.04623174 0.1202474 0.05298128 0.2319013
pv
         0     1
[1,] 0.874 0.126


[1] 206


 sim range :
 -----
ev
          mean         sd      50%       2.5%     97.5%
[1,] 0.1267375 0.04642299 0.120186 0.05275384 0.2324688
pv
         0     1
[1,] 0.854 0.146


[1] 207


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1266897 0.04661435 0.1201246 0.05252732 0.2330371
pv
         0     1
[1,] 0.875 0.125


[1] 208


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1266423 0.04680582 0.1200634 0.05230173 0.2336065
pv
         0     1
[1,] 0.868 0.132


[1] 209


 sim range :
 -----
ev
          mean        sd      50%       2.5%     97.5%
[1,] 0.1265951 0.0469974 0.119887 0.05207705 0.2341768
pv
        0    1
[1,] 0.87 0.13


[1] 210


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1265483 0.04718908 0.1197476 0.05185328 0.2347481
pv
         0     1
[1,] 0.872 0.128


[1] 211


 sim range :
 -----
ev
          mean         sd      50%       2.5%     97.5%
[1,] 0.1265017 0.04738087 0.119681 0.05163042 0.2353203
pv
         0     1
[1,] 0.863 0.137


[1] 212


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1264555 0.04757277 0.1196369 0.05140905 0.2358935
pv
         0     1
[1,] 0.875 0.125


[1] 213


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1264096 0.04776476 0.1195984 0.05118956 0.2364677
pv
        0    1
[1,] 0.87 0.13


[1] 214


 sim range :
 -----
ev
         mean         sd       50%       2.5%     97.5%
[1,] 0.126364 0.04795686 0.1195277 0.05097097 0.2370428
pv
         0     1
[1,] 0.873 0.127


[1] 215


 sim range :
 -----
ev
          mean         sd      50%       2.5%     97.5%
[1,] 0.1263188 0.04814906 0.119412 0.05071376 0.2376189
pv
        0    1
[1,] 0.89 0.11


[1] 216


 sim range :
 -----
ev
          mean         sd       50%      2.5%    97.5%
[1,] 0.1262738 0.04834137 0.1192599 0.0504481 0.238196
pv
         0     1
[1,] 0.875 0.125


[1] 217


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1262291 0.04853377 0.1191022 0.05018377 0.2387691
pv
         0     1
[1,] 0.862 0.138


[1] 218


 sim range :
 -----
ev
          mean         sd   50%       2.5%     97.5%
[1,] 0.1261848 0.04872627 0.119 0.04992075 0.2393161
pv
         0     1
[1,] 0.877 0.123


[1] 219


 sim range :
 -----
ev
          mean         sd       50%       2.5%    97.5%
[1,] 0.1261407 0.04891886 0.1188801 0.04965904 0.239864
pv
         0     1
[1,] 0.864 0.136


[1] 220


 sim range :
 -----
ev
         mean         sd       50%       2.5%     97.5%
[1,] 0.126097 0.04911156 0.1187211 0.04947034 0.2404128
pv
         0     1
[1,] 0.876 0.124


[1] 221


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1260536 0.04930435 0.1185292 0.04928727 0.2409624
pv
         0     1
[1,] 0.875 0.125


[1] 222


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1260105 0.04949723 0.1183855 0.04910485 0.2415129
pv
         0     1
[1,] 0.879 0.121


[1] 223


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1259676 0.04969021 0.1182754 0.04892307 0.2420642
pv
         0     1
[1,] 0.876 0.124


[1] 224


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1259251 0.04988328 0.1181654 0.04874194 0.2426162
pv
         0     1
[1,] 0.859 0.141


[1] 225


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1258829 0.05007645 0.1180555 0.04856144 0.2431687
pv
         0     1
[1,] 0.886 0.114


[1] 226


 sim range :
 -----
ev
         mean        sd       50%       2.5%     97.5%
[1,] 0.125841 0.0502697 0.1179457 0.04838158 0.2437221
pv
         0     1
[1,] 0.883 0.117


[1] 227


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1257994 0.05046305 0.1178695 0.04819301 0.2442734
pv
         0     1
[1,] 0.885 0.115


[1] 228


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1257581 0.05065649 0.1177855 0.04798463 0.2448253
pv
         0     1
[1,] 0.894 0.106


[1] 229


 sim range :
 -----
ev
          mean         sd       50%       2.5%   97.5%
[1,] 0.1257171 0.05085001 0.1176393 0.04777793 0.24531
pv
         0     1
[1,] 0.869 0.131


[1] 230


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1256764 0.05104363 0.1174941 0.04757223 0.2457525
pv
         0     1
[1,] 0.863 0.137


[1] 231


 sim range :
 -----
ev
          mean         sd      50%       2.5%     97.5%
[1,] 0.1256361 0.05123733 0.117349 0.04736737 0.2461955
pv
         0     1
[1,] 0.877 0.123


[1] 232


 sim range :
 -----
ev
         mean         sd       50%       2.5%     97.5%
[1,] 0.125596 0.05143112 0.1172042 0.04716335 0.2466391
pv
         0     1
[1,] 0.878 0.122


[1] 233


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1255562 0.05162499 0.1170595 0.04696017 0.2470832
pv
         0     1
[1,] 0.886 0.114


[1] 234


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1255167 0.05181896 0.1169149 0.04675781 0.2475279
pv
         0     1
[1,] 0.886 0.114


[1] 235


 sim range :
 -----
ev
          mean       sd       50%       2.5%     97.5%
[1,] 0.1254775 0.052013 0.1167964 0.04655629 0.2479699
pv
         0     1
[1,] 0.873 0.127


[1] 236


 sim range :
 -----
ev
          mean         sd       50%       2.5%   97.5%
[1,] 0.1254386 0.05220713 0.1167044 0.04635559 0.24841
pv
         0     1
[1,] 0.877 0.123


[1] 237


 sim range :
 -----
ev
       mean         sd       50%       2.5%     97.5%
[1,] 0.1254 0.05240135 0.1166126 0.04615572 0.2488506
pv
         0     1
[1,] 0.884 0.116


[1] 238


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1253617 0.05259564 0.1165207 0.04595666 0.2492917
pv
         0     1
[1,] 0.855 0.145


[1] 239


 sim range :
 -----
ev
          mean         sd      50%       2.5%     97.5%
[1,] 0.1253237 0.05279002 0.116429 0.04575842 0.2497334
pv
         0     1
[1,] 0.878 0.122


[1] 240


 sim range :
 -----
ev
          mean         sd       50%     2.5%     97.5%
[1,] 0.1252859 0.05298448 0.1163373 0.045561 0.2501756
pv
         0     1
[1,] 0.846 0.154


[1] 241


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1252485 0.05317902 0.1162193 0.04536439 0.2506183
pv
        0    1
[1,] 0.88 0.12


[1] 242


 sim range :
 -----
ev
          mean         sd      50%       2.5%     97.5%
[1,] 0.1252114 0.05337365 0.116113 0.04516859 0.2509668
pv
         0     1
[1,] 0.875 0.125


[1] 243


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1251745 0.05356835 0.1160315 0.04497359 0.2513081
pv
         0     1
[1,] 0.871 0.129


[1] 244


 sim range :
 -----
ev
         mean         sd     50%       2.5%     97.5%
[1,] 0.125138 0.05376313 0.11595 0.04477543 0.2518024
pv
         0     1
[1,] 0.882 0.118


[1] 245


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1251017 0.05395799 0.1158982 0.04457352 0.2524061
pv
         0     1
[1,] 0.889 0.111


[1] 246


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1250658 0.05415292 0.1158228 0.04437248 0.2528593
pv
         0     1
[1,] 0.858 0.142


[1] 247


 sim range :
 -----
ev
          mean         sd       50%      2.5%    97.5%
[1,] 0.1250301 0.05434794 0.1157441 0.0441723 0.253313
pv
         0     1
[1,] 0.875 0.125


[1] 248


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1249947 0.05454303 0.1156373 0.04397298 0.2537672
pv
        0    1
[1,] 0.88 0.12


[1] 249


 sim range :
 -----
ev
          mean         sd       50%       2.5%    97.5%
[1,] 0.1249596 0.05473819 0.1155305 0.04377453 0.254222
pv
         0     1
[1,] 0.841 0.159


[1] 250


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1249248 0.05493344 0.1154382 0.04357692 0.2546773
pv
         0     1
[1,] 0.874 0.126


[1] 251


 sim range :
 -----
ev
          mean         sd    50%       2.5%     97.5%
[1,] 0.1248903 0.05512875 0.1154 0.04338017 0.2551332
pv
         0     1
[1,] 0.884 0.116


[1] 252


 sim range :
 -----
ev
         mean         sd       50%       2.5%     97.5%
[1,] 0.124856 0.05532414 0.1153651 0.04318427 0.2555896
pv
         0     1
[1,] 0.866 0.134


[1] 253


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1248221 0.05551961 0.1152766 0.04298921 0.2560465
pv
        0    1
[1,] 0.88 0.12


[1] 254


 sim range :
 -----
ev
          mean         sd       50%     2.5%    97.5%
[1,] 0.1247884 0.05571514 0.1151941 0.042795 0.256504
pv
        0    1
[1,] 0.87 0.13


[1] 255


 sim range :
 -----
ev
         mean         sd       50%       2.5%    97.5%
[1,] 0.124755 0.05591075 0.1150625 0.04260162 0.256962
pv
        0    1
[1,] 0.88 0.12


[1] 256


 sim range :
 -----
ev
         mean         sd       50%       2.5%     97.5%
[1,] 0.124722 0.05610644 0.1149311 0.04240908 0.2574204
pv
         0     1
[1,] 0.877 0.123


[1] 257


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1246891 0.05630219 0.1147998 0.04221737 0.2578741
pv
         0     1
[1,] 0.886 0.114


[1] 258


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1246566 0.05649801 0.1146982 0.04202649 0.2583283
pv
         0     1
[1,] 0.876 0.124


[1] 259


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1246244 0.05669391 0.1145976 0.04183643 0.2587831
pv
         0     1
[1,] 0.874 0.126


[1] 260


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1245924 0.05688987 0.1144969 0.04164719 0.2592384
pv
         0     1
[1,] 0.865 0.135


[1] 261


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1245607 0.05708591 0.1144258 0.04145949 0.2596942
pv
         0     1
[1,] 0.885 0.115


[1] 262


 sim range :
 -----
ev
          mean         sd      50%       2.5%     97.5%
[1,] 0.1245293 0.05728201 0.114387 0.04127263 0.2601505
pv
         0     1
[1,] 0.879 0.121


[1] 263


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1244982 0.05747818 0.1143093 0.04108657 0.2606073
pv
         0     1
[1,] 0.873 0.127


[1] 264


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1244673 0.05767442 0.1142026 0.04090132 0.2610646
pv
         0     1
[1,] 0.879 0.121


[1] 265


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1244367 0.05787072 0.1141264 0.04071687 0.2615225
pv
         0     1
[1,] 0.892 0.108


[1] 266


 sim range :
 -----
ev
          mean        sd       50%       2.5%     97.5%
[1,] 0.1244064 0.0580671 0.1140523 0.04053321 0.2619809
pv
         0     1
[1,] 0.868 0.132


[1] 267


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1243764 0.05826353 0.1139594 0.04035035 0.2624398
pv
        0    1
[1,] 0.89 0.11


[1] 268


 sim range :
 -----
ev
          mean         sd      50%       2.5%    97.5%
[1,] 0.1243467 0.05846004 0.113863 0.04016828 0.262929
pv
         0     1
[1,] 0.874 0.126


[1] 269


 sim range :
 -----
ev
          mean         sd       50%     2.5%     97.5%
[1,] 0.1243172 0.05865661 0.1137189 0.039987 0.2634977
pv
         0     1
[1,] 0.877 0.123


[1] 270


 sim range :
 -----
ev
         mean         sd     50%      2.5%     97.5%
[1,] 0.124288 0.05885324 0.11357 0.0398065 0.2640672
pv
         0     1
[1,] 0.878 0.122


[1] 271


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1242591 0.05904994 0.1134634 0.03962679 0.2646374
pv
         0     1
[1,] 0.886 0.114


[1] 272


 sim range :
 -----
ev
          mean        sd       50%       2.5%     97.5%
[1,] 0.1242305 0.0592467 0.1134113 0.03944076 0.2652085
pv
         0     1
[1,] 0.892 0.108


[1] 273


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1242021 0.05944352 0.1133317 0.03922152 0.2657803
pv
         0     1
[1,] 0.885 0.115


[1] 274


 sim range :
 -----
ev
         mean        sd       50%       2.5%     97.5%
[1,] 0.124174 0.0596404 0.1132428 0.03900344 0.2663529
pv
         0     1
[1,] 0.901 0.099


[1] 275


 sim range :
 -----
ev
          mean         sd      50%       2.5%     97.5%
[1,] 0.1241462 0.05983735 0.113154 0.03878654 0.2669263
pv
         0     1
[1,] 0.872 0.128


[1] 276


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1241186 0.06003436 0.1130653 0.03858978 0.2675005
pv
         0     1
[1,] 0.891 0.109


[1] 277


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1240913 0.06023143 0.1129991 0.03844208 0.2680755
pv
         0     1
[1,] 0.886 0.114


[1] 278


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1240643 0.06042856 0.1129604 0.03829494 0.2686513
pv
         0     1
[1,] 0.862 0.138


[1] 279


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1240375 0.06062575 0.1128828 0.03814834 0.2692278
pv
         0     1
[1,] 0.876 0.124


[1] 280


 sim range :
 -----
ev
         mean       sd       50%       2.5%     97.5%
[1,] 0.124011 0.060823 0.1127715 0.03800228 0.2698051
pv
         0     1
[1,] 0.879 0.121


[1] 281


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1239848 0.06102031 0.1126721 0.03785676 0.2703832
pv
         0     1
[1,] 0.863 0.137


[1] 282


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1239589 0.06121768 0.1125402 0.0376898 0.2709621
pv
         0     1
[1,] 0.868 0.132


[1] 283


 sim range :
 -----
ev
          mean        sd       50%       2.5%     97.5%
[1,] 0.1239332 0.0614151 0.1124028 0.03751756 0.2715418
pv
         0     1
[1,] 0.872 0.128


[1] 284


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1239078 0.06161258 0.1122807 0.03734608 0.2721222
pv
         0     1
[1,] 0.874 0.126


[1] 285


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1238826 0.06181012 0.1121662 0.03717535 0.2727034
pv
         0     1
[1,] 0.876 0.124


[1] 286


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1238577 0.06200772 0.1120535 0.03700538 0.2732854
pv
         0     1
[1,] 0.876 0.124


[1] 287


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1238331 0.06220537 0.1119408 0.03683709 0.2737898
pv
         0     1
[1,] 0.883 0.117


[1] 288


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1238087 0.06240307 0.1118283 0.0366696 0.2741349
pv
         0     1
[1,] 0.867 0.133


[1] 289


 sim range :
 -----
ev
          mean         sd       50%       2.5%   97.5%
[1,] 0.1237846 0.06260083 0.1117515 0.03650285 0.27448
pv
         0     1
[1,] 0.884 0.116


[1] 290


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1237608 0.06279865 0.1116258 0.03633683 0.2748254
pv
         0     1
[1,] 0.897 0.103


[1] 291


 sim range :
 -----
ev
          mean         sd       50%       2.5%    97.5%
[1,] 0.1237372 0.06299652 0.1115262 0.03617153 0.275171
pv
         0     1
[1,] 0.877 0.123


[1] 292


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1237139 0.06319444 0.1114126 0.03600695 0.2755169
pv
        0    1
[1,] 0.88 0.12


[1] 293


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1236908 0.06339241 0.1112992 0.0358431 0.2760474
pv
         0     1
[1,] 0.891 0.109


[1] 294


 sim range :
 -----
ev
         mean         sd       50%       2.5%     97.5%
[1,] 0.123668 0.06359044 0.1111859 0.03567997 0.2766412
pv
         0     1
[1,] 0.865 0.135


[1] 295


 sim range :
 -----
ev
          mean         sd      50%       2.5%     97.5%
[1,] 0.1236455 0.06378852 0.111071 0.03551755 0.2772359
pv
         0     1
[1,] 0.883 0.117


[1] 296


 sim range :
 -----
ev
          mean         sd      50%       2.5%     97.5%
[1,] 0.1236232 0.06398665 0.110937 0.03535584 0.2778313
pv
         0     1
[1,] 0.877 0.123


[1] 297


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1236012 0.06418483 0.1108033 0.03519484 0.2784275
pv
         0     1
[1,] 0.867 0.133


[1] 298


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1235794 0.06438306 0.1106696 0.03503455 0.2790245
pv
         0     1
[1,] 0.864 0.136


[1] 299


 sim range :
 -----
ev
          mean         sd       50%       2.5%    97.5%
[1,] 0.1235579 0.06458134 0.1105361 0.03485331 0.279505
pv
         0     1
[1,] 0.881 0.119


[1] 300


 sim range :
 -----
ev
          mean         sd      50%       2.5%     97.5%
[1,] 0.1235366 0.06477967 0.110407 0.03464535 0.2799798
pv
         0     1
[1,] 0.857 0.143


[1] 301


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1235156 0.06497805 0.1102607 0.03443859 0.2804551
pv
         0     1
[1,] 0.873 0.127


[1] 302


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1234949 0.06517648 0.1101531 0.03423302 0.2809309
pv
         0     1
[1,] 0.887 0.113


[1] 303


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1234744 0.06537495 0.1100067 0.03402864 0.2814072
pv
         0     1
[1,] 0.879 0.121


[1] 304


 sim range :
 -----
ev
          mean         sd       50%       2.5%    97.5%
[1,] 0.1234541 0.06557347 0.1098924 0.03382543 0.281884
pv
         0     1
[1,] 0.874 0.126


[1] 305


 sim range :
 -----
ev
          mean         sd     50%      2.5%     97.5%
[1,] 0.1234341 0.06577204 0.10978 0.0336234 0.2823612
pv
         0     1
[1,] 0.879 0.121


[1] 306


 sim range :
 -----
ev
          mean         sd       50%       2.5%    97.5%
[1,] 0.1234144 0.06597065 0.1096858 0.03342253 0.282839
pv
         0     1
[1,] 0.855 0.145


[1] 307


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1233949 0.06616931 0.1095711 0.03322283 0.2833172
pv
         0     1
[1,] 0.888 0.112


[1] 308


 sim range :
 -----
ev
          mean         sd       50%       2.5%    97.5%
[1,] 0.1233757 0.06636801 0.1094631 0.03302428 0.283796
pv
         0     1
[1,] 0.891 0.109


[1] 309


 sim range :
 -----
ev
          mean         sd      50%       2.5%     97.5%
[1,] 0.1233567 0.06656676 0.109356 0.03284336 0.2842752
pv
         0     1
[1,] 0.862 0.138


[1] 310


 sim range :
 -----
ev
          mean         sd      50%       2.5%     97.5%
[1,] 0.1233379 0.06676556 0.109249 0.03269319 0.2847549
pv
         0     1
[1,] 0.857 0.143


[1] 311


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1233194 0.06696439 0.1091442 0.03254369 0.2852351
pv
         0     1
[1,] 0.875 0.125


[1] 312


 sim range :
 -----
ev
          mean         sd      50%       2.5%     97.5%
[1,] 0.1233012 0.06716327 0.109055 0.03239485 0.2857158
pv
        0    1
[1,] 0.89 0.11


[1] 313


 sim range :
 -----
ev
          mean         sd      50%       2.5%     97.5%
[1,] 0.1232832 0.06736219 0.108954 0.03224668 0.2861939
pv
         0     1
[1,] 0.869 0.131


[1] 314


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1232654 0.06756115 0.1088309 0.03209915 0.2866709
pv
         0     1
[1,] 0.893 0.107


[1] 315


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1232479 0.06776016 0.1087079 0.03195229 0.2871816
pv
         0     1
[1,] 0.859 0.141


[1] 316


 sim range :
 -----
ev
          mean        sd       50%       2.5%     97.5%
[1,] 0.1232306 0.0679592 0.1086012 0.03180607 0.2877323
pv
         0     1
[1,] 0.877 0.123


[1] 317


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1232136 0.06815829 0.1085181 0.03166051 0.2882835
pv
        0    1
[1,] 0.88 0.12


[1] 318


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1231968 0.06835741 0.1084607 0.03151559 0.2888355
pv
         0     1
[1,] 0.868 0.132


[1] 319


 sim range :
 -----
ev
          mean         sd       50%       2.5%    97.5%
[1,] 0.1231803 0.06855658 0.1084064 0.03137131 0.289388
pv
         0     1
[1,] 0.866 0.134


[1] 320


 sim range :
 -----
ev
         mean         sd       50%       2.5%     97.5%
[1,] 0.123164 0.06875578 0.1083431 0.03122767 0.2899412
pv
         0     1
[1,] 0.873 0.127


[1] 321


 sim range :
 -----
ev
          mean         sd       50%       2.5%    97.5%
[1,] 0.1231479 0.06895502 0.1082544 0.03108467 0.290495
pv
         0     1
[1,] 0.869 0.131


[1] 322


 sim range :
 -----
ev
          mean        sd       50%       2.5%     97.5%
[1,] 0.1231321 0.0691543 0.1081658 0.03094231 0.2910494
pv
         0     1
[1,] 0.879 0.121


[1] 323


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1231166 0.06935362 0.1080772 0.03080058 0.2916044
pv
         0     1
[1,] 0.878 0.122


[1] 324


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1231012 0.06955297 0.1079887 0.03065948 0.2921601
pv
         0     1
[1,] 0.897 0.103


[1] 325


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1230861 0.06975236 0.1079003 0.03051901 0.2927164
pv
         0     1
[1,] 0.875 0.125


[1] 326


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1230713 0.06995178 0.1078119 0.03037916 0.2932733
pv
         0     1
[1,] 0.859 0.141


[1] 327


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1230567 0.07015124 0.1077236 0.03023994 0.2938308
pv
         0     1
[1,] 0.871 0.129


[1] 328


 sim range :
 -----
ev
          mean         sd       50%       2.5%    97.5%
[1,] 0.1230423 0.07035073 0.1076354 0.03010134 0.294389
pv
         0     1
[1,] 0.877 0.123


[1] 329


 sim range :
 -----
ev
          mean         sd       50%       2.5%    97.5%
[1,] 0.1230281 0.07055026 0.1075472 0.02996335 0.294839
pv
        0    1
[1,] 0.88 0.12


[1] 330


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1230142 0.07074981 0.1074591 0.02982598 0.2952729
pv
         0     1
[1,] 0.865 0.135


[1] 331


 sim range :
 -----
ev
          mean         sd      50%       2.5%     97.5%
[1,] 0.1230006 0.07094941 0.107371 0.02968922 0.2957071
pv
         0     1
[1,] 0.869 0.131


[1] 332


 sim range :
 -----
ev
          mean         sd      50%       2.5%     97.5%
[1,] 0.1229871 0.07114903 0.107283 0.02955306 0.2961416
pv
         0     1
[1,] 0.865 0.135


[1] 333


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1229739 0.07134869 0.1071961 0.02941825 0.2965766
pv
         0     1
[1,] 0.899 0.101


[1] 334


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1229609 0.07154837 0.1071139 0.02928478 0.2970119
pv
        0    1
[1,] 0.87 0.13


[1] 335


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1229482 0.07174809 0.1070319 0.02915189 0.2974476
pv
         0     1
[1,] 0.877 0.123


[1] 336


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1229357 0.07194784 0.1069498 0.02901959 0.2978836
pv
         0     1
[1,] 0.866 0.134


[1] 337


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1229234 0.07214761 0.1068386 0.02888787 0.2983201
pv
         0     1
[1,] 0.883 0.117


[1] 338


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1229114 0.07234742 0.1067536 0.02875673 0.2987569
pv
         0     1
[1,] 0.873 0.127


[1] 339


 sim range :
 -----
ev
          mean         sd       50%       2.5%    97.5%
[1,] 0.1228996 0.07254725 0.1066875 0.02862617 0.299194
pv
         0     1
[1,] 0.862 0.138


[1] 340


 sim range :
 -----
ev
         mean         sd       50%       2.5%     97.5%
[1,] 0.122888 0.07274712 0.1066096 0.02849618 0.2996315
pv
        0    1
[1,] 0.89 0.11


[1] 341


 sim range :
 -----
ev
          mean       sd       50%       2.5%     97.5%
[1,] 0.1228766 0.072947 0.1065194 0.02836677 0.3000694
pv
         0     1
[1,] 0.881 0.119


[1] 342


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1228655 0.07314692 0.1064292 0.02823792 0.3005077
pv
         0     1
[1,] 0.863 0.137


[1] 343


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1228546 0.07334686 0.1063391 0.02810965 0.3009463
pv
         0     1
[1,] 0.874 0.126


[1] 344


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1228439 0.07354683 0.1062491 0.02798194 0.3013853
pv
         0     1
[1,] 0.883 0.117


[1] 345


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1228335 0.07374682 0.1061592 0.0278548 0.3018247
pv
         0     1
[1,] 0.897 0.103


[1] 346


 sim range :
 -----
ev
          mean         sd      50%       2.5%     97.5%
[1,] 0.1228233 0.07394683 0.106053 0.02772822 0.3022644
pv
         0     1
[1,] 0.878 0.122


[1] 347


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1228133 0.07414687 0.1059735 0.02760219 0.3027892
pv
         0     1
[1,] 0.894 0.106


[1] 348


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1228035 0.07434694 0.1058898 0.02747673 0.3033531
pv
         0     1
[1,] 0.886 0.114


[1] 349


 sim range :
 -----
ev
         mean         sd       50%       2.5%     97.5%
[1,] 0.122794 0.07454702 0.1057783 0.02735181 0.3039175
pv
        0    1
[1,] 0.88 0.12


[1] 350


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1227847 0.07474713 0.1056534 0.02722746 0.3044825
pv
         0     1
[1,] 0.876 0.124


[1] 351


 sim range :
 -----
ev
          mean         sd      50%       2.5%     97.5%
[1,] 0.1227756 0.07494726 0.105525 0.02710365 0.3050482
pv
         0     1
[1,] 0.882 0.118


[1] 352


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1227667 0.07514741 0.1054234 0.02698038 0.3056144
pv
         0     1
[1,] 0.876 0.124


[1] 353


 sim range :
 -----
ev
         mean         sd       50%       2.5%     97.5%
[1,] 0.122758 0.07534758 0.1053285 0.02685767 0.3061812
pv
         0     1
[1,] 0.862 0.138


[1] 354


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1227496 0.07554777 0.1052336 0.02672965 0.3067485
pv
         0     1
[1,] 0.874 0.126


[1] 355


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1227414 0.07574798 0.1051388 0.02658742 0.3073165
pv
         0     1
[1,] 0.887 0.113


[1] 356


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1227334 0.07594821 0.1050441 0.02644593 0.3078851
pv
         0     1
[1,] 0.874 0.126


[1] 357


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1227256 0.07614846 0.1049494 0.02630517 0.3084542
pv
        0    1
[1,] 0.87 0.13


[1] 358


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1227181 0.07634872 0.1048548 0.02616514 0.3090239
pv
         0     1
[1,] 0.872 0.128


[1] 359


 sim range :
 -----
ev
          mean       sd       50%       2.5%     97.5%
[1,] 0.1227107 0.076549 0.1047603 0.02602584 0.3095916
pv
         0     1
[1,] 0.887 0.113


[1] 360


 sim range :
 -----
ev
          mean        sd       50%       2.5%     97.5%
[1,] 0.1227036 0.0767493 0.1046659 0.02588726 0.3101598
pv
         0     1
[1,] 0.874 0.126


[1] 361


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1226967 0.07694961 0.1045715 0.0257494 0.3107286
pv
         0     1
[1,] 0.884 0.116


[1] 362


 sim range :
 -----
ev
        mean         sd       50%       2.5%     97.5%
[1,] 0.12269 0.07714994 0.1044773 0.02561226 0.3112979
pv
         0     1
[1,] 0.883 0.117


[1] 363


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1226836 0.07735029 0.1043831 0.02547583 0.3118678
pv
         0     1
[1,] 0.876 0.124


[1] 364


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1226773 0.07755064 0.1042889 0.0253401 0.3124383
pv
         0     1
[1,] 0.871 0.129


[1] 365


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1226713 0.07775101 0.1041949 0.02520509 0.3130094
pv
         0     1
[1,] 0.883 0.117


[1] 366


 sim range :
 -----
ev
          mean        sd       50%       2.5%    97.5%
[1,] 0.1226655 0.0779514 0.1041009 0.02507077 0.313581
pv
         0     1
[1,] 0.887 0.113


[1] 367


 sim range :
 -----
ev
          mean        sd      50%       2.5%     97.5%
[1,] 0.1226599 0.0781518 0.104007 0.02493716 0.3141532
pv
         0     1
[1,] 0.887 0.113


[1] 368


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1226545 0.07835221 0.1039132 0.02480424 0.3147259
pv
         0     1
[1,] 0.905 0.095


[1] 369


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1226493 0.07855263 0.1038194 0.02467201 0.3152993
pv
         0     1
[1,] 0.887 0.113


[1] 370


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1226443 0.07875306 0.1037063 0.02454047 0.3158731
pv
        0    1
[1,] 0.89 0.11


[1] 371


 sim range :
 -----
ev
          mean        sd       50%       2.5%     97.5%
[1,] 0.1226395 0.0789535 0.1035884 0.02441968 0.3164476
pv
         0     1
[1,] 0.874 0.126


[1] 372


 sim range :
 -----
ev
         mean         sd       50%       2.5%     97.5%
[1,] 0.122635 0.07915395 0.1034706 0.02432133 0.3170226
pv
         0     1
[1,] 0.909 0.091


[1] 373


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1226306 0.07935441 0.1033584 0.02422337 0.3175981
pv
         0     1
[1,] 0.866 0.134


[1] 374


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1226265 0.07955488 0.1032831 0.0241258 0.3181743
pv
         0     1
[1,] 0.878 0.122


[1] 375


 sim range :
 -----
ev
          mean         sd       50%      2.5%    97.5%
[1,] 0.1226226 0.07975535 0.1032078 0.0240286 0.318751
pv
         0     1
[1,] 0.881 0.119


[1] 376


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1226189 0.07995584 0.1031327 0.0239318 0.3193282
pv
         0     1
[1,] 0.886 0.114


[1] 377


 sim range :
 -----
ev
          mean         sd       50%       2.5%    97.5%
[1,] 0.1226153 0.08015633 0.1030575 0.02383537 0.319906
pv
         0     1
[1,] 0.883 0.117


[1] 378


 sim range :
 -----
ev
         mean         sd      50%       2.5%     97.5%
[1,] 0.122612 0.08035683 0.102979 0.02373932 0.3204843
pv
        0    1
[1,] 0.89 0.11


[1] 379


 sim range :
 -----
ev
          mean         sd      50%       2.5%     97.5%
[1,] 0.1226089 0.08055733 0.102886 0.02364366 0.3210632
pv
        0    1
[1,] 0.89 0.11


[1] 380


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1226061 0.08075784 0.1027802 0.02354837 0.3216427
pv
         0     1
[1,] 0.877 0.123


[1] 381


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1226034 0.08095835 0.1026684 0.02345345 0.3222227
pv
         0     1
[1,] 0.888 0.112


[1] 382


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1226009 0.08115887 0.1025567 0.02335891 0.3228032
pv
         0     1
[1,] 0.881 0.119


[1] 383


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1225986 0.08135939 0.1024921 0.02326475 0.3233843
pv
         0     1
[1,] 0.875 0.125


[1] 384


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1225965 0.08155991 0.1024221 0.02317095 0.3239468
pv
         0     1
[1,] 0.871 0.129


[1] 385


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1225947 0.08176044 0.1023295 0.02307753 0.3244324
pv
         0     1
[1,] 0.882 0.118


[1] 386


 sim range :
 -----
ev
         mean         sd       50%       2.5%     97.5%
[1,] 0.122593 0.08196096 0.1022371 0.02298447 0.3249183
pv
         0     1
[1,] 0.884 0.116


[1] 387


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1225915 0.08216149 0.1021446 0.02289178 0.3254047
pv
         0     1
[1,] 0.904 0.096


[1] 388


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1225903 0.08236202 0.1020523 0.02279946 0.3258889
pv
         0     1
[1,] 0.887 0.113


[1] 389


 sim range :
 -----
ev
          mean         sd     50%      2.5%    97.5%
[1,] 0.1225892 0.08256255 0.10196 0.0227075 0.326372
pv
         0     1
[1,] 0.883 0.117


[1] 390


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1225883 0.08276308 0.1018679 0.02261591 0.3268555
pv
         0     1
[1,] 0.876 0.124


[1] 391


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1225877 0.08296361 0.1017757 0.02252468 0.3273394
pv
         0     1
[1,] 0.863 0.137


[1] 392


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1225872 0.08316413 0.1016871 0.02243381 0.3278236
pv
        0    1
[1,] 0.87 0.13


[1] 393


 sim range :
 -----
ev
          mean         sd      50%       2.5%     97.5%
[1,] 0.1225869 0.08336466 0.101584 0.02234329 0.3283082
pv
         0     1
[1,] 0.872 0.128


[1] 394


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1225869 0.08356518 0.1014998 0.02224388 0.3287827
pv
         0     1
[1,] 0.871 0.129


[1] 395


 sim range :
 -----
ev
         mean         sd      50%       2.5%    97.5%
[1,] 0.122587 0.08376569 0.101408 0.02213731 0.329214
pv
         0     1
[1,] 0.886 0.114


[1] 396


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1225873 0.08396621 0.1013163 0.02203124 0.3296456
pv
         0     1
[1,] 0.868 0.132


[1] 397


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1225878 0.08416672 0.1012239 0.02192567 0.3300774
pv
         0     1
[1,] 0.864 0.136


[1] 398


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1225885 0.08436722 0.1011127 0.02182059 0.3305095
pv
         0     1
[1,] 0.883 0.117


[1] 399


 sim range :
 -----
ev
          mean         sd       50%     2.5%    97.5%
[1,] 0.1225894 0.08456772 0.1010079 0.021716 0.330942
pv
        0    1
[1,] 0.88 0.12


[1] 400


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1225905 0.08476822 0.1009229 0.02161191 0.3313747
pv
         0     1
[1,] 0.875 0.125


[1] 401


 sim range :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1225918 0.08496871 0.1008575 0.0215083 0.3318076
pv
         0     1
[1,] 0.862 0.138


[1] 402


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1225933 0.08516919 0.1007673 0.02140518 0.3322409
pv
         0     1
[1,] 0.869 0.131


[1] 403


 sim range :
 -----
ev
         mean         sd       50%       2.5%     97.5%
[1,] 0.122595 0.08536966 0.1006761 0.02130254 0.3326745
pv
         0     1
[1,] 0.871 0.129


[1] 404


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1225969 0.08557012 0.1005849 0.02120039 0.3331083
pv
         0     1
[1,] 0.869 0.131


[1] 405


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1225989 0.08577058 0.1004939 0.02109871 0.3335424
pv
         0     1
[1,] 0.891 0.109


[1] 406


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1226012 0.08597103 0.1004119 0.02099751 0.3339768
pv
         0     1
[1,] 0.861 0.139


[1] 407


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1226036 0.08617146 0.1003379 0.02089679 0.3344115
pv
         0     1
[1,] 0.898 0.102


[1] 408


 sim range :
 -----
ev
          mean         sd      50%       2.5%     97.5%
[1,] 0.1226063 0.08637189 0.100264 0.02079654 0.3348465
pv
         0     1
[1,] 0.878 0.122


[1] 409


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1226091 0.08657231 0.1001701 0.02069676 0.3352818
pv
         0     1
[1,] 0.887 0.113


[1] 410


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1226121 0.08677271 0.1000601 0.02059745 0.3357173
pv
         0     1
[1,] 0.884 0.116


[1] 411


 sim range :
 -----
ev
          mean        sd       50%       2.5%     97.5%
[1,] 0.1226153 0.0869731 0.0999502 0.02049861 0.3361531
pv
         0     1
[1,] 0.896 0.104


[1] 412


 sim range :
 -----
ev
          mean         sd        50%       2.5%     97.5%
[1,] 0.1226187 0.08717348 0.09986399 0.02040023 0.3365892
pv
         0     1
[1,] 0.884 0.116


[1] 413


 sim range :
 -----
ev
          mean         sd        50%       2.5%     97.5%
[1,] 0.1226223 0.08737385 0.09979677 0.02029959 0.3370256
pv
        0    1
[1,] 0.91 0.09


[1] 414


 sim range :
 -----
ev
         mean         sd        50%      2.5%     97.5%
[1,] 0.122626 0.08757421 0.09969595 0.0201949 0.3374943
pv
         0     1
[1,] 0.887 0.113


[1] 415


 sim range :
 -----
ev
        mean         sd        50%       2.5%     97.5%
[1,] 0.12263 0.08777455 0.09959521 0.02009073 0.3381891
pv
         0     1
[1,] 0.861 0.139


[1] 416


 sim range :
 -----
ev
          mean         sd        50%       2.5%     97.5%
[1,] 0.1226341 0.08797487 0.09949674 0.01998709 0.3388845
pv
        0    1
[1,] 0.87 0.13


[1] 417


 sim range :
 -----
ev
          mean         sd        50%       2.5%     97.5%
[1,] 0.1226385 0.08817518 0.09940653 0.01988398 0.3395806
pv
         0     1
[1,] 0.875 0.125


[1] 418


 sim range :
 -----
ev
         mean         sd       50%       2.5%     97.5%
[1,] 0.122643 0.08837548 0.0993164 0.01978139 0.3402774
pv
         0     1
[1,] 0.886 0.114


[1] 419


 sim range :
 -----
ev
          mean         sd        50%       2.5%     97.5%
[1,] 0.1226476 0.08857576 0.09922633 0.01967931 0.3409749
pv
         0     1
[1,] 0.882 0.118


[1] 420


 sim range :
 -----
ev
          mean         sd        50%       2.5%     97.5%
[1,] 0.1226525 0.08877602 0.09913634 0.01957775 0.3415689
pv
         0     1
[1,] 0.889 0.111


[1] 421


 sim range :
 -----
ev
          mean         sd        50%       2.5%     97.5%
[1,] 0.1226576 0.08897627 0.09904643 0.01947671 0.3420701
pv
         0     1
[1,] 0.884 0.116


[1] 422


 sim range :
 -----
ev
          mean        sd        50%       2.5%     97.5%
[1,] 0.1226628 0.0891765 0.09895658 0.01937618 0.3425717
pv
         0     1
[1,] 0.895 0.105


[1] 423


 sim range :
 -----
ev
          mean         sd        50%       2.5%     97.5%
[1,] 0.1226682 0.08937671 0.09886681 0.01927615 0.3430737
pv
         0     1
[1,] 0.881 0.119


[1] 424


 sim range :
 -----
ev
          mean        sd        50%       2.5%    97.5%
[1,] 0.1226738 0.0895769 0.09877711 0.01917664 0.343576
pv
         0     1
[1,] 0.845 0.155


[1] 425


 sim range :
 -----
ev
          mean         sd        50%       2.5%     97.5%
[1,] 0.1226796 0.08977707 0.09868748 0.01907762 0.3440786
pv
         0     1
[1,] 0.879 0.121


[1] 426


 sim range :
 -----
ev
          mean         sd        50%       2.5%     97.5%
[1,] 0.1226856 0.08997722 0.09859793 0.01897911 0.3445816
pv
         0     1
[1,] 0.876 0.124


[1] 427


 sim range :
 -----
ev
          mean         sd        50%       2.5%     97.5%
[1,] 0.1226917 0.09017735 0.09850844 0.01888109 0.3450849
pv
         0     1
[1,] 0.877 0.123


[1] 428


 sim range :
 -----
ev
         mean         sd        50%       2.5%    97.5%
[1,] 0.122698 0.09037747 0.09841903 0.01878358 0.345583
pv
         0     1
[1,] 0.868 0.132


[1] 429


 sim range :
 -----
ev
          mean         sd        50%       2.5%    97.5%
[1,] 0.1227045 0.09057756 0.09832969 0.01868655 0.346079
pv
         0     1
[1,] 0.877 0.123


[1] 430


 sim range :
 -----
ev
          mean         sd        50%       2.5%     97.5%
[1,] 0.1227112 0.09077763 0.09824043 0.01859002 0.3465753
pv
        0    1
[1,] 0.88 0.12


[1] 431


 sim range :
 -----
ev
         mean         sd        50%       2.5%    97.5%
[1,] 0.122718 0.09097768 0.09815124 0.01849398 0.347072
pv
         0     1
[1,] 0.898 0.102


[1] 432


 sim range :
 -----
ev
          mean        sd        50%       2.5%     97.5%
[1,] 0.1227251 0.0911777 0.09806211 0.01839842 0.3475689
pv
         0     1
[1,] 0.871 0.129


[1] 433


 sim range :
 -----
ev
          mean        sd        50%       2.5%     97.5%
[1,] 0.1227323 0.0913777 0.09797306 0.01830335 0.3480662
pv
         0     1
[1,] 0.882 0.118


[1] 434


 sim range :
 -----
ev
          mean         sd        50%       2.5%     97.5%
[1,] 0.1227396 0.09157768 0.09788409 0.01820876 0.3486971
pv
        0    1
[1,] 0.87 0.13


[1] 435


 sim range :
 -----
ev
          mean         sd        50%       2.5%     97.5%
[1,] 0.1227472 0.09177764 0.09779518 0.01811466 0.3493365
pv
         0     1
[1,] 0.864 0.136


[1] 436


 sim range :
 -----
ev
          mean         sd        50%       2.5%     97.5%
[1,] 0.1227549 0.09197757 0.09772468 0.01802102 0.3499763
pv
         0     1
[1,] 0.888 0.112


[1] 437


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1227628 0.09217747 0.0976484 0.01792787 0.3505939
pv
        0    1
[1,] 0.87 0.13


[1] 438


 sim range :
 -----
ev
          mean         sd        50%       2.5%    97.5%
[1,] 0.1227709 0.09237735 0.09754942 0.01783519 0.350989
pv
         0     1
[1,] 0.892 0.108


[1] 439


 sim range :
 -----
ev
          mean        sd        50%       2.5%     97.5%
[1,] 0.1227792 0.0925772 0.09745053 0.01774297 0.3513842
pv
         0     1
[1,] 0.895 0.105


[1] 440


 sim range :
 -----
ev
          mean         sd        50%       2.5%     97.5%
[1,] 0.1227876 0.09277703 0.09735174 0.01765123 0.3517797
pv
         0     1
[1,] 0.889 0.111


[1] 441


 sim range :
 -----
ev
          mean         sd        50%       2.5%     97.5%
[1,] 0.1227962 0.09297683 0.09727287 0.01755995 0.3521754
pv
         0     1
[1,] 0.878 0.122


[1] 442


 sim range :
 -----
ev
          mean        sd        50%       2.5%     97.5%
[1,] 0.1228049 0.0931766 0.09719694 0.01746913 0.3525898
pv
         0     1
[1,] 0.873 0.127


[1] 443


 sim range :
 -----
ev
          mean         sd        50%       2.5%     97.5%
[1,] 0.1228139 0.09337635 0.09712107 0.01737878 0.3530969
pv
         0     1
[1,] 0.863 0.137


[1] 444


 sim range :
 -----
ev
         mean         sd        50%       2.5%     97.5%
[1,] 0.122823 0.09357607 0.09704525 0.01728889 0.3536044
pv
         0     1
[1,] 0.868 0.132


[1] 445


 sim range :
 -----
ev
          mean         sd        50%       2.5%     97.5%
[1,] 0.1228323 0.09377575 0.09696948 0.01719945 0.3541121
pv
         0     1
[1,] 0.875 0.125


[1] 446


 sim range :
 -----
ev
          mean         sd        50%       2.5%     97.5%
[1,] 0.1228417 0.09397541 0.09688009 0.01711047 0.3546202
pv
        0    1
[1,] 0.88 0.12


[1] 447


 sim range :
 -----
ev
          mean         sd        50%       2.5%     97.5%
[1,] 0.1228513 0.09417504 0.09678739 0.01702193 0.3551286
pv
         0     1
[1,] 0.883 0.117


[1] 448


 sim range :
 -----
ev
          mean         sd        50%       2.5%     97.5%
[1,] 0.1228611 0.09437464 0.09669478 0.01693385 0.3556374
pv
         0     1
[1,] 0.852 0.148


[1] 449


 sim range :
 -----
ev
          mean         sd        50%       2.5%     97.5%
[1,] 0.1228711 0.09457421 0.09660224 0.01684622 0.3561464
pv
         0     1
[1,] 0.888 0.112


[1] 450


 sim range :
 -----
ev
          mean         sd        50%       2.5%     97.5%
[1,] 0.1228812 0.09477375 0.09650979 0.01675904 0.3566558
pv
         0     1
[1,] 0.871 0.129


[1] 451


 sim range :
 -----
ev
          mean         sd        50%       2.5%     97.5%
[1,] 0.1228915 0.09497326 0.09643261 0.01667229 0.3571656
pv
         0     1
[1,] 0.862 0.138


[1] 452


 sim range :
 -----
ev
         mean         sd        50%       2.5%     97.5%
[1,] 0.122902 0.09517273 0.09635872 0.01658599 0.3576756
pv
         0     1
[1,] 0.866 0.134


[1] 453


 sim range :
 -----
ev
          mean         sd        50%       2.5%    97.5%
[1,] 0.1229126 0.09537217 0.09628488 0.01650013 0.358186
pv
         0     1
[1,] 0.877 0.123


[1] 454


 sim range :
 -----
ev
          mean         sd        50%      2.5%     97.5%
[1,] 0.1229234 0.09557158 0.09621109 0.0164147 0.3586967
pv
         0     1
[1,] 0.892 0.108


[1] 455


 sim range :
 -----
ev
          mean         sd        50%       2.5%     97.5%
[1,] 0.1229343 0.09577096 0.09613736 0.01632971 0.3592077
pv
         0     1
[1,] 0.864 0.136


[1] 456


 sim range :
 -----
ev
          mean        sd        50%       2.5%    97.5%
[1,] 0.1229454 0.0959703 0.09606367 0.01624516 0.359719
pv
         0     1
[1,] 0.892 0.108


[1] 457


 sim range :
 -----
ev
          mean         sd        50%       2.5%     97.5%
[1,] 0.1229567 0.09616961 0.09599004 0.01616103 0.3602307
pv
         0     1
[1,] 0.881 0.119


[1] 458


 sim range :
 -----
ev
          mean         sd        50%       2.5%     97.5%
[1,] 0.1229682 0.09636888 0.09591646 0.01607733 0.3607427
pv
         0     1
[1,] 0.874 0.126


[1] 459


 sim range :
 -----
ev
          mean         sd        50%       2.5%     97.5%
[1,] 0.1229798 0.09656812 0.09584293 0.01599406 0.3612549
pv
         0     1
[1,] 0.876 0.124


[1] 460


 sim range :
 -----
ev
          mean         sd        50%       2.5%     97.5%
[1,] 0.1229915 0.09676733 0.09576946 0.01591121 0.3617675
pv
         0     1
[1,] 0.888 0.112


[1] 461


 sim range :
 -----
ev
          mean         sd        50%       2.5%     97.5%
[1,] 0.1230035 0.09696649 0.09569603 0.01582879 0.3622805
pv
         0     1
[1,] 0.885 0.115


[1] 462


 sim range :
 -----
ev
          mean         sd        50%       2.5%     97.5%
[1,] 0.1230156 0.09716563 0.09562266 0.01574678 0.3627937
pv
         0     1
[1,] 0.872 0.128


[1] 463


 sim range :
 -----
ev
          mean         sd        50%       2.5%     97.5%
[1,] 0.1230278 0.09736472 0.09554934 0.01566519 0.3633073
pv
         0     1
[1,] 0.888 0.112


[1] 464


 sim range :
 -----
ev
          mean         sd        50%       2.5%     97.5%
[1,] 0.1230402 0.09756378 0.09547607 0.01558402 0.3638211
pv
         0     1
[1,] 0.878 0.122


[1] 465


 sim range :
 -----
ev
          mean        sd        50%       2.5%     97.5%
[1,] 0.1230528 0.0977628 0.09540285 0.01550327 0.3643353
pv
         0     1
[1,] 0.865 0.135


[1] 466


 sim range :
 -----
ev
          mean         sd        50%       2.5%     97.5%
[1,] 0.1230655 0.09796178 0.09533108 0.01542292 0.3648498
pv
         0     1
[1,] 0.892 0.108


[1] 467


 sim range :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1230784 0.09816073 0.0952662 0.01534299 0.3653646
pv
         0     1
[1,] 0.875 0.125


[1] 468


 sim range :
 -----
ev
          mean         sd        50%       2.5%     97.5%
[1,] 0.1230915 0.09835963 0.09520136 0.01526346 0.3658797
pv
         0     1
[1,] 0.879 0.121


[1] 469


 sim range :
 -----
ev
          mean        sd        50%       2.5%     97.5%
[1,] 0.1231047 0.0985585 0.09512684 0.01518434 0.3663951
pv
        0    1
[1,] 0.89 0.11


[1] 470


 sim range :
 -----
ev
          mean         sd        50%       2.5%     97.5%
[1,] 0.1231181 0.09875733 0.09504608 0.01510562 0.3669108
pv
         0     1
[1,] 0.864 0.136


[1] 471


 sim range :
 -----
ev
          mean         sd        50%      2.5%     97.5%
[1,] 0.1231316 0.09895612 0.09496046 0.0150273 0.3674268
pv
         0     1
[1,] 0.898 0.102


[1] 472


 sim range :
 -----
ev
          mean         sd        50%       2.5%     97.5%
[1,] 0.1231453 0.09915486 0.09486944 0.01494939 0.3679431
pv
         0     1
[1,] 0.848 0.152


[1] 473


 sim range :
 -----
ev
          mean         sd        50%       2.5%     97.5%
[1,] 0.1231591 0.09935357 0.09477851 0.01487187 0.3684598
pv
         0     1
[1,] 0.881 0.119


[1] 474


 sim range :
 -----
ev
          mean         sd        50%       2.5%     97.5%
[1,] 0.1231731 0.09955224 0.09468766 0.01479475 0.3689767
pv
         0     1
[1,] 0.862 0.138


[1] 475


 sim range :
 -----
ev
          mean         sd       50%       2.5%    97.5%
[1,] 0.1231873 0.09975086 0.0945969 0.01471802 0.369555
pv
         0     1
[1,] 0.896 0.104


[1] 476


 sim range :
 -----
ev
          mean         sd        50%       2.5%     97.5%
[1,] 0.1232016 0.09994944 0.09450621 0.01464169 0.3702206
pv
         0     1
[1,] 0.875 0.125


[1] 477


 sim range :
 -----
ev
         mean       sd        50%       2.5%     97.5%
[1,] 0.123216 0.100148 0.09441561 0.01456574 0.3708866
pv
         0     1
[1,] 0.879 0.121


[1] 478


 sim range :
 -----
ev
          mean        sd       50%       2.5%     97.5%
[1,] 0.1232306 0.1003465 0.0943251 0.01449018 0.3715531
pv
         0     1
[1,] 0.887 0.113


[1] 479


 sim range :
 -----
ev
          mean        sd        50%       2.5%     97.5%
[1,] 0.1232454 0.1005449 0.09424271 0.01441501 0.3722201
pv
         0     1
[1,] 0.878 0.122


[1] 480


 sim range :
 -----
ev
          mean        sd        50%       2.5%     97.5%
[1,] 0.1232603 0.1007433 0.09417462 0.01434022 0.3728876
pv
         0     1
[1,] 0.879 0.121


[1] 481


 sim range :
 -----
ev
          mean        sd        50%       2.5%     97.5%
[1,] 0.1232754 0.1009417 0.09410657 0.01426582 0.3735556
pv
         0     1
[1,] 0.888 0.112


[1] 482


 sim range :
 -----
ev
          mean      sd        50%      2.5%     97.5%
[1,] 0.1232906 0.10114 0.09403856 0.0141918 0.3742241
pv
         0     1
[1,] 0.883 0.117


[1] 483


 sim range :
 -----
ev
         mean        sd       50%       2.5%    97.5%
[1,] 0.123306 0.1013383 0.0939706 0.01411815 0.374893
pv
         0     1
[1,] 0.887 0.113


[1] 484


 sim range :
 -----
ev
          mean        sd        50%       2.5%     97.5%
[1,] 0.1233215 0.1015365 0.09390268 0.01404488 0.3755624
pv
         0     1
[1,] 0.865 0.135


[1] 485


 sim range :
 -----
ev
          mean        sd        50%       2.5%     97.5%
[1,] 0.1233372 0.1017347 0.09382691 0.01397199 0.3762323
pv
         0     1
[1,] 0.894 0.106


[1] 486


 sim range :
 -----
ev
         mean        sd        50%       2.5%     97.5%
[1,] 0.123353 0.1019329 0.09374732 0.01389947 0.3769027
pv
         0     1
[1,] 0.884 0.116


[1] 487


 sim range :
 -----
ev
         mean       sd       50%       2.5%     97.5%
[1,] 0.123369 0.102131 0.0936678 0.01382732 0.3775736
pv
         0     1
[1,] 0.866 0.134


[1] 488


 sim range :
 -----
ev
          mean       sd        50%       2.5%     97.5%
[1,] 0.1233851 0.102329 0.09358833 0.01375554 0.3782449
pv
         0     1
[1,] 0.874 0.126


[1] 489


 sim range :
 -----
ev
          mean       sd        50%       2.5%     97.5%
[1,] 0.1234014 0.102527 0.09350893 0.01368412 0.3789167
pv
         0     1
[1,] 0.869 0.131


[1] 490


 sim range :
 -----
ev
          mean       sd        50%       2.5%     97.5%
[1,] 0.1234178 0.102725 0.09342959 0.01361308 0.3795889
pv
         0     1
[1,] 0.878 0.122


[1] 491


 sim range :
 -----
ev
          mean        sd        50%       2.5%     97.5%
[1,] 0.1234344 0.1029229 0.09335031 0.01354239 0.3802616
pv
         0     1
[1,] 0.871 0.129


[1] 492


 sim range :
 -----
ev
          mean        sd        50%       2.5%     97.5%
[1,] 0.1234511 0.1031208 0.09327109 0.01347207 0.3809348
pv
         0     1
[1,] 0.884 0.116


[1] 493


 sim range :
 -----
ev
         mean        sd        50%       2.5%     97.5%
[1,] 0.123468 0.1033186 0.09319194 0.01340211 0.3816084
pv
         0     1
[1,] 0.877 0.123


[1] 494


 sim range :
 -----
ev
         mean        sd       50%       2.5%     97.5%
[1,] 0.123485 0.1035163 0.0931135 0.01333251 0.3822825
pv
         0     1
[1,] 0.869 0.131


[1] 495


 sim range :
 -----
ev
          mean        sd        50%       2.5%    97.5%
[1,] 0.1235021 0.1037141 0.09302098 0.01326326 0.382957
pv
         0     1
[1,] 0.884 0.116


[1] 496


 sim range :
 -----
ev
          mean        sd        50%       2.5%    97.5%
[1,] 0.1235194 0.1039117 0.09292854 0.01319437 0.383632
pv
         0     1
[1,] 0.873 0.127


[1] 497


 sim range :
 -----
ev
          mean        sd        50%       2.5%     97.5%
[1,] 0.1235369 0.1041094 0.09283619 0.01312583 0.3843074
pv
         0     1
[1,] 0.877 0.123


[1] 498


 sim range :
 -----
ev
          mean        sd        50%       2.5%     97.5%
[1,] 0.1235545 0.1043069 0.09274392 0.01305765 0.3849833
pv
         0     1
[1,] 0.866 0.134


[1] 499


 sim range :
 -----
ev
          mean        sd        50%       2.5%     97.5%
[1,] 0.1235722 0.1045044 0.09265174 0.01298981 0.3856597
pv
         0     1
[1,] 0.873 0.127


[1] 500


 sim range :
 -----
ev
          mean        sd        50%       2.5%     97.5%
[1,] 0.1235901 0.1047019 0.09255965 0.01292232 0.3863364
pv
         0     1
[1,] 0.881 0.119


[1] 501


 sim range :
 -----
ev
          mean        sd        50%       2.5%     97.5%
[1,] 0.1236081 0.1048993 0.09247832 0.01285518 0.3870136
pv
         0     1
[1,] 0.881 0.119


[1] 502


 sim range :
 -----
ev
          mean        sd        50%       2.5%     97.5%
[1,] 0.1236263 0.1050967 0.09241622 0.01278838 0.3876913
pv
         0     1
[1,] 0.883 0.117


[1] 503


 sim range :
 -----
ev
          mean       sd        50%       2.5%    97.5%
[1,] 0.1236446 0.105294 0.09235416 0.01272192 0.388429
pv
         0     1
[1,] 0.865 0.135


[1] 504


 sim range :
 -----
ev
         mean        sd        50%       2.5%     97.5%
[1,] 0.123663 0.1054913 0.09229214 0.01265581 0.3892345
pv
         0     1
[1,] 0.883 0.117


[1] 505


 sim range :
 -----
ev
          mean        sd        50%       2.5%     97.5%
[1,] 0.1236816 0.1056885 0.09223016 0.01259008 0.3900407
pv
        0    1
[1,] 0.88 0.12


[1] 506


 sim range :
 -----
ev
          mean        sd        50%       2.5%     97.5%
[1,] 0.1237004 0.1058857 0.09216822 0.01252477 0.3908474
pv
        0    1
[1,] 0.88 0.12


[1] 507


 sim range :
 -----
ev
          mean        sd        50%       2.5%     97.5%
[1,] 0.1237192 0.1060828 0.09210632 0.01245979 0.3916548
pv
         0     1
[1,] 0.896 0.104


[1] 508


 sim range :
 -----
ev
          mean        sd        50%       2.5%     97.5%
[1,] 0.1237382 0.1062798 0.09204445 0.01239515 0.3924627
pv
         0     1
[1,] 0.865 0.135


[1] 509


 sim range :
 -----
ev
          mean        sd        50%       2.5%     97.5%
[1,] 0.1237574 0.1064768 0.09198262 0.01233084 0.3932713
pv
         0     1
[1,] 0.894 0.106


[1] 510


 sim range :
 -----
ev
          mean        sd        50%       2.5%     97.5%
[1,] 0.1237767 0.1066738 0.09190652 0.01226686 0.3940804
pv
         0     1
[1,] 0.877 0.123


[1] 511


 sim range :
 -----
ev
          mean        sd        50%      2.5%     97.5%
[1,] 0.1237961 0.1068707 0.09181724 0.0122032 0.3948901
pv
         0     1
[1,] 0.863 0.137


[1] 512


 sim range :
 -----
ev
          mean        sd        50%       2.5%     97.5%
[1,] 0.1238157 0.1070675 0.09172804 0.01213988 0.3957004
pv
         0     1
[1,] 0.876 0.124
ci.plot(zhw.tit)

In the chunk above, I estimated the the effect of the price of ticket on surviving the titanic using zelig 5. First, I estimated the model, then I use the command “set range” to set the maximum and minimum values for fare, and then ran the simulate that estimates the probability of surviving based on fare price.

The results show that fare is not a significant contributor to one’s probability of surviving the titanic. There is only a .021 (2.1%) difference in the mean survival rate of those who paid 0 dollars for a ticket and those who paid 512 dollars for a ticket.

Sex difference using Zelig 5:

fd <- zhw.tit$getqi(xvalue="x1", qi="fd")
summary(fd)
       V1        
 Min.   :0.1169  
 1st Qu.:0.2249  
 Median :0.2553  
 Mean   :0.2556  
 3rd Qu.:0.2858  
 Max.   :0.3994  

In the chunk above, I continued to extend my analyses to the probability of surviving the titanic based off one’s sex. First I estimated the model, then I set x and x1 as male and female, then proceeded to run the simulation. Last, I estimated the first difference of sex.

The results of the simulation show that females have a higher probability of surviving the titanic than their male counterparts, such that females have have a 39% chance of surviving the titanic and males have only a 14% chance of survival. The first difference tells us that the difference between survival for males and females is 25%.

Class differences

First class:

zhw.tit <- zlogit$new()
zhw.tit$zelig(survival ~ age + sex*pclass + fare, data = titanic)
zhw.tit$setx(sex = "male", pclass = "1st")
zhw.tit$setx1(sex = "female", pclass = "1st")
zhw.tit$sim()
summary(zhw.tit)

 sim x :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.4530581 0.05060707 0.4524582 0.3592681 0.5521067
pv
         0     1
[1,] 0.549 0.451

 sim x1 :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.9730613 0.01434672 0.9760337 0.9391917 0.9910279
pv
         0     1
[1,] 0.028 0.972
fd
          mean         sd       50%      2.5%     97.5%
[1,] 0.5200032 0.04949576 0.5218192 0.4260899 0.6155547

In the above chunk, I ran a simulation to estimate the probability of sex AND class on surviving the titanic. In this simulation, I estimated the probability of surviving the titanic in first class, and estimated the difference between males and females.

The results show that there is a significant difference between the chance of survival for men and women in first class. Females have a 97% chance of survival if they were seated in first class, as comapred to males who have a 45% chance of survival. Taken together, males have a better chance at survival in first class than anywhere else.

Second class:

zhw.tit <- zlogit$new()
zhw.tit$zelig(survival ~ age + sex*pclass + fare, data = titanic)
zhw.tit$setx(sex = "male", pclass = "2nd")
zhw.tit$setx1(sex = "female", pclass = "2nd")
zhw.tit$sim()
summary(zhw.tit)

 sim x :
 -----
ev
          mean         sd       50%       2.5%     97.5%
[1,] 0.1386963 0.02931324 0.1352395 0.09025659 0.2052227
pv
         0     1
[1,] 0.865 0.135

 sim x1 :
 -----
ev
         mean         sd       50%      2.5%     97.5%
[1,] 0.889027 0.03279302 0.8915889 0.8155068 0.9417413
pv
         0     1
[1,] 0.106 0.894
fd
          mean         sd       50%     2.5%     97.5%
[1,] 0.7503307 0.04351493 0.7529219 0.660709 0.8264288

In this simulation, I estimated the probability of surviving the titanic in second class, and estimated the difference between males and females.

The results show, again, that females have a higher chance of surviving the titanic than males do, such that their probability of surviving is .89 (89%), and the chance of male survival is .14 (14%).

Third class:

zhw.tit <- zlogit$new()
zhw.tit$zelig(survival ~ age + sex*pclass + fare, data = titanic)
zhw.tit$setx(sex = "male", pclass = "3rd")
zhw.tit$setx1(sex = "female", pclass = "3rd")
zhw.tit$sim()
summary(zhw.tit)

 sim x :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.1405366 0.02023464 0.1394446 0.1051416 0.1852217
pv
         0     1
[1,] 0.855 0.145

 sim x1 :
 -----
ev
          mean         sd       50%      2.5%     97.5%
[1,] 0.3961808 0.04340501 0.3943631 0.3154749 0.4807446
pv
         0     1
[1,] 0.623 0.377
fd
          mean         sd       50%      2.5%     97.5%
[1,] 0.2556442 0.04356916 0.2553489 0.1766789 0.3358499

In this simulation, I estimated the probability of surviving the titanic in third class class, and estimated the difference between males and females.

The results show, again, that females have a higher chance of surviving the titanic than males do, such that their probability of surviving is .40 (40%), and the chance of male survival is .14 (14%). Taken with the other results, females in third class have a lower chance of survival than females in any other cabin. It is also worth noting that males do not have a different chance of surviving in the second or third class.

Regression using own data

library(readr)
health <- read_csv("/Users/sophia.halkitis/Desktop/R/Datasets/health.csv", col_names = TRUE)
Parsed with column specification:
cols(
  .default = col_double(),
  Geo_FIPS = col_integer(),
  Geo_NAME = col_character(),
  Geo_QNAME = col_character(),
  Geo_STATE = col_integer(),
  Geo_COUNTY = col_character(),
  SE_NV002_001 = col_integer(),
  SE_T004_001 = col_integer(),
  SE_T004_002 = col_integer(),
  SE_T004_003 = col_integer(),
  SE_NV005_001 = col_integer(),
  SE_NV005_002 = col_integer(),
  SE_NV005_003 = col_integer(),
  SE_T007_001 = col_integer(),
  SE_T008_001 = col_integer(),
  SE_T008_002 = col_integer(),
  SE_T008_003 = col_integer(),
  SE_T008_004 = col_integer(),
  SE_NV007_001 = col_character(),
  SE_NV007_002 = col_integer(),
  SE_T010_001 = col_integer()
  # ... with 2 more columns
)
See spec(...) for full column specifications.
head(health)
healthsub <- subset(health, select=c(Geo_NAME, SE_T001_001, SE_T001_002, SE_NV003_001, SE_NV003_002, SE_NV005_002, SE_NV006_001, SE_NV006_002, SE_NV006_004, SE_NV008_001, SE_NV008_002, SE_NV008_003, SE_T011_001, SE_T011_002, SE_T012_003)) %>%
    rename(
      PhysUn = SE_T001_001,
      MentUn  = SE_T001_002,
      PCPRate = SE_NV003_001,
      MHPRate = SE_NV003_002,
      NoInsR = SE_NV005_002,
      InfMortR = SE_NV006_001,
      ChldMortR = SE_NV006_002,
      DrgMortR = SE_NV006_004,
      TeenBirthR = SE_NV008_001,
      ChylamRate = SE_NV008_002,
      HIVRate = SE_NV008_003, 
      CurrentSmokP = SE_T011_001,
      CurrentDrinkP = SE_T011_002,
      ObeseP = SE_T012_003)
head(healthsub)

sapply(healthsub, max, na.rm=TRUE)
sapply(healthsub, min, na.rm=TRUE)

mean(healthsub$TeenBirthR)
healthsub <- mutate(healthsub, meanteenbirth = ifelse(TeenBirthR <= "23",0,1))

Here, I created a subset of the data called “Health” which is county level health data that I obtained from social explorer. I then renamed the column headers to be more inutive, and used the sapply command to obtain the maximum and minimum values. Last, I obtained the mean value of the teen birth rate for each county, and created a new dichotmous variable that indicates whether a county had a teen birth rate either below or above the mean. I did this so that I could interepret the results of my simulation more clearly.

Estimating the model

zhealth <- zlogit$new()
zhealth$zelig(meanteenbirth ~ CurrentSmokP, data = healthsub)
summary(zhealth)
Model: 

Call:
zhealth$zelig(formula = meanteenbirth ~ CurrentSmokP, data = healthsub)

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-2.0312  -1.0522   0.5105   1.0110   1.9250  

Coefficients:
             Estimate Std. Error z value Pr(>|z|)
(Intercept)  -10.5287     3.8338  -2.746  0.00603
CurrentSmokP   0.7077     0.2536   2.791  0.00526

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 85.692  on 61  degrees of freedom
Residual deviance: 74.961  on 60  degrees of freedom
AIC: 78.961

Number of Fisher Scoring iterations: 4

Next step: Use 'setx' method

After the data was more tidy, I continued to estimate my linear regression using Zelig 5. I decided to investigate how the percentage of current smokers in a county effect a county’s teen birth rate (above or below average). The dependent variable is if the teen birth rate is above or below average, and the IV is the percent of current smokers in county. I had originally estimated this model with 3 other independent variables, but the only variable with significant results was the percentage of current smokers in the county.

Setting the parameters, running the simulation, and plotting

zhealth <- zls$new()
zhealth$zelig(meanteenbirth ~ CurrentSmokP, data = healthsub)
summary(zhealth)
Model: 

Call:
zhealth$zelig(formula = meanteenbirth ~ CurrentSmokP, data = healthsub)

Residuals:
     Min       1Q   Median       3Q      Max 
-0.88618 -0.44564  0.09984  0.41451  0.82707 

Coefficients:
             Estimate Std. Error t value Pr(>|t|)
(Intercept)  -1.57524    0.63793  -2.469  0.01640
CurrentSmokP  0.13985    0.04215   3.318  0.00154

Residual standard error: 0.4662 on 60 degrees of freedom
Multiple R-squared:  0.155, Adjusted R-squared:  0.141 
F-statistic: 11.01 on 1 and 60 DF,  p-value: 0.001544

Next step: Use 'setx' method
zhealth$setrange(CurrentSmokP = 11:19)
zhealth$sim()
summary(zhealth)

[1] 11


 sim range :
 -----
ev
         mean        sd         50%       2.5%     97.5%
1 -0.04952409 0.1785407 -0.04740265 -0.3974758 0.3066978
pv
            mean        sd         50%      2.5%     97.5%
[1,] -0.06137712 0.4904594 -0.05004501 -1.043983 0.9054191


[1] 12


 sim range :
 -----
ev
        mean        sd       50%       2.5%     97.5%
1 0.09256927 0.1403152 0.0917177 -0.1820235 0.3619543
pv
          mean        sd        50%       2.5%    97.5%
[1,] 0.1024087 0.4917095 0.09850477 -0.8932502 1.066323


[1] 13


 sim range :
 -----
ev
       mean        sd       50%       2.5%     97.5%
1 0.2346626 0.1043666 0.2312377 0.03283123 0.4352798
pv
          mean        sd       50%       2.5%   97.5%
[1,] 0.2435011 0.4689076 0.2177115 -0.6436969 1.18905


[1] 14


 sim range :
 -----
ev
      mean         sd       50%      2.5%     97.5%
1 0.376756 0.07408701 0.3751693 0.2349771 0.5215725
pv
          mean        sd       50%       2.5%    97.5%
[1,] 0.3737993 0.4768082 0.3804816 -0.5604263 1.295684


[1] 15


 sim range :
 -----
ev
       mean         sd       50%      2.5%     97.5%
1 0.5188494 0.05897324 0.5187428 0.4039926 0.6362471
pv
          mean        sd       50%       2.5%    97.5%
[1,] 0.5196384 0.4680793 0.5162193 -0.4063806 1.419462


[1] 16


 sim range :
 -----
ev
       mean        sd       50%      2.5%     97.5%
1 0.6609427 0.0697084 0.6617376 0.5241665 0.7974846
pv
          mean        sd       50%       2.5%    97.5%
[1,] 0.6421691 0.4742686 0.6455761 -0.2958969 1.553411


[1] 17


 sim range :
 -----
ev
       mean         sd       50%     2.5%     97.5%
1 0.8030361 0.09814856 0.8031085 0.600841 0.9811039
pv
          mean        sd       50%        2.5%    97.5%
[1,] 0.7948796 0.4711027 0.7823548 -0.08815457 1.742217


[1] 18


 sim range :
 -----
ev
       mean        sd       50%      2.5%    97.5%
1 0.9451295 0.1334147 0.9475112 0.6746808 1.187507
pv
         mean        sd      50%       2.5%    97.5%
[1,] 0.939379 0.4768482 0.922559 0.01115446 1.788113


[1] 19


 sim range :
 -----
ev
      mean        sd      50%      2.5%  97.5%
1 1.087223 0.1713426 1.092677 0.7535076 1.3997
pv
         mean        sd     50%      2.5%    97.5%
[1,] 1.110659 0.4951587 1.08434 0.1533872 2.094216
ci.plot(zhealth)

After I estimated the model, I set the range for the simulation as the maximum and minimum percent of smokers in the county. I then proceeded to run the simulation.

The results show a 1.04 change in means from the counties with the lowest percent of smokers to the counties with the most. This indicates that there is a higher chance of an above average teen birth rate in counties with a higher percent of current smokers.

LS0tCnRpdGxlOiAiU2ltdWxhdGlvbiBIb21ld29yayIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQoKYGBge3J9CmxpYnJhcnkobWFncml0dHIpCmxpYnJhcnkoWmVsaWcpCmxpYnJhcnkoZHBseXIpCmxpYnJhcnkodGV4cmVnKQpsaWJyYXJ5KHJhZGlhbnQuZGF0YSkKZGF0YSh0aXRhbmljKQpoZWFkKHRpdGFuaWMpCmBgYAoKCiMjVXNlIFplbGlnIDUgc3ludGF4IHRvIHJlY3JlYXRlIHNsaWRlcyA1LjIgLT4gNS4xOApgYGB7cn0KdGl0YW5pYyA8LSBtdXRhdGUodGl0YW5pYywgc3Vydml2YWwgPSBpZmVsc2Uoc3Vydml2ZWQgPT0gIk5vIiwwLDEpKQoKbWluKHRpdGFuaWMkYWdlKQptYXgodGl0YW5pYyRhZ2UpCgptaW4odGl0YW5pYyRmYXJlKQptYXgodGl0YW5pYyRmYXJlKQpgYGAKSW4gdGhpcyBjaHVuayBJIGNyZWF0ZWQgYSBuZXcgdmFyaWFibGUgY2FsbGVkIHN1cnZpdmFsLCB3aGljaCBpcyB0aGUgbnVtZXJpY2FsIHZlcnNpb24gb2YgdGhlIHZhcmlhYmxlICJzdXJ2aXZlZCIsIHdoZXJlIDAgaXMgZXF1aXZsZW5ldCB0byBub3Qgc3Vydml2aW5nLCBhbmQgMSBpcyBlcXVpdmVsZW50IHRvIHN1cnZpdmluZy4gVGhlbiwgSSBvYnRhaW5lZCB0aGUgbWF4aW11bSBhbmQgbWluaW11bSB2YWx1ZXMgZm9yIGFnZSBhbmQgZmFyZSBzbyB0aGF0IEkgY291bGQgdXNlIHRoZW0gaW4gbXkgcmFuZ2UgZm9yIHRoZSBmb2xsb3dpbmcgYW5hbHlzZXMuIAoKIyMjQWdlIGVmZmVjdCBVc2luZyBaZWxpZyA1OiAKYGBge3J9Cnpody50aXQgPC0gemxvZ2l0JG5ldygpCnpody50aXQkemVsaWcoc3Vydml2YWwgfiBhZ2UgKyBzZXgqcGNsYXNzICsgZmFyZSwgZGF0YSA9IHRpdGFuaWMpCnN1bW1hcnkoemh3LnRpdCkKCgp6aHcudGl0JHNldHJhbmdlKGFnZSA9IDA6ODApCnpody50aXQkc2ltKCkKI3N1bW1hcnkoemh3LnRpdCkKCmNpLnBsb3Qoemh3LnRpdCkKYGBgCkluIHRoZSBhYm92ZSBhbmFseXNpcywgSSBlc3RpbWF0ZWQgdGhlIGFnZSBlZmZlY3Qgb2Ygc3Vydml2aW5nIHRoZSB0aXRhbmljIHVzaW5nIHplbGlnIDUuIEZpcnN0LCBJIGVzdGltYXRlZCB0aGUgbW9kZWwsIHRoZW4gSSBzZXQgdXNlZCB0aGUgY29tbWFuZCAic2V0IHJhbmdlIiB0byBzZXQgdGhlIG1heGltdW0gYW5kIG1pbmltdW0gdmFsdWVzIGZvciBhZ2UsIGFuZCB1bHRpbWF0ZWx5IHJ1biB0aGUgc2ltdWxhdGlvbiB0aGF0IGVzdGltYXRlcyB0aGUgcHJvYmFiaWxpdHkgb2Ygc3Vydml2aW5nIGJhc2VkIG9uIGFnZS4KClRoZSByZXN1bHRzIHNob3cgdGhhdCB0aGVyZSBpcyBhIGRpZmZlcmVuY2UgaW4gdGhlIHByb2JhYmlsaXR5IG9mIHN1cnZpdmluZyBiYXNlZCBvbiBvbmVzIGFnZSwgc3VjaCB0aGF0IHRoZSBvbGRlciBvbmUgZ2V0cywgdGhlIGxlc3MgbGlrZWx5IHRoYXQgdGhleSBhcmUgdG8gc3Vydml2ZSB0aGUgdGl0YW5pYy4gVGhlcmUgaXMgYW4gYXBwcm94aW1hdGUgLjMxNiAoMzEuNiUpIGRpZmZlcmVuY2UgaW4gdGhlIG1lYW4gc3Vydml2YWwgcHJvYmFiaWxpdHkgaW4gMCB5ZWFyIG9sZHMgYXMgY29tcGFyZWQgdG8gODAgeWVhciBvbGRzLiAKCiMjI0ZhcmUgZWZmZWN0IHVzaW5nIHplbGlnIDU6IApgYGB7cn0Kemh3LnRpdCA8LSB6bG9naXQkbmV3KCkKemh3LnRpdCR6ZWxpZyhzdXJ2aXZhbCB+IGFnZSArIHNleCpwY2xhc3MgKyBmYXJlLCBkYXRhID0gdGl0YW5pYykKCnpody50aXQkc2V0cmFuZ2UoZmFyZSA9IDA6NTEyKQp6aHcudGl0JHNpbSgpCiNzdW1tYXJ5KHpody50aXQpCgpjaS5wbG90KHpody50aXQpCgpgYGAKSW4gdGhlIGNodW5rIGFib3ZlLCBJIGVzdGltYXRlZCB0aGUgdGhlIGVmZmVjdCBvZiB0aGUgcHJpY2Ugb2YgdGlja2V0IG9uIHN1cnZpdmluZyB0aGUgdGl0YW5pYyB1c2luZyB6ZWxpZyA1LiBGaXJzdCwgSSBlc3RpbWF0ZWQgdGhlIG1vZGVsLCB0aGVuIEkgdXNlIHRoZSBjb21tYW5kICJzZXQgcmFuZ2UiIHRvIHNldCB0aGUgbWF4aW11bSBhbmQgbWluaW11bSB2YWx1ZXMgZm9yIGZhcmUsIGFuZCB0aGVuIHJhbiB0aGUgc2ltdWxhdGUgdGhhdCBlc3RpbWF0ZXMgdGhlIHByb2JhYmlsaXR5IG9mIHN1cnZpdmluZyBiYXNlZCBvbiBmYXJlIHByaWNlLiAKClRoZSByZXN1bHRzIHNob3cgdGhhdCBmYXJlIGlzIG5vdCBhIHNpZ25pZmljYW50IGNvbnRyaWJ1dG9yIHRvIG9uZSdzIHByb2JhYmlsaXR5IG9mIHN1cnZpdmluZyB0aGUgdGl0YW5pYy4gVGhlcmUgaXMgb25seSBhIC4wMjEgKDIuMSUpIGRpZmZlcmVuY2UgaW4gdGhlIG1lYW4gc3Vydml2YWwgcmF0ZSBvZiB0aG9zZSB3aG8gcGFpZCAwIGRvbGxhcnMgZm9yIGEgdGlja2V0IGFuZCB0aG9zZSB3aG8gcGFpZCA1MTIgZG9sbGFycyBmb3IgYSB0aWNrZXQuIAoKCiMjI1NleCBkaWZmZXJlbmNlIHVzaW5nIFplbGlnIDU6IApgYGB7cn0Kemh3LnRpdCA8LSB6bG9naXQkbmV3KCkKemh3LnRpdCR6ZWxpZyhzdXJ2aXZhbCB+IGFnZSArIHNleCpwY2xhc3MgKyBmYXJlLCBkYXRhID0gdGl0YW5pYykKCnpody50aXQkc2V0eChzZXggPSAiZmVtYWxlIikKemh3LnRpdCRzZXR4MShzZXggPSAibWFsZSIpCnpody50aXQkc2ltKCkKc3VtbWFyeSh6aHcudGl0KQoKZmQgPC0gemh3LnRpdCRnZXRxaSh4dmFsdWU9IngxIiwgcWk9ImZkIikKc3VtbWFyeShmZCkKYGBgCkluIHRoZSBjaHVuayBhYm92ZSwgSSBjb250aW51ZWQgdG8gZXh0ZW5kIG15IGFuYWx5c2VzIHRvIHRoZSBwcm9iYWJpbGl0eSBvZiBzdXJ2aXZpbmcgdGhlIHRpdGFuaWMgYmFzZWQgb2ZmIG9uZSdzIHNleC4gRmlyc3QgSSBlc3RpbWF0ZWQgdGhlIG1vZGVsLCB0aGVuIEkgc2V0IHggYW5kIHgxIGFzIG1hbGUgYW5kIGZlbWFsZSwgdGhlbiBwcm9jZWVkZWQgdG8gcnVuIHRoZSBzaW11bGF0aW9uLiBMYXN0LCBJIGVzdGltYXRlZCB0aGUgZmlyc3QgZGlmZmVyZW5jZSAgb2Ygc2V4LiAKClRoZSByZXN1bHRzIG9mIHRoZSBzaW11bGF0aW9uIHNob3cgdGhhdCBmZW1hbGVzIGhhdmUgYSBoaWdoZXIgcHJvYmFiaWxpdHkgb2Ygc3Vydml2aW5nIHRoZSB0aXRhbmljIHRoYW4gdGhlaXIgbWFsZSBjb3VudGVycGFydHMsIHN1Y2ggdGhhdCBmZW1hbGVzIGhhdmUgaGF2ZSBhIDM5JSBjaGFuY2Ugb2Ygc3Vydml2aW5nIHRoZSB0aXRhbmljIGFuZCBtYWxlcyBoYXZlIG9ubHkgYSAxNCUgY2hhbmNlIG9mIHN1cnZpdmFsLiBUaGUgZmlyc3QgZGlmZmVyZW5jZSB0ZWxscyB1cyB0aGF0IHRoZSBkaWZmZXJlbmNlIGJldHdlZW4gc3Vydml2YWwgZm9yIG1hbGVzIGFuZCBmZW1hbGVzIGlzIDI1JS4gCgoKIyMjQ2xhc3MgZGlmZmVyZW5jZXMgCl9fRmlyc3QgY2xhc3M6X18KYGBge3J9Cnpody50aXQgPC0gemxvZ2l0JG5ldygpCnpody50aXQkemVsaWcoc3Vydml2YWwgfiBhZ2UgKyBzZXgqcGNsYXNzICsgZmFyZSwgZGF0YSA9IHRpdGFuaWMpCgp6aHcudGl0JHNldHgoc2V4ID0gIm1hbGUiLCBwY2xhc3MgPSAiMXN0IikKemh3LnRpdCRzZXR4MShzZXggPSAiZmVtYWxlIiwgcGNsYXNzID0gIjFzdCIpCnpody50aXQkc2ltKCkKCnN1bW1hcnkoemh3LnRpdCkKCmBgYApJbiB0aGUgYWJvdmUgY2h1bmssIEkgcmFuIGEgc2ltdWxhdGlvbiB0byBlc3RpbWF0ZSB0aGUgcHJvYmFiaWxpdHkgb2Ygc2V4IEFORCBjbGFzcyBvbiBzdXJ2aXZpbmcgdGhlIHRpdGFuaWMuIEluIHRoaXMgc2ltdWxhdGlvbiwgSSBlc3RpbWF0ZWQgdGhlIHByb2JhYmlsaXR5IG9mIHN1cnZpdmluZyB0aGUgdGl0YW5pYyBpbiBmaXJzdCBjbGFzcywgYW5kIGVzdGltYXRlZCB0aGUgZGlmZmVyZW5jZSBiZXR3ZWVuIG1hbGVzIGFuZCBmZW1hbGVzLiAKClRoZSByZXN1bHRzIHNob3cgdGhhdCB0aGVyZSBpcyBhIHNpZ25pZmljYW50IGRpZmZlcmVuY2UgYmV0d2VlbiB0aGUgY2hhbmNlIG9mIHN1cnZpdmFsIGZvciBtZW4gYW5kIHdvbWVuIGluIGZpcnN0IGNsYXNzLiBGZW1hbGVzIGhhdmUgYSA5NyUgY2hhbmNlIG9mIHN1cnZpdmFsIGlmIHRoZXkgd2VyZSBzZWF0ZWQgaW4gZmlyc3QgY2xhc3MsIGFzIGNvbWFwcmVkIHRvIG1hbGVzIHdobyBoYXZlIGEgNDUlIGNoYW5jZSBvZiBzdXJ2aXZhbC4gVGFrZW4gdG9nZXRoZXIsIG1hbGVzIGhhdmUgYSBiZXR0ZXIgY2hhbmNlIGF0IHN1cnZpdmFsIGluIGZpcnN0IGNsYXNzIHRoYW4gYW55d2hlcmUgZWxzZS4gCgpfX1NlY29uZCBjbGFzczpfXwpgYGB7cn0Kemh3LnRpdCA8LSB6bG9naXQkbmV3KCkKemh3LnRpdCR6ZWxpZyhzdXJ2aXZhbCB+IGFnZSArIHNleCpwY2xhc3MgKyBmYXJlLCBkYXRhID0gdGl0YW5pYykKCnpody50aXQkc2V0eChzZXggPSAibWFsZSIsIHBjbGFzcyA9ICIybmQiKQp6aHcudGl0JHNldHgxKHNleCA9ICJmZW1hbGUiLCBwY2xhc3MgPSAiMm5kIikKemh3LnRpdCRzaW0oKQoKc3VtbWFyeSh6aHcudGl0KQoKYGBgCiBJbiB0aGlzIHNpbXVsYXRpb24sIEkgZXN0aW1hdGVkIHRoZSBwcm9iYWJpbGl0eSBvZiBzdXJ2aXZpbmcgdGhlIHRpdGFuaWMgaW4gc2Vjb25kIGNsYXNzLCBhbmQgZXN0aW1hdGVkIHRoZSBkaWZmZXJlbmNlIGJldHdlZW4gbWFsZXMgYW5kIGZlbWFsZXMuIAogCiBUaGUgcmVzdWx0cyBzaG93LCBhZ2FpbiwgdGhhdCBmZW1hbGVzIGhhdmUgYSBoaWdoZXIgY2hhbmNlIG9mIHN1cnZpdmluZyB0aGUgdGl0YW5pYyB0aGFuIG1hbGVzIGRvLCBzdWNoIHRoYXQgdGhlaXIgcHJvYmFiaWxpdHkgb2Ygc3Vydml2aW5nIGlzIC44OSAoODklKSwgYW5kIHRoZSBjaGFuY2Ugb2YgbWFsZSBzdXJ2aXZhbCBpcyAuMTQgKDE0JSkuCgpfX1RoaXJkIGNsYXNzOl9fCmBgYHtyfQp6aHcudGl0IDwtIHpsb2dpdCRuZXcoKQp6aHcudGl0JHplbGlnKHN1cnZpdmFsIH4gYWdlICsgc2V4KnBjbGFzcyArIGZhcmUsIGRhdGEgPSB0aXRhbmljKQoKemh3LnRpdCRzZXR4KHNleCA9ICJtYWxlIiwgcGNsYXNzID0gIjNyZCIpCnpody50aXQkc2V0eDEoc2V4ID0gImZlbWFsZSIsIHBjbGFzcyA9ICIzcmQiKQp6aHcudGl0JHNpbSgpCgpzdW1tYXJ5KHpody50aXQpCgpgYGAKSW4gdGhpcyBzaW11bGF0aW9uLCBJIGVzdGltYXRlZCB0aGUgcHJvYmFiaWxpdHkgb2Ygc3Vydml2aW5nIHRoZSB0aXRhbmljIGluIHRoaXJkIGNsYXNzIGNsYXNzLCBhbmQgZXN0aW1hdGVkIHRoZSBkaWZmZXJlbmNlIGJldHdlZW4gbWFsZXMgYW5kIGZlbWFsZXMuIAoKVGhlIHJlc3VsdHMgc2hvdywgYWdhaW4sIHRoYXQgZmVtYWxlcyBoYXZlIGEgaGlnaGVyIGNoYW5jZSBvZiBzdXJ2aXZpbmcgdGhlIHRpdGFuaWMgdGhhbiBtYWxlcyBkbywgc3VjaCB0aGF0IHRoZWlyIHByb2JhYmlsaXR5IG9mIHN1cnZpdmluZyBpcyAuNDAgKDQwJSksIGFuZCB0aGUgY2hhbmNlIG9mIG1hbGUgc3Vydml2YWwgaXMgLjE0ICgxNCUpLiBUYWtlbiB3aXRoIHRoZSBvdGhlciByZXN1bHRzLCBmZW1hbGVzIGluIHRoaXJkIGNsYXNzIGhhdmUgYSBsb3dlciBjaGFuY2Ugb2Ygc3Vydml2YWwgdGhhbiBmZW1hbGVzIGluIGFueSBvdGhlciBjYWJpbi4gSXQgaXMgYWxzbyB3b3J0aCBub3RpbmcgdGhhdCBtYWxlcyBkbyBub3QgaGF2ZSBhIGRpZmZlcmVudCBjaGFuY2Ugb2Ygc3Vydml2aW5nIGluIHRoZSBzZWNvbmQgb3IgdGhpcmQgY2xhc3MuIAoKIyMjUmVncmVzc2lvbiB1c2luZyBvd24gZGF0YSAKYGBge3J9CmxpYnJhcnkocmVhZHIpCmhlYWx0aCA8LSByZWFkX2NzdigiL1VzZXJzL3NvcGhpYS5oYWxraXRpcy9EZXNrdG9wL1IvRGF0YXNldHMvaGVhbHRoLmNzdiIsIGNvbF9uYW1lcyA9IFRSVUUpCmhlYWQoaGVhbHRoKQpgYGAKCmBgYHtyfQpoZWFsdGhzdWIgPC0gc3Vic2V0KGhlYWx0aCwgc2VsZWN0PWMoR2VvX05BTUUsIFNFX1QwMDFfMDAxLCBTRV9UMDAxXzAwMiwgU0VfTlYwMDNfMDAxLCBTRV9OVjAwM18wMDIsIFNFX05WMDA1XzAwMiwgU0VfTlYwMDZfMDAxLCBTRV9OVjAwNl8wMDIsIFNFX05WMDA2XzAwNCwgU0VfTlYwMDhfMDAxLCBTRV9OVjAwOF8wMDIsIFNFX05WMDA4XzAwMywgU0VfVDAxMV8wMDEsIFNFX1QwMTFfMDAyLCBTRV9UMDEyXzAwMykpICU+JQogICAgcmVuYW1lKAogICAgICBQaHlzVW4gPSBTRV9UMDAxXzAwMSwKICAgICAgTWVudFVuICA9IFNFX1QwMDFfMDAyLAogICAgICBQQ1BSYXRlID0gU0VfTlYwMDNfMDAxLAogICAgICBNSFBSYXRlID0gU0VfTlYwMDNfMDAyLAogICAgICBOb0luc1IgPSBTRV9OVjAwNV8wMDIsCiAgICAgIEluZk1vcnRSID0gU0VfTlYwMDZfMDAxLAogICAgICBDaGxkTW9ydFIgPSBTRV9OVjAwNl8wMDIsCiAgICAgIERyZ01vcnRSID0gU0VfTlYwMDZfMDA0LAogICAgICBUZWVuQmlydGhSID0gU0VfTlYwMDhfMDAxLAogICAgICBDaHlsYW1SYXRlID0gU0VfTlYwMDhfMDAyLAogICAgICBISVZSYXRlID0gU0VfTlYwMDhfMDAzLCAKICAgICAgQ3VycmVudFNtb2tQID0gU0VfVDAxMV8wMDEsCiAgICAgIEN1cnJlbnREcmlua1AgPSBTRV9UMDExXzAwMiwKICAgICAgT2Jlc2VQID0gU0VfVDAxMl8wMDMpCmhlYWQoaGVhbHRoc3ViKQoKc2FwcGx5KGhlYWx0aHN1YiwgbWF4LCBuYS5ybT1UUlVFKQpzYXBwbHkoaGVhbHRoc3ViLCBtaW4sIG5hLnJtPVRSVUUpCgptZWFuKGhlYWx0aHN1YiRUZWVuQmlydGhSKQpoZWFsdGhzdWIgPC0gbXV0YXRlKGhlYWx0aHN1YiwgbWVhbnRlZW5iaXJ0aCA9IGlmZWxzZShUZWVuQmlydGhSIDw9ICIyMyIsMCwxKSkKYGBgCkhlcmUsIEkgY3JlYXRlZCBhIHN1YnNldCBvZiB0aGUgZGF0YSBjYWxsZWQgIkhlYWx0aCIgd2hpY2ggaXMgY291bnR5IGxldmVsIGhlYWx0aCBkYXRhIHRoYXQgSSBvYnRhaW5lZCBmcm9tIHNvY2lhbCBleHBsb3Jlci4gSSB0aGVuIHJlbmFtZWQgdGhlIGNvbHVtbiBoZWFkZXJzIHRvIGJlIG1vcmUgaW51dGl2ZSwgYW5kIHVzZWQgdGhlIHNhcHBseSBjb21tYW5kIHRvIG9idGFpbiB0aGUgbWF4aW11bSBhbmQgbWluaW11bSB2YWx1ZXMuIExhc3QsIEkgb2J0YWluZWQgdGhlIG1lYW4gdmFsdWUgb2YgdGhlIHRlZW4gYmlydGggcmF0ZSBmb3IgZWFjaCBjb3VudHksIGFuZCBjcmVhdGVkIGEgbmV3IGRpY2hvdG1vdXMgdmFyaWFibGUgdGhhdCBpbmRpY2F0ZXMgd2hldGhlciBhIGNvdW50eSBoYWQgYSB0ZWVuIGJpcnRoIHJhdGUgZWl0aGVyIGJlbG93IG9yIGFib3ZlIHRoZSBtZWFuLiBJIGRpZCB0aGlzIHNvIHRoYXQgSSBjb3VsZCBpbnRlcmVwcmV0IHRoZSByZXN1bHRzIG9mIG15IHNpbXVsYXRpb24gbW9yZSBjbGVhcmx5LiAKCiMjI0VzdGltYXRpbmcgdGhlIG1vZGVsIApgYGB7cn0KemhlYWx0aCA8LSB6bG9naXQkbmV3KCkKemhlYWx0aCR6ZWxpZyhtZWFudGVlbmJpcnRoIH4gQ3VycmVudFNtb2tQLCBkYXRhID0gaGVhbHRoc3ViKQpzdW1tYXJ5KHpoZWFsdGgpCmBgYApBZnRlciB0aGUgZGF0YSB3YXMgbW9yZSB0aWR5LCBJIGNvbnRpbnVlZCB0byBlc3RpbWF0ZSBteSBsaW5lYXIgcmVncmVzc2lvbiB1c2luZyBaZWxpZyA1LiBJIGRlY2lkZWQgdG8gaW52ZXN0aWdhdGUgaG93IHRoZSBwZXJjZW50YWdlIG9mIGN1cnJlbnQgc21va2VycyBpbiBhIGNvdW50eSBlZmZlY3QgYSBjb3VudHkncyB0ZWVuIGJpcnRoIHJhdGUgKGFib3ZlIG9yIGJlbG93IGF2ZXJhZ2UpLiBUaGUgZGVwZW5kZW50IHZhcmlhYmxlIGlzIGlmIHRoZSB0ZWVuIGJpcnRoIHJhdGUgaXMgYWJvdmUgb3IgYmVsb3cgYXZlcmFnZSwgYW5kIHRoZSBJViBpcyB0aGUgcGVyY2VudCBvZiBjdXJyZW50IHNtb2tlcnMgaW4gY291bnR5LiAgSSBoYWQgb3JpZ2luYWxseSBlc3RpbWF0ZWQgdGhpcyBtb2RlbCB3aXRoIDMgb3RoZXIgaW5kZXBlbmRlbnQgdmFyaWFibGVzLCBidXQgdGhlIG9ubHkgdmFyaWFibGUgd2l0aCBzaWduaWZpY2FudCByZXN1bHRzIHdhcyB0aGUgcGVyY2VudGFnZSBvZiBjdXJyZW50IHNtb2tlcnMgaW4gdGhlIGNvdW50eS4gCgojIyNTZXR0aW5nIHRoZSBwYXJhbWV0ZXJzLCBydW5uaW5nIHRoZSBzaW11bGF0aW9uLCBhbmQgcGxvdHRpbmcKYGBge3J9CnpoZWFsdGggPC0gemxzJG5ldygpCnpoZWFsdGgkemVsaWcobWVhbnRlZW5iaXJ0aCB+IEN1cnJlbnRTbW9rUCwgZGF0YSA9IGhlYWx0aHN1YikKc3VtbWFyeSh6aGVhbHRoKQoKemhlYWx0aCRzZXRyYW5nZShDdXJyZW50U21va1AgPSAxMToxOSkKemhlYWx0aCRzaW0oKQpzdW1tYXJ5KHpoZWFsdGgpCgpjaS5wbG90KHpoZWFsdGgpCmBgYApBZnRlciBJIGVzdGltYXRlZCB0aGUgbW9kZWwsIEkgc2V0IHRoZSByYW5nZSBmb3IgdGhlIHNpbXVsYXRpb24gYXMgdGhlIG1heGltdW0gYW5kIG1pbmltdW0gcGVyY2VudCBvZiBzbW9rZXJzIGluIHRoZSBjb3VudHkuIEkgdGhlbiBwcm9jZWVkZWQgdG8gcnVuIHRoZSBzaW11bGF0aW9uLgoKVGhlIHJlc3VsdHMgc2hvdyBhIDEuMDQgY2hhbmdlIGluIG1lYW5zIGZyb20gdGhlIGNvdW50aWVzIHdpdGggdGhlIGxvd2VzdCBwZXJjZW50IG9mIHNtb2tlcnMgdG8gdGhlIGNvdW50aWVzIHdpdGggdGhlIG1vc3QuIFRoaXMgaW5kaWNhdGVzIHRoYXQgdGhlcmUgaXMgYSBoaWdoZXIgY2hhbmNlIG9mIGFuIGFib3ZlIGF2ZXJhZ2UgdGVlbiBiaXJ0aCByYXRlIGluIGNvdW50aWVzIHdpdGggYSBoaWdoZXIgcGVyY2VudCBvZiBjdXJyZW50IHNtb2tlcnMuIAo=