R Under development (unstable) (2026-07-09 r90225 ucrt) -- "Unsuffered Consequences" Copyright (C) 2026 The R Foundation for Statistical Computing Platform: x86_64-w64-mingw32/x64 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > # This file is part of the standard setup for testthat. > # It is recommended that you do not modify it. > # > # Where should you do additional test configuration? > # Learn more about the roles of various files in: > # * https://r-pkgs.org/testing-design.html#sec-tests-files-overview > # * https://testthat.r-lib.org/articles/special-files.html > > library(testthat) > library(bayesics) > > test_check("bayesics") --- Bayes factor in favor of the full vs. null model: 4.81e+55; =>Level of evidence: Decisive --- Summary of factor level means --- # A tibble: 10 x 5 Variable `Post Mean` Lower Upper `Prob Dir` 1 Mean : x1 : a -0.911 -1.11 -0.716 1.000 2 Mean : x1 : b -1.09 -1.31 -0.882 1 3 Mean : x1 : c -0.915 -1.09 -0.737 1 4 Mean : x1 : d 0.905 0.707 1.10 1.000 5 Mean : x1 : e 1.05 0.854 1.25 1 6 Var : x1 : a 0.992 0.750 1.31 NA 7 Var : x1 : b 1.17 0.884 1.54 NA 8 Var : x1 : c 0.829 0.627 1.09 NA 9 Var : x1 : d 1.03 0.776 1.36 NA 10 Var : x1 : e 1.01 0.764 1.33 NA --- Summary of pairwise differences --- # A tibble: 10 x 9 Comparison `Post Mean` Lower Upper `Prob Dir` `ROPE (0.1)` EPR 1 a-b 0.182 -0.103 0.465 0.894 0.270 0.549 2 a-c 0.00419 -0.261 0.268 0.506 0.525 0.501 3 a-d -1.82 -2.09 -1.54 1 0 0.101 4 a-e -1.96 -2.24 -1.69 1 0 0.0836 5 b-c -0.178 -0.452 0.0985 0.900 0.261 0.450 6 b-d -2.00 -2.28 -1.71 1 0 0.0893 7 b-e -2.14 -2.43 -1.85 1 0 0.0740 8 c-d -1.82 -2.09 -1.55 1 0 0.0913 9 c-e -1.97 -2.23 -1.70 1 0 0.0744 10 d-e -0.145 -0.428 0.135 0.844 0.340 0.460 # i 2 more variables: `EPR Lower` , `EPR Upper` *Note: EPR (Exceedence in Pairs Rate) for a Comparison of g-h = Pr(Y_(gi) > Y_(hi)|parameters) --- Bayes factor in favor of the full vs. null model: 5.79e+65; =>Level of evidence: Decisive --- Summary of factor level means --- # A tibble: 6 x 5 Variable `Post Mean` Lower Upper `Prob Dir` 1 Mean : x1 : a -0.911 -1.11 -0.717 1 2 Mean : x1 : b -1.09 -1.29 -0.899 1 3 Mean : x1 : c -0.915 -1.11 -0.721 1 4 Mean : x1 : d 0.905 0.711 1.10 1 5 Mean : x1 : e 1.05 0.856 1.24 1 6 Var 0.989 0.874 1.12 NA --- Summary of pairwise differences --- # A tibble: 10 x 9 Comparison `Post Mean` Lower Upper `Prob Dir` `ROPE (0.1)` EPR 1 a-b 0.182 -0.0926 0.453 0.903 0.254 0.551 2 a-c 0.00419 -0.272 0.279 0.512 0.527 0.501 3 a-d -1.82 -2.09 -1.54 1 0 0.0990 4 a-e -1.96 -2.24 -1.69 1 0 0.0825 5 b-c -0.178 -0.450 0.0973 0.898 0.262 0.450 6 b-d -2.00 -2.27 -1.73 1 0 0.0784 7 b-e -2.14 -2.42 -1.87 1 0 0.0647 8 c-d -1.82 -2.10 -1.55 1 0 0.0986 9 c-e -1.97 -2.24 -1.69 1 0 0.0822 10 d-e -0.145 -0.422 0.131 0.850 0.339 0.460 # i 2 more variables: `EPR Lower` , `EPR Upper` *Note: EPR (Exceedence in Pairs Rate) for a Comparison of g-h = Pr(Y_(gi) > Y_(hi)|parameters) --- Summary of factor level means --- # A tibble: 10 x 5 Variable `Post Mean` Lower Upper `Prob Dir` 1 Mean : x1 : a -0.926 -1.12 -0.731 1.000 2 Mean : x1 : b -1.11 -1.32 -0.898 1 3 Mean : x1 : c -0.930 -1.11 -0.753 1 4 Mean : x1 : d 0.900 0.700 1.10 1.000 5 Mean : x1 : e 1.05 0.847 1.24 1 6 Var : x1 : a 0.980 0.740 1.30 NA 7 Var : x1 : b 1.16 0.873 1.53 NA 8 Var : x1 : c 0.815 0.616 1.08 NA 9 Var : x1 : d 1.03 0.781 1.37 NA 10 Var : x1 : e 1.02 0.769 1.35 NA --- Summary of pairwise differences --- # A tibble: 10 x 9 Comparison `Post Mean` Lower Upper `Prob Dir` `ROPE (0.1)` EPR 1 a-b 0.183 -0.101 0.467 0.895 0.267 0.550 2 a-c 0.00421 -0.260 0.266 0.507 0.523 0.501 3 a-d -1.83 -2.11 -1.55 1 0 0.0997 4 a-e -1.97 -2.25 -1.69 1 0 0.0823 5 b-c -0.179 -0.454 0.0971 0.901 0.257 0.449 6 b-d -2.01 -2.30 -1.72 1 0 0.0880 7 b-e -2.15 -2.44 -1.86 1 0 0.0728 8 c-d -1.83 -2.10 -1.56 1 0 0.0899 9 c-e -1.98 -2.24 -1.71 1 0 0.0731 10 d-e -0.146 -0.431 0.137 0.843 0.340 0.460 # i 2 more variables: `EPR Lower` , `EPR Upper` *Note: EPR (Exceedence in Pairs Rate) for a Comparison of g-h = Pr(Y_(gi) > Y_(hi)|parameters) --- Summary of factor level means --- # A tibble: 6 x 5 Variable `Post Mean` Lower Upper `Prob Dir` 1 Mean : x1 : a -0.926 -1.12 -0.732 1 2 Mean : x1 : b -1.11 -1.30 -0.915 1 3 Mean : x1 : c -0.930 -1.12 -0.736 1 4 Mean : x1 : d 0.900 0.706 1.09 1 5 Mean : x1 : e 1.05 0.852 1.24 1 6 Var 0.977 0.863 1.11 NA --- Summary of pairwise differences --- # A tibble: 10 x 9 Comparison `Post Mean` Lower Upper `Prob Dir` `ROPE (0.1)` EPR 1 a-b 0.183 -0.0907 0.452 0.905 0.251 0.552 2 a-c 0.00421 -0.271 0.278 0.512 0.526 0.501 3 a-d -1.83 -2.10 -1.55 1 0 0.0964 4 a-e -1.97 -2.25 -1.70 1 0 0.0801 5 b-c -0.179 -0.449 0.0954 0.900 0.258 0.449 6 b-d -2.01 -2.28 -1.74 1 0 0.0761 7 b-e -2.15 -2.43 -1.88 1 0 0.0625 8 c-d -1.83 -2.11 -1.56 1 0 0.0961 9 c-e -1.98 -2.25 -1.70 1 0 0.0798 10 d-e -0.146 -0.421 0.130 0.852 0.335 0.459 # i 2 more variables: `EPR Lower` , `EPR Upper` *Note: EPR (Exceedence in Pairs Rate) for a Comparison of g-h = Pr(Y_(gi) > Y_(hi)|parameters) # A tibble: 15 x 7 Variable `Post Mean` Lower Upper `Prob Dir` ROPE `ROPE bounds` 1 (Intercept) -1.08 -1.24 -0.929 1 NA (NA,NA) 2 x1 1.02 0.926 1.12 1 0 (-0.089,0.08~ 3 x2 0.0130 -0.0373 0.105 0.283 0.947 (-0.088,0.08~ 4 x3b -0.0166 -0.239 0.156 0.231 0.996 (-0.351,0.35~ 5 x3c 0.0388 -0.0923 0.298 0.291 0.991 (-0.351,0.35~ 6 x3d 2.01 1.76 2.25 1 0 (-0.351,0.35~ 7 x3e 2.14 1.89 2.38 1 0 (-0.351,0.35~ 8 x4 -0.0000604 -0.0697 0.0696 0.177 0.973 (-0.083,0.08~ 9 x5 -0.00599 -0.0917 0.0568 0.221 0.966 (-0.088,0.08~ 10 x6 -0.000371 -0.0731 0.0703 0.176 0.974 (-0.087,0.08~ 11 x7 0.0132 -0.0350 0.107 0.282 0.948 (-0.089,0.08~ 12 x8 0.00545 -0.0581 0.0867 0.224 0.970 (-0.088,0.08~ 13 x9 0.0571 0 0.192 0.503 0.648 (-0.089,0.08~ 14 x10 -0.0121 -0.105 0.0387 0.281 0.948 (-0.083,0.08~ 15 Residual variance 1.05 0.936 1.28 NA NA (NA,NA) Cell sizes were too small for large sample approximation. Instead, setting uniform prior on Pr(exposure|outcome) and making exact finite sample inference. Cell sizes were too small for large sample approximation. Instead, setting uniform prior on Pr(exposure|outcome) and making exact finite sample inference. Cell sizes were too small for large sample approximation. Instead, setting uniform prior on Pr(exposure|outcome) and making exact finite sample inference. Cell sizes were too small for large sample approximation. Instead, setting uniform prior on Pr(exposure|outcome) and making exact finite sample inference. Prior shape parameters were not supplied. Jeffrey's prior will be used. Prior shape parameters were not supplied. Jeffrey's prior will be used. Prior shape parameters were not supplied. A uniform prior will be used. Prior shape parameters were not supplied. Jeffrey's prior will be used. Prior shape parameters were not supplied. Jeffrey's prior will be used. Prior shape parameters were not supplied. A uniform prior will be used. Prior shape parameters were not supplied. Jeffrey's prior will be used. Prior shape parameters were not supplied. Jeffrey's prior will be used. Prior shape parameters were not supplied. A uniform prior will be used. Prior shape parameters were not supplied. Jeffrey's prior will be used. Prior shape parameters were not supplied. Jeffrey's prior will be used. Prior shape parameters were not supplied. Beta(2,2) prior will be used. Prior shape parameters were not supplied. Beta(2,2) prior will be used. Prior shape parameters were not supplied. Beta(2,2) prior will be used. Prior shape parameters were not supplied. Beta(2,2) prior will be used. Prior shape parameters were not supplied. Beta(2,2) prior will be used. Prior shape parameters were not supplied. Beta(2,2) prior will be used. Prior shape parameters were not supplied. A uniform prior will be used. Prior shape parameters were not supplied. Beta(2,3.9) prior will be used. Prior shape parameters were not supplied. Beta(3.9,2) prior will be used. The g hyperparameter in Zellner's g prior is not specified. It will be set automatically to n. Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. ---------- Values given in terms of odds ratios ---------- # A tibble: 6 x 7 Variable `Post Mean` Lower Upper `Prob Dir` ROPE `ROPE bounds` 1 x1 2.22 1.26 3.93 0.997 0.00180 (0.972,1.029) 2 x2 0.838 0.473 1.49 0.727 0.0694 (0.97,1.031) 3 x3b 0.429 0.0508 3.62 0.781 0.0638 (0.889,1.125) 4 x3c 0.911 0.126 6.58 0.537 0.0925 (0.889,1.125) 5 x3d 8.49 1.54 46.8 0.993 0.00535 (0.889,1.125) 6 x3e 21.5 3.91 118. 1.000 0.000223 (0.889,1.125) Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. The g hyperparameter in Zellner's g prior is not specified. It will be set automatically to n. Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. ---------- Values given in terms of odds ratios ---------- # A tibble: 6 x 7 Variable `Post Mean` Lower Upper `Prob Dir` ROPE `ROPE bounds` 1 x1 2.12 1.48 3.06 1.000 0.0000502 (0.972,1.029) 2 x2 0.848 0.588 1.22 0.822 0.0912 (0.97,1.031) 3 x3b 0.481 0.124 1.81 0.865 0.0752 (0.889,1.125) 4 x3c 0.910 0.265 3.13 0.561 0.155 (0.889,1.125) 5 x3d 7.32 2.58 21.8 1.000 0.000261 (0.889,1.125) 6 x3e 17.6 6.23 53.5 1.000 0.000000637 (0.889,1.125) Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. The g hyperparameter in Zellner's g prior is not specified. It will be set automatically to n. Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. The g hyperparameter in Zellner's g prior is not specified. It will be set automatically to n. Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. The g hyperparameter in Zellner's g prior is not specified. It will be set automatically to n. Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. ---------- Values given in terms of odds ratios ---------- # A tibble: 6 x 7 Variable `Post Mean` Lower Upper `Prob Dir` ROPE `ROPE bounds` 1 x1 2.09 1.21 3.62 0.996 0.00258 (0.972,1.029) 2 x2 0.850 0.490 1.47 0.719 0.0732 (0.97,1.031) 3 x3b 0.492 0.0630 3.84 0.751 0.0712 (0.889,1.125) 4 x3c 0.908 0.137 6.03 0.540 0.0966 (0.889,1.125) 5 x3d 7.08 1.38 36.3 0.991 0.00726 (0.889,1.125) 6 x3e 16.8 3.28 86.0 1.000 0.000381 (0.889,1.125) Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. The g hyperparameter in Zellner's g prior is not specified. It will be set automatically to n. The g hyperparameter in Zellner's g prior is not specified. It will be set automatically to n. The g hyperparameter in Zellner's g prior is not specified. It will be set automatically to n. The g hyperparameter in Zellner's g prior is not specified. It will be set automatically to n. The g hyperparameter in Zellner's g prior is not specified. It will be set automatically to n. The g hyperparameter in Zellner's g prior is not specified. It will be set automatically to n. The g hyperparameter in Zellner's g prior is not specified. It will be set automatically to n. ---------- Values given in terms of rate ratios ---------- # A tibble: 6 x 7 Variable `Post Mean` Lower Upper `Prob Dir` ROPE `ROPE bounds` 1 x1 2.53 2.04 3.13 1 1.59e-16 (0.972,1.029) 2 x2 0.989 0.751 1.30 0.531 1.71e- 1 (0.97,1.031) 3 x3b 1.69 0.375 7.63 0.753 9.66e- 2 (0.889,1.125) 4 x3c 1.35 0.300 6.11 0.653 1.13e- 1 (0.889,1.125) 5 x3d 8.89 2.43 32.6 1.000 6.52e- 4 (0.889,1.125) 6 x3e 10.0 2.73 36.8 1.000 3.57e- 4 (0.889,1.125) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) The g hyperparameter in Zellner's g prior is not specified. It will be set automatically to n. ---------- Values given in terms of rate ratios ---------- # A tibble: 6 x 7 Variable `Post Mean` Lower Upper `Prob Dir` ROPE `ROPE bounds` 1 x1 2.51 2.19 2.89 1 8.39e-29 (0.972,1.029) 2 x2 0.988 0.827 1.18 0.557 2.70e- 1 (0.97,1.031) 3 x3b 1.61 0.629 4.25 0.840 1.16e- 1 (0.889,1.125) 4 x3c 1.25 0.493 3.35 0.680 1.79e- 1 (0.889,1.125) 5 x3d 7.82 3.54 18.4 1.000 4.75e- 7 (0.889,1.125) 6 x3e 8.88 4.07 21.0 1.000 6.63e- 8 (0.889,1.125) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) The g hyperparameter in Zellner's g prior is not specified. It will be set automatically to n. ---------- Values given in terms of rate ratios ---------- # A tibble: 6 x 7 Variable `Post Mean` Lower Upper `Prob Dir` ROPE `ROPE bounds` 1 x1 2.51 2.02 3.11 1 2.55e-16 (0.972,1.029) 2 x2 0.987 0.750 1.30 0.536 1.71e- 1 (0.97,1.031) 3 x3b 1.58 0.360 6.95 0.728 1.03e- 1 (0.889,1.125) 4 x3c 1.23 0.281 5.40 0.609 1.19e- 1 (0.889,1.125) 5 x3d 7.58 2.13 27.0 0.999 1.15e- 3 (0.889,1.125) 6 x3e 8.62 2.42 30.8 1.000 6.15e- 4 (0.889,1.125) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) The g hyperparameter in Zellner's g prior is not specified. It will be set automatically to n. ---------- Values given in terms of rate ratios ---------- # A tibble: 7 x 7 Variable `Post Mean` Lower Upper `Prob Dir` ROPE `ROPE bounds` 1 x1 3.15 1.86 5.33 1.000 0.0000102 (0.972,1.029) 2 x2 0.737 0.428 1.27 0.864 0.0481 (0.97,1.031) 3 x3b 0.536 0.0796 3.61 0.739 0.0785 (0.889,1.125) 4 x3c 0.672 0.105 4.31 0.662 0.0906 (0.889,1.125) 5 x3d 3.95 0.770 20.3 0.950 0.0292 (0.889,1.125) 6 x3e 2.38 0.470 12.1 0.853 0.0655 (0.889,1.125) 7 phi 0.614 0.390 0.966 0.983 NA (NA,NA) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) The g hyperparameter in Zellner's g prior is not specified. It will be set automatically to n. ---------- Values given in terms of rate ratios ---------- # A tibble: 7 x 7 Variable `Post Mean` Lower Upper `Prob Dir` ROPE `ROPE bounds` 1 x1 3.02 2.19 4.20 1.000 1.16e-9 (0.972,1.029) 2 x2 0.721 0.521 1.02 0.972 2.03e-2 (0.97,1.031) 3 x3b 0.540 0.162 1.75 0.857 9.23e-2 (0.889,1.125) 4 x3c 0.694 0.210 2.20 0.754 1.39e-1 (0.889,1.125) 5 x3d 3.58 1.33 9.71 0.995 6.26e-3 (0.889,1.125) 6 x3e 2.25 0.846 6.27 0.952 4.67e-2 (0.889,1.125) 7 phi 0.730 0.434 1.20 0.898 NA (NA,NA) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) The g hyperparameter in Zellner's g prior is not specified. It will be set automatically to n. ---------- Values given in terms of rate ratios ---------- # A tibble: 7 x 7 Variable `Post Mean` Lower Upper `Prob Dir` ROPE `ROPE bounds` 1 x1 3.01 1.83 4.96 1.000 0.00000833 (0.972,1.029) 2 x2 0.722 0.433 1.20 0.894 0.0429 (0.97,1.031) 3 x3b 0.540 0.0902 3.23 0.750 0.0818 (0.889,1.125) 4 x3c 0.686 0.119 3.97 0.663 0.0958 (0.889,1.125) 5 x3d 3.54 0.770 16.3 0.948 0.0325 (0.889,1.125) 6 x3e 2.22 0.483 10.2 0.847 0.0714 (0.889,1.125) 7 phi 0.756 0.368 1.55 0.778 NA (NA,NA) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The hyperparameters for the residual variance were not provided. Instead, the prior will put 50% prior probability that R^2 is between 0.1^2 and 0.9^2. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The hyperparameters for the residual variance were not provided. Instead, the prior will put 50% prior probability that R^2 is between 0.1^2 and 0.9^2. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) The g hyperparameter in Zellner's g prior is not specified. It will be set automatically to n. The hyperparameters for the residual variance were not provided. Instead, the prior will put 50% prior probability that R^2 is between 0.1^2 and 0.9^2. The g hyperparameter in Zellner's g prior is not specified. It will be set automatically to n. The hyperparameters for the residual variance were not provided. Instead, the prior will put 50% prior probability that R^2 is between 0.1^2 and 0.9^2. ---------- Test for heteroscedasticity in 1-way ANOVA models. Bayes factor in favor of homoscedasticity = 592366584.582148 Level of evidence: Decisive in favor of homoscedasticity ---------- ---------- Test for heteroscedasticity in 1-way ANOVA models. Bayes factor in favor of homoscedasticity = 0.028467538606378 Level of evidence: Strong in favor of heteroscedasticity ---------- The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The hyperparameters for the residual variance were not provided. Instead, the prior will put 50% prior probability that R^2 is between 0.1^2 and 0.9^2. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The hyperparameters for the residual variance were not provided. Instead, the prior will put 50% prior probability that R^2 is between 0.1^2 and 0.9^2. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) # A tibble: 8 x 7 Variable `Post Mean` Lower Upper `Prob Dir` ROPE `ROPE bounds` 1 (Intercept) -0.996 -1.20 -0.790 1 NA (NA,NA) 2 x1 0.943 0.849 1.04 1 1.85e-56 (-0.085,0.08~ 3 x2 0.0897 -0.00309 0.182 0.971 4.56e- 1 (-0.084,0.08~ 4 x3b 0.0428 -0.249 0.335 0.613 9.71e- 1 (-0.337,0.33~ 5 x3c -0.0155 -0.307 0.276 0.542 9.76e- 1 (-0.337,0.33~ 6 x3d 1.90 1.61 2.19 1 8.75e-24 (-0.337,0.33~ 7 x3e 1.94 1.65 2.24 1 5.61e-25 (-0.337,0.33~ 8 Residual varia~ 1.10 0.978 1.25 NA NA (NA,NA) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The hyperparameters for the residual variance were not provided. Instead, the prior will put 50% prior probability that R^2 is between 0.1^2 and 0.9^2. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) By default, the fraction of data "used" is max(ncol(X) + 1,log(n)) / n. ---------- The fractional Bayes factor equaled 1.76e+60. Interpretation: Decisive (in favor of the first model) ---------- The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The hyperparameters for the residual variance were not provided. Instead, the prior will put 50% prior probability that R^2 is between 0.1^2 and 0.9^2. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The hyperparameters for the residual variance were not provided. Instead, the prior will put 50% prior probability that R^2 is between 0.1^2 and 0.9^2. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The hyperparameters for the residual variance were not provided. Instead, the prior will put 50% prior probability that R^2 is between 0.1^2 and 0.9^2. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The hyperparameters for the residual variance were not provided. Instead, the prior will put 50% prior probability that R^2 is between 0.1^2 and 0.9^2. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The hyperparameters for the residual variance were not provided. Instead, the prior will put 50% prior probability that R^2 is between 0.1^2 and 0.9^2. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) # A tibble: 8 x 7 Variable `Post Mean` Lower Upper `Prob Dir` ROPE `ROPE bounds` 1 (Intercept) -0.996 -1.20 -0.790 1 NA (NA,NA) 2 x1 0.943 0.849 1.04 1 1.85e-56 (-0.085,0.08~ 3 x2 0.0897 -0.00309 0.182 0.971 4.56e- 1 (-0.084,0.08~ 4 x3b 0.0428 -0.249 0.335 0.613 9.71e- 1 (-0.337,0.33~ 5 x3c -0.0155 -0.307 0.276 0.542 9.76e- 1 (-0.337,0.33~ 6 x3d 1.90 1.61 2.19 1 8.75e-24 (-0.337,0.33~ 7 x3e 1.94 1.65 2.24 1 5.61e-25 (-0.337,0.33~ 8 Residual varia~ 1.10 0.978 1.25 NA NA (NA,NA) Missing other covariate values in 'exemplar_covariates.' Using medoid observation instead. The g hyperparameter in Zellner's g prior is not specified. It will be set automatically to n. The hyperparameters for the residual variance were not provided. Instead, the prior will put 50% prior probability that R^2 is between 0.1^2 and 0.9^2. The g hyperparameter in Zellner's g prior is not specified. It will be set automatically to n. The hyperparameters for the residual variance were not provided. Instead, the prior will put 50% prior probability that R^2 is between 0.1^2 and 0.9^2. # A tibble: 8 x 7 Variable `Post Mean` Lower Upper `Prob Dir` ROPE `ROPE bounds` 1 (Intercept) -0.996 -1.20 -0.790 1 NA (NA,NA) 2 x1 0.941 0.848 1.04 1 2.84e-56 (-0.085,0.08~ 3 x2 0.0895 -0.00332 0.182 0.971 4.58e- 1 (-0.084,0.08~ 4 x3b 0.0450 -0.247 0.337 0.619 9.70e- 1 (-0.337,0.33~ 5 x3c -0.0132 -0.305 0.279 0.535 9.76e- 1 (-0.337,0.33~ 6 x3d 1.90 1.61 2.19 1 1.03e-23 (-0.337,0.33~ 7 x3e 1.94 1.65 2.23 1 6.66e-25 (-0.337,0.33~ 8 Residual varia~ 1.11 0.980 1.26 NA NA (NA,NA) The g hyperparameter in Zellner's g prior is not specified. It will be set automatically to n. The hyperparameters for the residual variance were not provided. Instead, the prior will put 50% prior probability that R^2 is between 0.1^2 and 0.9^2. By default, the fraction of data "used" is max(ncol(X) + 1,log(n)) / n. ---------- The fractional Bayes factor equaled 1.76e+60. Interpretation: Decisive (in favor of the first model) ---------- The g hyperparameter in Zellner's g prior is not specified. It will be set automatically to n. The hyperparameters for the residual variance were not provided. Instead, the prior will put 50% prior probability that R^2 is between 0.1^2 and 0.9^2. The g hyperparameter in Zellner's g prior is not specified. It will be set automatically to n. The hyperparameters for the residual variance were not provided. Instead, the prior will put 50% prior probability that R^2 is between 0.1^2 and 0.9^2. The g hyperparameter in Zellner's g prior is not specified. It will be set automatically to n. The hyperparameters for the residual variance were not provided. Instead, the prior will put 50% prior probability that R^2 is between 0.1^2 and 0.9^2. The g hyperparameter in Zellner's g prior is not specified. It will be set automatically to n. The hyperparameters for the residual variance were not provided. Instead, the prior will put 50% prior probability that R^2 is between 0.1^2 and 0.9^2. The g hyperparameter in Zellner's g prior is not specified. It will be set automatically to n. The hyperparameters for the residual variance were not provided. Instead, the prior will put 50% prior probability that R^2 is between 0.1^2 and 0.9^2. The hyperparameters for the residual variance were not provided. Instead, the prior will put 50% prior probability that R^2 is between 0.1^2 and 0.9^2. # A tibble: 8 x 7 Variable `Post Mean` Lower Upper `Prob Dir` ROPE `ROPE bounds` 1 (Intercept) -0.998 -1.21 -0.790 1 NA (NA,NA) 2 x1 0.943 0.849 1.04 1 1.17e-55 (-0.085,0.08~ 3 x2 0.0896 -0.00380 0.183 0.970 4.57e- 1 (-0.084,0.08~ 4 x3b 0.0451 -0.249 0.340 0.618 9.69e- 1 (-0.337,0.33~ 5 x3c -0.0132 -0.307 0.281 0.535 9.75e- 1 (-0.337,0.33~ 6 x3d 1.90 1.61 2.20 1 1.69e-23 (-0.337,0.33~ 7 x3e 1.95 1.65 2.24 1 1.13e-24 (-0.337,0.33~ 8 Residual varia~ 1.12 0.991 1.27 NA NA (NA,NA) By default, the fraction of data "used" is max(ncol(X) + 1,log(n)) / n. ---------- The fractional Bayes factor equaled 1.76e+60. Interpretation: Decisive (in favor of the first model) ---------- # A tibble: 8 x 7 Variable `Post Mean` Lower Upper `Prob Dir` ROPE `ROPE bounds` 1 (Intercept) -0.998 -1.21 -0.790 1 NA (NA,NA) 2 x1 0.943 0.849 1.04 1 1.17e-55 (-0.085,0.08~ 3 x2 0.0896 -0.00380 0.183 0.970 4.57e- 1 (-0.084,0.08~ 4 x3b 0.0451 -0.249 0.340 0.618 9.69e- 1 (-0.337,0.33~ 5 x3c -0.0132 -0.307 0.281 0.535 9.75e- 1 (-0.337,0.33~ 6 x3d 1.90 1.61 2.20 1 1.69e-23 (-0.337,0.33~ 7 x3e 1.95 1.65 2.24 1 1.13e-24 (-0.337,0.33~ 8 Residual varia~ 1.12 0.991 1.27 NA NA (NA,NA) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The hyperparameters for the residual variance were not provided. Instead, the prior will put 50% prior probability that R^2 is between 0.1^2 and 0.9^2. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The hyperparameters for the residual variance were not provided. Instead, the prior will put 50% prior probability that R^2 is between 0.1^2 and 0.9^2. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The hyperparameters for the residual variance were not provided. Instead, the prior will put 50% prior probability that R^2 is between 0.1^2 and 0.9^2. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) The g hyperparameter in Zellner's g prior is not specified. It will be set automatically to n. The hyperparameters for the residual variance were not provided. Instead, the prior will put 50% prior probability that R^2 is between 0.1^2 and 0.9^2. The g hyperparameter in Zellner's g prior is not specified. It will be set automatically to n. The hyperparameters for the residual variance were not provided. Instead, the prior will put 50% prior probability that R^2 is between 0.1^2 and 0.9^2. The g hyperparameter in Zellner's g prior is not specified. It will be set automatically to n. The hyperparameters for the residual variance were not provided. Instead, the prior will put 50% prior probability that R^2 is between 0.1^2 and 0.9^2. The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The hyperparameters for the residual variance were not provided. Instead, the prior will put 50% prior probability that R^2 is between 0.1^2 and 0.9^2. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) # A tibble: 10 x 7 Variable `Post Mean` Lower Upper `Prob Dir` ROPE `ROPE bounds` 1 (Intercept) -1.10 -1.31 -0.886 1 NA (NA,NA) 2 x1 0.947 0.855 1.04 1 3.79e-57 (-0.091,0.09~ 3 x2 0.0814 -0.0108 0.174 0.958 5.76e- 1 (-0.09,0.09) 4 x3b 0.0202 -0.269 0.310 0.555 9.85e- 1 (-0.361,0.36~ 5 x3c -0.000627 -0.290 0.288 0.502 9.86e- 1 (-0.361,0.36~ 6 x3d -2.10 -2.39 -1.81 1 0 (-0.361,0.36~ 7 x3e -2.06 -2.35 -1.77 1 0 (-0.361,0.36~ 8 I(x1^2) 0.361 0.294 0.428 1 3.27e-17 (-0.065,0.06~ 9 x1:x2 -0.0557 -0.150 0.0382 0.878 7.66e- 1 (-0.091,0.09~ 10 Residual vari~ 1.08 0.957 1.23 NA NA (NA,NA) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The hyperparameters for the residual variance were not provided. Instead, the prior will put 50% prior probability that R^2 is between 0.1^2 and 0.9^2. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) # A tibble: 8 x 7 Variable `Post Mean` Lower Upper `Prob Dir` ROPE `ROPE bounds` 1 (Intercept) 5.41 4.32 6.50 1 NA (NA,NA) 2 ns(x1, df = 5)1 -6.26 -7.31 -5.21 1 0 (-0.407,0.40~ 3 ns(x1, df = 5)2 -5.28 -6.50 -4.07 1 8.88e-15 (-0.398,0.39~ 4 ns(x1, df = 5)3 -1.31 -2.18 -0.433 0.998 3.25e- 2 (-0.484,0.48~ 5 ns(x1, df = 5)4 -3.53 -6.03 -1.04 0.997 1.41e- 2 (-0.756,0.75~ 6 ns(x1, df = 5)5 8.22 7.01 9.42 1 8.92e-30 (-0.84,0.84) 7 x2 0.119 -0.00581 0.244 0.969 4.69e- 1 (-0.114,0.11~ 8 Residual varia~ 1.97 1.74 2.23 NA NA (NA,NA) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The hyperparameters for the residual variance were not provided. Instead, the prior will put 50% prior probability that R^2 is between 0.1^2 and 0.9^2. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The hyperparameters for the residual variance were not provided. Instead, the prior will put 50% prior probability that R^2 is between 0.1^2 and 0.9^2. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) # A tibble: 10 x 7 Variable `Post Mean` Lower Upper `Prob Dir` ROPE `ROPE bounds` 1 (Intercept) -1.09 -1.31 -0.881 1 NA (NA,NA) 2 x1 0.947 0.854 1.04 1 1.27e- 54 (-0.115,0.11~ 3 x2 0.0815 -0.0109 0.174 0.958 7.54e- 1 (-0.114,0.11~ 4 x3b 0.0158 -0.274 0.306 0.543 9.98e- 1 (-0.455,0.45~ 5 x3c -0.00524 -0.295 0.284 0.514 9.98e- 1 (-0.455,0.45~ 6 x3d 1.89 1.60 2.18 1 5.41e- 21 (-0.455,0.45~ 7 x3e 1.94 1.65 2.23 1 3.47e- 22 (-0.455,0.45~ 8 I(x1^2) 1.11 1.04 1.18 1 1.27e-114 (-0.082,0.08~ 9 x1:x2 -0.0556 -0.150 0.0384 0.877 8.88e- 1 (-0.114,0.11~ 10 Residual var~ 1.09 0.961 1.23 NA NA (NA,NA) The mu hyperparameter in the normal prior is not specified. It will be set automatically to 0. The hyperparameters for the residual variance were not provided. Instead, the prior will put 50% prior probability that R^2 is between 0.1^2 and 0.9^2. The V hyperparameter in the normal prior is not specified. It will be set automatically to 4/25Diag(s^2_{X_j}) # A tibble: 8 x 7 Variable `Post Mean` Lower Upper `Prob Dir` ROPE `ROPE bounds` 1 (Intercept) -0.996 -1.20 -0.790 1 NA (NA,NA) 2 x1 0.943 0.849 1.04 1 1.85e-56 (-0.085,0.08~ 3 x2 0.0897 -0.00309 0.182 0.971 4.56e- 1 (-0.084,0.08~ 4 x3b 0.0428 -0.249 0.335 0.613 9.71e- 1 (-0.337,0.33~ 5 x3c -0.0155 -0.307 0.276 0.542 9.76e- 1 (-0.337,0.33~ 6 x3d 1.90 1.61 2.19 1 8.75e-24 (-0.337,0.33~ 7 x3e 1.94 1.65 2.24 1 5.61e-25 (-0.337,0.33~ 8 Residual varia~ 1.10 0.978 1.25 NA NA (NA,NA) Missing other covariate values in 'exemplar_covariates.' Using medoid observation instead. The g hyperparameter in Zellner's g prior is not specified. It will be set automatically to n. The hyperparameters for the residual variance were not provided. Instead, the prior will put 50% prior probability that R^2 is between 0.1^2 and 0.9^2. The g hyperparameter in Zellner's g prior is not specified. It will be set automatically to n. Finished with 500 preliminary posterior draws. # A tibble: 8 x 5 Estimand Estimate Lower Upper `Prob Dir` 1 ACME (Control) 1.08 0.486 1.80 1 2 ACME (Treatment) 1.90 0.952 3.30 1 3 ADE (Control) 0.865 0.642 1.15 1 4 ADE (Treatment) 1.68 1.28 2.23 1 5 Total Effect 2.76 1.91 4.05 1 6 ACME (Average) 1.49 0.741 2.53 1 7 ADE (Average) 1.27 1.02 1.57 1 8 Prop. Mediated (Average) 0.529 0.373 0.633 NA The g hyperparameter in Zellner's g prior is not specified. It will be set automatically to n. The hyperparameters for the residual variance were not provided. Instead, the prior will put 50% prior probability that R^2 is between 0.1^2 and 0.9^2. The g hyperparameter in Zellner's g prior is not specified. It will be set automatically to n. Finished with 500 preliminary posterior draws. control_value missing; set to be the 1st quintile of tr treat_value missing; set to be the 4th quintile of tr Finished with 500 preliminary posterior draws. Finished with 500 preliminary posterior draws. # A tibble: 8 x 5 Estimand Estimate Lower Upper `Prob Dir` 1 ACME (Control) 0.638 0.367 1.01 1 2 ACME (Treatment) 2.30 1.44 3.60 1 3 ADE (Control) 0.927 0.687 1.26 1 4 ADE (Treatment) 2.59 2.02 3.52 1 5 Total Effect 3.23 2.45 4.50 1 6 ACME (Average) 1.47 0.933 2.30 1 7 ADE (Average) 1.76 1.43 2.24 1 8 Prop. Mediated (Average) 0.450 0.359 0.527 NA The g hyperparameter in Zellner's g prior is not specified. It will be set automatically to n. The hyperparameters for the residual variance were not provided. Instead, the prior will put 50% prior probability that R^2 is between 0.1^2 and 0.9^2. The g hyperparameter in Zellner's g prior is not specified. It will be set automatically to n. The hyperparameters for the residual variance were not provided. Instead, the prior will put 50% prior probability that R^2 is between 0.1^2 and 0.9^2. # A tibble: 4 x 5 Estimand Estimate Lower Upper `Prob Dir` 1 ACME 2.04 1.48 2.55 1 2 ADE 2.27 1.89 2.63 1 3 Total Effect 4.31 3.76 4.88 1 4 Prop. Mediated 0.473 0.381 0.555 NA Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. Finished with 50 preliminary Bayesian bootstraps. Continuing on with 1 more Bayesian bootstraps. ---------- Values given in terms of odds ratios ---------- # A tibble: 6 x 7 Variable `Post Mean` Lower Upper `Prob Dir` ROPE `ROPE bounds` 1 x1 2.56 1.19 6.45 0.980 0 (0.972,1.029) 2 x2 0.842 0.502 1.30 0.686 0.0588 (0.97,1.031) 3 x3b 0.218 0.00459 6.24 0.843 0.0392 (0.889,1.125) 4 x3c 0.782 0.0543 5.87 0.569 0.0392 (0.889,1.125) 5 x3d 12.7 2.64 92.8 1 0 (0.889,1.125) 6 x3e 30.0 4.60 430. 1 0 (0.889,1.125) # A tibble: 7 x 7 Variable `Post Mean` Lower Upper `Prob Dir` ROPE `ROPE bounds` 1 (Intercept) -2.54 -4.62 -1.27 1 NA (NA,NA) 2 x1 0.940 0.174 1.86 0.980 0 (-0.029,0.029) 3 x2 -0.172 -0.690 0.264 0.686 0.0588 (-0.031,0.031) 4 x3b -1.52 -5.38 1.83 0.843 0.0392 (-0.118,0.118) 5 x3c -0.246 -2.91 1.77 0.569 0.0392 (-0.118,0.118) 6 x3d 2.54 0.972 4.53 1 0 (-0.118,0.118) 7 x3e 3.40 1.53 6.06 1 0 (-0.118,0.118) Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. Missing other covariate values in 'exemplar_covariates.' Using medoid observation instead. Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. Missing other covariate values in 'exemplar_covariates.' Using medoid observation instead. Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. Missing other covariate values in 'exemplar_covariates.' Using medoid observation instead. Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. Assuming all observations correspond to Bernoulli, i.e., Binomial with one trial. Finished with 50 preliminary Bayesian bootstraps. Continuing on with 1 more Bayesian bootstraps. Prior shape parameters were not supplied. Jeffrey's prior will be used. Prior shape parameters were not supplied. Jeffrey's prior will be used. Prior shape parameters were not supplied. Jeffrey's prior will be used. Prior shape parameters were not supplied. A flat Gamma(0.001,0.001) prior will be used. Prior shape parameters were not supplied. Jeffrey's prior will be used. Prior shape parameters were not supplied. Jeffrey's prior will be used. Prior shape parameters were not supplied. Jeffrey's prior will be used. Prior shape parameters were not supplied. A flat Gamma(0.001,0.001) prior will be used. Prior shape parameters were not supplied. Jeffrey's prior will be used. Prior shape parameters were not supplied. Jeffrey's prior will be used. Prior shape parameters were not supplied. Jeffrey's prior will be used. Prior shape parameters were not supplied. Jeffrey's prior will be used. Prior shape parameters were not supplied. Jeffrey's prior will be used. Prior shape parameters were not supplied. Jeffrey's prior will be used. Prior shape parameters were not supplied. Jeffrey's prior will be used. Prior shape parameters were not supplied. Jeffrey's prior will be used. Prior shape parameters were not supplied. A uniform prior will be used. Prior shape parameters were not supplied. Jeffrey's prior will be used. Prior shape parameters were not supplied. Jeffrey's prior will be used. Prior shape parameters were not supplied. Jeffrey's prior will be used. Prior shape parameters were not supplied. Jeffrey's prior will be used. Prior shape parameters were not supplied. Beta(2,2) prior will be used. Prior shape parameters were not supplied. Beta(2,2) prior will be used. Prior shape parameters were not supplied. A uniform prior will be used. Prior shape parameters were not supplied. Beta(2,2) prior will be used. Prior shape parameters were not supplied. Beta(2,2) prior will be used. Prior shape parameters were not supplied. Beta(2,2) prior will be used. Prior shape parameters were not supplied. A uniform prior will be used. Prior shape parameters were not supplied. Beta(2,2) prior will be used. Prior shape parameters were not supplied. Beta(2,2) prior will be used. [ FAIL 0 | WARN 5 | SKIP 0 | PASS 1335 ] [ FAIL 0 | WARN 5 | SKIP 0 | PASS 1335 ] > > proc.time() user system elapsed 1747.93 250.76 2029.20