Single Replication from the Simulation Study
Ivan Jacob Agaloos Pesigan
Source:vignettes/replication.Rmd
replication.Rmd
Data Generation
data <- RandomMeasurement(sim)
str(data)
#> 'data.frame': 9000 obs. of 5 variables:
#> $ id : num 1 1 1 1 1 1 1 1 1 1 ...
#> $ time: num 0 0.0333 0.0667 0.1 0.1333 ...
#> $ x : num -1.427 -0.67 -2.044 -0.977 -0.604 ...
#> $ m : num -0.972 -1.492 -1.343 -0.702 -0.842 ...
#> $ y : num -0.0563 0.4158 0.1407 -0.1219 -0.119 ...
Model Fitting
fit <- FitDynr(data = data, taskid = taskid)
#> [1] "Get ready!!!!"
#> using C compiler: ‘gcc (Ubuntu 13.3.0-6ubuntu2~24.04) 13.3.0’
#> Optimization function called.
#> Starting Hessian calculation ...
#> Finished Hessian calculation.
#> Original exit flag: 3
#> Modified exit flag: 3
#> Optimization terminated successfully: ftol_rel or ftol_abs was reached.
#> Original fitted parameters: -0.3346013 0.02166289 -0.01933279 0.7706831
#> -0.5507161 0.008072293 -0.5077725 0.7861072 -0.7079146 -1.646831 0.1363935
#> -0.1794374 -2.668642 0.1643978 -2.856534 -1.617031 -1.639487 -1.553734
#> 0.1037222 -0.02241764 0.1438094 0.02120675 0.3877461 0.3593095 -0.4531352
#> 0.1506019 -0.4459209
#>
#> Transformed fitted parameters: -0.3346013 0.02166289 -0.01933279 0.7706831
#> -0.5507161 0.008072293 -0.5077725 0.7861072 -0.7079146 0.1926596 0.02627752
#> -0.03457033 0.07293042 0.006685215 0.06554499 0.198487 0.1940796 0.2114569
#> 0.1037222 -0.02241764 0.1438094 1.021433 0.3960568 0.3670106 0.7892017
#> 0.2380344 0.7865216
#>
#> Doing end processing
#> Successful trial
#> Total Time: 9.705725
#> Backend Time: 9.692163
summary(fit)
#> Coefficients:
#> Estimate Std. Error t value ci.lower ci.upper Pr(>|t|)
#> phi_11 -0.334601 0.094078 -3.557 -0.518991 -0.150211 0.0002 ***
#> phi_12 0.021663 0.090875 0.238 -0.156448 0.199774 0.4058
#> phi_13 -0.019333 0.074466 -0.260 -0.165283 0.126617 0.3976
#> phi_21 0.770683 0.068892 11.187 0.635656 0.905710 <2e-16 ***
#> phi_22 -0.550716 0.067048 -8.214 -0.682128 -0.419305 <2e-16 ***
#> phi_23 0.008072 0.053580 0.151 -0.096942 0.113087 0.4401
#> phi_31 -0.507772 0.068188 -7.447 -0.641419 -0.374126 <2e-16 ***
#> phi_32 0.786107 0.066051 11.902 0.656650 0.915565 <2e-16 ***
#> phi_33 -0.707915 0.053656 -13.194 -0.813078 -0.602751 <2e-16 ***
#> sigma_11 0.192660 0.019483 9.888 0.154473 0.230846 <2e-16 ***
#> sigma_12 0.026278 0.009713 2.705 0.007239 0.045316 0.0034 **
#> sigma_13 -0.034570 0.009843 -3.512 -0.053861 -0.015279 0.0002 ***
#> sigma_22 0.072930 0.009634 7.570 0.054048 0.091813 <2e-16 ***
#> sigma_23 0.006685 0.006742 0.992 -0.006529 0.019900 0.1607
#> sigma_33 0.065545 0.009791 6.695 0.046356 0.084734 <2e-16 ***
#> theta_11 0.198487 0.004542 43.698 0.189584 0.207390 <2e-16 ***
#> theta_22 0.194080 0.004240 45.769 0.185769 0.202391 <2e-16 ***
#> theta_33 0.211457 0.004604 45.933 0.202434 0.220480 <2e-16 ***
#> mu0_1 0.103722 0.145704 0.712 -0.181853 0.389297 0.2383
#> mu0_2 -0.022418 0.128088 -0.175 -0.273465 0.228630 0.4305
#> mu0_3 0.143809 0.128032 1.123 -0.107129 0.394747 0.1307
#> sigma0_11 1.021433 0.211348 4.833 0.607199 1.435667 <2e-16 ***
#> sigma0_12 0.396057 0.143821 2.754 0.114173 0.677940 0.0030 **
#> sigma0_13 0.367011 0.143034 2.566 0.086669 0.647353 0.0052 **
#> sigma0_22 0.789202 0.165089 4.780 0.465632 1.112771 <2e-16 ***
#> sigma0_23 0.238034 0.122124 1.949 -0.001324 0.477393 0.0257 *
#> sigma0_33 0.786522 0.166425 4.726 0.460334 1.112709 <2e-16 ***
#> ---
#> Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
#>
#> -2 log-likelihood value at convergence = 21100.40
#> AIC = 21154.40
#> BIC = 21346.24
phi_hat <- PhiHat(fit)
theta_hat <- ThetaHat(fit)
Confidence Intervals (Unstandardized)
Model
Delta Method
summary(delta)
#> effect interval est se z p
#> 1 total 1 -0.1256906920 0.0308918759 -4.0687297 4.727015e-05
#> 2 direct 1 -0.3037507224 0.0419975788 -7.2325770 4.739138e-13
#> 3 indirect 1 0.1780600303 0.0200807224 8.8672124 7.499978e-19
#> 4 total 2 0.0547873238 0.0335737093 1.6318520 1.027107e-01
#> 5 direct 2 -0.3688024194 0.0596256657 -6.1852965 6.198587e-10
#> 6 indirect 2 0.4235897432 0.0505670167 8.3767992 5.438654e-17
#> 7 total 3 0.2322234227 0.0336008193 6.9112429 4.804257e-12
#> 8 direct 3 -0.3406933403 0.0677448446 -5.0290667 4.928730e-07
#> 9 indirect 3 0.5729167630 0.0760667805 7.5317604 5.006077e-14
#> 10 total 4 0.3345715006 0.0364851142 9.1700823 4.726579e-20
#> 11 direct 4 -0.2836293280 0.0696159090 -4.0742027 4.617224e-05
#> 12 indirect 4 0.6182008285 0.0919204611 6.7253887 1.751243e-11
#> 13 total 5 0.3671226357 0.0402454180 9.1220977 7.368405e-20
#> 14 direct 5 -0.2242547245 0.0670821528 -3.3429864 8.288196e-04
#> 15 indirect 5 0.5913773602 0.0979402040 6.0381471 1.558939e-09
#> 16 total 6 0.3530953009 0.0427896862 8.2518787 1.559236e-16
#> 17 direct 6 -0.1722832310 0.0617782365 -2.7887366 5.291408e-03
#> 18 indirect 6 0.5253785319 0.0960973509 5.4671490 4.573317e-08
#> 19 total 7 0.3140492805 0.0435718568 7.2076176 5.693927e-13
#> 20 direct 7 -0.1301149456 0.0550230922 -2.3647334 1.804305e-02
#> 21 indirect 7 0.4441642261 0.0889481801 4.9935167 5.928964e-07
#> 22 total 8 0.2652310175 0.0429392412 6.1768911 6.537613e-10
#> 23 direct 8 -0.0972392485 0.0477684656 -2.0356368 4.178684e-02
#> 24 indirect 8 0.3624702660 0.0789137467 4.5932462 4.364037e-06
#> 25 total 9 0.2159203686 0.0414880526 5.2043987 1.946257e-07
#> 26 direct 9 -0.0721899710 0.0406421182 -1.7762354 7.569412e-02
#> 27 indirect 9 0.2881103396 0.0679754135 4.2384492 2.250691e-05
#> 28 total 10 0.1710189987 0.0396745377 4.3105480 1.628505e-05
#> 29 direct 10 -0.0533669745 0.0340204734 -1.5686723 1.167243e-01
#> 30 indirect 10 0.2243859733 0.0575411679 3.8995728 9.636256e-05
#> 31 total 11 0.1325971732 0.0376886068 3.5182296 4.344363e-04
#> 32 direct 11 -0.0393443549 0.0280990247 -1.4002036 1.614524e-01
#> 33 indirect 11 0.1719415281 0.0484196097 3.5510722 3.836652e-04
#> 34 total 12 0.1010631543 0.0355250643 2.8448409 4.443362e-03
#> 35 direct 12 -0.0289549896 0.0229505143 -1.2616270 2.070830e-01
#> 36 indirect 12 0.1300181439 0.0408869663 3.1799411 1.473050e-03
#> 37 total 13 0.0759493417 0.0331202091 2.2931420 2.183983e-02
#> 38 direct 13 -0.0212845291 0.0185694334 -1.1462132 2.517070e-01
#> 39 indirect 13 0.0972338708 0.0348314757 2.7915519 5.245595e-03
#> 40 total 14 0.0564004968 0.0304478777 1.8523622 6.397383e-02
#> 41 direct 14 -0.0156342942 0.0149043087 -1.0489782 2.941882e-01
#> 42 indirect 14 0.0720347910 0.0299402769 2.4059494 1.613050e-02
#> 43 total 15 0.0414557535 0.0275480675 1.5048516 1.323622e-01
#> 44 direct 15 -0.0114783446 0.0118800638 -0.9661854 3.339514e-01
#> 45 indirect 15 0.0529340981 0.0258704656 2.0461208 4.074449e-02
#> 46 total 16 0.0301977504 0.0245117831 1.2319687 2.179608e-01
#> 47 direct 16 -0.0084244345 0.0094127759 -0.8950000 3.707871e-01
#> 48 indirect 16 0.0386221849 0.0223531169 1.7278210 8.402032e-02
#> 49 total 17 0.0218209161 0.0214519774 1.0171984 3.090591e-01
#> 50 direct 17 -0.0061817434 0.0074188690 -0.8332461 4.047060e-01
#> 51 indirect 17 0.0280026595 0.0192220522 1.4567986 1.451720e-01
#> 52 total 18 0.0156533681 0.0184781162 0.8471301 3.969226e-01
#> 53 direct 18 -0.0045354608 0.0058204007 -0.7792352 4.358412e-01
#> 54 indirect 18 0.0201888288 0.0163955957 1.2313568 2.181894e-01
#> 55 total 19 0.0111541299 0.0156802289 0.7113499 4.768674e-01
#> 56 direct 19 -0.0033273058 0.0045477173 -0.7316431 4.643865e-01
#> 57 indirect 19 0.0144814357 0.0138439814 1.0460456 2.955400e-01
#> 58 total 20 0.0078987343 0.0131220554 0.6019434 5.472118e-01
#> 59 direct 20 -0.0024408347 0.0035404148 -0.6894206 4.905587e-01
#> 60 indirect 20 0.0103395690 0.0115611877 0.8943345 3.711429e-01
#> 61 total 21 0.0055606962 0.0108407123 0.5129456 6.079893e-01
#> 62 direct 21 -0.0017904707 0.0027472658 -0.6517282 5.145765e-01
#> 63 indirect 21 0.0073511669 0.0095467789 0.7700154 4.412907e-01
#> 64 total 22 0.0038928671 0.0088501452 0.4398648 6.600351e-01
#> 65 direct 22 -0.0013133639 0.0021255685 -0.6178883 5.366490e-01
#> 66 indirect 22 0.0052062310 0.0077966995 0.6677481 5.042944e-01
#> 67 total 23 0.0027106046 0.0071462364 0.3793052 7.044612e-01
#> 68 direct 23 -0.0009633757 0.0016402113 -0.5873485 5.569697e-01
#> 69 indirect 23 0.0036739804 0.0063001179 0.5831606 5.597852e-01
#> 70 total 24 0.0018774952 0.0057121589 0.3286840 7.423946e-01
#> 71 direct 24 -0.0007066455 0.0012626439 -0.5596554 5.757145e-01
#> 72 indirect 24 0.0025841407 0.0050396514 0.5127618 6.081180e-01
#> 73 total 25 0.0012937263 0.0045231842 0.2860211 7.748619e-01
#> 74 direct 25 -0.0005183277 0.0009698639 -0.5344334 5.930417e-01
#> 75 indirect 25 0.0018120540 0.0039931039 0.4537959 6.499758e-01
#> 76 total 26 0.0008868868 0.0035505859 0.2497860 8.027528e-01
#> 77 direct 26 -0.0003801939 0.0007434821 -0.5113693 6.090925e-01
#> 78 indirect 26 0.0012670807 0.0031355954 0.4040957 6.861423e-01
#> 79 total 27 0.0006048452 0.0027645621 0.2187852 8.268174e-01
#> 80 direct 27 -0.0002788718 0.0005688940 -0.4902000 6.239924e-01
#> 81 indirect 27 0.0008837170 0.0024415006 0.3619565 7.173845e-01
#> 82 total 28 0.0004103317 0.0021362558 0.1920799 8.476797e-01
#> 83 direct 28 -0.0002045518 0.0004345666 -0.4707029 6.378529e-01
#> 84 indirect 28 0.0006148835 0.0018859565 0.3260327 7.443996e-01
#> 85 total 29 0.0002768726 0.0016390223 0.1689255 8.658553e-01
#> 86 direct 29 -0.0001500380 0.0003314373 -0.4526889 6.507728e-01
#> 87 indirect 29 0.0004269106 0.0014458914 0.2952577 7.677970e-01
#> 88 total 30 0.0001857764 0.0012491141 0.1487265 8.817694e-01
#> 89 direct 30 -0.0001100522 0.0002524158 -0.4359959 6.628397e-01
#> 90 indirect 30 0.0002958286 0.0011006275 0.2687818 7.880976e-01
#> 2.5% 97.5%
#> 1 -0.1862376561 -0.0651437279
#> 2 -0.3860644642 -0.2214369806
#> 3 0.1387025376 0.2174175231
#> 4 -0.0110159372 0.1205905847
#> 5 -0.4856665768 -0.2519382620
#> 6 0.3244802116 0.5226992748
#> 7 0.1663670270 0.2980798183
#> 8 -0.4734707959 -0.2079158848
#> 9 0.4238286128 0.7220049132
#> 10 0.2630619908 0.4060810103
#> 11 -0.4200740024 -0.1471846535
#> 12 0.4380400353 0.7983616218
#> 13 0.2882430658 0.4460022056
#> 14 -0.3557333280 -0.0927761209
#> 15 0.3994180877 0.7833366326
#> 16 0.2692290571 0.4369615447
#> 17 -0.2933663496 -0.0512001125
#> 18 0.3370311851 0.7137258787
#> 19 0.2286500105 0.3994485506
#> 20 -0.2379582246 -0.0222716666
#> 21 0.2698289966 0.6184994556
#> 22 0.1810716513 0.3493903837
#> 23 -0.1908637207 -0.0036147763
#> 24 0.2078021645 0.5171383675
#> 25 0.1346052798 0.2972354575
#> 26 -0.1518470588 0.0074671169
#> 27 0.1548809773 0.4213397018
#> 28 0.0932583337 0.2487796638
#> 29 -0.1200458770 0.0133119280
#> 30 0.1116073566 0.3371645899
#> 31 0.0587288612 0.2064654852
#> 32 -0.0944174314 0.0157287215
#> 33 0.0770408370 0.2668422193
#> 34 0.0314353077 0.1706910009
#> 35 -0.0739371711 0.0160271920
#> 36 0.0498811626 0.2101551252
#> 37 0.0110349248 0.1408637586
#> 38 -0.0576799497 0.0151108916
#> 39 0.0289654329 0.1655023086
#> 40 -0.0032762469 0.1160772405
#> 41 -0.0448462024 0.0135776140
#> 42 0.0133529265 0.1307166555
#> 43 -0.0125374666 0.0954489737
#> 44 -0.0347628417 0.0118061525
#> 45 0.0022289173 0.1036392789
#> 46 -0.0178444617 0.0782399625
#> 47 -0.0268731363 0.0100242673
#> 48 -0.0051891191 0.0824334889
#> 49 -0.0202241871 0.0638660192
#> 50 -0.0207224595 0.0083589727
#> 51 -0.0096718705 0.0656771894
#> 52 -0.0205630742 0.0518698103
#> 53 -0.0159432365 0.0068723149
#> 54 -0.0119459483 0.0523236060
#> 55 -0.0195785539 0.0418868137
#> 56 -0.0122406680 0.0055860564
#> 57 -0.0126522692 0.0416151406
#> 58 -0.0178200217 0.0336174903
#> 59 -0.0093799202 0.0044982507
#> 60 -0.0123199425 0.0329990806
#> 61 -0.0156867095 0.0268081019
#> 62 -0.0071750128 0.0035940714
#> 63 -0.0113601758 0.0260625096
#> 64 -0.0134530989 0.0212388330
#> 65 -0.0054794016 0.0028526739
#> 66 -0.0100750192 0.0204874811
#> 67 -0.0112957613 0.0167169706
#> 68 -0.0041781309 0.0022513794
#> 69 -0.0086740239 0.0160219846
#> 70 -0.0093181305 0.0130731209
#> 71 -0.0031813821 0.0017680911
#> 72 -0.0072933946 0.0124616760
#> 73 -0.0075715518 0.0101590045
#> 74 -0.0024192260 0.0013825706
#> 75 -0.0060142858 0.0096383939
#> 76 -0.0060721337 0.0078459072
#> 77 -0.0018373921 0.0010770042
#> 78 -0.0048785733 0.0074127347
#> 79 -0.0048135969 0.0060232874
#> 80 -0.0013938836 0.0008361399
#> 81 -0.0039015362 0.0056689703
#> 82 -0.0037766528 0.0045973162
#> 83 -0.0010562866 0.0006471831
#> 84 -0.0030815232 0.0043112902
#> 85 -0.0029355520 0.0034892973
#> 86 -0.0007996432 0.0004995672
#> 87 -0.0024069844 0.0032608057
#> 88 -0.0022624422 0.0026339950
#> 89 -0.0006047780 0.0003846736
#> 90 -0.0018613617 0.0024530190
Monte Carlo Method
summary(mc)
#> effect interval est se R 2.5%
#> 1 total 1 -0.1256906920 0.030918568 20000 -0.1869791011
#> 2 direct 1 -0.3037507224 0.041980003 20000 -0.3886010498
#> 3 indirect 1 0.1780600303 0.020059379 20000 0.1408047197
#> 4 total 2 0.0547873238 0.033876655 20000 -0.0142993962
#> 5 direct 2 -0.3688024194 0.060525370 20000 -0.4985977091
#> 6 indirect 2 0.4235897432 0.051025024 20000 0.3329876718
#> 7 total 3 0.2322234227 0.033838788 20000 0.1655282635
#> 8 direct 3 -0.3406933403 0.070165575 20000 -0.5008906509
#> 9 indirect 3 0.5729167630 0.077894454 20000 0.4426513533
#> 10 total 4 0.3345715006 0.036692167 20000 0.2663046973
#> 11 direct 4 -0.2836293280 0.073990654 20000 -0.4616496790
#> 12 indirect 4 0.6182008285 0.095851715 20000 0.4669791212
#> 13 total 5 0.3671226357 0.040667760 20000 0.2953307719
#> 14 direct 5 -0.2242547245 0.073644527 20000 -0.4096828571
#> 15 indirect 5 0.5913773602 0.104307653 20000 0.4363696913
#> 16 total 6 0.3530953009 0.043558468 20000 0.2796447652
#> 17 direct 6 -0.1722832310 0.070550034 20000 -0.3558646636
#> 18 indirect 6 0.5253785319 0.104828572 20000 0.3770350042
#> 19 total 7 0.3140492805 0.044699280 20000 0.2403623115
#> 20 direct 7 -0.1301149456 0.065853809 20000 -0.3067056768
#> 21 indirect 7 0.4441642261 0.099677470 20000 0.3107140101
#> 22 total 8 0.2652310175 0.044383404 20000 0.1943364233
#> 23 direct 8 -0.0972392485 0.060393567 20000 -0.2652269090
#> 24 indirect 8 0.3624702660 0.091106681 20000 0.2471933304
#> 25 total 9 0.2159203686 0.043196687 20000 0.1477593900
#> 26 direct 9 -0.0721899710 0.054737440 20000 -0.2280769702
#> 27 indirect 9 0.2881103396 0.081028641 20000 0.1900110111
#> 28 total 10 0.1710189987 0.041607039 20000 0.1060752703
#> 29 direct 10 -0.0533669745 0.049244682 20000 -0.1941477460
#> 30 indirect 10 0.2243859733 0.070859302 20000 0.1424503771
#> 31 total 11 0.1325971732 0.039820924 20000 0.0703529005
#> 32 direct 11 -0.0393443549 0.044123282 20000 -0.1650168399
#> 33 indirect 11 0.1719415281 0.061482017 20000 0.1033058716
#> 34 total 12 0.1010631543 0.037843040 20000 0.0407332048
#> 35 direct 12 -0.0289549896 0.039477163 20000 -0.1415685108
#> 36 indirect 12 0.1300181439 0.053309608 20000 0.0720594168
#> 37 total 13 0.0759493417 0.035609784 20000 0.0197460764
#> 38 direct 13 -0.0212845291 0.035342138 20000 -0.1212438690
#> 39 indirect 13 0.0972338708 0.046411235 20000 0.0477917921
#> 40 total 14 0.0564004968 0.033088958 20000 0.0045043655
#> 41 direct 14 -0.0156342942 0.031712003 20000 -0.1042632568
#> 42 indirect 14 0.0720347910 0.040660294 20000 0.0302870787
#> 43 total 15 0.0414557535 0.030313856 20000 -0.0055767877
#> 44 direct 15 -0.0114783446 0.028556633 20000 -0.0893071096
#> 45 indirect 15 0.0529340981 0.035861114 20000 0.0177769254
#> 46 total 16 0.0301977504 0.027371757 20000 -0.0117410907
#> 47 direct 16 -0.0084244345 0.025834034 20000 -0.0768759588
#> 48 indirect 16 0.0386221849 0.031828950 20000 0.0092249428
#> 49 total 17 0.0218209161 0.024375851 20000 -0.0143911975
#> 50 direct 17 -0.0061817434 0.023497967 20000 -0.0659026462
#> 51 indirect 17 0.0280026595 0.028421910 20000 0.0036407020
#> 52 total 18 0.0156533681 0.021438886 20000 -0.0147688384
#> 53 direct 18 -0.0045354608 0.021502508 20000 -0.0567494811
#> 54 indirect 18 0.0201888288 0.025540439 20000 0.0002997946
#> 55 total 19 0.0111541299 0.018655516 20000 -0.0140381821
#> 56 direct 19 -0.0033273058 0.019804554 20000 -0.0485077374
#> 57 indirect 19 0.0144814357 0.023113174 20000 -0.0015950536
#> 58 total 20 0.0078987343 0.016093820 20000 -0.0123533113
#> 59 direct 20 -0.0024408347 0.018365011 20000 -0.0417012187
#> 60 indirect 20 0.0103395690 0.021082378 20000 -0.0025279416
#> 61 total 21 0.0055606962 0.013793897 20000 -0.0104138547
#> 62 direct 21 -0.0017904707 0.017149186 20000 -0.0356180805
#> 63 indirect 21 0.0073511669 0.019394842 20000 -0.0027764958
#> 64 total 22 0.0038928671 0.011770954 20000 -0.0086014033
#> 65 direct 22 -0.0013133639 0.016126716 20000 -0.0305485402
#> 66 indirect 22 0.0052062310 0.017998912 20000 -0.0026262933
#> 67 total 23 0.0027106046 0.010020667 20000 -0.0069171159
#> 68 direct 23 -0.0009633757 0.015271258 20000 -0.0260906142
#> 69 indirect 23 0.0036739804 0.016845602 20000 -0.0022990453
#> 70 total 24 0.0018774952 0.008525186 20000 -0.0053513231
#> 71 direct 24 -0.0007066455 0.014560084 20000 -0.0224467596
#> 72 indirect 24 0.0025841407 0.015891172 20000 -0.0019372122
#> 73 total 25 0.0012937263 0.007258729 20000 -0.0041484145
#> 74 direct 25 -0.0005183277 0.013973644 20000 -0.0193071134
#> 75 indirect 25 0.0018120540 0.015099174 20000 -0.0016187147
#> 76 total 26 0.0008868868 0.006192160 20000 -0.0031654324
#> 77 direct 26 -0.0003801939 0.013495153 20000 -0.0165773195
#> 78 indirect 26 0.0012670807 0.014441113 20000 -0.0012513974
#> 79 total 27 0.0006048452 0.005296305 20000 -0.0023506567
#> 80 direct 27 -0.0002788718 0.013110213 20000 -0.0142498758
#> 81 indirect 27 0.0008837170 0.013895784 20000 -0.0009730558
#> 82 total 28 0.0004103317 0.004544056 20000 -0.0017255010
#> 83 direct 28 -0.0002045518 0.012806482 20000 -0.0122217295
#> 84 indirect 28 0.0006148835 0.013447819 20000 -0.0007458488
#> 85 total 29 0.0002768726 0.003911442 20000 -0.0012501279
#> 86 direct 29 -0.0001500380 0.012573387 20000 -0.0105101985
#> 87 indirect 29 0.0004269106 0.013085977 20000 -0.0005607599
#> 88 total 30 0.0001857764 0.003377964 20000 -0.0009091101
#> 89 direct 30 -0.0001100522 0.012401875 20000 -0.0090545986
#> 90 indirect 30 0.0002958286 0.012801618 20000 -0.0004141329
#> 97.5%
#> 1 -6.646446e-02
#> 2 -2.233577e-01
#> 3 2.199914e-01
#> 4 1.176580e-01
#> 5 -2.614643e-01
#> 6 5.333207e-01
#> 7 2.986343e-01
#> 8 -2.272604e-01
#> 9 7.478701e-01
#> 10 4.109161e-01
#> 11 -1.743764e-01
#> 12 8.395155e-01
#> 13 4.544647e-01
#> 14 -1.247865e-01
#> 15 8.428931e-01
#> 16 4.497320e-01
#> 17 -8.609432e-02
#> 18 7.846319e-01
#> 19 4.149341e-01
#> 20 -5.700074e-02
#> 21 6.964449e-01
#> 22 3.674155e-01
#> 23 -3.720898e-02
#> 24 5.988217e-01
#> 25 3.170367e-01
#> 26 -2.374214e-02
#> 27 5.010845e-01
#> 28 2.686865e-01
#> 29 -1.482227e-02
#> 30 4.118206e-01
#> 31 2.259415e-01
#> 32 -9.171390e-03
#> 33 3.354873e-01
#> 34 1.895172e-01
#> 35 -5.625974e-03
#> 36 2.720115e-01
#> 37 1.594271e-01
#> 38 -3.377325e-03
#> 39 2.209714e-01
#> 40 1.342035e-01
#> 41 -1.990189e-03
#> 42 1.806128e-01
#> 43 1.136657e-01
#> 44 -1.167150e-03
#> 45 1.481355e-01
#> 46 9.625064e-02
#> 47 -6.656567e-04
#> 48 1.227296e-01
#> 49 8.211953e-02
#> 50 -3.831208e-04
#> 51 1.024684e-01
#> 52 6.989265e-02
#> 53 -2.094467e-04
#> 54 8.650737e-02
#> 55 5.939192e-02
#> 56 -1.122974e-04
#> 57 7.354663e-02
#> 58 5.070657e-02
#> 59 -6.050309e-05
#> 60 6.317263e-02
#> 61 4.330184e-02
#> 62 -3.144843e-05
#> 63 5.404042e-02
#> 64 3.702237e-02
#> 65 -1.555783e-05
#> 66 4.669262e-02
#> 67 3.160801e-02
#> 68 -6.704140e-06
#> 69 4.051071e-02
#> 70 2.686572e-02
#> 71 -2.943748e-06
#> 72 3.532881e-02
#> 73 2.302549e-02
#> 74 -1.112480e-06
#> 75 3.081770e-02
#> 76 1.963603e-02
#> 77 -2.234064e-07
#> 78 2.673677e-02
#> 79 1.678864e-02
#> 80 9.555190e-08
#> 81 2.332481e-02
#> 82 1.438135e-02
#> 83 1.382529e-07
#> 84 2.041860e-02
#> 85 1.232948e-02
#> 86 1.276054e-07
#> 87 1.788115e-02
#> 88 1.055287e-02
#> 89 9.798301e-08
#> 90 1.550820e-02
Parametric Bootstrap Method
B <- 1000L
path <- root$find_file(
".setup",
"data-raw"
)
prefix <- "pb_replication"
boot <- BootPara(
fit = fit,
path = path,
prefix = prefix,
taskid = taskid,
B = B,
ncores = parallel::detectCores(),
seed = seed,
clean = FALSE
)
pb <- BootParaXMY(boot = boot, phi_hat = phi_hat)
plot(pb)
plot(pb, type = "bc")
summary(pb)
#> effect interval est se R 2.5% 97.5%
#> 1 total 1 -0.1256906920 0.028595885 1000 -0.1807940111 -7.056570e-02
#> 2 direct 1 -0.3037507224 0.038429144 1000 -0.3780736134 -2.285441e-01
#> 3 indirect 1 0.1780600303 0.018319975 1000 0.1441372577 2.148714e-01
#> 4 total 2 0.0547873238 0.031265118 1000 -0.0057722660 1.132452e-01
#> 5 direct 2 -0.3688024194 0.056117382 1000 -0.4865857228 -2.627405e-01
#> 6 indirect 2 0.4235897432 0.047433298 1000 0.3389181538 5.176740e-01
#> 7 total 3 0.2322234227 0.031077076 1000 0.1737472146 2.905343e-01
#> 8 direct 3 -0.3406933403 0.065191207 1000 -0.4772649743 -2.231497e-01
#> 9 indirect 3 0.5729167630 0.072897307 1000 0.4444458814 7.191563e-01
#> 10 total 4 0.3345715006 0.034198991 1000 0.2710307742 3.985612e-01
#> 11 direct 4 -0.2836293280 0.068324125 1000 -0.4308255198 -1.691142e-01
#> 12 indirect 4 0.6182008285 0.089533921 1000 0.4659064441 8.059694e-01
#> 13 total 5 0.3671226357 0.038364561 1000 0.2966544875 4.430815e-01
#> 14 direct 5 -0.2242547245 0.067255204 1000 -0.3793180016 -1.198886e-01
#> 15 indirect 5 0.5913773602 0.096695642 1000 0.4324062758 8.000867e-01
#> 16 total 6 0.3530953009 0.041074149 1000 0.2797465782 4.359995e-01
#> 17 direct 6 -0.1722832310 0.063523672 1000 -0.3275100846 -7.989026e-02
#> 18 indirect 6 0.5253785319 0.096018276 1000 0.3710340307 7.409865e-01
#> 19 total 7 0.3140492805 0.041770105 1000 0.2368084937 4.042859e-01
#> 20 direct 7 -0.1301149456 0.058336995 1000 -0.2781046097 -5.240774e-02
#> 21 indirect 7 0.4441642261 0.089874326 1000 0.3035324863 6.541076e-01
#> 22 total 8 0.2652310175 0.040971892 1000 0.1892984638 3.477480e-01
#> 23 direct 8 -0.0972392485 0.052547697 1000 -0.2345875915 -3.344158e-02
#> 24 indirect 8 0.3624702660 0.080623768 1000 0.2396918059 5.543130e-01
#> 25 total 9 0.2159203686 0.039443132 1000 0.1425766422 2.972618e-01
#> 26 direct 9 -0.0721899710 0.046711342 1000 -0.1970357712 -2.110724e-02
#> 27 indirect 9 0.2881103396 0.070249738 1000 0.1839000109 4.577474e-01
#> 28 total 10 0.1710189987 0.037718431 1000 0.0999825191 2.493221e-01
#> 29 direct 10 -0.0533669745 0.041161797 1000 -0.1662846507 -1.272641e-02
#> 30 indirect 10 0.2243859733 0.060182235 1000 0.1372070095 3.706810e-01
#> 31 total 11 0.1325971732 0.035973146 1000 0.0681688618 2.115010e-01
#> 32 direct 11 -0.0393443549 0.036078147 1000 -0.1395318915 -7.852610e-03
#> 33 indirect 11 0.1719415281 0.051258508 1000 0.1009940851 2.909715e-01
#> 34 total 12 0.1010631543 0.034139532 1000 0.0399098215 1.761026e-01
#> 35 direct 12 -0.0289549896 0.031536989 1000 -0.1173294331 -4.922461e-03
#> 36 indirect 12 0.1300181439 0.043793903 1000 0.0708004445 2.303240e-01
#> 37 total 13 0.0759493417 0.032091675 1000 0.0171257119 1.467951e-01
#> 38 direct 13 -0.0212845291 0.027550613 1000 -0.0986147630 -2.942374e-03
#> 39 indirect 13 0.0972338708 0.037729903 1000 0.0460181029 1.863131e-01
#> 40 total 14 0.0564004968 0.029764229 1000 0.0029638546 1.231369e-01
#> 41 direct 14 -0.0156342942 0.024093614 1000 -0.0828604852 -1.744862e-03
#> 42 indirect 14 0.0720347910 0.032813268 1000 0.0276702102 1.483906e-01
#> 43 total 15 0.0414557535 0.027180864 1000 -0.0065033007 1.036075e-01
#> 44 direct 15 -0.0114783446 0.021120660 1000 -0.0696024198 -1.029370e-03
#> 45 indirect 15 0.0529340981 0.028751481 1000 0.0155096651 1.221016e-01
#> 46 total 16 0.0301977504 0.024430563 1000 -0.0109195633 8.723559e-02
#> 47 direct 16 -0.0084244345 0.018577890 1000 -0.0584479657 -5.917611e-04
#> 48 indirect 16 0.0386221849 0.025305438 1000 0.0077752998 1.000256e-01
#> 49 total 17 0.0218209161 0.021630587 1000 -0.0136222947 7.395583e-02
#> 50 direct 17 -0.0061817434 0.016409883 1000 -0.0490786241 -3.204353e-04
#> 51 indirect 17 0.0280026595 0.022316992 1000 0.0021731778 8.334176e-02
#> 52 total 18 0.0156533681 0.018895948 1000 -0.0125409207 6.346468e-02
#> 53 direct 18 -0.0045354608 0.014563670 1000 -0.0413712057 -1.497913e-04
#> 54 indirect 18 0.0201888288 0.019696093 1000 -0.0011467484 7.127164e-02
#> 55 total 19 0.0111541299 0.016320824 1000 -0.0115574099 5.304095e-02
#> 56 direct 19 -0.0033273058 0.012990832 1000 -0.0346224618 -8.680201e-05
#> 57 indirect 19 0.0144814357 0.017394778 1000 -0.0024859468 5.954485e-02
#> 58 total 20 0.0078987343 0.013970621 1000 -0.0103768299 4.515121e-02
#> 59 direct 20 -0.0024408347 0.011648406 1000 -0.0290744621 -5.024221e-05
#> 60 indirect 20 0.0103395690 0.015384239 1000 -0.0028839802 4.965729e-02
#> 61 total 21 0.0055606962 0.011881511 1000 -0.0088196623 3.851327e-02
#> 62 direct 21 -0.0017904707 0.010499069 1000 -0.0244153449 -2.645781e-05
#> 63 indirect 21 0.0073511669 0.013640373 1000 -0.0027627996 4.213110e-02
#> 64 total 22 0.0038928671 0.010064353 1000 -0.0069973128 3.253887e-02
#> 65 direct 22 -0.0013133639 0.009510915 1000 -0.0205684930 -1.100752e-05
#> 66 indirect 22 0.0052062310 0.012137309 1000 -0.0027309361 3.603524e-02
#> 67 total 23 0.0027106046 0.008510664 1000 -0.0057509337 2.782735e-02
#> 68 direct 23 -0.0009633757 0.008657005 1000 -0.0173321162 -5.296160e-06
#> 69 indirect 23 0.0036739804 0.010846371 1000 -0.0023772351 3.047804e-02
#> 70 total 24 0.0018774952 0.007199029 1000 -0.0046277917 2.345767e-02
#> 71 direct 24 -0.0007066455 0.007914818 1000 -0.0145166022 -2.676330e-06
#> 72 indirect 24 0.0025841407 0.009737758 1000 -0.0020648640 2.630233e-02
#> 73 total 25 0.0012937263 0.006100968 1000 -0.0036572996 1.989734e-02
#> 74 direct 25 -0.0005183277 0.007265677 1000 -0.0122364583 -1.031786e-06
#> 75 indirect 25 0.0018120540 0.008782916 1000 -0.0016044027 2.249668e-02
#> 76 total 26 0.0008868868 0.005185686 1000 -0.0028829073 1.703081e-02
#> 77 direct 26 -0.0003801939 0.006694191 1000 -0.0103853405 -3.578850e-07
#> 78 indirect 26 0.0012670807 0.007956379 1000 -0.0011992584 1.960282e-02
#> 79 total 27 0.0006048452 0.004423461 1000 -0.0021262875 1.441835e-02
#> 80 direct 27 -0.0002788718 0.006187730 1000 -0.0088142578 1.074175e-08
#> 81 indirect 27 0.0008837170 0.007236659 1000 -0.0008620253 1.658593e-02
#> 82 total 28 0.0004103317 0.003787661 1000 -0.0015668789 1.220657e-02
#> 83 direct 28 -0.0002045518 0.005735969 1000 -0.0074805362 3.793469e-08
#> 84 indirect 28 0.0006148835 0.006606259 1000 -0.0006146157 1.412502e-02
#> 85 total 29 0.0002768726 0.003255614 1000 -0.0010595526 1.033390e-02
#> 86 direct 29 -0.0001500380 0.005330475 1000 -0.0063476816 2.310074e-08
#> 87 indirect 29 0.0004269106 0.006051160 1000 -0.0004865328 1.210174e-02
#> 88 total 30 0.0001857764 0.002808647 1000 -0.0007613906 8.747282e-03
#> 89 direct 30 -0.0001100522 0.004964368 1000 -0.0053863958 2.924745e-08
#> 90 indirect 30 0.0002958286 0.005560089 1000 -0.0003585573 1.038221e-02
summary(pb, type = "bc")
#> effect interval est se R 2.5% 97.5%
#> 1 total 1 -0.1256906920 0.028595885 1000 -0.1812175346 -7.274220e-02
#> 2 direct 1 -0.3037507224 0.038429144 1000 -0.3850346753 -2.322560e-01
#> 3 indirect 1 0.1780600303 0.018319975 1000 0.1469516736 2.189493e-01
#> 4 total 2 0.0547873238 0.031265118 1000 -0.0069143044 1.125144e-01
#> 5 direct 2 -0.3688024194 0.056117382 1000 -0.4941467468 -2.715494e-01
#> 6 indirect 2 0.4235897432 0.047433298 1000 0.3465092639 5.297533e-01
#> 7 total 3 0.2322234227 0.031077076 1000 0.1742399805 2.919720e-01
#> 8 direct 3 -0.3406933403 0.065191207 1000 -0.4944338430 -2.326845e-01
#> 9 indirect 3 0.5729167630 0.072897307 1000 0.4561683228 7.379779e-01
#> 10 total 4 0.3345715006 0.034198991 1000 0.2762625279 4.053200e-01
#> 11 direct 4 -0.2836293280 0.068324125 1000 -0.4462073273 -1.750791e-01
#> 12 indirect 4 0.6182008285 0.089533921 1000 0.4751433378 8.225573e-01
#> 13 total 5 0.3671226357 0.038364561 1000 0.3043964537 4.525063e-01
#> 14 direct 5 -0.2242547245 0.067255204 1000 -0.3908028416 -1.239319e-01
#> 15 indirect 5 0.5913773602 0.096695642 1000 0.4415918687 8.331003e-01
#> 16 total 6 0.3530953009 0.041074149 1000 0.2830390994 4.449133e-01
#> 17 direct 6 -0.1722832310 0.063523672 1000 -0.3365904230 -8.443203e-02
#> 18 indirect 6 0.5253785319 0.096018276 1000 0.3802736905 7.635681e-01
#> 19 total 7 0.3140492805 0.041770105 1000 0.2405254210 4.074022e-01
#> 20 direct 7 -0.1301149456 0.058336995 1000 -0.2869661022 -5.580684e-02
#> 21 indirect 7 0.4441642261 0.089874326 1000 0.3104722435 6.666143e-01
#> 22 total 8 0.2652310175 0.040971892 1000 0.1952500442 3.599167e-01
#> 23 direct 8 -0.0972392485 0.052547697 1000 -0.2403402229 -3.548973e-02
#> 24 indirect 8 0.3624702660 0.080623768 1000 0.2452996229 5.604305e-01
#> 25 total 9 0.2159203686 0.039443132 1000 0.1493125107 3.039196e-01
#> 26 direct 9 -0.0721899710 0.046711342 1000 -0.2012244132 -2.188136e-02
#> 27 indirect 9 0.2881103396 0.070249738 1000 0.1844222925 4.591420e-01
#> 28 total 10 0.1710189987 0.037718431 1000 0.1109292611 2.606071e-01
#> 29 direct 10 -0.0533669745 0.041161797 1000 -0.1674789292 -1.322698e-02
#> 30 indirect 10 0.2243859733 0.060182235 1000 0.1365770652 3.699140e-01
#> 31 total 11 0.1325971732 0.035973146 1000 0.0726647926 2.184899e-01
#> 32 direct 11 -0.0393443549 0.036078147 1000 -0.1407240860 -7.962796e-03
#> 33 indirect 11 0.1719415281 0.051258508 1000 0.1011234106 2.913109e-01
#> 34 total 12 0.1010631543 0.034139532 1000 0.0433791361 1.804081e-01
#> 35 direct 12 -0.0289549896 0.031536989 1000 -0.1201078240 -5.006815e-03
#> 36 indirect 12 0.1300181439 0.043793903 1000 0.0679509720 2.289250e-01
#> 37 total 13 0.0759493417 0.032091675 1000 0.0200524875 1.519540e-01
#> 38 direct 13 -0.0212845291 0.027550613 1000 -0.1005396347 -2.948745e-03
#> 39 indirect 13 0.0972338708 0.037729903 1000 0.0416574655 1.802592e-01
#> 40 total 14 0.0564004968 0.029764229 1000 0.0048287139 1.276160e-01
#> 41 direct 14 -0.0156342942 0.024093614 1000 -0.0828604852 -1.744862e-03
#> 42 indirect 14 0.0720347910 0.032813268 1000 0.0239998397 1.404121e-01
#> 43 total 15 0.0414557535 0.027180864 1000 -0.0034906274 1.076431e-01
#> 44 direct 15 -0.0114783446 0.021120660 1000 -0.0689281045 -1.021689e-03
#> 45 indirect 15 0.0529340981 0.028751481 1000 0.0109763926 1.108846e-01
#> 46 total 16 0.0301977504 0.024430563 1000 -0.0100029915 9.168399e-02
#> 47 direct 16 -0.0084244345 0.018577890 1000 -0.0573282279 -5.737258e-04
#> 48 indirect 16 0.0386221849 0.025305438 1000 0.0037398060 9.011877e-02
#> 49 total 17 0.0218209161 0.021630587 1000 -0.0111904103 7.668001e-02
#> 50 direct 17 -0.0061817434 0.016409883 1000 -0.0487631201 -2.937516e-04
#> 51 indirect 17 0.0280026595 0.022316992 1000 -0.0017413950 6.896862e-02
#> 52 total 18 0.0156533681 0.018895948 1000 -0.0118968157 6.524309e-02
#> 53 direct 18 -0.0045354608 0.014563670 1000 -0.0409595860 -1.399155e-04
#> 54 indirect 18 0.0201888288 0.019696093 1000 -0.0041287195 5.737093e-02
#> 55 total 19 0.0111541299 0.016320824 1000 -0.0111683019 5.518316e-02
#> 56 direct 19 -0.0033273058 0.012990832 1000 -0.0345659829 -8.435738e-05
#> 57 indirect 19 0.0144814357 0.017394778 1000 -0.0051315722 4.671783e-02
#> 58 total 20 0.0078987343 0.013970621 1000 -0.0099858689 4.644881e-02
#> 59 direct 20 -0.0024408347 0.011648406 1000 -0.0289663002 -4.894182e-05
#> 60 indirect 20 0.0103395690 0.015384239 1000 -0.0052009495 3.796733e-02
#> 61 total 21 0.0055606962 0.011881511 1000 -0.0084124261 4.060952e-02
#> 62 direct 21 -0.0017904707 0.010499069 1000 -0.0243646642 -2.575880e-05
#> 63 indirect 21 0.0073511669 0.013640373 1000 -0.0051512918 2.922659e-02
#> 64 total 22 0.0038928671 0.010064353 1000 -0.0067300505 3.425890e-02
#> 65 direct 22 -0.0013133639 0.009510915 1000 -0.0206825802 -1.286076e-05
#> 66 indirect 22 0.0052062310 0.012137309 1000 -0.0045500326 2.345156e-02
#> 67 total 23 0.0027106046 0.008510664 1000 -0.0055697567 2.884226e-02
#> 68 direct 23 -0.0009633757 0.008657005 1000 -0.0174292716 -6.049663e-06
#> 69 indirect 23 0.0036739804 0.010846371 1000 -0.0041204939 1.846173e-02
#> 70 total 24 0.0018774952 0.007199029 1000 -0.0044348104 2.460478e-02
#> 71 direct 24 -0.0007066455 0.007914818 1000 -0.0146805922 -2.731159e-06
#> 72 indirect 24 0.0025841407 0.009737758 1000 -0.0036932608 1.397821e-02
#> 73 total 25 0.0012937263 0.006100968 1000 -0.0034478702 2.111167e-02
#> 74 direct 25 -0.0005183277 0.007265677 1000 -0.0123288201 -1.080470e-06
#> 75 indirect 25 0.0018120540 0.008782916 1000 -0.0029394803 1.128154e-02
#> 76 total 26 0.0008868868 0.005185686 1000 -0.0028790674 1.708208e-02
#> 77 direct 26 -0.0003801939 0.006694191 1000 -0.0103853405 -3.578850e-07
#> 78 indirect 26 0.0012670807 0.007956379 1000 -0.0028498450 8.624406e-03
#> 79 total 27 0.0006048452 0.004423461 1000 -0.0021290434 1.419233e-02
#> 80 direct 27 -0.0002788718 0.006187730 1000 -0.0088142578 1.074175e-08
#> 81 indirect 27 0.0008837170 0.007236659 1000 -0.0021071653 6.857035e-03
#> 82 total 28 0.0004103317 0.003787661 1000 -0.0016338804 1.152157e-02
#> 83 direct 28 -0.0002045518 0.005735969 1000 -0.0073148602 4.872051e-08
#> 84 indirect 28 0.0006148835 0.006606259 1000 -0.0016910715 5.204838e-03
#> 85 total 29 0.0002768726 0.003255614 1000 -0.0012260393 8.697761e-03
#> 86 direct 29 -0.0001500380 0.005330475 1000 -0.0061720434 2.914898e-08
#> 87 indirect 29 0.0004269106 0.006051160 1000 -0.0013920173 3.996803e-03
#> 88 total 30 0.0001857764 0.002808647 1000 -0.0008942805 6.906815e-03
#> 89 direct 30 -0.0001100522 0.004964368 1000 -0.0052079151 3.009917e-08
#> 90 indirect 30 0.0002958286 0.005560089 1000 -0.0011656726 3.052632e-03
Model
Delta Method
summary(delta)
#> effect interval est se z p
#> 1 total 1 -1.155252e-02 4.446122e-02 -0.259833651 0.7949921
#> 2 direct 1 -1.156492e-02 4.435617e-02 -0.260728628 0.7943018
#> 3 indirect 1 1.240083e-05 2.871105e-04 0.043191826 0.9655486
#> 4 total 2 -1.408408e-02 5.438427e-02 -0.258973520 0.7956557
#> 5 direct 2 -1.404168e-02 5.375525e-02 -0.261215083 0.7939266
#> 6 indirect 2 -4.240204e-05 7.889860e-04 -0.053742453 0.9571404
#> 7 total 3 -1.309786e-02 5.082946e-02 -0.257682407 0.7966520
#> 8 direct 3 -1.297147e-02 4.968110e-02 -0.261094556 0.7940196
#> 9 indirect 3 -1.263934e-04 1.327836e-03 -0.095187493 0.9241659
#> 10 total 4 -1.098008e-02 4.280684e-02 -0.256502815 0.7975626
#> 11 direct 4 -1.079883e-02 4.146906e-02 -0.260406798 0.7945500
#> 12 indirect 4 -1.812500e-04 1.667480e-03 -0.108696964 0.9134428
#> 13 total 5 -8.727446e-03 3.411197e-02 -0.255847029 0.7980690
#> 14 direct 5 -8.538213e-03 3.293937e-02 -0.259209979 0.7954732
#> 15 indirect 5 -1.892332e-04 1.773180e-03 -0.106719681 0.9150114
#> 16 total 6 -6.718987e-03 2.624474e-02 -0.256012706 0.7979410
#> 17 direct 6 -6.559465e-03 2.546639e-02 -0.257573412 0.7967362
#> 18 indirect 6 -1.595224e-04 1.746925e-03 -0.091316124 0.9272414
#> 19 total 7 -5.063621e-03 1.968619e-02 -0.257216918 0.7970113
#> 20 direct 7 -4.953961e-03 1.938394e-02 -0.255570410 0.7982826
#> 21 indirect 7 -1.096601e-04 1.698556e-03 -0.064560757 0.9485237
#> 22 total 8 -3.757504e-03 1.447244e-02 -0.259631752 0.7951478
#> 23 direct 8 -3.702261e-03 1.461770e-02 -0.253272389 0.8000577
#> 24 indirect 8 -5.524344e-05 1.674953e-03 -0.032982092 0.9736889
#> 25 total 9 -2.755037e-03 1.045895e-02 -0.263414298 0.7922313
#> 26 direct 9 -2.748541e-03 1.096151e-02 -0.250744746 0.8020115
#> 27 indirect 9 -6.495647e-06 1.663364e-03 -0.003905127 0.9968842
#> 28 total 10 -2.000302e-03 7.443552e-03 -0.268729438 0.7881379
#> 29 direct 10 -2.031880e-03 8.191592e-03 -0.248044544 0.8040999
#> 30 indirect 10 3.157812e-05 1.634075e-03 0.019324772 0.9845820
#> 31 total 11 -1.440230e-03 5.222660e-03 -0.275765575 0.7827281
#> 32 direct 11 -1.497986e-03 6.108752e-03 -0.245219704 0.8062863
#> 33 indirect 11 5.775644e-05 1.569309e-03 0.036803733 0.9706415
#> 34 total 12 -1.029362e-03 3.615066e-03 -0.284742376 0.7758415
#> 35 direct 12 -1.102424e-03 4.549658e-03 -0.242309272 0.8085405
#> 36 indirect 12 7.306199e-05 1.467292e-03 0.049793757 0.9602867
#> 37 total 13 -7.308208e-04 2.469776e-03 -0.295905717 0.7673021
#> 38 direct 13 -8.103814e-04 3.385839e-03 -0.239344304 0.8108386
#> 39 indirect 13 7.956057e-05 1.336152e-03 0.059544528 0.9525184
#> 40 total 14 -5.156803e-04 1.666190e-03 -0.309496631 0.7569438
#> 41 direct 14 -5.952559e-04 2.518546e-03 -0.236349055 0.8131618
#> 42 indirect 14 7.957560e-05 1.187631e-03 0.067003615 0.9465788
#> 43 total 15 -3.617743e-04 1.110882e-03 -0.325664112 0.7446785
#> 44 direct 15 -4.370234e-04 1.872886e-03 -0.233342207 0.8154957
#> 45 indirect 15 7.524910e-05 1.033141e-03 0.072835299 0.9419372
#> 46 total 16 -2.524048e-04 7.331796e-04 -0.344260550 0.7306503
#> 47 direct 16 -3.207496e-04 1.392517e-03 -0.230338007 0.8178291
#> 48 indirect 16 6.834481e-05 8.819036e-04 0.077496918 0.9382282
#> 49 total 17 -1.751620e-04 4.806644e-04 -0.364416483 0.7155470
#> 50 direct 17 -2.353620e-04 1.035253e-03 -0.227347248 0.8201537
#> 51 indirect 17 6.020000e-05 7.404217e-04 0.081305013 0.9351994
#> 52 total 18 -1.209245e-04 3.150730e-04 -0.383798166 0.7011280
#> 53 direct 18 -1.726819e-04 7.696024e-04 -0.224378080 0.8224631
#> 54 indirect 18 5.175744e-05 6.126455e-04 0.084481879 0.9326733
#> 55 total 19 -8.305073e-05 2.088029e-04 -0.397746957 0.6908167
#> 56 direct 19 -1.266829e-04 5.720956e-04 -0.221436639 0.8247525
#> 57 indirect 19 4.363220e-05 5.004473e-04 0.087186399 0.9305234
#> 58 total 20 -5.674485e-05 1.420531e-04 -0.399462326 0.6895526
#> 59 direct 20 -9.293168e-05 4.252630e-04 -0.218527546 0.8270181
#> 60 indirect 20 3.618682e-05 4.041676e-04 0.089534194 0.9286574
#> 61 total 21 -3.856892e-05 1.006545e-04 -0.383181121 0.7015855
#> 62 direct 21 -6.816989e-05 3.161073e-04 -0.215654274 0.8292572
#> 63 indirect 21 2.960097e-05 3.231164e-04 0.091610852 0.9270072
#> 64 total 22 -2.607509e-05 7.465646e-05 -0.349267707 0.7268883
#> 65 direct 22 -5.000466e-05 2.349628e-04 -0.212819432 0.8314678
#> 66 indirect 22 2.392957e-05 2.559840e-04 0.093480722 0.9255217
#> 67 total 23 -1.753128e-05 5.749686e-05 -0.304908540 0.7604358
#> 68 direct 23 -3.667930e-05 1.746426e-04 -0.210024965 0.8336482
#> 69 indirect 23 1.914802e-05 2.011498e-04 0.095192864 0.9241616
#> 70 total 24 -1.171906e-05 4.528805e-05 -0.258767237 0.7958148
#> 71 direct 24 -2.690463e-05 1.298033e-04 -0.207272313 0.8357972
#> 72 indirect 24 1.518557e-05 1.568998e-04 0.096785125 0.9228970
#> 73 total 25 -7.786092e-06 3.597534e-05 -0.216428584 0.8286537
#> 74 direct 25 -1.973467e-05 9.647255e-05 -0.204562516 0.8379139
#> 75 indirect 25 1.194857e-05 1.215682e-04 0.098287000 0.9217044
#> 76 total 26 -5.139335e-06 2.855275e-05 -0.179994379 0.8571570
#> 77 direct 26 -1.447540e-05 7.169720e-05 -0.201896306 0.8399978
#> 78 indirect 26 9.336065e-06 9.362126e-05 0.099721642 0.9205653
#> 79 total 27 -3.368332e-06 2.252667e-05 -0.149526383 0.8811383
#> 80 direct 27 -1.061769e-05 5.328182e-05 -0.199274161 0.8420483
#> 81 indirect 27 7.249358e-06 7.169963e-05 0.101107326 0.9194653
#> 82 total 28 -2.190458e-06 1.762528e-05 -0.124279314 0.9010941
#> 83 direct 28 -7.788048e-06 3.959427e-05 -0.196696358 0.8440652
#> 84 indirect 28 5.597590e-06 5.463276e-05 0.102458493 0.9183928
#> 85 total 29 -1.412110e-06 1.366520e-05 -0.103336224 0.9176961
#> 86 direct 29 -5.712505e-06 2.942118e-05 -0.194163005 0.8460482
#> 87 indirect 29 4.300396e-06 4.143500e-05 0.103786537 0.9173388
#> 88 total 30 -9.013541e-07 1.049844e-05 -0.085856014 0.9315809
#> 89 direct 30 -4.190099e-06 2.186054e-05 -0.191674070 0.8479975
#> 90 indirect 30 3.288745e-06 3.129147e-05 0.105100370 0.9162962
#> 2.5% 97.5%
#> 1 -9.869491e-02 7.558987e-02
#> 2 -9.850141e-02 7.537156e-02
#> 3 -5.503254e-04 5.751270e-04
#> 4 -1.206753e-01 9.250712e-02
#> 5 -1.194000e-01 9.131668e-02
#> 6 -1.588786e-03 1.503982e-03
#> 7 -1.127218e-01 8.652606e-02
#> 8 -1.103446e-01 8.440171e-02
#> 9 -2.728905e-03 2.476118e-03
#> 10 -9.487995e-02 7.291979e-02
#> 11 -9.207669e-02 7.047904e-02
#> 12 -3.449450e-03 3.086950e-03
#> 13 -7.558568e-02 5.813079e-02
#> 14 -7.309819e-02 5.602176e-02
#> 15 -3.664602e-03 3.286135e-03
#> 16 -5.815773e-02 4.471976e-02
#> 17 -5.647267e-02 4.335374e-02
#> 18 -3.583432e-03 3.264387e-03
#> 19 -4.364784e-02 3.352060e-02
#> 20 -4.294578e-02 3.303786e-02
#> 21 -3.438769e-03 3.219449e-03
#> 22 -3.212296e-02 2.460795e-02
#> 23 -3.235243e-02 2.494791e-02
#> 24 -3.338090e-03 3.227603e-03
#> 25 -2.325420e-02 1.774413e-02
#> 26 -2.423271e-02 1.873563e-02
#> 27 -3.266629e-03 3.253638e-03
#> 28 -1.658940e-02 1.258879e-02
#> 29 -1.808710e-02 1.402335e-02
#> 30 -3.171150e-03 3.234306e-03
#> 31 -1.167646e-02 8.795995e-03
#> 32 -1.347092e-02 1.047495e-02
#> 33 -3.018033e-03 3.133546e-03
#> 34 -8.114761e-03 6.056036e-03
#> 35 -1.001959e-02 7.814742e-03
#> 36 -2.802778e-03 2.948902e-03
#> 37 -5.571492e-03 4.109851e-03
#> 38 -7.446504e-03 5.825742e-03
#> 39 -2.539250e-03 2.698371e-03
#> 40 -3.781353e-03 2.749993e-03
#> 41 -5.531514e-03 4.341003e-03
#> 42 -2.248139e-03 2.407290e-03
#> 43 -2.539062e-03 1.815514e-03
#> 44 -4.107813e-03 3.233766e-03
#> 45 -1.949669e-03 2.100167e-03
#> 46 -1.689410e-03 1.184601e-03
#> 47 -3.050033e-03 2.408534e-03
#> 48 -1.660154e-03 1.796844e-03
#> 49 -1.117247e-03 7.669228e-04
#> 50 -2.264421e-03 1.793697e-03
#> 51 -1.391000e-03 1.511400e-03
#> 52 -7.384563e-04 4.966074e-04
#> 53 -1.681075e-03 1.335711e-03
#> 54 -1.149006e-03 1.252521e-03
#> 55 -4.922970e-04 3.261955e-04
#> 56 -1.247970e-03 9.946039e-04
#> 57 -9.372265e-04 1.024491e-03
#> 58 -3.351638e-04 2.216741e-04
#> 59 -9.264318e-04 7.405684e-04
#> 60 -7.559671e-04 8.283408e-04
#> 61 -2.358482e-04 1.587104e-04
#> 62 -6.877288e-04 5.513891e-04
#> 63 -6.036956e-04 6.628975e-04
#> 64 -1.723991e-04 1.202489e-04
#> 65 -5.105234e-04 4.105140e-04
#> 66 -4.777898e-04 5.256490e-04
#> 67 -1.302231e-04 9.516049e-05
#> 68 -3.789725e-04 3.056139e-04
#> 69 -3.750983e-04 4.133943e-04
#> 70 -1.004820e-04 7.704388e-05
#> 71 -2.813144e-04 2.275051e-04
#> 72 -2.923324e-04 3.227035e-04
#> 73 -7.829647e-05 6.272428e-05
#> 74 -2.088174e-04 1.693481e-04
#> 75 -2.263207e-04 2.502179e-04
#> 76 -6.110170e-05 5.082303e-05
#> 77 -1.549993e-04 1.260485e-04
#> 78 -1.741582e-04 1.928304e-04
#> 79 -4.751980e-05 4.078313e-05
#> 80 -1.150481e-04 9.381276e-05
#> 81 -1.332793e-04 1.477781e-04
#> 82 -3.673537e-05 3.235445e-05
#> 83 -8.539138e-05 6.981529e-05
#> 84 -1.014807e-04 1.126758e-04
#> 85 -2.819540e-05 2.537118e-05
#> 86 -6.337697e-05 5.195195e-05
#> 87 -7.691072e-05 8.551151e-05
#> 88 -2.147791e-05 1.967521e-05
#> 89 -4.703597e-05 3.865577e-05
#> 90 -5.804140e-05 6.461889e-05
Monte Carlo Method
summary(mc)
#> effect interval est se R 2.5%
#> 1 total 1 -1.155252e-02 0.0450766621 20000 -0.0990919232
#> 2 direct 1 -1.156492e-02 0.0450087668 20000 -0.0972223947
#> 3 indirect 1 1.240083e-05 0.0014748635 20000 -0.0042250411
#> 4 total 2 -1.408408e-02 0.0554497249 20000 -0.1241844354
#> 5 direct 2 -1.404168e-02 0.0550795218 20000 -0.1180626904
#> 6 indirect 2 -4.240204e-05 0.0039405798 20000 -0.0118030092
#> 7 total 3 -1.309786e-02 0.0522240796 20000 -0.1202685772
#> 8 direct 3 -1.297147e-02 0.0517111445 20000 -0.1110756690
#> 9 indirect 3 -1.263934e-04 0.0061702704 20000 -0.0181241553
#> 10 total 4 -1.098008e-02 0.0444284706 20000 -0.1049248448
#> 11 direct 4 -1.079883e-02 0.0441313064 20000 -0.0969840784
#> 12 indirect 4 -1.812500e-04 0.0077625855 20000 -0.0212186277
#> 13 total 5 -8.727446e-03 0.0358651653 20000 -0.0856074927
#> 14 direct 5 -8.538213e-03 0.0361202221 20000 -0.0819461944
#> 15 indirect 5 -1.892332e-04 0.0087953252 20000 -0.0216337222
#> 16 total 6 -6.718987e-03 0.0280409654 20000 -0.0671326945
#> 17 direct 6 -6.559465e-03 0.0290529289 20000 -0.0682705760
#> 18 indirect 6 -1.595224e-04 0.0094852651 20000 -0.0198583996
#> 19 total 7 -5.063621e-03 0.0214687115 20000 -0.0503872324
#> 20 direct 7 -4.953961e-03 0.0232740415 20000 -0.0563467590
#> 21 indirect 7 -1.096601e-04 0.0099649702 20000 -0.0175140456
#> 22 total 8 -3.757504e-03 0.0162356052 20000 -0.0369771850
#> 23 direct 8 -3.702261e-03 0.0187195519 20000 -0.0469574072
#> 24 indirect 8 -5.524344e-05 0.0102546375 20000 -0.0150727514
#> 25 total 9 -2.755037e-03 0.0122447184 20000 -0.0263868686
#> 26 direct 9 -2.748541e-03 0.0151928219 20000 -0.0392840681
#> 27 indirect 9 -6.495647e-06 0.0103282732 20000 -0.0127701980
#> 28 total 10 -2.000302e-03 0.0093166897 20000 -0.0185455029
#> 29 direct 10 -2.031880e-03 0.0124807056 20000 -0.0327304252
#> 30 indirect 10 3.157812e-05 0.0101708508 20000 -0.0104025083
#> 31 total 11 -1.440230e-03 0.0072353332 20000 -0.0128948430
#> 32 direct 11 -1.497986e-03 0.0103964829 20000 -0.0273130747
#> 33 indirect 11 5.775644e-05 0.0097974716 20000 -0.0083829738
#> 34 total 12 -1.029362e-03 0.0057765204 20000 -0.0092412725
#> 35 direct 12 -1.102424e-03 0.0087905443 20000 -0.0227660937
#> 36 indirect 12 7.306199e-05 0.0092486001 20000 -0.0067880541
#> 37 total 13 -7.308208e-04 0.0047370087 20000 -0.0065721557
#> 38 direct 13 -8.103814e-04 0.0075481347 20000 -0.0191675413
#> 39 indirect 13 7.956057e-05 0.0085781132 20000 -0.0054257856
#> 40 total 14 -5.156803e-04 0.0039595266 20000 -0.0047768709
#> 41 direct 14 -5.952559e-04 0.0065830100 20000 -0.0161094720
#> 42 indirect 14 7.957560e-05 0.0078424310 20000 -0.0044426408
#> 43 total 15 -3.617743e-04 0.0033416507 20000 -0.0034921450
#> 44 direct 15 -4.370234e-04 0.0058307720 20000 -0.0134223334
#> 45 indirect 15 7.524910e-05 0.0070930749 20000 -0.0035742017
#> 46 total 16 -2.524048e-04 0.0028269755 20000 -0.0025637145
#> 47 direct 16 -3.207496e-04 0.0052432269 20000 -0.0113344913
#> 48 indirect 16 6.834481e-05 0.0063725246 20000 -0.0028811297
#> 49 total 17 -1.751620e-04 0.0023887627 20000 -0.0018800630
#> 50 direct 17 -2.353620e-04 0.0047840422 20000 -0.0095094545
#> 51 indirect 17 6.020000e-05 0.0057125325 20000 -0.0023315354
#> 52 total 18 -1.209245e-04 0.0020153851 20000 -0.0014105919
#> 53 direct 18 -1.726819e-04 0.0044255601 20000 -0.0080434717
#> 54 indirect 18 5.175744e-05 0.0051340204 20000 -0.0018930980
#> 55 total 19 -8.305073e-05 0.0017008255 20000 -0.0010579172
#> 56 direct 19 -1.266829e-04 0.0041465177 20000 -0.0068076944
#> 57 indirect 19 4.363220e-05 0.0046478780 20000 -0.0015452118
#> 58 total 20 -5.674485e-05 0.0014397439 20000 -0.0008172077
#> 59 direct 20 -9.293168e-05 0.0039304279 20000 -0.0057692385
#> 60 indirect 20 3.618682e-05 0.0042562636 20000 -0.0012565600
#> 61 total 21 -3.856892e-05 0.0012256634 20000 -0.0006459969
#> 62 direct 21 -6.816989e-05 0.0037644241 20000 -0.0048783215
#> 63 indirect 21 2.960097e-05 0.0039542656 20000 -0.0010291880
#> 64 total 22 -2.607509e-05 0.0010509857 20000 -0.0005237013
#> 65 direct 22 -5.000466e-05 0.0036384236 20000 -0.0041143651
#> 66 indirect 22 2.392957e-05 0.0037319004 20000 -0.0008331698
#> 67 total 23 -1.753128e-05 0.0009078544 20000 -0.0004306387
#> 68 direct 23 -3.667930e-05 0.0035445070 20000 -0.0034576963
#> 69 indirect 23 1.914802e-05 0.0035763224 20000 -0.0006658227
#> 70 total 24 -1.171906e-05 0.0007891357 20000 -0.0003690701
#> 71 direct 24 -2.690463e-05 0.0034764487 20000 -0.0029511770
#> 72 indirect 24 1.518557e-05 0.0034739374 20000 -0.0005438558
#> 73 total 25 -7.786092e-06 0.0006890611 20000 -0.0003209581
#> 74 direct 25 -1.973467e-05 0.0034293530 20000 -0.0024954583
#> 75 indirect 25 1.194857e-05 0.0034120502 20000 -0.0004439975
#> 76 total 26 -5.139335e-06 0.0006033943 20000 -0.0002881791
#> 77 direct 26 -1.447540e-05 0.0033993696 20000 -0.0021049203
#> 78 indirect 26 9.336065e-06 0.0033798471 20000 -0.0003639180
#> 79 total 27 -3.368332e-06 0.0005292303 20000 -0.0002547013
#> 80 direct 27 -1.061769e-05 0.0033834699 20000 -0.0017748022
#> 81 indirect 27 7.249358e-06 0.0033687536 20000 -0.0002975633
#> 82 total 28 -2.190458e-06 0.0004646283 20000 -0.0002170266
#> 83 direct 28 -7.788048e-06 0.0033792697 20000 -0.0015091745
#> 84 indirect 28 5.597590e-06 0.0033723593 20000 -0.0002436937
#> 85 total 29 -1.412110e-06 0.0004082454 20000 -0.0001834578
#> 86 direct 29 -5.712505e-06 0.0033848902 20000 -0.0012811507
#> 87 indirect 29 4.300396e-06 0.0033861192 20000 -0.0002009515
#> 88 total 30 -9.013541e-07 0.0003590597 20000 -0.0001521658
#> 89 direct 30 -4.190099e-06 0.0033988487 20000 -0.0010833878
#> 90 indirect 30 3.288745e-06 0.0034069817 20000 -0.0001663632
#> 97.5%
#> 1 0.0781222619
#> 2 0.0792662806
#> 3 0.0020032310
#> 4 0.0935027598
#> 5 0.0988202133
#> 6 0.0044038829
#> 7 0.0851064892
#> 8 0.0932953589
#> 9 0.0074423653
#> 10 0.0702879215
#> 11 0.0780846237
#> 12 0.0112608150
#> 13 0.0558335100
#> 14 0.0622351221
#> 15 0.0156243688
#> 16 0.0446650471
#> 17 0.0483389167
#> 18 0.0197675729
#> 19 0.0361614651
#> 20 0.0368825424
#> 21 0.0233798935
#> 22 0.0296534937
#> 23 0.0281567512
#> 24 0.0259947236
#> 25 0.0248617014
#> 26 0.0214321751
#> 27 0.0271836759
#> 28 0.0209040548
#> 29 0.0162568638
#> 30 0.0272756323
#> 31 0.0177299868
#> 32 0.0122934400
#> 33 0.0265848171
#> 34 0.0152823455
#> 35 0.0094009156
#> 36 0.0250429032
#> 37 0.0132099375
#> 38 0.0071365593
#> 39 0.0230729814
#> 40 0.0115810109
#> 41 0.0055340868
#> 42 0.0209143961
#> 43 0.0100240967
#> 44 0.0042938045
#> 45 0.0183039579
#> 46 0.0086229910
#> 47 0.0033747394
#> 48 0.0159047147
#> 49 0.0073369732
#> 50 0.0026024189
#> 51 0.0135570373
#> 52 0.0062047950
#> 53 0.0020602735
#> 54 0.0114081591
#> 55 0.0052249292
#> 56 0.0016315341
#> 57 0.0096441589
#> 58 0.0044055635
#> 59 0.0012972294
#> 60 0.0080979347
#> 61 0.0037098481
#> 62 0.0010180311
#> 63 0.0067877136
#> 64 0.0031044754
#> 65 0.0007992242
#> 66 0.0057146570
#> 67 0.0026008702
#> 68 0.0006361584
#> 69 0.0047989242
#> 70 0.0021906199
#> 71 0.0005016164
#> 72 0.0040748018
#> 73 0.0018568413
#> 74 0.0003975743
#> 75 0.0034172942
#> 76 0.0015695605
#> 77 0.0003136096
#> 78 0.0028906852
#> 79 0.0013335047
#> 80 0.0002482165
#> 81 0.0024810970
#> 82 0.0011336929
#> 83 0.0001970045
#> 84 0.0021199181
#> 85 0.0009619083
#> 86 0.0001590837
#> 87 0.0018210976
#> 88 0.0008282391
#> 89 0.0001267157
#> 90 0.0015695901
Parametric Bootstrap Method
pb <- BootParaYMX(boot = boot, phi_hat = phi_hat)
plot(pb)
plot(pb, type = "bc")
summary(pb)
#> effect interval est se R 2.5% 97.5%
#> 1 total 1 -1.155252e-02 0.0435677096 1000 -0.0958694237 7.346596e-02
#> 2 direct 1 -1.156492e-02 0.0434687647 1000 -0.0950729391 7.705411e-02
#> 3 indirect 1 1.240083e-05 0.0013224681 1000 -0.0038573273 1.693653e-03
#> 4 total 2 -1.408408e-02 0.0531265157 1000 -0.1194266017 8.772910e-02
#> 5 direct 2 -1.404168e-02 0.0526928335 1000 -0.1144580750 9.089847e-02
#> 6 indirect 2 -4.240204e-05 0.0035270982 1000 -0.0102120881 3.549647e-03
#> 7 total 3 -1.309786e-02 0.0495544514 1000 -0.1131581667 8.042671e-02
#> 8 direct 3 -1.297147e-02 0.0489538940 1000 -0.1063503147 8.670756e-02
#> 9 indirect 3 -1.263934e-04 0.0054698544 1000 -0.0156504915 6.245830e-03
#> 10 total 4 -1.098008e-02 0.0417049535 1000 -0.0984022455 6.795841e-02
#> 11 direct 4 -1.079883e-02 0.0412818967 1000 -0.0930176426 7.163681e-02
#> 12 indirect 4 -1.812500e-04 0.0068142776 1000 -0.0188837687 9.347666e-03
#> 13 total 5 -8.727446e-03 0.0332707863 1000 -0.0784429454 5.436900e-02
#> 14 direct 5 -8.538213e-03 0.0333166061 1000 -0.0780512300 5.651611e-02
#> 15 indirect 5 -1.892332e-04 0.0076738109 1000 -0.0200737221 1.380427e-02
#> 16 total 6 -6.718987e-03 0.0256960696 1000 -0.0605558795 4.274293e-02
#> 17 direct 6 -6.559465e-03 0.0263473196 1000 -0.0620341401 4.108887e-02
#> 18 indirect 6 -1.595224e-04 0.0082490818 1000 -0.0181912372 1.826480e-02
#> 19 total 7 -5.063621e-03 0.0194520377 1000 -0.0445088899 3.476500e-02
#> 20 direct 7 -4.953961e-03 0.0206734044 1000 -0.0526431837 3.107916e-02
#> 21 indirect 7 -1.096601e-04 0.0086302586 1000 -0.0163192412 2.205054e-02
#> 22 total 8 -3.757504e-03 0.0145889407 1000 -0.0319983285 2.724719e-02
#> 23 direct 8 -3.702261e-03 0.0162111342 1000 -0.0445066832 2.308570e-02
#> 24 indirect 8 -5.524344e-05 0.0088049718 1000 -0.0126173251 2.473030e-02
#> 25 total 9 -2.755037e-03 0.0109708958 1000 -0.0226298542 2.290323e-02
#> 26 direct 9 -2.748541e-03 0.0127579017 1000 -0.0362339722 1.716800e-02
#> 27 indirect 9 -6.495647e-06 0.0087387707 1000 -0.0101717593 2.538861e-02
#> 28 total 10 -2.000302e-03 0.0083779895 1000 -0.0165304791 1.872097e-02
#> 29 direct 10 -2.031880e-03 0.0101006151 1000 -0.0303313779 1.269067e-02
#> 30 indirect 10 3.157812e-05 0.0084259066 1000 -0.0080874469 2.602083e-02
#> 31 total 11 -1.440230e-03 0.0065594039 1000 -0.0112915494 1.586347e-02
#> 32 direct 11 -1.497986e-03 0.0080544226 1000 -0.0243573951 9.554255e-03
#> 33 indirect 11 5.775644e-05 0.0078974317 1000 -0.0062253944 2.481158e-02
#> 34 total 12 -1.029362e-03 0.0052748341 1000 -0.0082995669 1.401353e-02
#> 35 direct 12 -1.102424e-03 0.0064716180 1000 -0.0197692937 7.312133e-03
#> 36 indirect 12 7.306199e-05 0.0072087990 1000 -0.0050329983 2.322720e-02
#> 37 total 13 -7.308208e-04 0.0043300562 1000 -0.0058138126 1.177777e-02
#> 38 direct 13 -8.103814e-04 0.0052389512 1000 -0.0161940143 5.353742e-03
#> 39 indirect 13 7.956057e-05 0.0064242522 1000 -0.0040552829 2.089176e-02
#> 40 total 14 -5.156803e-04 0.0035941115 1000 -0.0043734107 1.037065e-02
#> 41 direct 14 -5.952559e-04 0.0042714027 1000 -0.0133913382 4.012763e-03
#> 42 indirect 14 7.957560e-05 0.0056049010 1000 -0.0031815096 1.898019e-02
#> 43 total 15 -3.617743e-04 0.0029930764 1000 -0.0031750490 8.605725e-03
#> 44 direct 15 -4.370234e-04 0.0035057371 1000 -0.0110723934 3.017559e-03
#> 45 indirect 15 7.524910e-05 0.0048016806 1000 -0.0025852968 1.675766e-02
#> 46 total 16 -2.524048e-04 0.0024913121 1000 -0.0023446740 7.136111e-03
#> 47 direct 16 -3.207496e-04 0.0028950065 1000 -0.0094582021 2.329541e-03
#> 48 indirect 16 6.834481e-05 0.0040523384 1000 -0.0021177097 1.439528e-02
#> 49 total 17 -1.751620e-04 0.0020730217 1000 -0.0017714858 5.916599e-03
#> 50 direct 17 -2.353620e-04 0.0024042430 1000 -0.0079046802 1.726695e-03
#> 51 indirect 17 6.020000e-05 0.0033811787 1000 -0.0017455430 1.178027e-02
#> 52 total 18 -1.209245e-04 0.0017296453 1000 -0.0012863367 4.937751e-03
#> 53 direct 18 -1.726819e-04 0.0020072307 1000 -0.0066526410 1.299868e-03
#> 54 indirect 18 5.175744e-05 0.0028004486 1000 -0.0013848248 9.738431e-03
#> 55 total 19 -8.305073e-05 0.0014531512 1000 -0.0009154029 4.069133e-03
#> 56 direct 19 -1.266829e-04 0.0016841453 1000 -0.0056687191 9.708863e-04
#> 57 indirect 19 4.363220e-05 0.0023125462 1000 -0.0010685087 7.971441e-03
#> 58 total 20 -5.674485e-05 0.0012336288 1000 -0.0007566349 3.247148e-03
#> 59 direct 20 -9.293168e-05 0.0014198486 1000 -0.0047196615 7.212730e-04
#> 60 indirect 20 3.618682e-05 0.0019125253 1000 -0.0008495774 6.660593e-03
#> 61 total 21 -3.856892e-05 0.0010596076 1000 -0.0006270144 2.762743e-03
#> 62 direct 21 -6.816989e-05 0.0012026619 1000 -0.0038688826 5.425605e-04
#> 63 indirect 21 2.960097e-05 0.0015906008 1000 -0.0006614730 5.476547e-03
#> 64 total 22 -2.607509e-05 0.0009196900 1000 -0.0004918352 2.296093e-03
#> 65 direct 22 -5.000466e-05 0.0010234834 1000 -0.0031746316 4.387517e-04
#> 66 indirect 22 2.392957e-05 0.0013344858 1000 -0.0005348933 4.232860e-03
#> 67 total 23 -1.753128e-05 0.0008042152 1000 -0.0004471469 1.914074e-03
#> 68 direct 23 -3.667930e-05 0.0008751510 1000 -0.0026000383 3.397428e-04
#> 69 indirect 23 1.914802e-05 0.0011314110 1000 -0.0004436033 3.462451e-03
#> 70 total 24 -1.171906e-05 0.0007061295 1000 -0.0003726414 1.610965e-03
#> 71 direct 24 -2.690463e-05 0.0007519768 1000 -0.0021327436 2.589997e-04
#> 72 indirect 24 1.518557e-05 0.0009696484 1000 -0.0003660600 2.838899e-03
#> 73 total 25 -7.786092e-06 0.0006209368 1000 -0.0003094031 1.355799e-03
#> 74 direct 25 -1.973467e-05 0.0006494067 1000 -0.0017828050 1.977075e-04
#> 75 indirect 25 1.194857e-05 0.0008393982 1000 -0.0002981540 2.368067e-03
#> 76 total 26 -5.139335e-06 0.0005460821 1000 -0.0002461929 1.145400e-03
#> 77 direct 26 -1.447540e-05 0.0005637667 1000 -0.0014900844 1.504807e-04
#> 78 indirect 26 9.336065e-06 0.0007330438 1000 -0.0002409071 1.993166e-03
#> 79 total 27 -3.368332e-06 0.0004801948 1000 -0.0002145393 9.859990e-04
#> 80 direct 27 -1.061769e-05 0.0004920735 1000 -0.0012454274 1.164825e-04
#> 81 indirect 27 7.249358e-06 0.0006449436 1000 -0.0001921843 1.684577e-03
#> 82 total 28 -2.190458e-06 0.0004224571 1000 -0.0001783019 8.487947e-04
#> 83 direct 28 -7.788048e-06 0.0004318911 1000 -0.0010409420 9.095779e-05
#> 84 indirect 28 5.597590e-06 0.0005710051 1000 -0.0001533979 1.438988e-03
#> 85 total 29 -1.412110e-06 0.0003721956 1000 -0.0001438534 7.281713e-04
#> 86 direct 29 -5.712505e-06 0.0003812212 1000 -0.0008700319 6.982327e-05
#> 87 indirect 29 4.300396e-06 0.0005082427 1000 -0.0001229109 1.231467e-03
#> 88 total 30 -9.013541e-07 0.0003286844 1000 -0.0001121523 6.122969e-04
#> 89 direct 30 -4.190099e-06 0.0003384189 1000 -0.0007271842 5.309047e-05
#> 90 indirect 30 3.288745e-06 0.0004544295 1000 -0.0001012807 1.055336e-03
summary(pb, type = "bc")
#> effect interval est se R 2.5% 97.5%
#> 1 total 1 -1.155252e-02 0.0435677096 1000 -0.0945674875 7.991905e-02
#> 2 direct 1 -1.156492e-02 0.0434687647 1000 -0.0933288314 8.005903e-02
#> 3 indirect 1 1.240083e-05 0.0013224681 1000 -0.0017710080 2.637875e-03
#> 4 total 2 -1.408408e-02 0.0531265157 1000 -0.1142102974 9.627751e-02
#> 5 direct 2 -1.404168e-02 0.0526928335 1000 -0.1087763078 1.017215e-01
#> 6 indirect 2 -4.240204e-05 0.0035270982 1000 -0.0045007812 7.707056e-03
#> 7 total 3 -1.309786e-02 0.0495544514 1000 -0.1076135554 8.893898e-02
#> 8 direct 3 -1.297147e-02 0.0489538940 1000 -0.1022827509 9.105333e-02
#> 9 indirect 3 -1.263934e-04 0.0054698544 1000 -0.0074808609 1.185692e-02
#> 10 total 4 -1.098008e-02 0.0417049535 1000 -0.0928004368 7.321774e-02
#> 11 direct 4 -1.079883e-02 0.0412818967 1000 -0.0888062348 7.686182e-02
#> 12 indirect 4 -1.812500e-04 0.0068142776 1000 -0.0104727376 1.616952e-02
#> 13 total 5 -8.727446e-03 0.0332707863 1000 -0.0748180129 5.855024e-02
#> 14 direct 5 -8.538213e-03 0.0333166061 1000 -0.0724495767 5.876262e-02
#> 15 indirect 5 -1.892332e-04 0.0076738109 1000 -0.0127158138 1.841111e-02
#> 16 total 6 -6.718987e-03 0.0256960696 1000 -0.0576835151 4.623283e-02
#> 17 direct 6 -6.559465e-03 0.0263473196 1000 -0.0605242201 4.551359e-02
#> 18 indirect 6 -1.595224e-04 0.0082490818 1000 -0.0145940684 2.007294e-02
#> 19 total 7 -5.063621e-03 0.0194520377 1000 -0.0425532161 3.692599e-02
#> 20 direct 7 -4.953961e-03 0.0206734044 1000 -0.0497989739 3.470518e-02
#> 21 indirect 7 -1.096601e-04 0.0086302586 1000 -0.0169357762 2.053978e-02
#> 22 total 8 -3.757504e-03 0.0145889407 1000 -0.0303476778 3.044197e-02
#> 23 direct 8 -3.702261e-03 0.0162111342 1000 -0.0415090939 2.483510e-02
#> 24 indirect 8 -5.524344e-05 0.0088049718 1000 -0.0154303830 2.044716e-02
#> 25 total 9 -2.755037e-03 0.0109708958 1000 -0.0226238779 2.295718e-02
#> 26 direct 9 -2.748541e-03 0.0127579017 1000 -0.0353411548 1.752614e-02
#> 27 indirect 9 -6.495647e-06 0.0087387707 1000 -0.0135020218 1.883801e-02
#> 28 total 10 -2.000302e-03 0.0083779895 1000 -0.0172391859 1.616464e-02
#> 29 direct 10 -2.031880e-03 0.0101006151 1000 -0.0287355075 1.277679e-02
#> 30 indirect 10 3.157812e-05 0.0084259066 1000 -0.0130150657 1.552266e-02
#> 31 total 11 -1.440230e-03 0.0065594039 1000 -0.0147332722 1.052396e-02
#> 32 direct 11 -1.497986e-03 0.0080544226 1000 -0.0241178631 9.761737e-03
#> 33 indirect 11 5.775644e-05 0.0078974317 1000 -0.0135436867 1.326376e-02
#> 34 total 12 -1.029362e-03 0.0052748341 1000 -0.0139279295 7.132784e-03
#> 35 direct 12 -1.102424e-03 0.0064716180 1000 -0.0196745130 7.371820e-03
#> 36 indirect 12 7.306199e-05 0.0072087990 1000 -0.0127396841 1.045825e-02
#> 37 total 13 -7.308208e-04 0.0043300562 1000 -0.0129464330 4.603915e-03
#> 38 direct 13 -8.103814e-04 0.0052389512 1000 -0.0161587150 5.382953e-03
#> 39 indirect 13 7.956057e-05 0.0064242522 1000 -0.0106439985 8.931038e-03
#> 40 total 14 -5.156803e-04 0.0035941115 1000 -0.0104799258 2.897505e-03
#> 41 direct 14 -5.952559e-04 0.0042714027 1000 -0.0136528224 3.982367e-03
#> 42 indirect 14 7.957560e-05 0.0056049010 1000 -0.0086808382 7.223920e-03
#> 43 total 15 -3.617743e-04 0.0029930764 1000 -0.0086906140 1.745253e-03
#> 44 direct 15 -4.370234e-04 0.0035057371 1000 -0.0125858849 2.890202e-03
#> 45 indirect 15 7.524910e-05 0.0048016806 1000 -0.0073899432 6.067789e-03
#> 46 total 16 -2.524048e-04 0.0024913121 1000 -0.0067701164 1.057699e-03
#> 47 direct 16 -3.207496e-04 0.0028950065 1000 -0.0106103975 2.066445e-03
#> 48 indirect 16 6.834481e-05 0.0040523384 1000 -0.0068672773 4.899244e-03
#> 49 total 17 -1.751620e-04 0.0020730217 1000 -0.0053247072 6.817321e-04
#> 50 direct 17 -2.353620e-04 0.0024042430 1000 -0.0087475319 1.444537e-03
#> 51 indirect 17 6.020000e-05 0.0033811787 1000 -0.0062800881 3.810206e-03
#> 52 total 18 -1.209245e-04 0.0017296453 1000 -0.0041508596 5.106974e-04
#> 53 direct 18 -1.726819e-04 0.0020072307 1000 -0.0073140670 1.077852e-03
#> 54 indirect 18 5.175744e-05 0.0028004486 1000 -0.0056228239 3.094852e-03
#> 55 total 19 -8.305073e-05 0.0014531512 1000 -0.0032971061 3.559549e-04
#> 56 direct 19 -1.266829e-04 0.0016841453 1000 -0.0061669304 8.373588e-04
#> 57 indirect 19 4.363220e-05 0.0023125462 1000 -0.0053312951 2.412601e-03
#> 58 total 20 -5.674485e-05 0.0012336288 1000 -0.0026152580 2.453708e-04
#> 59 direct 20 -9.293168e-05 0.0014198486 1000 -0.0052340210 6.394442e-04
#> 60 indirect 20 3.618682e-05 0.0019125253 1000 -0.0047984985 1.851953e-03
#> 61 total 21 -3.856892e-05 0.0010596076 1000 -0.0020489041 1.991637e-04
#> 62 direct 21 -6.816989e-05 0.0012026619 1000 -0.0046881060 4.693067e-04
#> 63 indirect 21 2.960097e-05 0.0015906008 1000 -0.0043235816 1.449086e-03
#> 64 total 22 -2.607509e-05 0.0009196900 1000 -0.0020632670 1.303672e-04
#> 65 direct 22 -5.000466e-05 0.0010234834 1000 -0.0039693706 3.358860e-04
#> 66 indirect 22 2.392957e-05 0.0013344858 1000 -0.0038982810 1.127460e-03
#> 67 total 23 -1.753128e-05 0.0008042152 1000 -0.0027339083 9.674486e-05
#> 68 direct 23 -3.667930e-05 0.0008751510 1000 -0.0033401476 2.432599e-04
#> 69 indirect 23 1.914802e-05 0.0011314110 1000 -0.0034374437 8.816023e-04
#> 70 total 24 -1.171906e-05 0.0007061295 1000 -0.0028382057 8.366180e-05
#> 71 direct 24 -2.690463e-05 0.0007519768 1000 -0.0028083935 1.767116e-04
#> 72 indirect 24 1.518557e-05 0.0009696484 1000 -0.0030886792 6.951400e-04
#> 73 total 25 -7.786092e-06 0.0006209368 1000 -0.0025452293 7.201053e-05
#> 74 direct 25 -1.973467e-05 0.0006494067 1000 -0.0023644636 1.214232e-04
#> 75 indirect 25 1.194857e-05 0.0008393982 1000 -0.0026953268 5.366695e-04
#> 76 total 26 -5.139335e-06 0.0005460821 1000 -0.0020137171 5.937005e-05
#> 77 direct 26 -1.447540e-05 0.0005637667 1000 -0.0019907280 9.197116e-05
#> 78 indirect 26 9.336065e-06 0.0007330438 1000 -0.0023830746 4.260677e-04
#> 79 total 27 -3.368332e-06 0.0004801948 1000 -0.0013146685 6.261941e-05
#> 80 direct 27 -1.061769e-05 0.0004920735 1000 -0.0016785760 6.846242e-05
#> 81 indirect 27 7.249358e-06 0.0006449436 1000 -0.0020240338 3.317069e-04
#> 82 total 28 -2.190458e-06 0.0004224571 1000 -0.0008581960 5.566068e-05
#> 83 direct 28 -7.788048e-06 0.0004318911 1000 -0.0014576212 5.367506e-05
#> 84 indirect 28 5.597590e-06 0.0005710051 1000 -0.0017176265 2.675025e-04
#> 85 total 29 -1.412110e-06 0.0003721956 1000 -0.0005407394 4.966752e-05
#> 86 direct 29 -5.712505e-06 0.0003812212 1000 -0.0012547387 4.143059e-05
#> 87 indirect 29 4.300396e-06 0.0005082427 1000 -0.0014564035 2.219397e-04
#> 88 total 30 -9.013541e-07 0.0003286844 1000 -0.0003856950 4.775617e-05
#> 89 direct 30 -4.190099e-06 0.0003384189 1000 -0.0010697714 3.084611e-05
#> 90 indirect 30 3.288745e-06 0.0004544295 1000 -0.0012339169 1.804109e-04
Confidence Intervals (Standardized)
Model
Delta Method
delta <- DeltaStdXMY(theta_hat = theta_hat)
plot(delta)
summary(delta)
#> effect interval est se z p
#> 1 total 1 -0.1364455732 0.0344870006 -3.9564349 7.607661e-05
#> 2 direct 1 -0.3297415327 0.0484631056 -6.8039703 1.017745e-11
#> 3 indirect 1 0.1932959594 0.0232769906 8.3041645 1.005243e-16
#> 4 total 2 0.0594752696 0.0360580517 1.6494310 9.905936e-02
#> 5 direct 2 -0.4003594595 0.0681111976 -5.8780270 4.151852e-09
#> 6 indirect 2 0.4598347291 0.0575485228 7.9903828 1.345204e-15
#> 7 total 3 0.2520939102 0.0346405996 7.2774118 3.402852e-13
#> 8 direct 3 -0.3698451919 0.0765398379 -4.8320613 1.351266e-06
#> 9 indirect 3 0.6219391021 0.0851670460 7.3025792 2.823029e-13
#> 10 total 4 0.3631995292 0.0367070364 9.8945479 4.395921e-23
#> 11 direct 4 -0.3078984260 0.0779971204 -3.9475615 7.895126e-05
#> 12 indirect 4 0.6710979552 0.1017830464 6.5934159 4.298203e-11
#> 13 total 5 0.3985359430 0.0403896602 9.8672765 5.771027e-23
#> 14 direct 5 -0.2434433603 0.0747101015 -3.2585066 1.120003e-03
#> 15 indirect 5 0.6419793033 0.1076648874 5.9627546 2.480208e-09
#> 16 total 6 0.3833083417 0.0432383710 8.8650042 7.650117e-19
#> 17 direct 6 -0.1870248610 0.0685062714 -2.7300400 6.332664e-03
#> 18 indirect 6 0.5703332028 0.1051192355 5.4255836 5.776548e-08
#> 19 total 7 0.3409212998 0.0444314312 7.6729759 1.680511e-14
#> 20 direct 7 -0.1412483935 0.0608193039 -2.3224270 2.020995e-02
#> 21 indirect 7 0.4821696934 0.0969497603 4.9733975 6.578957e-07
#> 22 total 8 0.2879258411 0.0441985492 6.5143731 7.299385e-11
#> 23 direct 8 -0.1055596464 0.0526702795 -2.0041596 4.505297e-02
#> 24 indirect 8 0.3934854874 0.0857689935 4.5877359 4.480791e-06
#> 25 total 9 0.2343958649 0.0431007993 5.4383183 5.378583e-08
#> 26 direct 9 -0.0783669961 0.0447258418 -1.7521637 7.974566e-02
#> 27 indirect 9 0.3127628609 0.0737042989 4.2434819 2.200781e-05
#> 28 total 10 0.1856524532 0.0415794961 4.4650001 8.006873e-06
#> 29 direct 10 -0.0579333864 0.0373807283 -1.5498196 1.211848e-01
#> 30 indirect 10 0.2435858396 0.0622635259 3.9121755 9.146841e-05
#> 31 total 11 0.1439430161 0.0398077916 3.6159508 2.992472e-04
#> 32 direct 11 -0.0427109039 0.0308353509 -1.3851279 1.660134e-01
#> 33 indirect 11 0.1866539200 0.0523077123 3.5683824 3.591920e-04
#> 34 total 12 0.1097107495 0.0377664676 2.9049778 3.672792e-03
#> 35 direct 12 -0.0314325594 0.0251591354 -1.2493497 2.115372e-01
#> 36 indirect 12 0.1411433089 0.0441217139 3.1989535 1.379274e-03
#> 37 total 13 0.0824480421 0.0353888033 2.3297776 1.981791e-02
#> 38 direct 13 -0.0231057664 0.0203386285 -1.1360533 2.559343e-01
#> 39 indirect 13 0.1055538085 0.0375692693 2.8095784 4.960643e-03
#> 40 total 14 0.0612264758 0.0326578636 1.8747851 6.082229e-02
#> 41 direct 14 -0.0169720621 0.0163122160 -1.0404510 2.981304e-01
#> 42 indirect 14 0.0781985380 0.0322964873 2.4212707 1.546636e-02
#> 43 total 15 0.0450029669 0.0296310645 1.5187766 1.288187e-01
#> 44 direct 15 -0.0124605035 0.0129940534 -0.9589389 3.375895e-01
#> 45 indirect 15 0.0574634704 0.0279196916 2.0581700 3.957382e-02
#> 46 total 16 0.0327816587 0.0264199814 1.2407904 2.146832e-01
#> 47 direct 16 -0.0091452818 0.0102897788 -0.8887734 3.741249e-01
#> 48 indirect 16 0.0419269404 0.0241397628 1.7368414 8.241516e-02
#> 49 total 17 0.0236880501 0.0231574722 1.0229117 3.063496e-01
#> 50 direct 17 -0.0067106920 0.0081062363 -0.8278431 4.077593e-01
#> 51 indirect 17 0.0303987420 0.0207730349 1.4633751 1.433648e-01
#> 52 total 18 0.0169927681 0.0199700303 0.8509135 3.948174e-01
#> 53 direct 18 -0.0049235432 0.0063570128 -0.7745058 4.386317e-01
#> 54 indirect 18 0.0219163112 0.0177301749 1.2361024 2.164205e-01
#> 55 total 19 0.0121085470 0.0169609111 0.7139090 4.752835e-01
#> 56 direct 19 -0.0036120109 0.0049651626 -0.7274708 4.669376e-01
#> 57 indirect 19 0.0157205578 0.0149795125 1.0494706 2.939616e-01
#> 58 total 20 0.0085745994 0.0142032098 0.6037086 5.460374e-01
#> 59 direct 20 -0.0026496878 0.0038641302 -0.6857139 4.928935e-01
#> 60 indirect 20 0.0112242872 0.0125156183 0.8968224 3.698137e-01
#> 61 total 21 0.0060365041 0.0117399440 0.5141851 6.071226e-01
#> 62 direct 21 -0.0019436746 0.0029975812 -0.6484143 5.167170e-01
#> 63 indirect 21 0.0079801787 0.0103391536 0.7718406 4.402088e-01
#> 64 total 22 0.0042259651 0.0095881447 0.4407490 6.593947e-01
#> 65 direct 22 -0.0014257435 0.0023186269 -0.6149085 5.386151e-01
#> 66 indirect 22 0.0056517086 0.0084466991 0.6691026 5.034300e-01
#> 67 total 23 0.0029425409 0.0077446509 0.3799449 7.039863e-01
#> 68 direct 23 -0.0010458082 0.0017887602 -0.5846553 5.587795e-01
#> 69 indirect 23 0.0039883491 0.0068272821 0.5841782 5.591005e-01
#> 70 total 24 0.0020381454 0.0061920995 0.3291526 7.420404e-01
#> 71 direct 24 -0.0007671105 0.0013766999 -0.5572097 5.773842e-01
#> 72 indirect 24 0.0028052560 0.0054626322 0.5135356 6.075767e-01
#> 73 total 25 0.0014044256 0.0049042685 0.2863680 7.745963e-01
#> 74 direct 25 -0.0005626791 0.0010572644 -0.5322028 5.945855e-01
#> 75 indirect 25 0.0019671047 0.0043290993 0.4543912 6.495473e-01
#> 76 total 26 0.0009627744 0.0038504020 0.2500452 8.025524e-01
#> 77 direct 26 -0.0004127257 0.0008103357 -0.5093268 6.105232e-01
#> 78 indirect 26 0.0013755001 0.0033999988 0.4045590 6.858017e-01
#> 79 total 27 0.0006565995 0.0029984424 0.2189802 8.266655e-01
#> 80 direct 27 -0.0003027338 0.0006199460 -0.4883229 6.253211e-01
#> 81 indirect 27 0.0009593334 0.0026477455 0.3623208 7.171123e-01
#> 82 total 28 0.0004454423 0.0023172654 0.1922276 8.475640e-01
#> 83 direct 28 -0.0002220545 0.0004734919 -0.4689721 6.390896e-01
#> 84 indirect 28 0.0006674968 0.0020455148 0.3263221 7.441806e-01
#> 85 total 29 0.0003005636 0.0017780844 0.1690379 8.657669e-01
#> 86 direct 29 -0.0001628762 0.0003610741 -0.4510879 6.519262e-01
#> 87 indirect 29 0.0004634398 0.0015683779 0.2954899 7.676197e-01
#> 88 total 30 0.0002016726 0.0013552143 0.1488123 8.817017e-01
#> 89 direct 30 -0.0001194690 0.0002749506 -0.4345108 6.639176e-01
#> 90 indirect 30 0.0003211416 0.0011939697 0.2689696 7.879531e-01
#> 2.5% 97.5%
#> 1 -0.2040388524 -0.0688522940
#> 2 -0.4247274742 -0.2347555911
#> 3 0.1476738961 0.2389180227
#> 4 -0.0111972131 0.1301477524
#> 5 -0.5338549537 -0.2668639653
#> 6 0.3470416970 0.5726277612
#> 7 0.1841995825 0.3199882379
#> 8 -0.5198605175 -0.2198298663
#> 9 0.4550147593 0.7888634449
#> 10 0.2912550600 0.4351439985
#> 11 -0.4607699729 -0.1550268791
#> 12 0.4716068501 0.8705890603
#> 13 0.3193736637 0.4776982223
#> 14 -0.3898724686 -0.0970142520
#> 15 0.4309600015 0.8529986051
#> 16 0.2985626918 0.4680539917
#> 17 -0.3212946857 -0.0527550364
#> 18 0.3643032870 0.7763631185
#> 19 0.2538372949 0.4280053048
#> 20 -0.2604520388 -0.0220447483
#> 21 0.2921516549 0.6721877318
#> 22 0.2012982765 0.3745534057
#> 23 -0.2087914972 -0.0023277955
#> 24 0.2253813491 0.5615896258
#> 25 0.1499198505 0.3188718793
#> 26 -0.1660280351 0.0092940430
#> 27 0.1683050896 0.4572206323
#> 28 0.1041581383 0.2671467681
#> 29 -0.1311982676 0.0153314948
#> 30 0.1215515712 0.3656201080
#> 31 0.0659211782 0.2219648540
#> 32 -0.1031470812 0.0177252734
#> 33 0.0841326877 0.2891751522
#> 34 0.0356898332 0.1837316658
#> 35 -0.0807435588 0.0178784399
#> 36 0.0546663388 0.2276202790
#> 37 0.0130872622 0.1518088219
#> 38 -0.0629687458 0.0167572130
#> 39 0.0319193937 0.1791882232
#> 40 -0.0027817606 0.1252347123
#> 41 -0.0489434179 0.0149992937
#> 42 0.0148985861 0.1414984898
#> 43 -0.0130728524 0.1030787862
#> 44 -0.0379283802 0.0130073733
#> 45 0.0027418805 0.1121850603
#> 46 -0.0190005533 0.0845638707
#> 47 -0.0293128776 0.0110223141
#> 48 -0.0053861252 0.0892400061
#> 49 -0.0216997615 0.0690758616
#> 50 -0.0225986232 0.0091772392
#> 51 -0.0103156582 0.0711131423
#> 52 -0.0221477722 0.0561333083
#> 53 -0.0173830592 0.0075359729
#> 54 -0.0128341930 0.0566668154
#> 55 -0.0211342280 0.0453513219
#> 56 -0.0133435508 0.0061195290
#> 57 -0.0136387472 0.0450798628
#> 58 -0.0192631802 0.0364123790
#> 59 -0.0102232438 0.0049238683
#> 60 -0.0133058739 0.0357544483
#> 61 -0.0169733632 0.0290463714
#> 62 -0.0078188257 0.0039314765
#> 63 -0.0122841900 0.0282445473
#> 64 -0.0145664532 0.0230183835
#> 65 -0.0059701687 0.0031186817
#> 66 -0.0109035173 0.0222069346
#> 67 -0.0122366960 0.0181217778
#> 68 -0.0045517137 0.0024600973
#> 69 -0.0093928779 0.0173695760
#> 70 -0.0100981465 0.0141744374
#> 71 -0.0034653928 0.0019311717
#> 72 -0.0079013064 0.0135118183
#> 73 -0.0082077640 0.0110166153
#> 74 -0.0026348791 0.0015095210
#> 75 -0.0065177741 0.0104519835
#> 76 -0.0065838749 0.0085094236
#> 77 -0.0020009546 0.0011755031
#> 78 -0.0052883750 0.0080393752
#> 79 -0.0052202397 0.0065334388
#> 80 -0.0015178057 0.0009123381
#> 81 -0.0042301525 0.0061488193
#> 82 -0.0040963145 0.0049871990
#> 83 -0.0011500816 0.0007059726
#> 84 -0.0033416385 0.0046766321
#> 85 -0.0031844177 0.0037855449
#> 86 -0.0008705684 0.0005448161
#> 87 -0.0026105244 0.0035374039
#> 88 -0.0024544986 0.0028578438
#> 89 -0.0006583622 0.0004194242
#> 90 -0.0020189961 0.0026612792
Monte Carlo Method
summary(mc)
#> effect interval est se R 2.5%
#> 1 total 1 -0.1364455732 0.034991098 20000 -0.2078373529
#> 2 direct 1 -0.3297415327 0.049422454 20000 -0.4315346086
#> 3 indirect 1 0.1932959594 0.023566517 20000 0.1493138688
#> 4 total 2 0.0594752696 0.036488158 20000 -0.0153256108
#> 5 direct 2 -0.4003594595 0.070137974 20000 -0.5525144959
#> 6 indirect 2 0.4598347291 0.058698512 20000 0.3551175854
#> 7 total 3 0.2520939102 0.034864640 20000 0.1822741532
#> 8 direct 3 -0.3698451919 0.080015618 20000 -0.5540782771
#> 9 indirect 3 0.6219391021 0.087834889 20000 0.4742520519
#> 10 total 4 0.3631995292 0.036849678 20000 0.2924783162
#> 11 direct 4 -0.3078984260 0.083406346 20000 -0.5108772738
#> 12 indirect 4 0.6710979552 0.106589297 20000 0.5024493987
#> 13 total 5 0.3985359430 0.040652282 20000 0.3247080507
#> 14 direct 5 -0.2434433603 0.082374417 20000 -0.4536821951
#> 15 indirect 5 0.6419793033 0.114947664 20000 0.4702363144
#> 16 total 6 0.3833083417 0.043754650 20000 0.3070340142
#> 17 direct 6 -0.1870248610 0.078504778 20000 -0.3954255311
#> 18 indirect 6 0.5703332028 0.114848328 20000 0.4095235989
#> 19 total 7 0.3409212998 0.045252584 20000 0.2648408177
#> 20 direct 7 -0.1412483935 0.073025383 20000 -0.3407462892
#> 21 indirect 7 0.4821696934 0.108788088 20000 0.3394230846
#> 22 total 8 0.2879258411 0.045315000 20000 0.2142290317
#> 23 direct 8 -0.1055596464 0.066820704 20000 -0.2931256135
#> 24 indirect 8 0.3934854874 0.099179019 20000 0.2690709067
#> 25 total 9 0.2343958649 0.044471555 20000 0.1625752809
#> 26 direct 9 -0.0783669961 0.060485100 20000 -0.2510833388
#> 27 indirect 9 0.3127628609 0.088048927 20000 0.2079578166
#> 28 total 10 0.1856524532 0.043162042 20000 0.1162086466
#> 29 direct 10 -0.0579333864 0.054390154 20000 -0.2153625086
#> 30 indirect 10 0.2435858396 0.076894558 20000 0.1556326919
#> 31 total 11 0.1439430161 0.041576146 20000 0.0770384052
#> 32 direct 11 -0.0427109039 0.048746778 20000 -0.1838620074
#> 33 indirect 11 0.1866539200 0.066649100 20000 0.1130759641
#> 34 total 12 0.1097107495 0.039711641 20000 0.0452800874
#> 35 direct 12 -0.0314325594 0.043655742 20000 -0.1566442465
#> 36 indirect 12 0.1411433089 0.057748028 20000 0.0788667560
#> 37 total 13 0.0824480421 0.037508416 20000 0.0219176985
#> 38 direct 13 -0.0231057664 0.039146047 20000 -0.1345384480
#> 39 indirect 13 0.1055538085 0.050261409 20000 0.0524446228
#> 40 total 14 0.0612264758 0.034945748 20000 0.0056040374
#> 41 direct 14 -0.0169720621 0.035202665 20000 -0.1150580032
#> 42 indirect 14 0.0781985380 0.044048621 20000 0.0327609539
#> 43 total 15 0.0450029669 0.032072607 20000 -0.0058788156
#> 44 direct 15 -0.0124605035 0.031785753 20000 -0.0982194090
#> 45 indirect 15 0.0574634704 0.038892270 20000 0.0190239591
#> 46 total 16 0.0327816587 0.028993027 20000 -0.0120634312
#> 47 direct 16 -0.0091452818 0.028843412 20000 -0.0841504476
#> 48 indirect 16 0.0419269404 0.034584302 20000 0.0097667293
#> 49 total 17 0.0236880501 0.025836064 20000 -0.0148009847
#> 50 direct 17 -0.0067106920 0.026319797 20000 -0.0722251981
#> 51 indirect 17 0.0303987420 0.030961751 20000 0.0036101288
#> 52 total 18 0.0169927681 0.022728366 20000 -0.0154860134
#> 53 direct 18 -0.0049235432 0.024159997 20000 -0.0620420696
#> 54 indirect 18 0.0219163112 0.027907637 20000 -0.0001943472
#> 55 total 19 0.0121085470 0.019775979 20000 -0.0146231162
#> 56 direct 19 -0.0036120109 0.022312757 20000 -0.0532257463
#> 57 indirect 19 0.0157205578 0.025336669 20000 -0.0020962834
#> 58 total 20 0.0085745994 0.017055638 20000 -0.0133194565
#> 59 direct 20 -0.0026496878 0.020731814 20000 -0.0455406810
#> 60 indirect 20 0.0112242872 0.023179959 20000 -0.0029276890
#> 61 total 21 0.0060365041 0.014613316 20000 -0.0113574615
#> 62 direct 21 -0.0019436746 0.019376365 20000 -0.0390713256
#> 63 indirect 21 0.0079801787 0.021375283 20000 -0.0030459149
#> 64 total 22 0.0042259651 0.012467437 20000 -0.0093618627
#> 65 direct 22 -0.0014257435 0.018211019 20000 -0.0335902220
#> 66 indirect 22 0.0056517086 0.019863724 20000 -0.0028534100
#> 67 total 23 0.0029425409 0.010614515 20000 -0.0074313849
#> 68 direct 23 -0.0010458082 0.017205466 20000 -0.0288661941
#> 69 indirect 23 0.0039883491 0.018590563 20000 -0.0025035814
#> 70 total 24 0.0020381454 0.009035612 20000 -0.0057934226
#> 71 direct 24 -0.0007671105 0.016333991 20000 -0.0245503216
#> 72 indirect 24 0.0028052560 0.017507676 20000 -0.0020929306
#> 73 total 25 0.0014044256 0.007702525 20000 -0.0044948837
#> 74 direct 25 -0.0005626791 0.015574944 20000 -0.0210464235
#> 75 indirect 25 0.0019671047 0.016575358 20000 -0.0017225641
#> 76 total 26 0.0009627744 0.006583013 20000 -0.0033840531
#> 77 direct 26 -0.0004127257 0.014910197 20000 -0.0179831388
#> 78 indirect 26 0.0013755001 0.015762759 20000 -0.0013477127
#> 79 total 27 0.0006565995 0.005644745 20000 -0.0025487413
#> 80 direct 27 -0.0003027338 0.014324638 20000 -0.0154060886
#> 81 indirect 27 0.0009593334 0.015047022 20000 -0.0010834225
#> 82 total 28 0.0004454423 0.004857855 20000 -0.0019041610
#> 83 direct 28 -0.0002220545 0.013805703 20000 -0.0131734129
#> 84 indirect 28 0.0006674968 0.014411665 20000 -0.0008264481
#> 85 total 29 0.0003005636 0.004196261 20000 -0.0013960430
#> 86 direct 29 -0.0001628762 0.013342970 20000 -0.0113119172
#> 87 indirect 29 0.0004634398 0.013844790 20000 -0.0006215314
#> 88 total 30 0.0002016726 0.003638034 20000 -0.0010159622
#> 89 direct 30 -0.0001194690 0.012927801 20000 -0.0096748257
#> 90 indirect 30 0.0003211416 0.013337503 20000 -0.0004614661
#> 97.5%
#> 1 -7.022729e-02
#> 2 -2.367113e-01
#> 3 2.418280e-01
#> 4 1.267040e-01
#> 5 -2.770975e-01
#> 6 5.854145e-01
#> 7 3.185299e-01
#> 8 -2.412815e-01
#> 9 8.195278e-01
#> 10 4.373799e-01
#> 11 -1.850781e-01
#> 12 9.181567e-01
#> 13 4.843004e-01
#> 14 -1.329373e-01
#> 15 9.205153e-01
#> 16 4.784531e-01
#> 17 -9.125614e-02
#> 18 8.588052e-01
#> 19 4.415058e-01
#> 20 -6.090528e-02
#> 21 7.654345e-01
#> 22 3.909504e-01
#> 23 -3.962517e-02
#> 24 6.551271e-01
#> 25 3.373563e-01
#> 26 -2.523033e-02
#> 27 5.468551e-01
#> 28 2.866119e-01
#> 29 -1.595525e-02
#> 30 4.494650e-01
#> 31 2.410517e-01
#> 32 -9.917266e-03
#> 33 3.661131e-01
#> 34 2.012933e-01
#> 35 -6.114857e-03
#> 36 2.960791e-01
#> 37 1.692625e-01
#> 38 -3.657192e-03
#> 39 2.394135e-01
#> 40 1.427016e-01
#> 41 -2.167826e-03
#> 42 1.953933e-01
#> 43 1.198719e-01
#> 44 -1.267497e-03
#> 45 1.597108e-01
#> 46 1.017272e-01
#> 47 -7.194679e-04
#> 48 1.325475e-01
#> 49 8.645041e-02
#> 50 -4.046531e-04
#> 51 1.101881e-01
#> 52 7.351062e-02
#> 53 -2.214863e-04
#> 54 9.311858e-02
#> 55 6.228402e-02
#> 56 -1.205276e-04
#> 57 7.896397e-02
#> 58 5.288233e-02
#> 59 -6.364049e-05
#> 60 6.724753e-02
#> 61 4.505398e-02
#> 62 -3.232021e-05
#> 63 5.782323e-02
#> 64 3.825172e-02
#> 65 -1.643241e-05
#> 66 4.980618e-02
#> 67 3.265107e-02
#> 68 -8.024309e-06
#> 69 4.328525e-02
#> 70 2.783807e-02
#> 71 -3.361541e-06
#> 72 3.744309e-02
#> 73 2.363482e-02
#> 74 -1.307609e-06
#> 75 3.257179e-02
#> 76 2.024984e-02
#> 77 -4.019225e-07
#> 78 2.849622e-02
#> 79 1.727727e-02
#> 80 3.986087e-09
#> 81 2.495068e-02
#> 82 1.472298e-02
#> 83 1.247700e-07
#> 84 2.178099e-02
#> 85 1.250671e-02
#> 86 1.059859e-07
#> 87 1.890333e-02
#> 88 1.065633e-02
#> 89 8.426307e-08
#> 90 1.652079e-02
Parametric Bootstrap Method
pb <- BootParaStdXMY(boot = boot, theta_hat = theta_hat)
plot(pb)
plot(pb, type = "bc")
summary(pb)
#> effect interval est se R 2.5% 97.5%
#> 1 total 1 -0.1364455732 0.031593961 1000 -0.1975991938 -7.233299e-02
#> 2 direct 1 -0.3297415327 0.043923747 1000 -0.4199197586 -2.463023e-01
#> 3 indirect 1 0.1932959594 0.021492058 1000 0.1546151472 2.375674e-01
#> 4 total 2 0.0594752696 0.033873719 1000 -0.0065023692 1.232890e-01
#> 5 direct 2 -0.4003594595 0.063984852 1000 -0.5290582788 -2.827258e-01
#> 6 indirect 2 0.4598347291 0.054673025 1000 0.3611191462 5.728244e-01
#> 7 total 3 0.2520939102 0.032834865 1000 0.1881148988 3.169078e-01
#> 8 direct 3 -0.3698451919 0.073904772 1000 -0.5236261765 -2.413211e-01
#> 9 indirect 3 0.6219391021 0.082658130 1000 0.4768460307 7.870564e-01
#> 10 total 4 0.3631995292 0.035259574 1000 0.2944976476 4.312574e-01
#> 11 direct 4 -0.3078984260 0.077074879 1000 -0.4714310913 -1.799793e-01
#> 12 indirect 4 0.6710979552 0.100368965 1000 0.5004946245 8.752099e-01
#> 13 total 5 0.3985359430 0.039167706 1000 0.3248311338 4.739805e-01
#> 14 direct 5 -0.2434433603 0.075614135 1000 -0.4158594068 -1.275165e-01
#> 15 indirect 5 0.6419793033 0.107585694 1000 0.4705548047 8.677319e-01
#> 16 total 6 0.3833083417 0.041911551 1000 0.3036323433 4.644002e-01
#> 17 direct 6 -0.1870248610 0.071281595 1000 -0.3540514453 -8.422741e-02
#> 18 indirect 6 0.5703332028 0.106313483 1000 0.4010701742 8.088294e-01
#> 19 total 7 0.3409212998 0.042766092 1000 0.2588435117 4.273353e-01
#> 20 direct 7 -0.1412483935 0.065414729 1000 -0.2992163752 -5.506163e-02
#> 21 indirect 7 0.4821696934 0.099205133 1000 0.3310924345 7.116195e-01
#> 22 total 8 0.2879258411 0.042184471 1000 0.2085719389 3.750422e-01
#> 23 direct 8 -0.1055596464 0.058941305 1000 -0.2530451988 -3.639010e-02
#> 24 indirect 8 0.3934854874 0.088833087 1000 0.2568498691 5.986736e-01
#> 25 total 9 0.2343958649 0.040911280 1000 0.1592620387 3.183185e-01
#> 26 direct 9 -0.0783669961 0.052458201 1000 -0.2114715093 -2.285889e-02
#> 27 indirect 9 0.3127628609 0.077337777 1000 0.1967533119 4.931251e-01
#> 28 total 10 0.1856524532 0.039454001 1000 0.1139613970 2.690459e-01
#> 29 direct 10 -0.0579333864 0.046319663 1000 -0.1766021103 -1.411563e-02
#> 30 indirect 10 0.2435858396 0.066253196 1000 0.1486876245 3.976850e-01
#> 31 total 11 0.1439430161 0.037942222 1000 0.0757162544 2.246408e-01
#> 32 direct 11 -0.0427109039 0.040712399 1000 -0.1484065627 -8.325197e-03
#> 33 indirect 11 0.1866539200 0.056468830 1000 0.1095885319 3.171905e-01
#> 34 total 12 0.1097107495 0.036265358 1000 0.0428822948 1.891057e-01
#> 35 direct 12 -0.0314325594 0.035713060 1000 -0.1250197781 -5.072628e-03
#> 36 indirect 12 0.1411433089 0.048308370 1000 0.0753551435 2.508105e-01
#> 37 total 13 0.0824480421 0.034276544 1000 0.0195708061 1.570603e-01
#> 38 direct 13 -0.0231057664 0.031329730 1000 -0.1052893685 -3.131193e-03
#> 39 indirect 13 0.1055538085 0.041692532 1000 0.0500576170 2.025977e-01
#> 40 total 14 0.0612264758 0.031913041 1000 0.0032380723 1.308252e-01
#> 41 direct 14 -0.0169720621 0.027530637 1000 -0.0890803605 -1.884744e-03
#> 42 indirect 14 0.0781985380 0.036335553 1000 0.0308671522 1.625826e-01
#> 43 total 15 0.0450029669 0.029216311 1000 -0.0072562513 1.105741e-01
#> 44 direct 15 -0.0124605035 0.024263263 1000 -0.0751949316 -1.109461e-03
#> 45 indirect 15 0.0574634704 0.031914301 1000 0.0168796180 1.326180e-01
#> 46 total 16 0.0327816587 0.026299147 1000 -0.0119773600 9.454183e-02
#> 47 direct 16 -0.0091452818 0.021466560 1000 -0.0629136689 -6.261175e-04
#> 48 indirect 16 0.0419269404 0.028167351 1000 0.0082746581 1.076031e-01
#> 49 total 17 0.0236880501 0.023302301 1000 -0.0147545553 8.065505e-02
#> 50 direct 17 -0.0067106920 0.019078405 1000 -0.0531577346 -3.158708e-04
#> 51 indirect 17 0.0303987420 0.024922700 1000 0.0023039160 8.938988e-02
#> 52 total 18 0.0169927681 0.020360570 1000 -0.0143612238 6.754034e-02
#> 53 direct 18 -0.0049235432 0.017039869 1000 -0.0449141198 -1.584464e-04
#> 54 indirect 18 0.0219163112 0.022082178 1000 -0.0012670151 7.684385e-02
#> 55 total 19 0.0121085470 0.017582899 1000 -0.0127964584 5.687202e-02
#> 56 direct 19 -0.0036120109 0.015297431 1000 -0.0378556576 -9.214649e-05
#> 57 indirect 19 0.0157205578 0.019592504 1000 -0.0028310622 6.462320e-02
#> 58 total 20 0.0085745994 0.015044445 1000 -0.0117764447 4.860998e-02
#> 59 direct 20 -0.0026496878 0.013803885 1000 -0.0317038782 -5.377416e-05
#> 60 indirect 20 0.0112242872 0.017420414 1000 -0.0030884104 5.446372e-02
#> 61 total 21 0.0060365041 0.012786798 1000 -0.0096565280 4.138777e-02
#> 62 direct 21 -0.0019436746 0.012518480 1000 -0.0265603817 -2.818491e-05
#> 63 indirect 21 0.0079801787 0.015537397 1000 -0.0030077251 4.597434e-02
#> 64 total 22 0.0042259651 0.010822886 1000 -0.0075348315 3.537892e-02
#> 65 direct 22 -0.0014257435 0.011406595 1000 -0.0224531960 -1.126313e-05
#> 66 indirect 22 0.0056517086 0.013913206 1000 -0.0029677099 3.856225e-02
#> 67 total 23 0.0029425409 0.009143959 1000 -0.0061634235 2.956641e-02
#> 68 direct 23 -0.0010458082 0.010439184 1000 -0.0190411051 -5.892516e-06
#> 69 indirect 23 0.0039883491 0.012515148 1000 -0.0025657418 3.314683e-02
#> 70 total 24 0.0020381454 0.007726920 1000 -0.0049703741 2.511935e-02
#> 71 direct 24 -0.0007671105 0.009592114 1000 -0.0161475337 -2.839721e-06
#> 72 indirect 24 0.0028052560 0.011310100 1000 -0.0021866391 2.883148e-02
#> 73 total 25 0.0014044256 0.006540906 1000 -0.0038881738 2.137321e-02
#> 74 direct 25 -0.0005626791 0.008845483 1000 -0.0136936856 -1.094864e-06
#> 75 indirect 25 0.0019671047 0.010267020 1000 -0.0016842904 2.452729e-02
#> 76 total 26 0.0009627744 0.005552560 1000 -0.0030291954 1.818521e-02
#> 77 direct 26 -0.0004127257 0.008182965 1000 -0.0116127378 -3.908304e-07
#> 78 indirect 26 0.0013755001 0.009358728 1000 -0.0012854013 2.068919e-02
#> 79 total 27 0.0006565995 0.004729719 1000 -0.0023054117 1.546882e-02
#> 80 direct 27 -0.0003027338 0.007591221 1000 -0.0098480221 1.133570e-08
#> 81 indirect 27 0.0009593334 0.008562613 1000 -0.0009485646 1.755764e-02
#> 82 total 28 0.0004454423 0.004043602 1000 -0.0016447241 1.315819e-02
#> 83 direct 28 -0.0002220545 0.007059370 1000 -0.0083068821 4.048862e-08
#> 84 indirect 28 0.0006674968 0.007860453 1000 -0.0006677754 1.505159e-02
#> 85 total 29 0.0003005636 0.003469717 1000 -0.0011515120 1.119270e-02
#> 86 direct 29 -0.0001628762 0.006578538 1000 -0.0070162085 2.527131e-08
#> 87 indirect 29 0.0004634398 0.007237723 1000 -0.0005486341 1.294006e-02
#> 88 total 30 0.0002016726 0.002987889 1000 -0.0008296492 9.520806e-03
#> 89 direct 30 -0.0001194690 0.006141484 1000 -0.0059542765 3.198615e-08
#> 90 indirect 30 0.0003211416 0.006682773 1000 -0.0003969616 1.107222e-02
summary(pb, type = "bc")
#> effect interval est se R 2.5% 97.5%
#> 1 total 1 -0.1364455732 0.031593961 1000 -0.2001150637 -7.417674e-02
#> 2 direct 1 -0.3297415327 0.043923747 1000 -0.4212635669 -2.474501e-01
#> 3 indirect 1 0.1932959594 0.021492058 1000 0.1546769534 2.377262e-01
#> 4 total 2 0.0594752696 0.033873719 1000 -0.0066916122 1.229905e-01
#> 5 direct 2 -0.4003594595 0.063984852 1000 -0.5430702065 -2.873910e-01
#> 6 indirect 2 0.4598347291 0.054673025 1000 0.3689392176 5.800168e-01
#> 7 total 3 0.2520939102 0.032834865 1000 0.1887390795 3.170547e-01
#> 8 direct 3 -0.3698451919 0.073904772 1000 -0.5293064214 -2.472064e-01
#> 9 indirect 3 0.6219391021 0.082658130 1000 0.4890206086 8.034230e-01
#> 10 total 4 0.3631995292 0.035259574 1000 0.2987371222 4.361692e-01
#> 11 direct 4 -0.3078984260 0.077074879 1000 -0.4836615828 -1.866888e-01
#> 12 indirect 4 0.6710979552 0.100368965 1000 0.5098083893 8.975428e-01
#> 13 total 5 0.3985359430 0.039167706 1000 0.3296063352 4.822814e-01
#> 14 direct 5 -0.2434433603 0.075614135 1000 -0.4214048093 -1.342398e-01
#> 15 indirect 5 0.6419793033 0.107585694 1000 0.4737499172 8.911837e-01
#> 16 total 6 0.3833083417 0.041911551 1000 0.3075847717 4.677228e-01
#> 17 direct 6 -0.1870248610 0.071281595 1000 -0.3668076272 -9.056218e-02
#> 18 indirect 6 0.5703332028 0.106313483 1000 0.4102759267 8.201475e-01
#> 19 total 7 0.3409212998 0.042766092 1000 0.2619826331 4.287256e-01
#> 20 direct 7 -0.1412483935 0.065414729 1000 -0.3103752165 -6.074976e-02
#> 21 indirect 7 0.4821696934 0.099205133 1000 0.3314423328 7.120999e-01
#> 22 total 8 0.2879258411 0.042184471 1000 0.2122553502 3.767076e-01
#> 23 direct 8 -0.1055596464 0.058941305 1000 -0.2593218933 -3.814254e-02
#> 24 indirect 8 0.3934854874 0.088833087 1000 0.2571208817 6.003871e-01
#> 25 total 9 0.2343958649 0.040911280 1000 0.1643500295 3.235900e-01
#> 26 direct 9 -0.0783669961 0.052458201 1000 -0.2145050657 -2.360276e-02
#> 27 indirect 9 0.3127628609 0.077337777 1000 0.1959935561 4.926099e-01
#> 28 total 10 0.1856524532 0.039454001 1000 0.1198589655 2.721142e-01
#> 29 direct 10 -0.0579333864 0.046319663 1000 -0.1816691948 -1.454534e-02
#> 30 indirect 10 0.2435858396 0.066253196 1000 0.1469159015 3.946416e-01
#> 31 total 11 0.1439430161 0.037942222 1000 0.0786564328 2.282601e-01
#> 32 direct 11 -0.0427109039 0.040712399 1000 -0.1535192299 -8.944975e-03
#> 33 indirect 11 0.1866539200 0.056468830 1000 0.1079382021 3.123179e-01
#> 34 total 12 0.1097107495 0.036265358 1000 0.0466845312 1.913531e-01
#> 35 direct 12 -0.0314325594 0.035713060 1000 -0.1270246988 -5.124805e-03
#> 36 indirect 12 0.1411433089 0.048308370 1000 0.0743690342 2.446450e-01
#> 37 total 13 0.0824480421 0.034276544 1000 0.0202184483 1.608152e-01
#> 38 direct 13 -0.0231057664 0.031329730 1000 -0.1049456585 -3.129980e-03
#> 39 indirect 13 0.1055538085 0.041692532 1000 0.0461553286 1.911369e-01
#> 40 total 14 0.0612264758 0.031913041 1000 0.0046813775 1.324776e-01
#> 41 direct 14 -0.0169720621 0.027530637 1000 -0.0889292953 -1.866745e-03
#> 42 indirect 14 0.0781985380 0.036335553 1000 0.0267189120 1.529708e-01
#> 43 total 15 0.0450029669 0.029216311 1000 -0.0063278742 1.124642e-01
#> 44 direct 15 -0.0124605035 0.024263263 1000 -0.0748429446 -1.092589e-03
#> 45 indirect 15 0.0574634704 0.031914301 1000 0.0122747636 1.207636e-01
#> 46 total 16 0.0327816587 0.026299147 1000 -0.0110534977 9.675331e-02
#> 47 direct 16 -0.0091452818 0.021466560 1000 -0.0628320003 -5.810800e-04
#> 48 indirect 16 0.0419269404 0.028167351 1000 0.0038667721 9.576623e-02
#> 49 total 17 0.0236880501 0.023302301 1000 -0.0118557897 8.231364e-02
#> 50 direct 17 -0.0067106920 0.019078405 1000 -0.0527265203 -3.026780e-04
#> 51 indirect 17 0.0303987420 0.024922700 1000 -0.0022331747 7.430412e-02
#> 52 total 18 0.0169927681 0.020360570 1000 -0.0130203262 6.973183e-02
#> 53 direct 18 -0.0049235432 0.017039869 1000 -0.0442103494 -1.558855e-04
#> 54 indirect 18 0.0219163112 0.022082178 1000 -0.0057285190 5.934170e-02
#> 55 total 19 0.0121085470 0.017582899 1000 -0.0126691297 5.750018e-02
#> 56 direct 19 -0.0036120109 0.015297431 1000 -0.0375147139 -8.938810e-05
#> 57 indirect 19 0.0157205578 0.019592504 1000 -0.0059520447 4.771417e-02
#> 58 total 20 0.0085745994 0.015044445 1000 -0.0112613240 5.035818e-02
#> 59 direct 20 -0.0026496878 0.013803885 1000 -0.0315568686 -5.177949e-05
#> 60 indirect 20 0.0112242872 0.017420414 1000 -0.0057110714 3.880493e-02
#> 61 total 21 0.0060365041 0.012786798 1000 -0.0096026927 4.341762e-02
#> 62 direct 21 -0.0019436746 0.012518480 1000 -0.0265413366 -2.623903e-05
#> 63 indirect 21 0.0079801787 0.015537397 1000 -0.0052101888 3.103701e-02
#> 64 total 22 0.0042259651 0.010822886 1000 -0.0075037606 3.717111e-02
#> 65 direct 22 -0.0014257435 0.011406595 1000 -0.0224159855 -1.126054e-05
#> 66 indirect 22 0.0056517086 0.013913206 1000 -0.0045452416 2.470047e-02
#> 67 total 23 0.0029425409 0.009143959 1000 -0.0057979963 3.190040e-02
#> 68 direct 23 -0.0010458082 0.010439184 1000 -0.0191359574 -6.083640e-06
#> 69 indirect 23 0.0039883491 0.012515148 1000 -0.0044108811 1.965244e-02
#> 70 total 24 0.0020381454 0.007726920 1000 -0.0047911922 2.642192e-02
#> 71 direct 24 -0.0007671105 0.009592114 1000 -0.0160846373 -2.675314e-06
#> 72 indirect 24 0.0028052560 0.011310100 1000 -0.0038603617 1.506565e-02
#> 73 total 25 0.0014044256 0.006540906 1000 -0.0037098687 2.295770e-02
#> 74 direct 25 -0.0005626791 0.008845483 1000 -0.0136936856 -1.094864e-06
#> 75 indirect 25 0.0019671047 0.010267020 1000 -0.0031879860 1.204300e-02
#> 76 total 26 0.0009627744 0.005552560 1000 -0.0030266436 1.840484e-02
#> 77 direct 26 -0.0004127257 0.008182965 1000 -0.0116127378 -3.908304e-07
#> 78 indirect 26 0.0013755001 0.009358728 1000 -0.0028239921 9.402269e-03
#> 79 total 27 0.0006565995 0.004729719 1000 -0.0023280478 1.531298e-02
#> 80 direct 27 -0.0003027338 0.007591221 1000 -0.0098559634 -1.790600e-08
#> 81 indirect 27 0.0009593334 0.008562613 1000 -0.0022528378 7.480420e-03
#> 82 total 28 0.0004454423 0.004043602 1000 -0.0017457588 1.243853e-02
#> 83 direct 28 -0.0002220545 0.007059370 1000 -0.0083182963 3.493473e-08
#> 84 indirect 28 0.0006674968 0.007860453 1000 -0.0017702733 5.664130e-03
#> 85 total 29 0.0003005636 0.003469717 1000 -0.0013394338 9.153354e-03
#> 86 direct 29 -0.0001628762 0.006578538 1000 -0.0070162085 2.527131e-08
#> 87 indirect 29 0.0004634398 0.007237723 1000 -0.0014185891 4.204394e-03
#> 88 total 30 0.0002016726 0.002987889 1000 -0.0009785316 7.252425e-03
#> 89 direct 30 -0.0001194690 0.006141484 1000 -0.0059282504 3.276935e-08
#> 90 indirect 30 0.0003211416 0.006682773 1000 -0.0011736425 3.239673e-03
Model
Delta Method
delta <- DeltaStdYMX(theta_hat = theta_hat)
plot(delta)
summary(delta)
#> effect interval est se z p
#> 1 total 1 3.909316e-01 0.0301423089 12.96952993 1.821439e-38
#> 2 direct 1 3.910746e-01 0.0318130176 12.29291079 9.888030e-35
#> 3 indirect 1 -1.430098e-04 0.0043397564 -0.03295341 9.737117e-01
#> 4 total 2 5.103374e-01 0.0403560142 12.64588097 1.178956e-36
#> 5 direct 2 5.095388e-01 0.0412695904 12.34659237 5.081871e-35
#> 6 indirect 2 7.985419e-04 0.0062991797 0.12676919 8.991231e-01
#> 7 total 3 5.040020e-01 0.0449964720 11.20092190 4.035123e-29
#> 8 direct 3 5.025374e-01 0.0450474596 11.15573193 6.713803e-29
#> 9 indirect 3 1.464586e-03 0.0045528069 0.32168852 7.476887e-01
#> 10 total 4 4.452000e-01 0.0462614445 9.62356431 6.358868e-22
#> 11 direct 4 4.445387e-01 0.0491049274 9.05283363 1.393054e-19
#> 12 indirect 4 6.612487e-04 0.0123124680 0.05370562 9.571697e-01
#> 13 total 5 3.702766e-01 0.0458508047 8.07568317 6.709990e-16
#> 14 direct 5 3.718598e-01 0.0552352316 6.73229415 1.670086e-11
#> 15 indirect 5 -1.583255e-03 0.0258309865 -0.06129288 9.511260e-01
#> 16 total 6 2.964921e-01 0.0450235860 6.58526094 4.540872e-11
#> 17 direct 6 3.010943e-01 0.0616475611 4.88412354 1.038899e-06
#> 18 indirect 6 -4.602242e-03 0.0398509218 -0.11548646 9.080596e-01
#> 19 total 7 2.312213e-01 0.0443158266 5.21757945 1.812765e-07
#> 20 direct 7 2.388786e-01 0.0661984874 3.60852057 3.079481e-04
#> 21 indirect 7 -7.657257e-03 0.0518503297 -0.14768001 8.825953e-01
#> 22 total 8 1.768004e-01 0.0436268372 4.05256013 5.066020e-05
#> 23 direct 8 1.870132e-01 0.0679235731 2.75328865 5.899985e-03
#> 24 indirect 8 -1.021282e-02 0.0604819729 -0.16885729 8.659089e-01
#> 25 total 9 1.331107e-01 0.0425855439 3.12572592 1.773668e-03
#> 26 direct 9 1.451059e-01 0.0668570618 2.17039003 2.997731e-02
#> 27 indirect 9 -1.199516e-02 0.0653408771 -0.18357822 8.543444e-01
#> 28 total 10 9.895651e-02 0.0408921644 2.41993810 1.552315e-02
#> 29 direct 10 1.119028e-01 0.0635407232 1.76112002 7.821809e-02
#> 30 indirect 10 -1.294633e-02 0.0666768162 -0.19416544 8.460463e-01
#> 31 total 11 7.278402e-02 0.0384544412 1.89273378 5.839328e-02
#> 32 direct 11 8.593217e-02 0.0586709527 1.46464582 1.430176e-01
#> 33 indirect 11 -1.314815e-02 0.0651056435 -0.20195094 8.399551e-01
#> 34 total 12 5.304095e-02 0.0353645644 1.49983331 1.336576e-01
#> 35 direct 12 6.579349e-02 0.0529069663 1.24356942 2.136581e-01
#> 36 indirect 12 -1.275253e-02 0.0613812738 -0.20775935 8.354169e-01
#> 37 total 13 3.833844e-02 0.0318201245 1.20484876 2.282617e-01
#> 38 direct 13 5.026946e-02 0.0467880151 1.07440901 2.826394e-01
#> 39 indirect 13 -1.193103e-02 0.0562441758 -0.21212912 8.320063e-01
#> 40 total 14 2.750793e-02 0.0280502897 0.98066458 3.267582e-01
#> 41 direct 14 3.835179e-02 0.0407129157 0.94200555 3.461898e-01
#> 42 indirect 14 -1.084387e-02 0.0503369598 -0.21542555 8.294356e-01
#> 43 total 15 1.960435e-02 0.0242666300 0.80787293 4.191637e-01
#> 44 direct 15 2.922897e-02 0.0349506911 0.83629168 4.029908e-01
#> 45 indirect 15 -9.624619e-03 0.0441691930 -0.21790343 8.275044e-01
#> 46 total 16 1.388442e-02 0.0206380241 0.67275906 5.011006e-01
#> 47 direct 16 2.225970e-02 0.0296637718 0.75040019 4.530137e-01
#> 48 indirect 16 -8.375283e-03 0.0381138236 -0.21974396 8.260706e-01
#> 49 total 17 9.775619e-03 0.0172829780 0.56562120 5.716513e-01
#> 50 direct 17 1.694322e-02 0.0249337092 0.67953069 4.968017e-01
#> 51 indirect 17 -7.167602e-03 0.0324212072 -0.22107757 8.250320e-01
#> 52 total 18 6.844225e-03 0.0142724461 0.47954114 6.315537e-01
#> 53 direct 18 1.289168e-02 0.0207845734 0.62025226 5.350917e-01
#> 54 indirect 18 -6.047454e-03 0.0272409363 -0.22199874 8.243149e-01
#> 55 total 19 4.766034e-03 0.0116378395 0.40952912 6.821514e-01
#> 56 direct 19 9.806335e-03 0.0172021853 0.57006331 5.686348e-01
#> 57 indirect 19 -5.040301e-03 0.0226453118 -0.22257590 8.238656e-01
#> 58 total 20 3.301445e-03 0.0093807066 0.35193986 7.248834e-01
#> 59 direct 20 7.457977e-03 0.0141488971 0.52710662 5.981195e-01
#> 60 indirect 20 -4.156533e-03 0.0186510096 -0.22285833 8.236458e-01
#> 61 total 21 2.275101e-03 0.0074820391 0.30407494 7.610708e-01
#> 62 direct 21 5.671217e-03 0.0115743694 0.48998066 6.241476e-01
#> 63 indirect 21 -3.396117e-03 0.0152373467 -0.22288110 8.236280e-01
#> 64 total 22 1.559764e-03 0.0059101963 0.26391068 7.918487e-01
#> 65 direct 22 4.312105e-03 0.0094230587 0.45761201 6.472312e-01
#> 66 indirect 22 -2.752341e-03 0.0123606989 -0.22266871 8.237934e-01
#> 67 total 23 1.063818e-03 0.0046271022 0.22991026 8.181615e-01
#> 68 direct 23 3.278478e-03 0.0076391595 0.42916732 6.678015e-01
#> 69 indirect 23 -2.214659e-03 0.0099652651 -0.22223788 8.241287e-01
#> 70 total 24 7.217575e-04 0.0035927523 0.20089263 8.407825e-01
#> 71 direct 24 2.492491e-03 0.0061696609 0.40399162 6.862189e-01
#> 72 indirect 24 -1.770734e-03 0.0079906878 -0.22159968 8.246255e-01
#> 73 total 25 4.870462e-04 0.0027682541 0.17593983 8.603412e-01
#> 74 direct 25 1.894871e-03 0.0049660618 0.38156421 7.027846e-01
#> 75 indirect 25 -1.407825e-03 0.0063771379 -0.22076130 8.252783e-01
#> 76 total 26 3.268255e-04 0.0021176940 0.15433083 8.773489e-01
#> 77 direct 26 1.440505e-03 0.0039851668 0.36146677 7.177505e-01
#> 78 indirect 26 -1.113680e-03 0.0050684617 -0.21972740 8.260835e-01
#> 79 total 27 2.180262e-04 0.0016091255 0.13549360 8.922216e-01
#> 80 direct 27 1.095071e-03 0.0031892790 0.34335999 7.313277e-01
#> 81 indirect 27 -8.770446e-04 0.0040139109 -0.21850127 8.270386e-01
#> 82 total 28 1.445409e-04 0.0012149357 0.11897003 9.052991e-01
#> 83 direct 28 8.324611e-04 0.0025460163 0.32696613 7.436935e-01
#> 84 indirect 28 -6.879202e-04 0.0031688870 -0.21708574 8.281415e-01
#> 85 total 29 9.518284e-05 0.0009118003 0.10439001 9.168599e-01
#> 86 direct 29 6.328222e-04 0.0020279126 0.31205595 7.549980e-01
#> 87 indirect 29 -5.376394e-04 0.0024950322 -0.21548394 8.293900e-01
#> 88 total 30 6.222295e-05 0.0006803892 0.09145200 9.271334e-01
#> 89 direct 30 4.810571e-04 0.0016119123 0.29843875 7.653683e-01
#> 90 indirect 30 -4.188341e-04 0.0019599180 -0.21369984 8.307811e-01
#> 2.5% 97.5%
#> 1 0.331853738 0.450009417
#> 2 0.328722218 0.453426956
#> 3 -0.008648776 0.008362756
#> 4 0.431241018 0.589433686
#> 5 0.428651899 0.590425721
#> 6 -0.011547623 0.013144707
#> 7 0.415810504 0.592193434
#> 8 0.414245985 0.590828782
#> 9 -0.007458752 0.010387923
#> 10 0.354529221 0.535870752
#> 11 0.348294849 0.540782627
#> 12 -0.023470745 0.024793243
#> 13 0.280410646 0.460142497
#> 14 0.263600762 0.480118892
#> 15 -0.052211059 0.049044548
#> 16 0.208247455 0.384736670
#> 17 0.180267305 0.421921304
#> 18 -0.082708613 0.073504130
#> 19 0.144363922 0.318078770
#> 20 0.109131952 0.368625255
#> 21 -0.109282036 0.093967522
#> 22 0.091293351 0.262307411
#> 23 0.053885446 0.320140960
#> 24 -0.128755311 0.108329667
#> 25 0.049644606 0.216576871
#> 26 0.014068467 0.276143333
#> 27 -0.140060927 0.116070604
#> 28 0.018809337 0.179103676
#> 29 -0.012634689 0.236440369
#> 30 -0.143630492 0.117737825
#> 31 -0.002585300 0.148153340
#> 32 -0.029060789 0.200925120
#> 33 -0.140752862 0.114456571
#> 34 -0.016272321 0.122354224
#> 35 -0.037902263 0.169489234
#> 36 -0.133057620 0.107552552
#> 37 -0.024027861 0.100704735
#> 38 -0.041433360 0.141972289
#> 39 -0.122167586 0.098305532
#> 40 -0.027469632 0.082485483
#> 41 -0.041444056 0.118147641
#> 42 -0.109502495 0.087814761
#> 43 -0.027957367 0.067166074
#> 44 -0.039273124 0.097731068
#> 45 -0.096194646 0.076945409
#> 46 -0.026565366 0.054334202
#> 47 -0.035880224 0.080399625
#> 48 -0.083077004 0.066326439
#> 49 -0.024098396 0.043649633
#> 50 -0.031925951 0.065812393
#> 51 -0.070712000 0.056376797
#> 52 -0.021129255 0.034817705
#> 53 -0.027845337 0.053628694
#> 54 -0.059438707 0.047343800
#> 55 -0.018043712 0.027575780
#> 56 -0.023909329 0.043521998
#> 57 -0.049424296 0.039343695
#> 58 -0.015084403 0.021687292
#> 59 -0.020273351 0.035189306
#> 60 -0.040711840 0.032398774
#> 61 -0.012389427 0.016939628
#> 62 -0.017014130 0.028356564
#> 63 -0.033260767 0.026468534
#> 64 -0.010024008 0.013143536
#> 65 -0.014156751 0.022780960
#> 66 -0.026978866 0.021474184
#> 67 -0.008005135 0.010132772
#> 68 -0.011694000 0.018250955
#> 69 -0.021746220 0.017316901
#> 70 -0.006319908 0.007763423
#> 71 -0.009599822 0.014584805
#> 72 -0.017432194 0.013890726
#> 73 -0.004938632 0.005912724
#> 74 -0.007838431 0.011628174
#> 75 -0.013906786 0.011091135
#> 76 -0.003823778 0.004477429
#> 77 -0.006370278 0.009251289
#> 78 -0.011047682 0.008820322
#> 79 -0.002935802 0.003371854
#> 80 -0.005155801 0.007345943
#> 81 -0.008744165 0.006990076
#> 82 -0.002236689 0.002525771
#> 83 -0.004157639 0.005822561
#> 84 -0.006898825 0.005522984
#> 85 -0.001691913 0.001882279
#> 86 -0.003341814 0.004607458
#> 87 -0.005427813 0.004352534
#> 88 -0.001271315 0.001395761
#> 89 -0.002678233 0.003640347
#> 90 -0.004260203 0.003422534
Monte Carlo Method
summary(mc)
#> effect interval est se R 2.5%
#> 1 total 1 -1.064193e-02 0.0412137844 20000 -0.0903676659
#> 2 direct 1 -1.065335e-02 0.0411676271 20000 -0.0886639271
#> 3 indirect 1 1.142337e-05 0.0013552845 20000 -0.0038547767
#> 4 total 2 -1.297395e-02 0.0506955441 20000 -0.1121523880
#> 5 direct 2 -1.293489e-02 0.0503988964 20000 -0.1069204842
#> 6 indirect 2 -3.905984e-05 0.0036015956 20000 -0.0106742220
#> 7 total 3 -1.206546e-02 0.0477294762 20000 -0.1089461488
#> 8 direct 3 -1.194903e-02 0.0473168354 20000 -0.1011747434
#> 9 indirect 3 -1.164308e-04 0.0056226255 20000 -0.0164375096
#> 10 total 4 -1.011461e-02 0.0405829646 20000 -0.0944772416
#> 11 direct 4 -9.947643e-03 0.0403656546 20000 -0.0890951901
#> 12 indirect 4 -1.669635e-04 0.0070611404 20000 -0.0195525770
#> 13 total 5 -8.039534e-03 0.0327413523 20000 -0.0774853334
#> 14 direct 5 -7.865216e-03 0.0330106571 20000 -0.0748434397
#> 15 indirect 5 -1.743175e-04 0.0079835555 20000 -0.0199958362
#> 16 total 6 -6.189385e-03 0.0255856629 20000 -0.0603950764
#> 17 direct 6 -6.042436e-03 0.0265146119 20000 -0.0621751947
#> 18 indirect 6 -1.469486e-04 0.0085841696 20000 -0.0183318432
#> 19 total 7 -4.664497e-03 0.0195843461 20000 -0.0459345077
#> 20 direct 7 -4.563481e-03 0.0211957578 20000 -0.0516036859
#> 21 indirect 7 -1.010165e-04 0.0089872904 20000 -0.0159089006
#> 22 total 8 -3.461331e-03 0.0148148513 20000 -0.0335972130
#> 23 direct 8 -3.410442e-03 0.0169973413 20000 -0.0427460222
#> 24 indirect 8 -5.088905e-05 0.0092168553 20000 -0.0137015807
#> 25 total 9 -2.537880e-03 0.0111853811 20000 -0.0237539629
#> 26 direct 9 -2.531896e-03 0.0137403338 20000 -0.0354016089
#> 27 indirect 9 -5.983648e-06 0.0092521121 20000 -0.0112892474
#> 28 total 10 -1.842634e-03 0.0085280660 20000 -0.0168116371
#> 29 direct 10 -1.871723e-03 0.0112297523 20000 -0.0294640367
#> 30 indirect 10 2.908908e-05 0.0090794736 20000 -0.0093091299
#> 31 total 11 -1.326708e-03 0.0066410874 20000 -0.0117506764
#> 32 direct 11 -1.379912e-03 0.0092941227 20000 -0.0246439930
#> 33 indirect 11 5.320397e-05 0.0087116678 20000 -0.0075873777
#> 34 total 12 -9.482262e-04 0.0053163947 20000 -0.0082978141
#> 35 direct 12 -1.015529e-03 0.0077955442 20000 -0.0205115596
#> 36 indirect 12 6.730311e-05 0.0081846080 20000 -0.0061212314
#> 37 total 13 -6.732162e-04 0.0043675009 20000 -0.0059045950
#> 38 direct 13 -7.465057e-04 0.0066279241 20000 -0.0172132410
#> 39 indirect 13 7.328946e-05 0.0075468245 20000 -0.0049032742
#> 40 total 14 -4.750334e-04 0.0036524313 20000 -0.0042613667
#> 41 direct 14 -5.483367e-04 0.0057114007 20000 -0.0142678921
#> 42 indirect 14 7.330331e-05 0.0068493957 20000 -0.0039145875
#> 43 total 15 -3.332586e-04 0.0030804456 20000 -0.0030808573
#> 44 direct 15 -4.025764e-04 0.0049863730 20000 -0.0120290643
#> 45 indirect 15 6.931784e-05 0.0061388844 20000 -0.0031270968
#> 46 total 16 -2.325098e-04 0.0026026091 20000 -0.0022727717
#> 47 direct 16 -2.954676e-04 0.0044083816 20000 -0.0100134115
#> 48 indirect 16 6.295775e-05 0.0054533451 20000 -0.0025492219
#> 49 total 17 -1.613554e-04 0.0021962623 20000 -0.0017078171
#> 50 direct 17 -2.168104e-04 0.0039441154 20000 -0.0083397222
#> 51 indirect 17 5.545493e-05 0.0048207016 20000 -0.0020566611
#> 52 total 18 -1.113930e-04 0.0018516713 20000 -0.0012870824
#> 53 direct 18 -1.590708e-04 0.0035684405 20000 -0.0070420534
#> 54 indirect 18 4.767783e-05 0.0042587200 20000 -0.0016658853
#> 55 total 19 -7.650453e-05 0.0015634709 20000 -0.0009790171
#> 56 direct 19 -1.166976e-04 0.0032622395 20000 -0.0060311066
#> 57 indirect 19 4.019303e-05 0.0037759560 20000 -0.0013472070
#> 58 total 20 -5.227212e-05 0.0013263176 20000 -0.0007464966
#> 59 direct 20 -8.560664e-05 0.0030108529 20000 -0.0050952335
#> 60 indirect 20 3.333451e-05 0.0033732619 20000 -0.0010870276
#> 61 total 21 -3.552885e-05 0.0011334368 20000 -0.0005847391
#> 62 direct 21 -6.279662e-05 0.0028029550 20000 -0.0042973714
#> 63 indirect 21 2.726777e-05 0.0030456143 20000 -0.0008798495
#> 64 total 22 -2.401981e-05 0.0009769039 20000 -0.0004787708
#> 65 direct 22 -4.606320e-05 0.0026297365 20000 -0.0036119949
#> 66 indirect 22 2.204340e-05 0.0027841007 20000 -0.0007013338
#> 67 total 23 -1.614944e-05 0.0008487117 20000 -0.0004033328
#> 68 direct 23 -3.378818e-05 0.0024843063 20000 -0.0030373359
#> 69 indirect 23 1.763874e-05 0.0025778753 20000 -0.0005659080
#> 70 total 24 -1.079535e-05 0.0007418526 20000 -0.0003437172
#> 71 direct 24 -2.478396e-05 0.0023612465 20000 -0.0025618773
#> 72 indirect 24 1.398861e-05 0.0024158437 20000 -0.0004591081
#> 73 total 25 -7.172379e-06 0.0006509408 20000 -0.0002972903
#> 74 direct 25 -1.817915e-05 0.0022562795 20000 -0.0021516428
#> 75 indirect 25 1.100677e-05 0.0022878664 20000 -0.0003754139
#> 76 total 26 -4.734244e-06 0.0005722740 20000 -0.0002621608
#> 77 direct 26 -1.333442e-05 0.0021660150 20000 -0.0018099702
#> 78 indirect 26 8.600180e-06 0.0021854063 20000 -0.0003045437
#> 79 total 27 -3.102834e-06 0.0005035115 20000 -0.0002304673
#> 80 direct 27 -9.780785e-06 0.0020877564 20000 -0.0015395206
#> 81 indirect 27 6.677951e-06 0.0021017049 20000 -0.0002481945
#> 82 total 28 -2.017802e-06 0.0004432130 20000 -0.0001965841
#> 83 direct 28 -7.174180e-06 0.0020193513 20000 -0.0013085775
#> 84 indirect 28 5.156378e-06 0.0020316545 20000 -0.0002023299
#> 85 total 29 -1.300805e-06 0.0003904205 20000 -0.0001654409
#> 86 direct 29 -5.262236e-06 0.0019590761 20000 -0.0011067193
#> 87 indirect 29 3.961431e-06 0.0019715300 20000 -0.0001660443
#> 88 total 30 -8.303077e-07 0.0003443666 20000 -0.0001374621
#> 89 direct 30 -3.859828e-06 0.0019055453 20000 -0.0009356131
#> 90 indirect 30 3.029520e-06 0.0019186925 20000 -0.0001364863
#> 97.5%
#> 1 0.0716936576
#> 2 0.0725556062
#> 3 0.0018833820
#> 4 0.0868975506
#> 5 0.0914420726
#> 6 0.0041574414
#> 7 0.0790351637
#> 8 0.0862204337
#> 9 0.0069647533
#> 10 0.0651437649
#> 11 0.0725500321
#> 12 0.0105079074
#> 13 0.0520155584
#> 14 0.0581300647
#> 15 0.0142524015
#> 16 0.0411661842
#> 17 0.0448837448
#> 18 0.0178502886
#> 19 0.0330759397
#> 20 0.0338430776
#> 21 0.0210506262
#> 22 0.0271153845
#> 23 0.0258851313
#> 24 0.0230923656
#> 25 0.0227584355
#> 26 0.0193633096
#> 27 0.0243335374
#> 28 0.0192129880
#> 29 0.0146499461
#> 30 0.0249381998
#> 31 0.0163012096
#> 32 0.0112807826
#> 33 0.0240962714
#> 34 0.0140818491
#> 35 0.0086048229
#> 36 0.0227460115
#> 37 0.0121531258
#> 38 0.0065638815
#> 39 0.0209823728
#> 40 0.0105563145
#> 41 0.0050700600
#> 42 0.0188747317
#> 43 0.0091099995
#> 44 0.0039208683
#> 45 0.0166912393
#> 46 0.0078141701
#> 47 0.0030378293
#> 48 0.0144978438
#> 49 0.0066197523
#> 50 0.0023531274
#> 51 0.0124115825
#> 52 0.0055655390
#> 53 0.0018227801
#> 54 0.0105730424
#> 55 0.0046348005
#> 56 0.0014347368
#> 57 0.0089327377
#> 58 0.0038418634
#> 59 0.0011201634
#> 60 0.0074776630
#> 61 0.0032244350
#> 62 0.0008699452
#> 63 0.0061766235
#> 64 0.0027305025
#> 65 0.0006798552
#> 66 0.0051462621
#> 67 0.0023223153
#> 68 0.0005324945
#> 69 0.0043672076
#> 70 0.0019589586
#> 71 0.0004251609
#> 72 0.0036608553
#> 73 0.0016517524
#> 74 0.0003361540
#> 75 0.0030707785
#> 76 0.0013927032
#> 77 0.0002702123
#> 78 0.0025989900
#> 79 0.0011803600
#> 80 0.0002109471
#> 81 0.0022160093
#> 82 0.0010027695
#> 83 0.0001636514
#> 84 0.0019064453
#> 85 0.0008542762
#> 86 0.0001296969
#> 87 0.0016490144
#> 88 0.0007243895
#> 89 0.0001032633
#> 90 0.0014177742
Parametric Bootstrap Method
pb <- BootParaStdYMX(boot = boot, theta_hat = theta_hat)
plot(pb)
plot(pb, type = "bc")
summary(pb)
#> effect interval est se R 2.5% 97.5%
#> 1 total 1 -1.064193e-02 0.0399988609 1000 -8.757136e-02 7.095323e-02
#> 2 direct 1 -1.065335e-02 0.0399120927 1000 -8.573149e-02 7.267064e-02
#> 3 indirect 1 1.142337e-05 0.0012236390 1000 -3.486168e-03 1.588027e-03
#> 4 total 2 -1.297395e-02 0.0487535123 1000 -1.084955e-01 8.339292e-02
#> 5 direct 2 -1.293489e-02 0.0483519751 1000 -1.015839e-01 8.717543e-02
#> 6 indirect 2 -3.905984e-05 0.0032660976 1000 -9.233322e-03 3.075105e-03
#> 7 total 3 -1.206546e-02 0.0454670117 1000 -1.037020e-01 7.582257e-02
#> 8 direct 3 -1.194903e-02 0.0448864644 1000 -9.723642e-02 7.723472e-02
#> 9 indirect 3 -1.164308e-04 0.0050558376 1000 -1.456335e-02 6.090961e-03
#> 10 total 4 -1.011461e-02 0.0382693562 1000 -8.747563e-02 6.298688e-02
#> 11 direct 4 -9.947643e-03 0.0378173207 1000 -8.500670e-02 6.582357e-02
#> 12 indirect 4 -1.669635e-04 0.0062778559 1000 -1.736755e-02 8.559676e-03
#> 13 total 5 -8.039534e-03 0.0305422624 1000 -7.122720e-02 5.066253e-02
#> 14 direct 5 -7.865216e-03 0.0304887246 1000 -6.941307e-02 5.266500e-02
#> 15 indirect 5 -1.743175e-04 0.0070405267 1000 -1.857199e-02 1.265136e-02
#> 16 total 6 -6.189385e-03 0.0236035793 1000 -5.496501e-02 4.100156e-02
#> 17 direct 6 -6.042436e-03 0.0240831561 1000 -5.697437e-02 3.834032e-02
#> 18 indirect 6 -1.469486e-04 0.0075376993 1000 -1.705068e-02 1.671566e-02
#> 19 total 7 -4.664497e-03 0.0178816011 1000 -3.978584e-02 3.153211e-02
#> 20 direct 7 -4.563481e-03 0.0188731813 1000 -4.671912e-02 2.849569e-02
#> 21 indirect 7 -1.010165e-04 0.0078619277 1000 -1.525101e-02 1.929381e-02
#> 22 total 8 -3.461331e-03 0.0134216080 1000 -2.990367e-02 2.616844e-02
#> 23 direct 8 -3.410442e-03 0.0147797530 1000 -3.965005e-02 2.034796e-02
#> 24 indirect 8 -5.088905e-05 0.0080071393 1000 -1.115899e-02 2.156680e-02
#> 25 total 9 -2.537880e-03 0.0101003728 1000 -2.152307e-02 2.161133e-02
#> 26 direct 9 -2.531896e-03 0.0116151545 1000 -3.284854e-02 1.547690e-02
#> 27 indirect 9 -5.983648e-06 0.0079420617 1000 -9.179985e-03 2.305419e-02
#> 28 total 10 -1.842634e-03 0.0077185775 1000 -1.491802e-02 1.764929e-02
#> 29 direct 10 -1.871723e-03 0.0091825414 1000 -2.626458e-02 1.159364e-02
#> 30 indirect 10 2.908908e-05 0.0076589800 1000 -7.408961e-03 2.384183e-02
#> 31 total 11 -1.326708e-03 0.0060482375 1000 -1.025240e-02 1.459010e-02
#> 32 direct 11 -1.379912e-03 0.0073113996 1000 -2.169451e-02 8.672373e-03
#> 33 indirect 11 5.320397e-05 0.0071831649 1000 -6.055820e-03 2.286862e-02
#> 34 total 12 -9.482262e-04 0.0048697265 1000 -7.728072e-03 1.213238e-02
#> 35 direct 12 -1.015529e-03 0.0058656130 1000 -1.795202e-02 6.466236e-03
#> 36 indirect 12 6.730311e-05 0.0065625860 1000 -4.674126e-03 2.172118e-02
#> 37 total 13 -6.732162e-04 0.0040043510 1000 -5.441561e-03 1.004973e-02
#> 38 direct 13 -7.465057e-04 0.0047409033 1000 -1.491806e-02 4.741098e-03
#> 39 indirect 13 7.328946e-05 0.0058540199 1000 -3.673339e-03 1.987903e-02
#> 40 total 14 -4.750334e-04 0.0033308783 1000 -4.041515e-03 9.175988e-03
#> 41 direct 14 -5.483367e-04 0.0038590480 1000 -1.227222e-02 3.663181e-03
#> 42 indirect 14 7.330331e-05 0.0051121862 1000 -2.823236e-03 1.726394e-02
#> 43 total 15 -3.332586e-04 0.0027805955 1000 -3.063103e-03 7.708074e-03
#> 44 direct 15 -4.025764e-04 0.0031619161 1000 -1.013825e-02 2.715757e-03
#> 45 indirect 15 6.931784e-05 0.0043831867 1000 -2.392881e-03 1.458324e-02
#> 46 total 16 -2.325098e-04 0.0023204147 1000 -2.188322e-03 6.469946e-03
#> 47 direct 16 -2.954676e-04 0.0026063882 1000 -8.442392e-03 2.073302e-03
#> 48 indirect 16 6.295775e-05 0.0037015657 1000 -1.919050e-03 1.260366e-02
#> 49 total 17 -1.613554e-04 0.0019358900 1000 -1.552895e-03 5.436758e-03
#> 50 direct 17 -2.168104e-04 0.0021603798 1000 -7.234708e-03 1.583462e-03
#> 51 indirect 17 5.545493e-05 0.0030898810 1000 -1.549831e-03 1.070740e-02
#> 52 total 18 -1.113930e-04 0.0016194598 1000 -1.156927e-03 4.564530e-03
#> 53 direct 18 -1.590708e-04 0.0017998630 1000 -6.122202e-03 1.144526e-03
#> 54 indirect 18 4.767783e-05 0.0025597976 1000 -1.250870e-03 8.977121e-03
#> 55 total 19 -7.650453e-05 0.0013641051 1000 -8.591791e-04 3.805311e-03
#> 56 direct 19 -1.166976e-04 0.0015066886 1000 -4.993066e-03 8.535852e-04
#> 57 indirect 19 4.019303e-05 0.0021139697 1000 -9.863422e-04 7.333612e-03
#> 58 total 20 -5.227212e-05 0.0011610097 1000 -6.768359e-04 2.961810e-03
#> 59 direct 20 -8.560664e-05 0.0012670143 1000 -4.181330e-03 6.446479e-04
#> 60 indirect 20 3.333451e-05 0.0017482372 1000 -8.042840e-04 6.051531e-03
#> 61 total 21 -3.552885e-05 0.0009997955 1000 -5.711536e-04 2.506548e-03
#> 62 direct 21 -6.279662e-05 0.0010701744 1000 -3.429701e-03 5.015311e-04
#> 63 indirect 21 2.726777e-05 0.0014538717 1000 -6.458498e-04 4.856658e-03
#> 64 total 22 -2.401981e-05 0.0008700292 1000 -4.381894e-04 2.078884e-03
#> 65 direct 22 -4.606320e-05 0.0009078678 1000 -2.887124e-03 3.891989e-04
#> 66 indirect 22 2.204340e-05 0.0012197162 1000 -5.187564e-04 3.882149e-03
#> 67 total 23 -1.614944e-05 0.0007628003 1000 -4.026765e-04 1.746854e-03
#> 68 direct 23 -3.378818e-05 0.0007735711 1000 -2.420568e-03 3.102516e-04
#> 69 indirect 23 1.763874e-05 0.0010340746 1000 -4.145167e-04 3.104441e-03
#> 70 total 24 -1.079535e-05 0.0006715757 1000 -3.376401e-04 1.467868e-03
#> 71 direct 24 -2.478396e-05 0.0006621109 1000 -1.981076e-03 2.345646e-04
#> 72 indirect 24 1.398861e-05 0.0008861637 1000 -3.329160e-04 2.652949e-03
#> 73 total 25 -7.172379e-06 0.0005921888 1000 -2.889059e-04 1.233447e-03
#> 74 direct 25 -1.817915e-05 0.0005693511 1000 -1.634928e-03 1.808937e-04
#> 75 indirect 25 1.100677e-05 0.0007669758 1000 -2.768083e-04 2.186495e-03
#> 76 total 26 -4.734244e-06 0.0005222834 1000 -2.346370e-04 1.050077e-03
#> 77 direct 26 -1.333442e-05 0.0004919595 1000 -1.366500e-03 1.379408e-04
#> 78 indirect 26 8.600180e-06 0.0006695342 1000 -2.128582e-04 1.809283e-03
#> 79 total 27 -3.102834e-06 0.0004606122 1000 -2.038816e-04 9.039369e-04
#> 80 direct 27 -9.780785e-06 0.0004272339 1000 -1.141911e-03 1.061544e-04
#> 81 indirect 27 6.677951e-06 0.0005887002 1000 -1.773865e-04 1.507722e-03
#> 82 total 28 -2.017802e-06 0.0004064440 1000 -1.706861e-04 7.781474e-04
#> 83 direct 28 -7.174180e-06 0.0003729703 1000 -9.540211e-04 8.546195e-05
#> 84 indirect 28 5.156378e-06 0.0005207669 1000 -1.491985e-04 1.283343e-03
#> 85 total 29 -1.300805e-06 0.0003591783 1000 -1.342998e-04 6.698727e-04
#> 86 direct 29 -5.262236e-06 0.0003273617 1000 -7.970483e-04 6.670611e-05
#> 87 indirect 29 3.961431e-06 0.0004630382 1000 -1.232282e-04 1.102376e-03
#> 88 total 30 -8.303077e-07 0.0003181590 1000 -1.063455e-04 5.646820e-04
#> 89 direct 30 -3.859828e-06 0.0002889209 1000 -6.659048e-04 5.081528e-05
#> 90 indirect 30 3.029520e-06 0.0004134999 1000 -9.818338e-05 9.506165e-04
summary(pb, type = "bc")
#> effect interval est se R 2.5% 97.5%
#> 1 total 1 -1.064193e-02 0.0399988609 1000 -0.0850940927 7.387274e-02
#> 2 direct 1 -1.065335e-02 0.0399120927 1000 -0.0841193412 7.567809e-02
#> 3 indirect 1 1.142337e-05 0.0012236390 1000 -0.0016389619 2.515644e-03
#> 4 total 2 -1.297395e-02 0.0487535123 1000 -0.1036390561 8.914127e-02
#> 5 direct 2 -1.293489e-02 0.0483519751 1000 -0.0990265132 9.112423e-02
#> 6 indirect 2 -3.905984e-05 0.0032660976 1000 -0.0041499477 7.270168e-03
#> 7 total 3 -1.206546e-02 0.0454670117 1000 -0.0981818469 8.296048e-02
#> 8 direct 3 -1.194903e-02 0.0448864644 1000 -0.0911667193 8.428349e-02
#> 9 indirect 3 -1.164308e-04 0.0050558376 1000 -0.0069018931 1.100015e-02
#> 10 total 4 -1.011461e-02 0.0382693562 1000 -0.0839149911 6.666206e-02
#> 11 direct 4 -9.947643e-03 0.0378173207 1000 -0.0816621641 6.947993e-02
#> 12 indirect 4 -1.669635e-04 0.0062778559 1000 -0.0095942453 1.523312e-02
#> 13 total 5 -8.039534e-03 0.0305422624 1000 -0.0691346244 5.357418e-02
#> 14 direct 5 -7.865216e-03 0.0304887246 1000 -0.0684462987 5.392747e-02
#> 15 indirect 5 -1.743175e-04 0.0070405267 1000 -0.0117841709 1.654568e-02
#> 16 total 6 -6.189385e-03 0.0236035793 1000 -0.0533874516 4.237968e-02
#> 17 direct 6 -6.042436e-03 0.0240831561 1000 -0.0543793629 4.083813e-02
#> 18 indirect 6 -1.469486e-04 0.0075376993 1000 -0.0139509737 1.823683e-02
#> 19 total 7 -4.664497e-03 0.0178816011 1000 -0.0384090648 3.386214e-02
#> 20 direct 7 -4.563481e-03 0.0188731813 1000 -0.0455286180 2.978084e-02
#> 21 indirect 7 -1.010165e-04 0.0078619277 1000 -0.0158869959 1.884033e-02
#> 22 total 8 -3.461331e-03 0.0134216080 1000 -0.0279040930 2.765859e-02
#> 23 direct 8 -3.410442e-03 0.0147797530 1000 -0.0378593996 2.163207e-02
#> 24 indirect 8 -5.088905e-05 0.0080071393 1000 -0.0143150551 1.877408e-02
#> 25 total 9 -2.537880e-03 0.0101003728 1000 -0.0214925444 2.184835e-02
#> 26 direct 9 -2.531896e-03 0.0116151545 1000 -0.0317689258 1.582699e-02
#> 27 indirect 9 -5.983648e-06 0.0079420617 1000 -0.0121672008 1.656888e-02
#> 28 total 10 -1.842634e-03 0.0077185775 1000 -0.0166303186 1.484984e-02
#> 29 direct 10 -1.871723e-03 0.0091825414 1000 -0.0261913900 1.164824e-02
#> 30 indirect 10 2.908908e-05 0.0076589800 1000 -0.0112799332 1.457964e-02
#> 31 total 11 -1.326708e-03 0.0060482375 1000 -0.0130163925 9.763080e-03
#> 32 direct 11 -1.379912e-03 0.0073113996 1000 -0.0209742211 8.697316e-03
#> 33 indirect 11 5.320397e-05 0.0071831649 1000 -0.0119257489 1.218782e-02
#> 34 total 12 -9.482262e-04 0.0048697265 1000 -0.0120190743 6.616435e-03
#> 35 direct 12 -1.015529e-03 0.0058656130 1000 -0.0185622374 6.287367e-03
#> 36 indirect 12 6.730311e-05 0.0065625860 1000 -0.0121010596 9.930568e-03
#> 37 total 13 -6.732162e-04 0.0040043510 1000 -0.0112494165 4.253606e-03
#> 38 direct 13 -7.465057e-04 0.0047409033 1000 -0.0149746772 4.720977e-03
#> 39 indirect 13 7.328946e-05 0.0058540199 1000 -0.0097182863 8.139744e-03
#> 40 total 14 -4.750334e-04 0.0033308783 1000 -0.0096598491 2.681612e-03
#> 41 direct 14 -5.483367e-04 0.0038590480 1000 -0.0129518986 3.613477e-03
#> 42 indirect 14 7.330331e-05 0.0051121862 1000 -0.0084125444 6.767638e-03
#> 43 total 15 -3.332586e-04 0.0027805955 1000 -0.0082849823 1.627667e-03
#> 44 direct 15 -4.025764e-04 0.0031619161 1000 -0.0113740940 2.622420e-03
#> 45 indirect 15 6.931784e-05 0.0043831867 1000 -0.0070908449 5.737716e-03
#> 46 total 16 -2.325098e-04 0.0023204147 1000 -0.0065053537 1.042185e-03
#> 47 direct 16 -2.954676e-04 0.0026063882 1000 -0.0095704872 1.780139e-03
#> 48 indirect 16 6.295775e-05 0.0037015657 1000 -0.0061333548 4.628230e-03
#> 49 total 17 -1.613554e-04 0.0019358900 1000 -0.0051754289 6.097761e-04
#> 50 direct 17 -2.168104e-04 0.0021603798 1000 -0.0080346023 1.305081e-03
#> 51 indirect 17 5.545493e-05 0.0030898810 1000 -0.0055750614 3.486811e-03
#> 52 total 18 -1.113930e-04 0.0016194598 1000 -0.0040363947 4.752519e-04
#> 53 direct 18 -1.590708e-04 0.0017998630 1000 -0.0066391400 9.968676e-04
#> 54 indirect 18 4.767783e-05 0.0025597976 1000 -0.0051043393 2.755816e-03
#> 55 total 19 -7.650453e-05 0.0013641051 1000 -0.0031682928 3.193049e-04
#> 56 direct 19 -1.166976e-04 0.0015066886 1000 -0.0057387799 7.555583e-04
#> 57 indirect 19 4.019303e-05 0.0021139697 1000 -0.0047168904 2.137021e-03
#> 58 total 20 -5.227212e-05 0.0011610097 1000 -0.0025175973 2.239409e-04
#> 59 direct 20 -8.560664e-05 0.0012670143 1000 -0.0048314773 5.772942e-04
#> 60 indirect 20 3.333451e-05 0.0017482372 1000 -0.0041673903 1.716129e-03
#> 61 total 21 -3.552885e-05 0.0009997955 1000 -0.0019656372 1.833369e-04
#> 62 direct 21 -6.279662e-05 0.0010701744 1000 -0.0040996609 4.053364e-04
#> 63 indirect 21 2.726777e-05 0.0014538717 1000 -0.0038000260 1.340027e-03
#> 64 total 22 -2.401981e-05 0.0008700292 1000 -0.0019962529 1.199377e-04
#> 65 direct 22 -4.606320e-05 0.0009078678 1000 -0.0034804750 3.050543e-04
#> 66 indirect 22 2.204340e-05 0.0012197162 1000 -0.0033755847 1.041580e-03
#> 67 total 23 -1.614944e-05 0.0007628003 1000 -0.0026103493 9.427607e-05
#> 68 direct 23 -3.378818e-05 0.0007735711 1000 -0.0029755872 2.253302e-04
#> 69 indirect 23 1.763874e-05 0.0010340746 1000 -0.0029495930 8.250661e-04
#> 70 total 24 -1.079535e-05 0.0006715757 1000 -0.0027129511 7.640367e-05
#> 71 direct 24 -2.478396e-05 0.0006621109 1000 -0.0025769067 1.530000e-04
#> 72 indirect 24 1.398861e-05 0.0008861637 1000 -0.0026813033 6.545636e-04
#> 73 total 25 -7.172379e-06 0.0005921888 1000 -0.0024425171 6.565484e-05
#> 74 direct 25 -1.817915e-05 0.0005693511 1000 -0.0021698257 1.103666e-04
#> 75 indirect 25 1.100677e-05 0.0007669758 1000 -0.0023691438 4.931002e-04
#> 76 total 26 -4.734244e-06 0.0005222834 1000 -0.0019290509 5.493112e-05
#> 77 direct 26 -1.333442e-05 0.0004919595 1000 -0.0018244634 8.230734e-05
#> 78 indirect 26 8.600180e-06 0.0006695342 1000 -0.0020418637 3.925042e-04
#> 79 total 27 -3.102834e-06 0.0004606122 1000 -0.0012339619 6.047519e-05
#> 80 direct 27 -9.780785e-06 0.0004272339 1000 -0.0015402728 6.367310e-05
#> 81 indirect 27 6.677951e-06 0.0005887002 1000 -0.0018239310 3.034498e-04
#> 82 total 28 -2.017802e-06 0.0004064440 1000 -0.0007937599 5.234278e-05
#> 83 direct 28 -7.174180e-06 0.0003729703 1000 -0.0013021625 4.884146e-05
#> 84 indirect 28 5.156378e-06 0.0005207669 1000 -0.0015164592 2.454520e-04
#> 85 total 29 -1.300805e-06 0.0003591783 1000 -0.0005039784 4.596516e-05
#> 86 direct 29 -5.262236e-06 0.0003273617 1000 -0.0011145562 3.902827e-05
#> 87 indirect 29 3.961431e-06 0.0004630382 1000 -0.0013069141 2.028427e-04
#> 88 total 30 -8.303077e-07 0.0003181590 1000 -0.0003754469 4.552413e-05
#> 89 direct 30 -3.859828e-06 0.0002889209 1000 -0.0009465892 3.026967e-05
#> 90 indirect 30 3.029520e-06 0.0004134999 1000 -0.0010970773 1.686548e-04