Rename Covariates in draws Object — rename_draws

Rename Covariates in draws Object

Usage

rename_draws_covariates(draws, analysis)

Arguments

draws: draws created from sampled analysis object. See example.
analysis: Analysis as created by create_analysis_obj().

Value

A draws[posterior::draws] object with covariate names.

Examples

if (check_cmdstan()) {
  analysis_object <- create_analysis_obj(
    data_matrix = example_matrix,
    covariates = add_covariates(
      covariates = c("cov1", "cov2"),
      priors = prior_normal(0, 1000)
    ),
    outcome = outcome_surv_exponential(
      "time",
      "cnsr",
      baseline_prior = prior_normal(0, 1000)
    ),
    borrowing = borrowing_hierarchical_commensurate(
      "ext",
      prior_exponential(.001)
    ),
    treatment = treatment_details(
      "trt",
      prior_normal(0, 1000)
    )
  )
  samples <- mcmc_sample(analysis_object, 200, 400, 1)
  draws <- samples$draws()
  renamed_draws <- rename_draws_covariates(draws, analysis_object)
  summary(renamed_draws)
}
#> Inputs look good.
#> Stan program compiled successfully!
#> Ready to go! Now call `mcmc_sample()`.
#> Running MCMC with 1 chain...
#> 
#> Chain 1 finished in 0.4 seconds.
#> # A tibble: 8 × 10
#>   variable           mean   median     sd   mad       q5      q95  rhat ess_bulk
#>   <chr>             <dbl>    <dbl>  <dbl> <dbl>    <dbl>    <dbl> <dbl>    <dbl>
#> 1 lp__           -1.57e+3 -1.57e+3  1.89  1.64  -1.57e+3 -1.56e+3 0.997     175.
#> 2 treatment log… -5.73e-1 -5.68e-1  0.227 0.229 -9.83e-1 -2.01e-1 1.00      169.
#> 3 commensurabil…  3.18e+1  9.73e+0 75.8   9.81   1.50e+0  1.23e+2 1.00      144.
#> 4 baseline log … -3.73e+0 -3.72e+0  0.185 0.172 -4.08e+0 -3.44e+0 0.999     206.
#> 5 baseline log … -3.23e+0 -3.23e+0  0.130 0.117 -3.46e+0 -3.03e+0 1.00      135.
#> 6 cov1            8.32e-1  8.29e-1  0.115 0.116  6.49e-1  1.03e+0 1.00      168.
#> 7 cov2            5.88e-1  5.88e-1  0.106 0.107  4.28e-1  7.69e-1 1.02      210.
#> 8 treatment HR    5.78e-1  5.67e-1  0.133 0.126  3.74e-1  8.18e-1 1.00      169.
#> # ℹ 1 more variable: ess_tail <dbl>