Next Article in Journal
Gender Marginalization in Sports Participation through Advertising: The Case of Nike
Next Article in Special Issue
Bayesian Sequential Monitoring of Single-Arm Trials: A Comparison of Futility Rules Based on Binary Data
Previous Article in Journal
Roles of Knowledge and Attitude in the Willingness of Nursing Students to Care for Older Adults in Hong Kong
Previous Article in Special Issue
Study Protocol: Phase I Dose Escalation Study of Oxaliplatin, Cisplatin and Doxorubicin Applied as PIPAC in Patients with Peritoneal Metastases
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Adjustment for Baseline Covariates to Increase Efficiency in RCTs with Binary Endpoint: A Comparison of Bayesian and Frequentist Approaches

1
Department of Clinical and Biological Sciences, University of Torino, 10100 Torino, Italy
2
Unit of Biostatistics, Epidemiology and Public Health, Department of Cardiac, Thoracic, Vascular Sciences and Public Health, University of Padova, 35121 Padova, Italy
3
Department of Medical Sciences, University of Ferrara, 44121 Ferrara, Italy
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2021, 18(15), 7758; https://doi.org/10.3390/ijerph18157758
Submission received: 11 June 2021 / Revised: 20 July 2021 / Accepted: 21 July 2021 / Published: 22 July 2021
(This article belongs to the Special Issue Bayesian Design in Clinical Trials)

Abstract

:
Background: In a randomized controlled trial (RCT) with binary outcome the estimate of the marginal treatment effect can be biased by prognostic baseline covariates adjustment. Methods that target the marginal odds ratio, allowing for improved precision and power, have been developed. Methods: The performance of different estimators for the treatment effect in the frequentist (targeted maximum likelihood estimator, inverse-probability-of-treatment weighting, parametric G-computation, and the semiparametric locally efficient estimator) and Bayesian (model averaging), adjustment for confounding, and generalized Bayesian causal effect estimation frameworks are assessed and compared in a simulation study under different scenarios. The use of these estimators is illustrated on an RCT in type II diabetes. Results: Model mis-specification does not increase the bias. The approaches that are not doubly robust have increased standard error (SE) under the scenario of mis-specification of the treatment model. The Bayesian estimators showed a higher type II error than frequentist estimators if noisy covariates are included in the treatment model. Conclusions: Adjusting for prognostic baseline covariates in the analysis of RCTs can have more power than intention-to-treat based tests. However, for some classes of model, when the regression model is mis-specified, inflated type I error and potential bias on treatment effect estimate may arise.

1. Introduction

Baseline covariates impact the outcome in many randomized controlled trials, and a recent systematic review reported that 84% of the trials present adjusted analysis. Among them, 91% pre-specified in the protocol such adjusted analysis [1]. It has been shown that models that adjust for baseline covariates can substantially improve the statistical power of the analysis when the covariates are moderately to strongly prognostic.
While this is justified for continuous outcomes, for binary outcomes, which require non-linear models, covariate adjustment may change the magnitude of the treatment effect, and thus the situation is subtler [2]. Due to the non-collapsibility of odds ratios, the non-adjusted and adjusted analyses estimate the marginal and the conditional treatment effect, respectively. However, the overall effect of adjusting for baseline covariates in logistic regression is still an increase in power. This is because the marginal estimate is always closer to the null effect than the conditional one, and the impact of the loss of precision on the power of the conditional estimate is offset by the larger effect size, leading to a net increase in power for the adjusted analyses [2,3].
The large amount of baseline covariates collected in an RCT opens the possibility to select the combination of covariates that results in the most favorable treatment effect estimate and/or the lowest p-value [4].
This is well-recognized in the “Guideline on adjustment for baseline covariates in clinical trials”, issued by the European Medicines Agency (EMA) in 2015, which requires pre-specification in the protocol of the variables to be included in the primary analysis for preventing the potential selection of the combination of covariates that may influence the treatment effect, especially in non-linear models [5].
However, pre-specification of the variables to be adjusted for is not always feasible as all prognostic variables may be not known in advance.
Under the frequentist approach, doubly robust and semi-parametric efficient estimators allow for the separation of treatment effect estimation from baseline covariate adjustment [6,7,8]. This is achieved by the inverse-probability-of-treatment weighting (IPTW) estimator, the parametric G-computation, the semiparametric locally efficient (SLE) estimator and the more recent targeted maximum likelihood estimator (TMLE). Under the Bayesian framework, model averaging is an alternative to the more common approach of model selection [9], which relies on estimation from a single model. While Bayesian model averaging (BMA) successfully accounts for model uncertainty in making a prediction, its advantages are less straightforward when used within the causal inference framework. In the context of causal treatment effect estimation, BMA tends to assign large posterior probabilities to models that may not accurately adjust for confounding. To overcome this drawback, the Bayesian adjustment for confounding (BAC) algorithm has been proposed as an alternative approach based on the specification of both an outcome and a treatment model, as in the propensity score framework [10].
However, since BMA and BAC are based on models that likely contain noisy prognostic covariates, they lose precision in estimating the treatment effect. To overcome this limitation, the generalized Bayesian causal effect estimation (GBCEE) has been proposed as a further unbiased and efficient estimator [11].
This study investigates which methods of adjusting for baseline covariate in the analysis of RCTs with binary endpoint maximize the statistical power while retaining the type I error rate and unbiased estimate of treatment effect. Such comparison is justified because type I error and power are still the study operating characteristics of concern healthcare regulators require when appraising the results of confirmatory clinical trials [12].
In the following, in Section 2.1, the motivation example is introduced. Then in Section 2.2, the simulation study is explained, and the frequentists and Bayesian estimators are briefly presented. Results of both simulation and the illustrative study will be reported in Section 3 and finally discussed in Section 4.

2. Materials and Methods

2.1. Illustrative Study and Simulated Data

Our simulation study was based on the motivating example of re-analyzing the PROLOGUE RCT [13]. The PROLOGUE study is among the largest trials investigating whether DPP-4 inhibitors provide cardiovascular protective effects to patients with type 2 diabetes by slowing carotid stiffness progression associated with conventional diabetes treatment.
The study participants were either allocated to add-on DPP-4 inhibitor (sitagliptin) treatment or to continue therapy with conventional anti-diabetic agents. The primary endpoint was the arterial stiffness of annual changes, which resulted in being not significantly different between the two groups. However, the study showed that the decrease in glycated haemoglobin (HbA1c) in patients treated with sitagliptin was superior to conventional therapy, proving a better glycemic control. As a re-analysis of the PROLOGUE study, we then investigated a potential sitagliptin effect on the improvement of HbA1c.

2.2. Simulation Study

The simulation study was carried out to compare the performance of several estimators applied to obtain a marginal treatment effect estimate in the case of a binomial outcome and was based on the same scheme adopted in Zhang et al. [14].
There was a 50% chance of being assigned to either the treatment or the control group. The treatment assignment variable (Z) was generated as Bernoulli with P(Z = 1) = P(Z = 0) = 0.5. The assignment Z = 1 corresponds to the treatment group. The baseline covariates were generated as follows:
  • X 1 , X 3 ,   X 8 follow a Normal (0,1) distribution;
  • X 4 follows a Bernoulli (0.3) distribution;
  • X 6 follows a Bernoulli (0.5) distribution;
  • X 2 was generated as 0.2 × X 1 + 0.98   U 1
  • X 5 was generated as 0.1 × X 1 + 0.2 × X 2 + 0.97   U 2
  • X 7 was generated as 0.1 × X 3 + 0.99   U 3
where U 1 , U 2 ,   U 3 are Normal (0,1) variables.
Finally, Y was generated as Y = logit   ( P ( Y = 1 | Z ,   X ) ) = α + β Z + γ X , where X = (X1, …, X8) is the matrix of covariates, α = 0.9 ,   β = 1.3 ,   γ = ( 0.5 ,   1.3 ,   0.5 ,   1.5 , 0 , 0 , 0 , 0 ). The parameter β is the conditional treatment effect; α is the intercept and γ is the vector of the coefficients of covariates X 1 , , X 8 . Thus, X 1 , , X 4 variables represent prognostic patient features for treatment effect.
To find the marginal treatment effect, one million individuals were simulated, and 30 repetitions were performed. The marginal treatment effect was then calculated as the mean of the treatment effects as the log odds ratios using unadjusted logistic regression. The true marginal treatment effect resulted in being equal to −0.871 ± 0.004.
For the simulation study, 5000 datasets of sample size n = 200 were generated. For the frequentist estimators, several scenarios were defined to evaluate the effects of model selection and are reported in Table 1.
The model estimated under the correct scenario is the same used to generate the outcome data when all prognostic variables are known. The model estimated under the mis-specification scenario includes only one prognostic variable and an additional noisy variable. Finally, the model estimated under the all-variables scenario includes all the prognostic variables as well as non-prognostic variables and mimics the situation of using all patient features for the treatment effect estimation in the case of uncertainty about knowledge on prognostic variables.

2.2.1. Frequentist Estimators

The frequentist estimators employed for the estimation of the treatment effect are briefly presented. In describing the estimators, we will refer to the treatment model as the conditional probability (likelihood) of being treated given the covariates, i.e., P ( Z | X ) , and to the outcome model as the probability, i.e., likelihood, of the outcome given the treatment and the prognostic covariates, i.e., P(Y|Z, X).
G-computation. To address confounding, G-computation relies on the estimation of the outcome model, i.e., the conditional expectation of the outcome given the treatment and the prognostic covariates. Contrary to the propensity score methods, it does not require estimating the exposure mechanism or treatment model, i.e., the conditional probability of being treated given the observed confounders [15].
Doubly Robust Estimation. Doubly robust (DR) estimation combines the outcome model and treatment model. Both the models are unbiased only if they are correctly specified. The DR estimation ensures that when combining the two models for the treatment effect estimation, only one of them must be correctly specified to obtain an unbiased estimate. The estimates of the parameters of interest of the outcome model and the treatment model are used to predict patients’ responses under the treatment condition and the treatment assignment (propensity score), respectively. Inversely weighting the expected response under treatment condition by the propensity score allows us to represent the estimator of the quantity of interest as an unbiased estimate plus a second term. This term reduces to 0 if either the treatment model or the outcome model are correctly specified and if, and only if, the possibly incorrect conditional density has the same support as the true conditional density [6].
Semi-Parametric Locally Efficient Estimator. It uses a semi-parametric model for the outcome model, which is used to generate predicted values separately from the treatment model. Finally, it computes the average treatment effect as the mean difference in predicted outcome pair across individuals [16].
Targeted Maximum Likelihood Estimator. TMLE is a doubly robust, maximum-likelihood–based estimation method that includes a secondary targeting step that optimizes the bias-variance tradeoff for the estimation of the parameter of interest. TMLE is particularly attractive for causal effect estimation in RCT analysis. First, it is a doubly robust method, which yields unbiased estimates if the treatment model is correctly specified, which is the case of RCT setup [17].
Augmented Inverse Probability Weighting. Propensity scores are estimated and used to create inverse probability weights; all observations are weighted. Finally, it computes the average treatment effect as the mean difference between weighted outcomes among exposed and unexposed [18].

2.2.2. Bayesian Estimators

Bayesian Model Averaging. BMA is an extension of the Bayesian inference methods. In addition to the usual modelling of parameter uncertainty through the prior distribution, it models the uncertainty of the model selection process, obtaining a posterior parameter and posterior probability model through Bayes’ theorem. In the present work, we considered Zellner’s g distribution as a-priori distribution on coefficients for the variable selection [19] and the Bayesian adaptive sampling algorithm for the model selection [9].
Bayesian Adjustment for Confounding. As in the propensity score framework, BAC requires the definition of the outcome model, which is a function of the treatment and potential confounders, and the treatment model, which is a function of the potential confounders to treatment assignment. Then it applies a Bayesian variable selection process in both models to select covariates and introduces a dependence parameter between the outcome and treatment model, ω, which denotes the prior odds of including a confounder in the outcome model, given that the same confounder is in the exposure model. In the special case of dependence parameter ω = 1, BAC reduces to BMA [10].
Generalized Bayesian Causal Effect Estimation. The generalized Bayesian causal effect estimation (GBCEE) algorithm performs variable selection and delivers doubly robust estimates. It employs a prior distribution that targets the selection of true confounders and predictors of the outcome. It thus takes advantage of the Bayesian framework to account for uncertainty in the model selection process. It is different from BMA in building the prior distribution. It uses a prior distribution tailored to identify the potential confounders, which uses information from the data and thus relies on the empirical Bayes approach. Finally, it adds a doubly robust estimation, employing the posterior distribution of the parameters and adopting the TMLE framework to estimate the causal effect that protects against model mis-specification [11].

3. Results

3.1. Simulation Study

For each method, we computed the bias as the difference between the average of the estimates and the true effect. The standard error (SE) of the estimates, the Monte Carlo standard error for the standard deviation (MC SD) and the coverage probability (CP), i.e., the proportion of simulation replicates in which the 95% confidence intervals included the true effect. For the frequentist estimators, 95% confidence intervals were computed. For the Bayesian estimator GBCEE, 95% confidence intervals were computed as well, using 50 non-parametric bootstrap replicates with the percentile method. For the BAC and BMA approach, the 95% credible intervals were computed.
Type I error and power were also calculated. For both frequentist and Bayesian estimators, type I error was computed simulating 5000 datasets under the null hypothesis that the treatment is not effective. For BAC and BMA, type I error was estimated by the proportion of the simulations incorrectly declared the treatment effective, based on the posterior probability P ( β   <   0 | Y ,   X 1 , ,   X 8 ) 0.95 .
Similarly, the power was calculated as the proportion of simulations that declare the trial successful based on the given decision criteria when the target treatment effect is assumed to be the true value. This approach has been recommended by the FDA [20].
The performance of the frequentist estimators was assessed under three scenarios: the ideal case of the correct model specification (correct scenario); the case of important prognostic variable not identified in the model specification (mis-specification scenario); finally, the case when noisy prognostic variables are introduced in the model (all-variables scenario). In the all-variables scenario, for SLE estimator, a model selection process was foreseen based on backward and forward stepwise techniques.
Bayesian estimators’ performance was assessed on the all-variables scenario only since they do not require selecting a final model but allow for averaging over the space of potential models that could have generated the data.
Results of the simulation study are reported in Table 2. The bias is similar across all methods, while more variation is observed in the power of the estimators, ranging between 84.6% and 94.9%. For Bayesian estimators, the power is given by the posterior probability of observing a treatment marginal effect greater than zero. A slight inflation of type I error is observed, except for BMA, however, it is not greater than 6.5%.

3.2. Illustrative Study

To illustrate the effect of baseline adjustment on the treatment effect estimation, we applied the introduced methods to re-analyze the PROLOGUE RCT [13].
The PROLOGUE RCT aimed to investigate whether DPP-4 inhibitors provide cardiovascular protective effects to patients with type 2 diabetes.
The study participants were either allocated to add-on DPP-4 inhibitor (sitagliptin) treatment or to continue therapy with conventional anti-diabetic agents. The study showed that the decrease in glycated haemoglobin (HbA1c) in patients treated with sitagliptin was superior to conventional therapy, proving a better glycemic control.
We set as outcome an improvement of at least 1% in HbA1c, obtaining a dichotomised outcome. This choice is motivated by the observation that two large-scale studies—the UK Prospective Diabetes Study (UKPDS) and the Diabetes Control and Complications Trial (DCCT)—demonstrated that improving HbA1c by 1% (or 11 mmol/mol) for people with type 1 diabetes or type 2 diabetes cuts the risk of microvascular complications by 25%.
As prognostic covariates, we used age (years), gender (female, male), body mass index (BMI, kg/cm2), systolic blood pressure (SBP, mmHg), low-density lipoprotein (LDL, mg/dL), high-density lipoprotein (HDL, mg/dL), HbA1c (%), fasting plasma glucose (FPG, mmol/L), dyslipidemia (LDL ≥ 130 mg/dL odds ratio (OR) HDL < 35 mg/dL OR triglyceride ≥ 150 mg/dL OR total cholesterol (=LDL + HDL + (Triglyceride/5)) ≥ 200 mg/dL).
In Figure 1, the unadjusted OR of improving HbA1c by 1% is reported along with 95% confidence interval. Frequentist estimates are reported with 95% confidence intervals, and finally, Bayesian estimates with 95% credible intervals are listed.

4. Discussion

We have presented a study to compare different approaches to address covariate adjustment to estimate treatment effect in RCTs.
Baseline covariates adjustment impacts the outcome in many RCTs in terms of power, type I error and bias of the marginal effect estimation.
In fact, variable selection methods based on p-values, and large observed differences between arms and stepwise approaches, have increased type I error rates [4]. The guideline on adjustment for baseline covariates in clinical trials issued by the EMA in 2015 strongly recommends pre-specifying the variables to be included in the primary analysis in the protocol to avoid bias and potential selection of the combination of covariates that may favour the treatment [4]. Moreover, the Consolidated Standards of Reporting Trials (CONSORT) [21] and the International Conference on Harmonization [22] recommend to pre-specify the potential prognostic variables to employ in adjusted analysis. However, there is still debate on how to identify prognostic covariates correctly.
Several approaches have been proposed to estimate a marginal causal effect, which is the standard measure of treatment effect reported when analyzing RCTs [23]. Although it is known from the literature that adjustment for prognostic covariates can increase the efficiency of the analysis, there is still a lack of attempts to assess comparatively the performance of all the methods under real scenarios of analysis, including adjustment for non-prognostic variables and model mis-specification.
We compared several frequentist and Bayesian estimators under different scenarios. The selected estimators were: SLE estimation, TMLE, G-computation, AIPTW, DR estimator, GBCEE, BMA and the Bayesian adjustment for confounding algorithm. We assessed their performance under three scenarios: the ideal case of the correct model specification; the case of important prognostic variable not identified in the model specification (model mis-specification); finally, the case when noisy prognostic variables are introduced in the model (all variables selected in the adjusted analysis). Since the Bayesian estimators can handle the uncertainty of the model selection process assigning a posterior probability to each set of covariates [24], they were assessed only under the scenario of all variables included in the adjusted analysis.
Our results from the simulation study showed that model mis-specification does not increase the bias. This holds also for the G-computation estimator, which is not theoretically guaranteed to be a consistent estimator under model mis-specification.
The approaches that are not doubly robust have increased MC SD. They also showed increased SE under the scenario of mis-specification.
Across different scenarios, frequentist estimators showed a similar precision (SE ranges between 0.244 and 0.298). This observation is particularly interesting since the correct specification of a parametric model with many covariates is nearly impossible [17]. Bayesian estimators behave differently, showing a high precision for GBCEE (SE = 0.147). The uncertainty of BMA is not directly comparable with SE, since it is the standard deviation computed on the posterior distribution. Thus it showed a larger uncertainty (0.451), which is expected since it embeds the uncertainty of estimates in the posterior probability function.
Covariate adjustment reduced type II error but under the scenario of mis-specification of the outcome model. The Bayesian estimators showed a higher type II error than frequentist estimators under the scenario, including all prognostic variables and noisy covariates in the model specification. On the other hand, BMA showed the largest bias, even if offset by the smallest type I error, which is not surprising since it has been shown that the bias can be relevant when covariates are only slightly associated with the outcome [11].
In the re-analysis of PROLOGUE RCT, we estimated the odds ratio of improving HbA1c using the frequentist and the Bayesian estimators introduced for avoiding confounding. Among the Bayesian estimators, GBCEE resulted similar to other frequentist estimators due to its doubly robust property. In contrast, BMA and BAC showed a smaller treatment effect, compared to the unadjusted estimate.
Adjusting for prognostic covariates leads to an increase in power, as seen by observing that the adjusted estimate is farther from the null value of 1 than the unadjusted estimate (odds ratio equal to 1 indicates no treatment effect). Compared to the unadjusted analysis, we did not observe a dramatic increase of SE, thus a loss of precision, except for BMA and BAC. However, the GBCEE Bayesian estimators showed performances comparable to other frequentist estimators.

5. Conclusions

Adjusting for baseline covariates predictive of outcome in the analysis of RCTs can have more power than intention-to-treat based tests. However, for some classes of model, when the regression model is mis-specified, inflated type I error and potential bias on treatment effect estimate may arise. Estimators that allow for separating the baseline covariate adjustment from the treatment effect estimation can avoid potential bias for covariates’ post hoc selection retaining the focus on objective inference on treatment effect. Among Bayesian estimators, BMA presents the largest bias.
Our simulations were carried out in the context of a binary outcome. Similar conclusions are likely to be applied to the hazard ratio since odds-ratio and hazard ratio showed the same non-collapsibility issue.
Limitations of this study rely on the assumption of independent, identical distribution of data, which is not necessarily the case in RCTs. Patients in RCTs often have wide variability in their response to treatment resulting in heterogeneity of treatment effect. Further research should include realistic synthetic datasets, which capture the relationships across clinical features among patients. Probabilistic models, classification-based imputation models, and generative adversarial neural networks are an example of data-driven approaches of synthetic data generation methods [25].

Author Contributions

Conceptualization, P.B. and I.B.; formal analysis, P.B.; data curation, C.L.; writing—original draft preparation, P.B.; writing—review and editing, S.U., V.S., D.A. and D.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the University of Torino, grant number BERP_RILO_17_01.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable. This study used secondary data already collected and publicly available online.

Data Availability Statement

The data and code for the analysis are available upon request to the authors.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ciolino, J.D.; Palac, H.L.; Yang, A.; Vaca, M.; Belli, H.M. Ideal vs. Real: A Systematic Review on Handling Covariates in Randomized Controlled Trials. BMC Med. Res. Methodol. 2019, 19, 136. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Robinson, L.D.; Jewell, N.P. Some Surprising Results about Covariate Adjustment in Logistic Regression Models. Int. Stat. Rev. Rev. Int. Stat. 1991, 59, 227. [Google Scholar] [CrossRef]
  3. Kahan, B.C.; Jairath, V.; Doré, C.J.; Morris, T.P. The Risks and Rewards of Covariate Adjustment in Randomized Trials: An Assessment of 12 Outcomes from 8 Studies. Trials 2014, 15, 139. [Google Scholar] [CrossRef] [PubMed]
  4. Raab, G.M.; Day, S.; Sales, J. How to Select Covariates to Include in the Analysis of a Clinical Trial. Control. Clin. Trials 2000, 21, 330–342. [Google Scholar] [CrossRef]
  5. Committee for Medicinal Products for Human Use (CHMP) Guideline on Adjustment for Baseline Covariates in Clinical Trials. 2015. Available online: www.Ema.Europa.Eu/Contact (accessed on 20 May 2021).
  6. Tsiatis, A.A. Locally Efficient Semiparametric Estimators for Functional Measurement Error Models. Biometrika 2004, 91, 835–848. [Google Scholar] [CrossRef]
  7. Funk, M.J.; Westreich, D.; Wiesen, C.; Stürmer, T.; Brookhart, M.A.; Davidian, M. Doubly Robust Estimation of Causal Effects. Am. J. Epidemiol. 2011, 173, 761–767. [Google Scholar] [CrossRef] [PubMed]
  8. Zhang, M.; Gilbert, P.B. Increasing the Efficiency of Prevention Trials by Incorporating Baseline Covariates. Stat. Commun. Infect. Dis. 2010, 2. [Google Scholar] [CrossRef]
  9. Clyde, M.A.; Ghosh, J.; Littman, M.L. Bayesian Adaptive Sampling for Variable Selection and Model Averaging. J. Comput. Graph. Stat. 2011, 20, 80–101. [Google Scholar] [CrossRef]
  10. Wang, C.; Parmigiani, G.; Dominici, F. Bayesian Effect Estimation Accounting for Adjustment Uncertainty. Biometrics 2012, 68, 661–671. [Google Scholar] [CrossRef] [Green Version]
  11. Talbot, D.; Lefebvre, G.; Atherton, J. The Bayesian Causal Effect Estimation Algorithm. J. Causal Inference 2015, 3, 207–236. [Google Scholar] [CrossRef]
  12. Ryan, E.G.; Brock, K.; Gates, S.; Slade, D. Do We Need to Adjust for Interim Analyses in a Bayesian Adaptive Trial Design? BMC Med. Res. Methodol. 2020, 20, 150. [Google Scholar] [CrossRef]
  13. Oyama, J.; Murohara, T.; Kitakaze, M.; Ishizu, T.; Sato, Y.; Kitagawa, K.; Kamiya, H.; Ajioka, M.; Ishihara, M.; Dai, K.; et al. The Effect of Sitagliptin on Carotid Artery Atherosclerosis in Type 2 Diabetes: The PROLOGUE Randomized Controlled Trial. PLoS Med. 2016, 13, e1002051. [Google Scholar] [CrossRef] [PubMed]
  14. Zhang, M.; Tsiatis, A.A.; Davidian, M. Improving Efficiency of Inferences in Randomized Clinical Trials Using Auxiliary Covariates. Biometrics 2008, 64, 707–715. [Google Scholar] [CrossRef] [PubMed]
  15. Wang, A.; Nianogo, R.A.; Arah, O.A. G-Computation of Average Treatment Effects on the Treated and the Untreated. BMC Med. Res. Methodol. 2017, 17, 3. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Robins, J.M.; Rotnitzky, A.; Zhao, L.P. Estimation of Regression Coefficients When Some Regressors Are Not Always Observed. J. Am. Stat. Assoc. 1994, 89, 846–866. [Google Scholar] [CrossRef]
  17. Laan, M.V.D.; Rose, S. Targeted Learning: Causal Inference for Observational and Experimental Data; Series in Statistics; Springer: New York, NY, USA, 2011. [Google Scholar]
  18. Glynn, A.N.; Quinn, K.M. An Introduction to the Augmented Inverse Propensity Weighted Estimator. Polit. Anal. 2010, 18, 36–56. [Google Scholar] [CrossRef] [Green Version]
  19. Liang, F.; Paulo, R.; Molina, G.; Clyde, M.A.; Berger, J.O. Mixtures of g Priors for Bayesian Variable Selection. J. Am. Stat. Assoc. 2008, 103, 410–423. [Google Scholar] [CrossRef]
  20. U.S. Food and Drug Administration Adaptive Designs for Clinical Trials of Drugs and Biologics: Guidance for Industry. 2019. Available online: https://www.fda.gov/regulatory-information/search-fda-guidance-documents/adaptive-design-clinical-trials-drugs-and-biologics-guidance-industry (accessed on 18 July 2021).
  21. Moher, D.; Hopewell, S.; Schulz, K.F.; Montori, V.; Gøtzsche, P.C.; Devereaux, P.J.; Elbourne, D.; Egger, M.; Altman, D.G. CONSORT 2010 Explanation and Elaboration: Updated Guidelines for Reporting Parallel Group Randomised Trials. J. Clin. Epidemiol. 2010, 63, e1–e37. [Google Scholar] [CrossRef] [Green Version]
  22. ICH Official Web Site: ICH. Available online: https://www.ich.org/ (accessed on 9 July 2021).
  23. Martens, E.P.; Pestman, W.R.; Klungel, O.H. Conditioning on the Propensity Score Can Result in Biased Estimation of Common Measures of Treatment Effect: A Monte Carlo Study (p n/a) by Peter C. Austin, Paul Grootendorst, Sharon-Lise T. Normand, Geoffrey M. Anderson, Statistics in Medicine. Stat. Med. 2007, 26, 3208–3210, author reply 3210–3212, doi:10.1002/sim.2878. [Google Scholar] [CrossRef]
  24. Volinsky, C.T.; Madigan, D.; Raftery, A.E.; Kronmal, R.A. Bayesian Model Averaging in Proportional Hazard Models: Assessing the Risk of a Stroke. J. R. Stat. Soc. Ser. C 1997, 46, 433–448. [Google Scholar] [CrossRef] [Green Version]
  25. Goncalves, A.; Ray, P.; Soper, B.; Stevens, J.; Coyle, L.; Sales, A.P. Generation and Evaluation of Synthetic Patient Data. BMC Med. Res. Methodol. 2020, 20, 108. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Odds ratio 1% change in Hba1c at 24 months.
Figure 1. Odds ratio 1% change in Hba1c at 24 months.
Ijerph 18 07758 g001
Table 1. Scenarios under which the estimators were compared. The model for outcome generation is Y = logit   P ( Y | Z ,   X ) = 0.9 + 1.3 Z + 0.5 X 1 + 1.3 X 2 + 0.5 X 3 + 1.5 X 4 .
Table 1. Scenarios under which the estimators were compared. The model for outcome generation is Y = logit   P ( Y | Z ,   X ) = 0.9 + 1.3 Z + 0.5 X 1 + 1.3 X 2 + 0.5 X 3 + 1.5 X 4 .
ScenarioOutcome Model EstimatedPrognostic Variables in the Outcome Model EstimatedNon-Prognostic Variables in the Outcome Model Estimated
Correct E ( Y | Z ,   X ) = α + β Z +
γ 1 X 1 + γ 2 X 2 + γ 3 X 3 + γ 4 X 4
X 1 , X 2 ,   X 3 ,   X 4 none
Misspecification E ( Y | Z ,   X ) = α + β Z + γ 3 X 3 + γ 5 X 5 X 3 X 8
All-variables E ( Y | Z ,   X ) = α + β Z +
γ 1 X 1 + γ 2 X 2 + γ 3 X 3 + γ 4 X 4 + γ 5 X 5 + γ 6 X 6 + γ 7 X 7 + γ 8 X 8
X 1 , X 2 ,   X 3 ,   X 4 X 5 , X 6 ,   X 7 ,   X 8
Frequentist estimators were compared on all the three scenarios. Bayesian estimators were compared on all-variables scenario, only.
Table 2. Results of the simulation study. In semiparametric locally efficient (SLE)/all-variables scenario, a model selection process based on backward (SLE backward) and forward (SLE forward) stepwise techniques were foreseen.
Table 2. Results of the simulation study. In semiparametric locally efficient (SLE)/all-variables scenario, a model selection process based on backward (SLE backward) and forward (SLE forward) stepwise techniques were foreseen.
MethodScenarioBIASSEMC SDPowerType II ErrorCPType I
Error
SLEAll variables−0.0100.2460.0150.9470.0530.9410.061
SLE BackwardAll variables−0.0110.2470.0150.9470.0530.9410.061
SLE ForwardAll variables−0.0100.2460.0150.9470.0530.9410.061
TMLEAll variables−0.0110.2480.0150.9470.0530.9390.065
G-CompAll variables−0.0110.2480.0150.9430.0570.9440.060
AIPTWAll variables−0.0110.2480.0150.9420.0580.9420.058
DRAll variables−0.0130.2980.0100.9430.0570.9400.059
SLEMis-specification−0.0100.2460.0150.8560.1440.9510.051
TMLEMis-specification−0.0110.2480.0160.8520.1480.9470.056
G-CompMis-specification−0.0110.2480.0160.8460.1540.9490.054
AIPTWMis-specification−0.0080.2440.0150.8470.1530.9480.053
DRMis-specification−0.0080.2440.0150.8470.1530.9480.052
SLECorrect−0.0080.2440.0150.9490.0510.9460.056
TMLECorrect−0.0120.2960.0100.9490.0510.9460.054
G-CompCorrect−0.0120.2960.0100.9490.0510.9470.052
AIPTWCorrect−0.0120.2980.0110.9490.0510.9480.052
DRCorrect−0.0130.2950.0110.9490.0510.9470.053
GBCEEAll variables−0.0140.1470.0220.9160.0840.9520.063
BACAll variables−0.0150.299 10.0460.9020.0980.9450.045
BMAAll variables−0.0510.451 10.0900.9220.0780.9420.02
1 The value is the standard deviation of the posterior distribution. Semi-parametric Locally Efficient (SLE) Estimator: Targeted Maximum Likelihood Estimator (TMLE); G-Computation (G-Comp); Augmented Inverse Probability Weighting (AIPTW); Doubly Robust (DR); Generalized Bayesian Causal Effect Estimation (GBCEE); Bayesian Adjustment for Confounding (BAC); Bayesian Model Average (BMA).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Berchialla, P.; Sciannameo, V.; Urru, S.; Lanera, C.; Azzolina, D.; Gregori, D.; Baldi, I. Adjustment for Baseline Covariates to Increase Efficiency in RCTs with Binary Endpoint: A Comparison of Bayesian and Frequentist Approaches. Int. J. Environ. Res. Public Health 2021, 18, 7758. https://doi.org/10.3390/ijerph18157758

AMA Style

Berchialla P, Sciannameo V, Urru S, Lanera C, Azzolina D, Gregori D, Baldi I. Adjustment for Baseline Covariates to Increase Efficiency in RCTs with Binary Endpoint: A Comparison of Bayesian and Frequentist Approaches. International Journal of Environmental Research and Public Health. 2021; 18(15):7758. https://doi.org/10.3390/ijerph18157758

Chicago/Turabian Style

Berchialla, Paola, Veronica Sciannameo, Sara Urru, Corrado Lanera, Danila Azzolina, Dario Gregori, and Ileana Baldi. 2021. "Adjustment for Baseline Covariates to Increase Efficiency in RCTs with Binary Endpoint: A Comparison of Bayesian and Frequentist Approaches" International Journal of Environmental Research and Public Health 18, no. 15: 7758. https://doi.org/10.3390/ijerph18157758

APA Style

Berchialla, P., Sciannameo, V., Urru, S., Lanera, C., Azzolina, D., Gregori, D., & Baldi, I. (2021). Adjustment for Baseline Covariates to Increase Efficiency in RCTs with Binary Endpoint: A Comparison of Bayesian and Frequentist Approaches. International Journal of Environmental Research and Public Health, 18(15), 7758. https://doi.org/10.3390/ijerph18157758

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop