 Research article
 Open Access
 Published:
Inverse probability weighting to estimate causal effect of a singular phase in a multiphase randomized clinical trial for multiple myeloma
BMC Medical Research Methodology volume 16, Article number: 150 (2016)
Abstract
Background
Randomization procedure in randomized controlled trials (RCTs) permits an unbiased estimation of causal effects. However, in clinical practice, differential compliance between arms may cause a strong violation of randomization balance and biased treatment effect among those who comply. We evaluated the effect of the consolidation phase on diseasefree survival of patients with multiple myeloma in an RCT designed for another purpose, adjusting for potential selection bias due to different compliance to previous treatment phases.
Methods
We computed two propensity scores (PS) to model two different selection processes: the first to undergo autologous stem cell transplantation, the second to begin consolidation therapy. Combined stabilized inverse probability treatment weights were then introduced in the Cox model to estimate the causal effect of consolidation therapy miming an ad hoc RCT protocol.
Results
We found that the effect of consolidation therapy was restricted to the first 18 months of the phase (HR: 0.40, robust 95 % CI: 0.170.96), after which it disappeared.
Conclusions
PSbased methods could be a complementary approach within an RCT context to evaluate the effect of the last phase of a complex therapeutic strategy, adjusting for potential selection bias caused by different compliance to the previous phases of the therapeutic scheme, in order to simulate an ad hoc randomization procedure.
Trial registration
ClinicalTrials.gov: NCT01134484 May 28, 2010 (retrospectively registered)
EudraCT: 200500372339 December 17, 2008 (retrospectively registered)
Background
Randomized controlled trials (RCTs) are the gold standard for clinical research because the randomization procedure is expected to achieve a balanced distribution of both known and unknown baseline characteristics between experimental and control treatment arms [1]. Unfortunately, selection mechanisms may occur after randomization, i.e. some patients are lost to follow up or drop out early or occurrence of competing risk events such as death from other diseases. Reasons for early drop out include informed consent withdrawal or toxicities. When a large number of patients do not comply with the therapeutic protocol, the effect of treatment may be difficult or impossible to estimate. In fact, the consequence of selection bias is that the association between treatment and outcome among those selected for analysis (because entirely compliant with treatment protocol) may differ from that among those who are eligible. Patients may decide not to follow, partially follow, or entirely follow a complex therapeutic program for reasons related to the outcome [2]. The intentiontotreat (ITT) analysis provides an unbiased estimate of gross treatment effect but is shrunk toward the null value on the basis of the percentage of compliers. Conversely, the perprotocol (PP) analysis may introduce bias because the groups of patients being compared no longer have similar characteristics. Thus, it is possible that some patients are more likely to be treated than others (possibly depending on prognostic factors) in which case treatment may erroneously appear more efficacious [3–5].
Another important issue is that RCTs usually compare the entire therapeutic strategy without any detailed evaluation of each singular therapeutic phase. The evaluation of a specific phase is complicated because it may be biased by different compliance to the previous phases. In fact, during the course of a trial, especially when the study protocol is complex and the disease is severe, randomization balance may be lost right from the earliest phases due to selective patient withdrawal between one phase and another.
Several methods have been proposed to tackle these problems in an RCT setting where postrandomization imbalance occurs, and also within the context of observational studies where randomization is not possible. One popular approach is the use of the propensity score (PS), a method developed to reduce the effect of observable confounding factors [6].
The aim of this article was to evaluate the effect of the consolidation therapy, the last phase of a complex therapeutic strategy composed by two arms of induction therapy randomly administered to prepare patients to receive at least one autologous stem cell transplantation (ASCT), in patients with previously untreated multiple myeloma. Only patients that received at least one ASCT can proceed to receive consolidation therapy. Those receiving consolidation therapy receive the same therapy to which they were randomized during the induction phase. For these reasons, we needed to adjust for potential selection bias due to different compliance to previous phases of therapeutic protocol in order to approximate an RCT with randomization procedure (consolidation vs not consolidation phase) after ASCT (with time 0 being the end of the last ASCT received). We did it by modelling both the probability to receive at least 1 ASCT and the probability to receive consolidation phase through a combined PS approach.
This study is justified by the increasing interest being shown in consolidation therapy on progression free survival endpoint [7]. Consolidation treatment is one of the latest phases of MM therapeutic protocols and aims to further increase the frequency and depth of clinical response obtained with previous treatments phases.
The paper is structured as follows. In section 2, we describe the statistical methods used, in particular, the definition of the PS, the weighting procedure in the survival analysis, and the Aalen’s additive hazards model to investigate timevarying covariates. In section 3, we report the results of the role of the consolidation phase in the treatment of MM. In section 4 and 5, we discuss the results and draw our conclusions.
Methods
In this section we describe the statistical methods used to restore randomization balance to obtain an unbiased estimate of the causal effect of a treatment. In particular, we review the definition of the PS and its possible uses, the weighting procedure in survival analysis, and the Aalen’s additive hazards model to investigate timevarying effects.
Propensity score
The propensity score (PS) is the conditional probability of being treated given a set of observed potential confounders. In this way all the information from a large number of potential confounders is summarized into a unique balancing score variable (the socalled propensity score). The PS may warrant that the distribution of measured baseline covariates is the same in treated and untreated subjects. Rosenbaum and Rubin said that PS can account for imbalance in treatment groups and reduce bias by simulating a sort of “virtual randomization” of subjects into treatment groups (conditional exchangeability) [6].
The PS is the conditional probability of receiving a treatment given pretreatment characteristics:
where Z = {0, 1} is the indicator of exposure to treatment and X is the multidimensional vector of pretreatment characteristics. Thus, if exposure to treatment is random within cells defined by X, it is also random within cells defined by the values of the onedimensional variable p(X). Biasremoving adjustments can therefore be made using the PS alone rather than modeling all of the background covariates individually.
The expected difference in observed responses at p(X) is equal to the average treatment effect at p(X), E{r_{1}p(X), Z = 1}  E{r_{0}p(X), Z = 0} = E{r_{1}  r_{0}p(X)} where r = {0, 1} is the indicator of the response that would have resulted if a patient had or had not received treatment, respectively; these are the potential outcomes in the two counterfactual situations of treatment and no treatment ([6], Theorem 4 in Rosenbaum and Rubin); [8, 9].
There is still some debate about how to determine a sufficient set of covariates in the PS model. Brookhart states that variables that are unrelated to the exposure but related to the outcome should always be included in a PS model. The inclusion of these variables decreases the variance of an estimated treatment effect without increasing bias. In contrast, including variables that are related to the exposure but not to the outcome increases the variance of the estimated exposure effect without decreasing bias [10]. Of course if a variable is related to both exposure and outcome, it must be included. Conversely, Pearl argues that, although treated and untreated units are balanced in each stratum of the PS, the balance only holds relative to the covariates measured; unobserved confounders may be highly unbalanced in each stratum. “Such imbalance may be dormant in the crude estimate and awakened through the use of PS methods” ([11], page 1415). The effectiveness of PS methods rests critically on the choice of a sufficient set of covariates, which requires in depth knowledge about the causal relationships among both observed and unobserved covariates [12].
Inverse probability of treatment weight
Austin explains that there are mainly four ways of using the PS to reduce or minimize the effects of confounding when estimating the effects of treatments on outcomes: matching on the PS, stratification on the PS, inverse probability of treatment weighting (IPTW) using the PS, and covariate adjustment using the PS [13]. We chose the IPTW method because we were mainly interested in population effect, i.e. the average treatment effect over the marginal distribution of observed covariates in the study sample. Moreover, Austin’s results showed that, while all ways of using PS allow for the estimation of marginal hazard ratios with minimal bias, IPTW provides estimates with lower mean squared error; it is also less subject to loss of information due to lack of matching.
IPTW is calculated by the inverse of the conditional probability of receiving the exposure that a patient indeed received:
where Z indicates whether or not the subject was treated while p(X) represents the conditional probability for the subject to be treated. Both treated subjects with a very low PS and untreated subjects with a high PS have large IPTWs to account for unequal probability of receiving treatment in the original sample.
The aim of IPTW is to reduce selection bias by creating a “pseudopopulation” in which the exposure is independent of the measured confounders so that the treatment effect estimate in a sample thus weighted will be less biased [14–17].
Stabilized inverse probability of treatment weight
It may happen that treated subjects have a PS near 0 or that untreated subjects have a PS near 1, making the relative IPTW excessively high and unstable. Computationally, Xu and Ross noticed that, as in any weighted regression, unstabilized IPTW changes the sample size of the original sample, generating an underestimate of the variance of the estimate of the effect, producing inappropriately narrow confidence intervals and leading to the lack of control of the probability of a type I error [15].
Stabilized inverse probability of treatment weight (SIPTW) can be obtained by multiplying the IPTW by the marginal probability of receiving the actual treatment received. Moreover, it preserves the sample size of the original data, produces appropriate estimation of the variance of the main effect, and adequately controls the type I error rate [18, 19].
Adjusted KaplanMeier survival curves
Cole et al. demonstrated that, under the assumption of no unmeasured confounding, adjusted KaplanMeier estimates for survival curves of treated and untreated  where each subject is weighted by its IPTW  represent the survival curve of the entire sample had none been exposed, and the survival curve of the entire sample had everyone been exposed, respectively [19]. Xie et al. showed that the adjusted KaplanMeier estimator is consistent even if weights came from the PS computed by a logistic regression model [20]. This is also true for stabilized inverse probability weighted estimates [21].
Adjusted Cox proportional hazard model
Cole et al. demonstrated that the stabilized inverse probability of treatment weighting (SIPTW) Cox regression model provides unbiased estimates, while robust variance estimation, such as those suggested by Lin and Wei, can be used to account for the weighting procedure. These results are also valid in the presence of timevarying confounding [18, 19, 22].
Aalen’s additive hazards model
Exploration of SIPTW weighted KaplanMeier survival curves is recommended in the presence of violation of the proportional hazard (PH) assumption [21]. In this case, Aalen’s additive hazards model represents a valid exploratory graphical method to detect and describe the nature of timevarying covariate effects [23]. The hazard function at time t for a model containing p+1 covariates is: h{t,X,β(t)} = β_{0}(t) + β_{1}(t)x_{1} + β_{2}(t)x_{2} + … β_{p}(t)x_{p}.
The coefficient β_{k}(t) provides the change in hazard at time t from the baseline hazard function β_{0}(t) for a oneunit change in the respective covariate x_{k}, holding all other covariates constant. A weighted version of the Aalen’s model was discussed by Huffer and McKeague [24].
The principal difference between Aalen’s and Cox models is that Aalen’s model allows the effect of the k^{th} covariate to change continuously over time and represents a valid method to describe the function. The Aalen estimator is based on the cumulative hazard function (obtained by integrating the hazard function at time t for a model containing p+1 covariates) where β_{k}(t) represents the cumulative regression coefficient for the k^{th} covariate over time. A plot of \( {\displaystyle {\widehat{\beta}}_k}(t) \) versus t with its 95 % confidence interval is a worthwhile method to identify and describe the possible interaction between the coefficient with time and also to evaluate its relevance.
Once a significant timevarying effect has been found for a given covariate, a simple approach would be to introduce it into the Cox proportional hazard model as a timevarying covariate and to estimate a specific hazard ratio for each particular timeinterval identified. An appropriate modeling of the timevarying covariates makes the Cox model suitable for satisfying the proportional hazard assumption.
Results
Study population characteristics and results of the role of the consolidation phase in the treatment of MM are described in this section.
Study design
The original study consisted of 480 patients enrolled between May 2006 and April 2008 from 73 Italian hospitals in a phase III openlabel randomized clinical trial (RCT). Eligible patients were aged 18–65 years and had previously untreated symptomatic and measurable multiple myeloma [25].
Enrolled patients were randomized (1:1 ratio) into two treatment arms: the experimental arm (Arm A) versus the standard arm (Arm B). Six patients withdrew consent immediately after randomization without starting therapy. In both arms, patients received induction therapy before undergoing up to two planned autologous stem cell transplantation (ASCT) procedures followed by a consolidation phase. The experimental/standard therapy was administered in both the induction and consolidation phases.
The study was approved by an independent ethics committee or by the institutional review board of each participating institution, and was performed in accordance with International Conference on Harmonisation guidelines on Good Clinical Practice and with the principles laid down in the Declaration of Helsinki. All patients provided written informed consent prior to enrollment.
Complete information on pretreatment baseline characteristics was available for 414 (87.34 %) of the 474 randomized patients (Table 1), which represents the core of the analysis.
At a median follow up of 54 months from randomization, 85/200 and 121/214 events/patient (progression, relapse, or death from any cause) had been recorded in the experimental and control arms, respectively.
Overall, 52 patients in arm A and 65 in arm B discontinued treatment before undergoing the consolidation phase, mainly because of toxicity in Arm A and disease progression in Arm B (Fig. 1a). To evaluate the specific role of the consolidation phase miming an ad hoc RCT (Fig. 1b) it was therefore necessary to consider the selection process in a differential way on the basis of baseline characteristics, treatment arm and clinical course. In doing so, we were able to restore a balanced comparison for the consolidation phase either conditionally or marginally by weighting for the conditional probability of receiving consolidation therapy.
Propensity score, inverse probability of treatment weight and stabilized inverse probability of treatment weight estimation
The consolidation phase was administered to patients who received at least one of the two planned ASCT procedures. Two selection processes were then modeled, the first by computing the conditional probability of receiving at least one ASCT and the second by computing the conditional probability of undergoing the consolidation phase.
The first PS – of undergoing at least one of the two planned ASCTs – was calculated on the basis of pretreatment baseline characteristics (age, sex, hemoglobin, platelets, creatinine, LDH, ISS stage, isotype of disease and cytogenetic abnormalities) and treatment arm. The second PS – of undergoing the consolidation phase – was calculated using pretreatment baseline characteristics, treatment arm and clinical response to the induction phase (categorized as complete or partial) (Table 2). We showed detailed models because there is still some debate about the choice of the proper covariates in the PS model. To note that, in our case, no variable in the logistic models resulted particularly important. Nevertheless, it may be possible to proceed with the reweighting procedure.
From each PS, the relative inverse probability of treatment weight (IPTW) and the stabilized inverse probability of treatment weight (SIPTW) were calculated as previously described.
Total stabilized inverse probability of treatment weight (TSIPTW) for each patient was obtained by multiplying the SIPTW of undergoing the 1^{st} ASCT by the SIPTW of undergoing the consolidation phase (Table 3). It may be noted that the stabilization procedure made every inverse probability weight more stable than their respective non stabilized weights.
Three hundred and sixtytwo (87.44 %) patients started the 1^{st} ASCT with a median PS (PS1) of 0.89 (range: 0.56–0.97; interquartile range (IQR): 0.84–0.93). The SIPTW of starting the 1^{st} ASCT was obtained as [0.8744/PS1] or [(10.8744)/(1 PS1)] in the event of not starting it. The median value was 0.97 (range: 0.29–2.62; IQR: 0.93–1.03).
Of the 362 who started the 1^{st} ASCT, 81.77 % also started the consolidation phase with a median PS (PS2) of 0.80 (range: 0.58–1.00; IQR: 0.75–0.88).
The SIPTW of starting consolidation therapy was obtained as [0.8177/PS2] or [(10.8177)/(1PS2)] in the event of not starting it. The median value was 1.00 (range: 0.43–5.79; IQR: 0.87–1.06).
Finally, the TSIPTW was obtained as the product of SIPTW of starting the 1^{st} ASCT by the SIPTW of starting consolidation therapy. The median value was 0.96 (range: 0.29–5.53; IQR: 0.86–1.07). A TSIPTW > 2 was observed in a minority of patients (n = 8) who had not started consolidation therapy. Five of these had not finished induction therapy due to serious toxicities (n = 3), progression (n = 1) or death (n = 1) and were removed from the study; the remaining 3, who had not reached complete clinical response (CR), dropped out of the study after the ASCT because of severe toxicity, progression or protocol violation. For those who did not start the 1^{st} ASCT (n = 52), the TSIPTW corresponded to the SIPTW of undergoing the 1^{st} ASCT because the SIPTW of starting consolidation therapy could not be computed, obviously.
Weighted KaplanMeier curves
After having estimated the TSIPTW to undergo the consolidation phase, 346 patients who could potentially receive consolidation therapy as they had undergone at least one ASCT were evaluated in terms of progressionfree survival (PFS), defined as the time elapsed from the last ASCT evaluation to the date of progression or death (events), or last follow up (censored).
The consolidation phase did not have an effect on PFS: a gradual decrease was seen in the survival curves of patients who had undergone consolidation therapy, while those who had not started this last phase of therapy rapidly relapsed till about 12 months from the last ASCT evaluation. After the first year, the two PFS survival curves became similar (after TSIPTW 1year PFS: 93 % vs 75 %, weighted Logrank test p = 0.9554) (Fig. 2). This happened because TSIPTWs gave more weight to patients who did not receive consolidation therapy and who were more likely to fail during the first year after ASCT.
Aalen’s additive hazards model to investigate timevarying effect
As the PH assumption did not hold (Schoenfeld’s test, p < 0.0001), Aalen’s additive hazards models were fit to investigate the timevarying effect of the consolidation phase. We fitted a model containing all the prognostic factors – i.e. age, sex, hemoglobin, platelets, creatinine, ISS stage, cytogenetic alterations – and treatment arm and an indicator for the consolidation phase (yes/no).
The plot of the cumulative regression coefficient for the last covariate (consolidation phase) over time showed a decreasing pattern up to the 18^{th} month and then a roughly zero slope. The upper confidence band excluded the null value in the first 18 months (Fig. 3). This plot suggests that the consolidation phase may have had an effect on early followup times up to 18 months, but no late effect. The consolidation phase appears to have lost its efficacy after one and half years.
Weighted Cox proportional hazard models with timedependent covariate
We performed a weighted Cox proportional hazard model analysis including the consolidation phase as covariate where all measured confounding factors were controlled by weighting. The consolidation phase was specified as an interaction term with the followup time. In particular, we specified 2 dummies to denote early and late benefit, one for the first 18 months and the other for the period >18 months, respectively. These dummies were set to zero for patients who did not receive the consolidation phase treatment. The first dummy (early effect of consolidation) was set at 1 for patients who received the consolidation treatment and followup time from t_{0} to t_{18} (or t_{end of fup}, if previous) and was set to 0 after t_{18}. The second dummy (late consolidation effect) was set at 1 for patients who received consolidation treatment and followup time from t_{18} to t_{end of fup}, and was set to 0 before t_{18}.
Consolidation therapy resulted in an early effect with weighted HR = 0.40 (robust 95 % confidence interval 0.17–0.96; p = 0.040), whereas it did not show any impact at later followup times (weighted HR = 1.98; robust 95 % confidence interval 0.93–4.20; p = 0.074). The proportional hazard assumption was satisfied (p = 0.9062).
Discussion
The novelty of this article was to use a propensity scorebased approach in an RCT context designed for another purpose to evaluate the effect of the last phase of a complex therapeutic strategy, adjusting for potential selection bias caused by different compliance to the previous phases of the therapeutic scheme. Differential compliance to earlier therapeutic phases of the protocol and the consequent failure of the randomization balance may have caused biased results.
Over the last few decades in different research settings, when the randomization procedure was absent by definition (e.g. observational studies) or was simply lost (e.g. RCTs with low compliance rate), propensity scorebased approaches were proposed and became popular. The aim was to recreate an artificial population in which treatment assignment could be ignored, and to compare outcomes in treated and untreated subjects, mimicking randomization when it was absent or when the RCT was designed with another purpose.
Using combined stabilized inverse probability treatment weights, we estimated that consolidation therapy in newly diagnosed multiple myeloma patients had an effect on PFS that was restricted to the first 18 months after the start of therapy.
We applied two sets of weights derived from two distinct PS equations to try to restore balance to receiving autologous stem cell transplantation (ASCT) and to ensure ignorability of the last treatment phase (the socalled consolidation phase). We introduced the PS into the model by inverse probability weighting (IPTW). Stabilization of weights was sufficient to gain robustness.
We used the same approach (IPTW) to calculate weighted KaplanMeier survival curves. Adjusted survival curves are very useful in exploratory and descriptive phases and comparing them to unweighted curves permits a simple inspection of the selection forces in a study. In particular, considering survival curves for patients receiving or not receiving the consolidation phase, a selection of poor prognosis patients was evident in the first followup period, after which weighted and unweighted curves become similar. The weighted analysis of the consolidation phase showed a timedependent benefit which was evident for a short time span after ASCT. Subsequently, any advantage disappeared. There was no evidence of interaction between consolidation and randomization arms and the sample size was too small for subgroup analysis.
In a setting such as that of our case study, where the proportional hazards assumption was violated, robustness of weighting procedures is critical. When the effects are timevarying, the gap between survival curves may depend on few risk sets and the weights attributed to a small number of subjects in those risk sets becomes very influential. An exploratory and descriptive solution is to investigate the presence of the timevarying effect through Aalen’s additive hazards models. A weighted version of this model was previously evaluated by Huffer and McKeague [24]. We included timevarying variables in the weighted Cox models as a set of dichotomous covariates after defining time intervals of different lengths based on Aalen’s analysis, so that proportional hazard assumption be satisfied.
Conclusion
Estimating stabilized inverse probability weights by PS logistic models, combining them to consider the sequence of therapeutic phases, carrying out sensitivity analysis to evaluate the proportional hazard assumption and timevarying effects, fitting Aalen’s additive hazard models to investigate timevarying effects and, finally, estimating a weighted Cox proportional hazard model with timevarying covariates undoubtedly appears a complex strategy. However, it is a more appropriate one to use in the event of lack of compliance to a series of treatment phases and when there is a timevarying effect. Both situations are common in cancer trials.
Notably, PS methods can adjust for imbalance caused by observable covariates. Unobservable covariates can still exert an effect and distort estimates. Thus, analysis using PS methods cannot be conclusive in this context. Ad hoc designed RCTs with randomization procedures planned after ASCT are important to understand what kind of patient really benefits from consolidation treatment.
Abbreviations
 ASCT:

Autologous stem cell transplantation
 CI:

Confidence interval
 CR:

Complete response
 HR:

Hazard ratio
 IPTW:

Inverse probability treatment weight
 IQR:

Interquartile range
 ITT:

Intention to treat
 MM:

Multiple myeloma
 Mos:

Months
 PD:

Progression disease
 PFS:

Progression free survival
 PH:

Proportional hazard
 PP:

Perprotocol
 PS:

Propensity score
 RCT:

Randomized clinical trial
 SIPTW:

Stabilized inverse probability treatment weight
 TSIPTW:

Total stabilized inverse probability treatment weight
References
 1.
Armitage P. The role of randomization in clinical trials. Stat Med. 1982;4:345–52.
 2.
Hernan MA, HernandezDiaz S, Robins JM. A structural approach to selection bias. Epidemiology. 2004;15(5):615–25.
 3.
Frangakis CE, Rubin DB. Addressing complications of intentionto treat analysis in the combined presence of allornone treatmentnoncompliance and subsequent missing outcomes. Biometrika. 1999;86:365–79.
 4.
Mealli F, Rubin D. Assumptions when Analyzing Randomized Experiments with Noncompliance and missing Outcomes. Health Serv Outcome Res Methodol. 2002;3:225–32.
 5.
Mealli F, Imbens GW, Ferro S, Biggeri A. Analyzing a randomized trial on breast selfexamination with noncompliance and missing outcomes. Biostatistics. 2004;5:2,207–22.
 6.
Rosenbaum P, Rubin D. The central role of the propensity score in observational studies for causal effects. Biometrika. 1983;70:41–55.
 7.
Cavo M, Brioli A, Tacchetti P, Zannetti BA, Mancuso K, Zamagni E. Role of consolidation therapy in transplant eligible multiple myeloma patients. Semin Oncol. 2013;40(5):610–7.
 8.
D’Agostino RB. Tutorial in biostatistics. Propensity score methods for bias reduction in the comparison of a treatment to a nonrandomized control group. Statist Med. 1998;17:2265–81.
 9.
Becker SO, Ichino A. Estimation of average treatment effects based on propensity scores. Stata J. 2002;2(4):358–77.
 10.
Brookhart MA, Schneeweiss S, Rothman KJ, Glynn RJ, Avorn J, Sturmer T. Variable selection for propensity score models. Am J Epidemiol. 2006;163(12):1149–56.
 11.
Pearl J. Letter to the editor. Remarks on the method of propensity score. Statist Med. 2009;28:1415–24.
 12.
Guo SY, Fraser MW. Propensity Score Analysis: Statistical Methods and Applications, Second Edition. Sage; 2015. https://0uksagepubcom.brum.beds.ac.uk/engb/eur/propensityscoreanalysis/book238151#description.
 13.
Austin P. The performance of different propensity score methods for estimating marginal hazard ratios. Statist Med. 2013;32(16):2837–49.
 14.
Austin P. The use of propensity score methods with survival or timetoevent outcomes: reporting measures of effect similar to those used in randomized experiments. Statist Med. 2014;33(7):1242–58.
 15.
Xu S, Ross C, Raebel MA, Shetterly S, Blanchette C, Smith D. Use of stabilized inverse propensity score as weight to directly estimate relative risk and its confidence intervals. Value Health. 2010;2(13):273–7.
 16.
Hernan MA, Robins J. Estimating causal effects from epidemiological data. J Epidemiol Community Health. 2006;60:578–86.
 17.
Cole SR, Hernan MA. Constructing Inverse Probability Weights for Marginal Structural models. Am J Epidemiol. 2008;168(6):656–64.
 18.
Westreich D, Cole SR, Schisterman EF, Platt RW. A simulation study of finitesample properties of marginal structural Cox proportional hazards models. Statist Med. 2012;31(19):2098–109.
 19.
Cole SR, Hernan MA. Adjusted survival curves with inverse probability weights. Comput Methods Prog Biomed. 2004;75:45–9.
 20.
Xie J, Liu C. Adjusted KaplanMeier Estimator and Logrank Test with Inverse Probability of Treatment Weighting for Survival Data. Statist. Med. 2005;24(20):3089–110.
 21.
Xu S, Shetterly S, Powers D, Raebel MA, Tsai T, Ho PM, Magic D. Extension of KaplanMeier methods in observational studies with timevarying treatment. Value Health. 2012;15(1):167–74.
 22.
Lin DY, Wei LJ. The robust inference for the proportional hazards model. J Am Stat Assoc. 1989;84:1074–8.
 23.
Hosmer DW, Royston P. Using Aalen’s linear hazards model to investigate timevarying effects in the proportional hazards regression model. Stata J. 2002;2(4):331–50.
 24.
Huffer FW, McKeague IW. Weighted least squares estimation of Aalen’s additive risk model. J Am Stat Assoc. 1991;86:114–29.
 25.
Cavo M, Tacchetti P, Patriarca F, Petrucci MT, Pantani L, Galli M, et al. Bortezomib with thalidomide plus dexamethasone compared with thalidomide plus dexamethasone as induction therapy before, and consolidation therapy after, double autologous stemcell transplantation in newly diagnosed multiple myeloma: a randomised phase 3 study. Lancet. 2010;376(9758):2075–85. Erratum in: Lancet. 2011;378(9806):1846.
Acknowledgements
Not applicable.
Funding
No specific funding.
Availability of data and materials
The datasets generated during and/or analyzed during the current study are not publicly available due to legal restrictions in Italy regarding patientlevel data sharing.
Authors’ contributions
AP performed the analysis, interpreted results, and wrote the first and final drafts of the paper. ON and AB oversaw the analysis, interpreted results, and critically revised the manuscript. MC and EZ critically revised the manuscript. All authors read and approved the final manuscript.
Competing interests
The authors declare that they have no competing interests.
Consent for publication
Not applicable.
Ethics approval and consent to participate
The study was approved by an independent ethics committee or by the institutional review board of each participating institution, and was performed in accordance with International Conference on Harmonisation guidelines on Good Clinical Practice and with the principles laid down in the Declaration of Helsinki. All patients provided written informed consent prior to enrollment.
Author information
Affiliations
Corresponding author
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Pezzi, A., Cavo, M., Biggeri, A. et al. Inverse probability weighting to estimate causal effect of a singular phase in a multiphase randomized clinical trial for multiple myeloma. BMC Med Res Methodol 16, 150 (2016). https://0doiorg.brum.beds.ac.uk/10.1186/s1287401602539
Received:
Accepted:
Published:
DOI: https://0doiorg.brum.beds.ac.uk/10.1186/s1287401602539
Keywords
 Selection bias
 Weighting sample
 Propensity score
 Causal effect
 RCT
 Compliance