Run Gibbs sampler for the random moderator placement design
Usage
prepost_gibbs(
formula,
data,
prepost,
moderator,
covariates,
iter = 1000,
thin = 1,
burn = 0,
offset = 0,
monotonicity = TRUE,
stable = TRUE,
saturated = TRUE,
priors
)
Arguments
- formula
A formula with syntax
y ~ t
, wherey
is the (unquoted) name of the outcome andt
is the (unquoted) name of the treatment.- data
A data.frame containing variables in the formula, moderator, and covariates arguments.
- prepost
A one-sided formula with syntax
~ z
, wherez
is the indicator variable for whether the moderator was measured pre- or post-treatment.- moderator
A one-sided formuala with syntax
~ d
, whered
is the (unquoted) name of the moderator variable for the CATE.- covariates
A one-sided formula with syntax
~ x1 + x2
, where the right-hand side variables signify which covariates the Gibbs will use to try and narrow the bounds.- iter
Integer indicating the number of iterations for the Gibbs sampler.
- thin
Integer indicating how often the Gibbs sampler should save a draw.
- burn
Integer indicating how many iterations should be performed before saving draws in the Gibbs sampler.
- offset
A numeric value indicating the center of the prior distribution for the covariate coefficients.
- monotonicity
A logical signifying whether the model assumes monotonicity.
- stable
A logical signifying whether the model assumes that the the pre vs post indicator does not affect the moderator under the control condition for treatment.
- saturated
A logical indicating whether the coefficients on the covariates are allowed to vary by the principal strata.
- priors
A list object containing the priors for the Gibbs sampler. Priors include beta.precision, psi.precision, alpha, y.alpha, and y.beta.
Examples
data(land_experiment)
prepost_gibbs(
support ~ treat_comb,
data = land_experiment,
moderator = ~ land_insecure,
prepost = ~ prepost,
covariates = ~ educ + age,
iter = 50
)
#> $betas
#> (Intercept) educ age treat_comb prepost
#> [1,] 30.212575 -11.8537972 2.7309365459 4.823811749 5.29340999
#> [2,] 17.347981 -11.6205965 1.3750416798 2.376185810 6.01355064
#> [3,] 17.063923 -8.2752723 0.7386389708 -0.668864652 -2.91078855
#> [4,] 12.183307 -5.5925127 0.3632198760 0.454410167 1.80916241
#> [5,] 13.688792 -4.2170128 0.2033732662 -0.861457855 0.37006097
#> [6,] 10.922087 -2.7178304 0.0746654762 -0.175505732 -0.10901944
#> [7,] 11.006588 -2.0502704 0.0147585403 -1.086761237 -0.51287137
#> [8,] 8.950237 -1.5977228 0.0033455239 -0.106795931 0.37992285
#> [9,] 7.330855 -1.0669566 -0.0254633815 0.399247503 0.30056232
#> [10,] 7.297227 -0.9357218 -0.0142156622 -0.476618773 -1.02979146
#> [11,] 4.912926 -0.6553579 0.0083146234 0.242060485 -0.12586741
#> [12,] 4.901685 -0.5969041 -0.0121492009 0.271876410 -1.02519872
#> [13,] 4.040699 -0.3545045 -0.0161197523 0.119241875 -0.54934397
#> [14,] 4.145066 -0.4392895 -0.0270533618 0.445127049 0.83497893
#> [15,] 3.599293 -0.4526549 -0.0064483567 0.512666844 0.57357107
#> [16,] 4.694732 -0.4779392 -0.0233520915 0.175285515 0.27727913
#> [17,] 3.428479 -0.3548127 -0.0059827652 0.497414121 -0.73498836
#> [18,] 3.800330 -0.2594181 -0.0192530758 -0.043071991 -0.70506497
#> [19,] 4.003899 -0.3124689 -0.0152006947 0.008227955 -0.24163257
#> [20,] 2.804720 -0.2600400 -0.0049623161 0.441984856 0.03918917
#> [21,] 2.747150 -0.3556296 0.0058852814 0.786562649 0.15980025
#> [22,] 2.985998 -0.3338013 0.0005566419 0.281011379 0.05687172
#> [23,] 2.791271 -0.3408166 0.0086912037 0.208047340 -0.36534928
#> [24,] 4.073964 -0.4512847 -0.0109389694 -0.362750012 -0.23847874
#> [25,] 4.325852 -0.3941545 -0.0291004076 0.409961672 0.54395510
#> [26,] 3.707239 -0.3387944 -0.0212457590 -0.018005373 0.73021423
#> [27,] 4.029704 -0.3263130 -0.0256834692 0.058645520 -0.54139032
#> [28,] 4.249006 -0.3347559 -0.0236237130 -0.305848090 -0.38014175
#> [29,] 4.974450 -0.3493922 -0.0403770422 -0.318752854 -0.64272830
#> [30,] 3.732185 -0.1944656 -0.0148292743 -0.381111574 -0.62270158
#> [31,] 4.846566 -0.4276565 -0.0249360265 -0.089333521 -0.52596283
#> [32,] 4.292095 -0.4004673 -0.0168863283 0.213068250 -0.27712678
#> [33,] 4.895765 -0.4489806 -0.0336228212 0.239584188 -0.40395924
#> [34,] 3.681832 -0.3828890 -0.0172059098 0.376950840 0.19585001
#> [35,] 4.188628 -0.3916016 -0.0205702133 0.046475689 0.28684475
#> [36,] 5.324942 -0.3947830 -0.0309248254 -0.361720247 -1.12355534
#> [37,] 4.944411 -0.3729182 -0.0306921257 0.150722021 -0.62928035
#> [38,] 4.486559 -0.4623836 -0.0287314194 0.554606711 -0.29102027
#> [39,] 4.364822 -0.3101202 -0.0278291419 -0.030334676 -0.06091105
#> [40,] 4.506959 -0.4318025 -0.0180695375 0.533079710 -0.31007643
#> [41,] 4.836545 -0.4973782 -0.0284422474 0.411153898 0.36543160
#> [42,] 4.614945 -0.4270191 -0.0353621807 0.246796047 0.74254676
#> [43,] 5.240896 -0.5073752 -0.0344816041 0.172709983 -0.13730172
#> [44,] 4.534536 -0.3864195 -0.0270793686 0.165161959 -0.05417443
#> [45,] 4.712136 -0.4405812 -0.0330813841 0.458692730 -0.61694161
#> [46,] 3.471315 -0.4166839 -0.0059984988 0.857158838 0.92978014
#> [47,] 3.759170 -0.3206376 -0.0063121806 0.471160481 -0.60242632
#> [48,] 3.090188 -0.4698050 0.0042258959 1.130620193 -0.20015615
#> [49,] 2.306876 -0.2739481 0.0077575500 0.708580982 -0.02784993
#> [50,] 3.509441 -0.3476037 -0.0058468577 -0.222844043 -0.88406421
#> treat_comb:prepost s111 s100 treat_comb:s111 treat_comb:s100
#> [1,] 12.38508036 -16.0854471 34.8900396 -1.79562827 2.7953828
#> [2,] 28.77298651 6.1564731 66.9314577 6.75905661 14.3515350
#> [3,] 26.71017604 5.0895341 62.4176858 -0.09275494 20.8968924
#> [4,] 15.27585025 0.8378360 60.7336731 3.73337903 21.2961738
#> [5,] 9.84313207 -1.6704474 54.5631320 1.09425860 17.8944962
#> [6,] 5.88857726 -1.7982097 51.0449948 2.09550292 16.9745082
#> [7,] 3.60591852 0.2266369 44.0807209 -0.27072828 17.8993154
#> [8,] 1.01684010 -2.0359767 37.0692414 2.11208211 21.6668082
#> [9,] 0.72957869 -1.4214501 33.9543130 0.96716921 20.0359952
#> [10,] 0.90020319 -2.3939860 30.5886513 0.90565259 20.6824489
#> [11,] -0.54519334 -1.2987377 26.8148044 0.53032047 18.6326663
#> [12,] -0.13449216 -1.3871678 25.1130590 1.11029700 16.5151741
#> [13,] 0.01966229 -2.9061950 25.9338554 1.73804931 11.1488347
#> [14,] -0.88912130 -1.6731489 20.5733907 1.19108026 11.1639277
#> [15,] -0.77906435 -1.7343319 17.0809838 0.11113084 12.5230598
#> [16,] -0.84890599 -0.5362716 19.2482317 -0.32387545 10.0894386
#> [17,] -0.12872987 -0.2569527 17.3164174 0.89757804 13.4261073
#> [18,] 0.02809284 -2.4145838 14.5268147 1.92596028 10.3813563
#> [19,] -0.69823663 -2.2667494 14.9392755 1.44796965 12.2762492
#> [20,] -0.72669377 -1.1833967 13.4602731 -0.64486009 14.9684720
#> [21,] -1.20576174 -0.8309939 15.3468947 -0.99437209 11.7998424
#> [22,] -0.64293179 -3.3399519 17.1281254 1.60335217 7.1953350
#> [23,] -0.13304321 -0.8635163 18.2163531 0.53298409 4.8727568
#> [24,] 0.42955891 -2.8840753 18.4868750 0.74379101 4.8801985
#> [25,] -1.37236594 -1.2790325 12.8530781 0.26864462 8.5853668
#> [26,] -1.13946345 -1.9378596 13.9491071 1.08065176 10.6051773
#> [27,] -0.46392840 -3.1875464 9.8933042 3.69908846 12.7339749
#> [28,] -0.13817933 -1.1550524 7.8714792 0.10462730 11.7705331
#> [29,] 0.31080921 -1.8624195 10.5050470 1.55705477 5.6645653
#> [30,] 0.25821147 -2.7913742 8.8717427 1.80600523 5.9997861
#> [31,] 0.20251689 -3.5938461 7.1395115 2.56021898 5.7203134
#> [32,] -0.70491284 -2.5842031 5.9880668 1.52305017 8.3868046
#> [33,] -0.13888481 -1.2428148 5.7400437 2.51895996 3.2404098
#> [34,] -0.78688316 -0.5314306 1.9704547 -0.77282832 4.3935361
#> [35,] -0.84177185 -1.5809878 2.0634683 2.03331287 3.9465056
#> [36,] 0.41488758 -3.1805812 0.4524743 3.10301080 5.7966402
#> [37,] -0.56582182 -1.1518807 -0.2815877 0.63092544 5.0518705
#> [38,] -0.68149696 -0.9446442 0.3661174 -0.89006578 4.2246814
#> [39,] -0.46278836 -1.6685370 -1.0881035 0.96888117 5.5890004
#> [40,] -0.93161729 -3.0423034 0.3737793 1.56632219 5.5037868
#> [41,] -1.39521079 -2.0732773 2.0801083 1.48960692 2.8872792
#> [42,] -1.00882729 -0.1659696 2.7554272 -0.04853148 2.9687928
#> [43,] -0.45263780 -2.1811078 3.0003752 2.98964639 1.1511771
#> [44,] -0.25057656 -1.8456675 2.5776401 2.08528871 2.4649329
#> [45,] -0.28530712 -1.6430893 0.1855946 2.20070839 4.2198279
#> [46,] -2.43454411 -1.2166178 1.3038905 -0.63598625 2.8369106
#> [47,] -1.00295627 -2.7132943 4.2654718 1.05444450 0.8263301
#> [48,] -0.94581353 -1.5824133 3.8272797 3.30484152 3.4094391
#> [49,] -1.14871852 -0.1807384 3.6623857 0.94106057 2.2826558
#> [50,] 1.02159834 -2.7651585 2.2549343 1.84652234 3.7990027
#> prepost:s111 prepost:s100 treat_comb:prepost:s111 treat_comb:prepost:s100
#> [1,] -9.08902544 2.5065562 -1.9102114 -1.6357955
#> [2,] -6.04116931 -4.5800302 3.2972060 -7.8315901
#> [3,] -1.83960951 9.9069841 2.5480292 -10.0029850
#> [4,] 0.76515878 8.5341838 -1.7267658 -4.0362529
#> [5,] 3.98573927 7.3881394 0.2924717 -3.8505013
#> [6,] 1.99432789 8.6845133 -0.7771206 -1.1736642
#> [7,] -2.05745721 9.0644806 -6.4011523 2.0080721
#> [8,] 0.80188072 13.9227472 -6.4055757 -0.5635995
#> [9,] -1.51562755 12.4858726 -7.1106378 4.6673503
#> [10,] -0.84422276 13.2238419 -0.8068188 6.5682080
#> [11,] 0.57702465 16.2733360 -5.5206109 5.6827149
#> [12,] 0.44903866 15.9719396 -2.9385794 8.0775324
#> [13,] 1.42879060 11.7598871 -1.9865231 8.4678067
#> [14,] -2.25551543 9.3170962 -0.4088152 7.5948510
#> [15,] -0.73866062 12.4435483 -2.9354140 7.9100637
#> [16,] 0.04766461 5.9542653 -2.1396561 10.6128300
#> [17,] -0.07435571 7.5887424 -2.1384005 5.3791589
#> [18,] 0.94717612 12.4675361 -3.7307751 10.3642018
#> [19,] 0.86540919 8.5894971 -6.3453442 8.4548159
#> [20,] -0.73225538 6.7724095 -6.9372600 7.2637186
#> [21,] 0.99141474 7.2876928 -3.6139709 1.3687238
#> [22,] 0.86498527 9.1330914 -1.9725577 1.7890724
#> [23,] 0.94614757 10.1367176 -1.0640929 4.0307486
#> [24,] 2.92521809 3.4780080 -1.9440249 6.8571470
#> [25,] -0.90721860 4.9293471 2.6082751 8.2009915
#> [26,] -0.56022435 4.4157635 0.8284288 3.2464224
#> [27,] 3.11704872 7.6102408 0.5458453 1.8141630
#> [28,] -0.19028395 8.8938214 -0.4457905 2.0348497
#> [29,] 1.45686097 4.5905190 -3.9542227 4.2876636
#> [30,] 1.13208443 4.8415439 -2.2252209 3.3745594
#> [31,] 1.30319413 4.2380672 -1.8316654 0.2516719
#> [32,] 1.21546897 3.6964359 -3.7863391 0.2422055
#> [33,] 0.83148810 3.2634909 -1.0104577 4.0925002
#> [34,] -0.42768511 7.5488986 3.6093205 -1.5474775
#> [35,] 0.02999038 0.3913469 3.5460432 3.7848681
#> [36,] 2.75499679 3.5256066 -0.2035966 -0.2483785
#> [37,] 0.84763188 -0.6312945 3.1880497 5.2939406
#> [38,] -0.19976060 -0.8324727 3.6544257 4.6877014
#> [39,] 0.85203154 0.5797567 -1.6306955 1.4190388
#> [40,] 0.90376297 -0.4233247 -0.2661265 0.5591580
#> [41,] 0.04960359 -0.6288505 -0.8975335 3.8465255
#> [42,] -2.30668322 0.3872884 2.5790673 2.0355833
#> [43,] 0.77114443 1.9363781 -3.0080613 1.0464749
#> [44,] -0.61724503 -0.3787547 -0.4978001 3.9426253
#> [45,] 1.16810714 4.8876608 1.6015822 -3.5990264
#> [46,] 0.69610217 2.0866492 2.5055598 -1.7740828
#> [47,] 2.14275554 3.9417728 -0.6139898 -3.5016447
#> [48,] 0.54297338 3.1377323 -1.1425364 -4.8285481
#> [49,] 0.39578283 2.1694092 -0.2807254 0.8186198
#> [50,] 0.47377589 4.6096484 -0.3773665 -1.0267630
#>
#> $psis
#> $psis$s111
#> (Intercept) educ age
#> [1,] -10.105394 -12.256262889 -2.60521207
#> [2,] -16.325836 -10.183464326 -2.02589106
#> [3,] -19.803014 -8.209267560 -1.67744206
#> [4,] -20.013498 -7.463966892 -1.36330815
#> [5,] -17.587259 -6.895760593 -1.15466787
#> [6,] -18.610069 -6.201205263 -0.90459214
#> [7,] -19.329208 -5.496692955 -0.71912031
#> [8,] -17.701598 -5.305527140 -0.55344911
#> [9,] -19.293968 -4.606317402 -0.36022775
#> [10,] -15.294976 -4.755113866 -0.29588882
#> [11,] -15.797084 -4.004791102 -0.19816940
#> [12,] -12.355360 -4.049429067 -0.16425935
#> [13,] -14.002349 -3.623110012 -0.06220627
#> [14,] -19.142667 -2.625839880 0.04555993
#> [15,] -18.478821 -2.320909190 0.07761438
#> [16,] -19.540601 -2.128627443 0.15049983
#> [17,] -18.556733 -1.814596711 0.15553779
#> [18,] -15.500877 -1.504175072 0.10797882
#> [19,] -13.602548 -1.469152059 0.12648841
#> [20,] -14.302768 -1.206809069 0.16860515
#> [21,] -12.299225 -1.391049074 0.19379851
#> [22,] -10.560982 -1.073460760 0.13909171
#> [23,] -8.465825 -0.988820745 0.11400070
#> [24,] -7.913076 -0.642180681 0.10117375
#> [25,] -7.839936 -0.484262163 0.10023441
#> [26,] -7.622826 -0.426068975 0.10307026
#> [27,] -7.278693 -0.271652505 0.09103297
#> [28,] -6.684070 -0.223556029 0.08211267
#> [29,] -6.255361 -0.184722456 0.07539422
#> [30,] -6.804764 0.072830310 0.07833816
#> [31,] -4.648826 0.042312260 0.03415097
#> [32,] -3.381602 -0.216536542 0.03296244
#> [33,] -3.958780 -0.121135648 0.03398095
#> [34,] -3.483508 -0.100245245 0.02332267
#> [35,] -3.511647 -0.034093286 0.01703376
#> [36,] -2.731719 -0.081011303 0.01126985
#> [37,] -4.407543 0.068783145 0.03412813
#> [38,] -5.053451 0.167632691 0.03895857
#> [39,] -3.983048 0.037658323 0.02428336
#> [40,] -3.894394 0.039493014 0.03363367
#> [41,] -4.718815 0.249970415 0.03842147
#> [42,] -4.550900 0.205779079 0.03532320
#> [43,] -3.483999 -0.004229671 0.02043677
#> [44,] -3.528087 0.003051194 0.03259845
#> [45,] -4.897411 0.166237230 0.04568311
#> [46,] -3.048005 -0.090378755 0.02248442
#> [47,] -2.810369 -0.075584238 0.01873685
#> [48,] -2.093205 -0.178230683 0.01683502
#> [49,] -3.454613 0.013795855 0.03620723
#> [50,] -3.116492 0.068775719 0.01763898
#>
#> $psis$s100
#> (Intercept) educ age
#> [1,] 10.24032930 -2.2688133135 0.905900143
#> [2,] 1.23668924 -0.4577817304 0.189169307
#> [3,] 2.04646291 -0.1837466112 -0.010533239
#> [4,] 1.67621042 -0.0008583885 -0.042080855
#> [5,] 0.62189508 0.0822501170 -0.034596345
#> [6,] 0.08458480 0.1269302837 -0.027650169
#> [7,] -1.50697937 0.3041969811 -0.014646723
#> [8,] -1.98500929 0.3633519137 -0.022155410
#> [9,] -2.15050329 0.2916179711 -0.012465025
#> [10,] -1.71235255 0.1926721399 -0.017941133
#> [11,] -2.60336906 0.1580859266 0.006293140
#> [12,] -2.11752724 0.0581695793 0.001305941
#> [13,] -2.04043520 0.0861976829 -0.010260214
#> [14,] -1.52875357 -0.0455236688 -0.011018042
#> [15,] -1.18135831 -0.0721741629 -0.016716285
#> [16,] -1.17373089 -0.0979851258 -0.016548651
#> [17,] -0.46667103 -0.1901438585 -0.024902937
#> [18,] -1.19726562 -0.0379825160 -0.030034280
#> [19,] 0.02820969 -0.1416997852 -0.050892071
#> [20,] 0.79128925 0.0039325418 -0.092894075
#> [21,] 1.07349875 -0.0721584424 -0.092045750
#> [22,] 2.62566542 -0.1896334036 -0.126093761
#> [23,] 1.36037652 0.0083853688 -0.114728687
#> [24,] 0.70840947 0.1350459231 -0.116116377
#> [25,] -0.23910774 0.2158799332 -0.093694689
#> [26,] -1.30699195 0.1446618996 -0.056515201
#> [27,] -0.94277957 0.1149267454 -0.063534257
#> [28,] -1.53946985 0.1063154863 -0.046094360
#> [29,] -0.13550161 -0.0181781155 -0.073067012
#> [30,] -1.35590640 0.2701932464 -0.074075767
#> [31,] -1.06710780 0.2986756712 -0.087654512
#> [32,] -0.87012736 0.2071931054 -0.084569887
#> [33,] -0.55977454 0.0969729272 -0.076292372
#> [34,] -0.48029519 0.1771468657 -0.090169013
#> [35,] -0.37043512 0.1987206816 -0.086859513
#> [36,] 1.03758501 0.0877936487 -0.115724128
#> [37,] 0.56581414 0.1647514821 -0.118135561
#> [38,] 0.80345289 0.1546164210 -0.123354537
#> [39,] 0.54744773 0.2407466095 -0.131216988
#> [40,] 0.10530164 0.1644955563 -0.104061979
#> [41,] -0.01818890 0.2351722279 -0.112720416
#> [42,] -0.56555960 0.3792883456 -0.108947813
#> [43,] 0.19067976 0.2961149962 -0.128897922
#> [44,] 1.10043473 0.2284136894 -0.147534033
#> [45,] 0.57296300 0.3265533034 -0.150391562
#> [46,] -0.01646022 0.3850728993 -0.141191549
#> [47,] 0.57190701 0.2402334306 -0.149223616
#> [48,] -0.07219741 0.3433111777 -0.139688669
#> [49,] 0.80769578 0.2127852104 -0.146543467
#> [50,] 0.96981691 0.0867615717 -0.141532260
#>
#> $psis$s000
#> (Intercept) educ age
#> [1,] 0 0 0
#> [2,] 0 0 0
#> [3,] 0 0 0
#> [4,] 0 0 0
#> [5,] 0 0 0
#> [6,] 0 0 0
#> [7,] 0 0 0
#> [8,] 0 0 0
#> [9,] 0 0 0
#> [10,] 0 0 0
#> [11,] 0 0 0
#> [12,] 0 0 0
#> [13,] 0 0 0
#> [14,] 0 0 0
#> [15,] 0 0 0
#> [16,] 0 0 0
#> [17,] 0 0 0
#> [18,] 0 0 0
#> [19,] 0 0 0
#> [20,] 0 0 0
#> [21,] 0 0 0
#> [22,] 0 0 0
#> [23,] 0 0 0
#> [24,] 0 0 0
#> [25,] 0 0 0
#> [26,] 0 0 0
#> [27,] 0 0 0
#> [28,] 0 0 0
#> [29,] 0 0 0
#> [30,] 0 0 0
#> [31,] 0 0 0
#> [32,] 0 0 0
#> [33,] 0 0 0
#> [34,] 0 0 0
#> [35,] 0 0 0
#> [36,] 0 0 0
#> [37,] 0 0 0
#> [38,] 0 0 0
#> [39,] 0 0 0
#> [40,] 0 0 0
#> [41,] 0 0 0
#> [42,] 0 0 0
#> [43,] 0 0 0
#> [44,] 0 0 0
#> [45,] 0 0 0
#> [46,] 0 0 0
#> [47,] 0 0 0
#> [48,] 0 0 0
#> [49,] 0 0 0
#> [50,] 0 0 0
#>
#>
#> $s.tab
#> s111 s100 s000
#> [1,] 0.05333333 0.65333333 0.2933333
#> [2,] 0.05333333 0.60000000 0.3466667
#> [3,] 0.05333333 0.57333333 0.3733333
#> [4,] 0.05333333 0.51200000 0.4346667
#> [5,] 0.05333333 0.40266667 0.5440000
#> [6,] 0.05333333 0.38666667 0.5600000
#> [7,] 0.05600000 0.31733333 0.6266667
#> [8,] 0.05600000 0.26133333 0.6826667
#> [9,] 0.05600000 0.19466667 0.7493333
#> [10,] 0.05600000 0.17066667 0.7733333
#> [11,] 0.05600000 0.14666667 0.7973333
#> [12,] 0.05600000 0.13066667 0.8133333
#> [13,] 0.05600000 0.12266667 0.8213333
#> [14,] 0.05600000 0.10400000 0.8400000
#> [15,] 0.05600000 0.12800000 0.8160000
#> [16,] 0.05600000 0.11733333 0.8266667
#> [17,] 0.05600000 0.12266667 0.8213333
#> [18,] 0.05600000 0.09600000 0.8480000
#> [19,] 0.05600000 0.10666667 0.8373333
#> [20,] 0.05600000 0.10666667 0.8373333
#> [21,] 0.05600000 0.11733333 0.8266667
#> [22,] 0.05866667 0.12533333 0.8160000
#> [23,] 0.05866667 0.15200000 0.7893333
#> [24,] 0.05866667 0.12533333 0.8160000
#> [25,] 0.05866667 0.11200000 0.8293333
#> [26,] 0.05866667 0.11466667 0.8266667
#> [27,] 0.05866667 0.08533333 0.8560000
#> [28,] 0.05866667 0.07200000 0.8693333
#> [29,] 0.05866667 0.08800000 0.8533333
#> [30,] 0.06133333 0.08266667 0.8560000
#> [31,] 0.06666667 0.08266667 0.8506667
#> [32,] 0.05866667 0.09066667 0.8506667
#> [33,] 0.05866667 0.08000000 0.8613333
#> [34,] 0.06933333 0.07733333 0.8533333
#> [35,] 0.06933333 0.08533333 0.8453333
#> [36,] 0.07200000 0.09600000 0.8320000
#> [37,] 0.06666667 0.10133333 0.8320000
#> [38,] 0.06933333 0.11466667 0.8160000
#> [39,] 0.06933333 0.08266667 0.8480000
#> [40,] 0.06400000 0.11200000 0.8240000
#> [41,] 0.07200000 0.09600000 0.8320000
#> [42,] 0.07466667 0.08533333 0.8400000
#> [43,] 0.07200000 0.11466667 0.8133333
#> [44,] 0.06666667 0.10933333 0.8240000
#> [45,] 0.06933333 0.09600000 0.8346667
#> [46,] 0.06666667 0.08800000 0.8453333
#> [47,] 0.06933333 0.09066667 0.8400000
#> [48,] 0.07466667 0.06933333 0.8560000
#> [49,] 0.06933333 0.09066667 0.8400000
#> [50,] 0.07466667 0.08000000 0.8453333
#>
#> $delta.1
#> [1] 0.24225352 0.26619718 0.30140845 0.28450704 0.11830986 0.21549296
#> [7] -0.49354318 -0.43422115 -0.50161421 -0.05891848 -0.68038741 -0.19047619
#> [13] 0.05649718 -0.15899919 -0.32728006 -0.61864407 -0.33373688 -0.35593220
#> [19] -0.53712672 -0.38337369 -0.61057304 -0.09940767 0.07391192 -0.15336080
#> [25] 0.16739634 0.34058717 0.20718517 -0.27852176 -0.34097347 0.03693182
#> [31] 0.25714286 -0.32977080 0.17048674 0.32367203 0.29953714 0.32247765
#> [37] 0.19428571 0.38461538 -0.37028874 0.11930199 -0.04214559 0.32204611
#> [43] 0.07151980 0.19142857 0.43001984 0.27428571 0.01311439 0.12958007
#> [49] 0.07857615 0.33923425
#>
#> $delta.2
#> [1] 0.04476431 0.21658958 0.32885535 0.22086946 0.10568854 0.15163813
#> [7] -0.45339832 -0.40968473 -0.47973660 0.05223910 -0.72598984 -0.33607924
#> [13] -0.04233769 0.09431371 -0.35681143 -0.52546449 -0.25204547 -0.34234290
#> [19] -0.53898528 -0.39003500 -0.74212221 -0.08471972 -0.09053455 -0.20657136
#> [25] 0.44336195 0.30036799 0.30001225 -0.11913143 -0.49566291 -0.11634907
#> [31] 0.17245026 -0.43453250 0.17687101 0.33500970 0.44100699 0.23426002
#> [37] 0.21084131 0.36148367 -0.20413154 0.22443539 0.02707730 0.40214294
#> [43] -0.03085452 0.32834225 0.30612169 0.21361178 0.05798255 0.35175680
#> [49] 0.09801565 0.37053218
#>