Simulation study of instrumental variable approaches with an application to a study of the antidiabetic effect of bezafibrate

Abstract

Purpose: We studied the application of the generalized structural mean model (GSMM) of instrumental variable (IV) methods in estimating treatment odds ratios (ORs) for binary outcomes in pharmacoepidemiologic studies and evaluated the bias of GSMM compared to other IV methods.

Methods: Because of the bias of standard IV methods, including two-stage predictor substitution (2SPS) and two-stage residual inclusion (2SRI) with binary outcomes, we implemented another IV approach based on the GSMM of Vansteelandt and Goetghebeur. We performed simulations under the principal stratification setting and evaluated whether GSMM provides approximately unbiased estimates of the causal OR and compared its bias and mean squared error to that of 2SPS and 2SRI. We then applied different IV methods to a study comparing bezafibrate versus other fibrates on the risk of diabetes.

Results: Our simulations showed that unlike the standard logistic, 2SPS, and 2SRI procedures, our implementation of GSMM provides an approximately unbiased estimate of the causal OR even under unmeasured confounding. However, for the effect of bezafibrate versus other fibrates on the risk of diabetes, the GSMM and two-stage approaches yielded similarly attenuated and statistically non-significant OR estimates. The attenuation of the OR by the two-stage and GSMM IV approaches suggests unmeasured confounding, although violations of the IV assumptions or differences in the parameters estimated could be playing a role.

Conclusion: The GSMM IV approach provides approximately unbiased adjustment for unmeasured confounding on binary outcomes when a valid IV is available.