tests/testthat/_snaps/predict_counterfactual.md

In RobinCar2: ROBust INference for Covariate Adjustment in Randomized Clinical Trials

predict_counterfactual works for guassian

Code
  predict_counterfactual(fit_glm, treatment ~ 1)
Output
  Model        :  y ~ treatment * s1 + covar 
  Randomization:  treatment ~ 1  ( Simple )
  Variance Type:  vcovG 
  Marginal Mean: 
       Estimate  Std.Err    2.5 % 97.5 %
  pbo  0.200321 0.067690 0.067651 0.3330
  trt1 0.763971 0.075929 0.615152 0.9128
  trt2 0.971250 0.076543 0.821228 1.1213

predict_counterfactual works for guassian with lm

Code
  predict_counterfactual(fit_lm, treatment ~ 1, data = glm_data)
Output
  Model        :  y ~ treatment * s1 + covar 
  Randomization:  treatment ~ 1  ( Simple )
  Variance Type:  vcovG 
  Marginal Mean: 
       Estimate  Std.Err    2.5 % 97.5 %
  pbo  0.200321 0.067690 0.067651 0.3330
  trt1 0.763971 0.075929 0.615152 0.9128
  trt2 0.971250 0.076543 0.821228 1.1213

predict_counterfactual works for binomial

Code
  predict_counterfactual(fit_binom, treatment ~ 1)
Output
  Model        :  y_b ~ treatment * s1 + covar 
  Randomization:  treatment ~ 1  ( Simple )
  Variance Type:  vcovG 
  Marginal Mean: 
       Estimate  Std.Err    2.5 % 97.5 %
  pbo  0.356097 0.033599 0.290243 0.4219
  trt1 0.580696 0.034418 0.513238 0.6482
  trt2 0.621386 0.034019 0.554711 0.6881

predict_counterfactual works if contrast are non-standard

Code
  pc
Output
  Model        :  y_b ~ treatment * s1 
  Randomization:  treatment ~ 1  ( Simple )
  Variance Type:  vcovG 
  Marginal Mean: 
       Estimate  Std.Err    2.5 % 97.5 %
  pbo  0.366007 0.033864 0.299634 0.4324
  trt1 0.580984 0.035027 0.512332 0.6496
  trt2 0.610154 0.034472 0.542589 0.6777

predict_counterfactual works with provided vcov function

Code
  predict_counterfactual(fit_lm, treatment ~ 1, data = glm_data, vcov = vcov_dummy)
Output
  Model        :  y ~ treatment * s1 + covar 
  Randomization:  treatment ~ 1  ( Simple )
  Variance Type:  vcov_dummy  
  Marginal Mean: 
       Estimate  Std.Err    2.5 % 97.5 %
  pbo   0.20032  0.31623 -0.41947 0.8201
  trt1  0.76397  0.31623  0.14418 1.3838
  trt2  0.97125  0.31623  0.35145 1.5910

Code
  predict_counterfactual(fit_lm, treatment ~ 1, data = glm_data, vcov = NULL)
Output
  Model        :  y ~ treatment * s1 + covar 
  Randomization:  treatment ~ 1  ( Simple )
  Variance Type:  
  Marginal Mean: 
       Estimate Std.Err 2.5 % 97.5 %
  pbo   0.20032      NA    NA     NA
  trt1  0.76397      NA    NA     NA
  trt2  0.97125      NA    NA     NA

confint method for prediction_cf works as expected

Code
  confint(predict_counterfactual(fit_binom, treatment = treatment ~ s1,
  eff_measure = h_diff))
Output
        Estimate     2.5 %    97.5 %
  pbo  0.3560965 0.2902435 0.4219496
  trt1 0.5806957 0.5132377 0.6481538
  trt2 0.6213865 0.5547112 0.6880618

RobinCar2 documentation built on Sept. 9, 2025, 5:28 p.m.