1
|
Matsouaka RA. Instrumental variable estimation of causal odds ratios using structural nested mean models. Biostatistics 2017; 18:465-476. [PMID: 28334061 PMCID: PMC5862265 DOI: 10.1093/biostatistics/kxw059] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2015] [Revised: 12/03/2016] [Accepted: 12/15/2016] [Indexed: 11/14/2022] Open
Abstract
We consider estimating causal odds ratios using an instrumental variable under a logistic structural nested mean model (LSNMM). Current methods for LSNMMs either rely heavily on possible "uncongenial" modeling assumptions or involve intricate numerical challenges, which have impeded their use. In this article, we present an alternative method that ensures a congenial parametrization, circumvents computational complexity of existing methods, and is easy to implement. We illustrate the proposed method to (1) estimate the causal effect of years of education on earnings using data from the NLSYM and (2) assess the impact of moving families from high to low-poverty neighborhoods had on lifetime major depressive disorder among adolescents in the "Moving to Opportunity (MTO) for Fair Housing Demonstration Project" from the Department of Housing and Urban Development.
Collapse
Affiliation(s)
- Roland A. Matsouaka
- Department of Biostatistics and Bioinformatics & Duke Clinical Research Institute Duke University, Durham, NC 27705, USA
| |
Collapse
|
2
|
Abstract
Background Randomized trials evaluating new cancer screening technologies may underestimate the efficacy of screening to reduce cancer mortality if study participants are noncompliant. Participants may fail to comply with the screening itself or fail to obtain appropriate diagnostic follow-up and treatment. Noncompliance with screening has drawn wide attention, but little attention has been paid to noncompliance with diagnostic follow-up and treatment. Purpose To examine the importance of noncompliance with screening, follow-up, and treatment in cancer screening trials. Methods The unique problems associated with noncompliance in screening trials are described and provide an example illustrating the potential impact of noncompliance in a screening trial. I discuss issues that arise with measurement of follow-up and therapeutic noncompliance, and the benefit of collecting information on health system and participant characteristics associated with noncompliance. Results The estimate of the efficacy of a screening program on cancer mortality can be adjusted for screening, follow-up, and treatment noncompliance. Noncompliance needs to be measured in a rigorous, systematic manner across all arms of the trial. Information on health system and participant characteristics associated with compliance may also be incorporated into statistical models to estimate screening effects with full compliance, plan interventions to increase compliance, and extrapolate results of screening trials from one population to another. Limitations Measuring compliance with follow-up and treatment can be difficult when these occur outside the trial, and when there is variation among providers in follow-up and treatment practices. Conclusions Noncompliance may alter the estimate of a screening effect on cancer mortality in clinical trials. It is possible to adjust screening efficacy estimates for noncompliance using existing statistical techniques. It is important that data describing compliance with screening, follow-up, and treatment are collected as part of standard data collection in cancer screening trials. Clinical Trials 2007; 4: 341—349. http://ctj.sagepub.com
Collapse
Affiliation(s)
- Ilana F Gareen
- Center for Statistical Sciences and the Department of Community Health, Brown University School of Medicine, Providence, RI 02912, USA.
| |
Collapse
|
3
|
Clarke PS, Palmer TM, Windmeijer F. Estimating Structural Mean Models with Multiple Instrumental Variables Using the Generalised Method of Moments. Stat Sci 2015. [DOI: 10.1214/14-sts503] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
4
|
Taguri M, Matsuyama Y, Ohashi Y. Model selection criterion for causal parameters in structural mean models based on a quasi-likelihood. Biometrics 2014; 70:721-30. [PMID: 24621405 DOI: 10.1111/biom.12165] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2012] [Revised: 01/01/2014] [Accepted: 01/01/2014] [Indexed: 11/27/2022]
Abstract
Structural mean models (SMMs) have been proposed for estimating causal parameters in the presence of non-ignorable non-compliance in clinical trials. To obtain a valid causal estimate, we must impose several assumptions. One of these is the correct specification of the structural model. Building on Pan's work (2001, Biometrics 57, 120-125) on developing a model selection criterion for generalized estimating equations, we propose a new approach for model selection of SMMs based on a quasi-likelihood. We provide a formal model selection criterion that is an extension of Akaike's information criterion. Using subset selection of baseline covariates, our method allows us to understand whether the treatment effect varies across the available baseline covariate levels, and/or to quantify the treatment effect on a specific covariates level to target specific individuals to maximize treatment benefit. We present simulation results in which our method performs reasonably well compared to other testing methods in terms of both the probability of selecting the correct model and the predictive performances of the individual treatment effects. We use a large randomized clinical trial of pravastatin as a motivation.
Collapse
Affiliation(s)
- Masataka Taguri
- Department of Biostatistics and Epidemiology, Graduate School of Medicine, Yokohama City University, 3-9 Fukuura, Kanazawa-ku, Yokohama, Kanagawa 236-0004, Japan
| | - Yutaka Matsuyama
- Department of Biostatistics, School of Public Health, Graduate School of Medicine, University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan
| | - Yasuo Ohashi
- Department of Biostatistics, School of Public Health, Graduate School of Medicine, University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan
| |
Collapse
|
5
|
Brumback BA, He Z, Prasad M, Freeman MC, Rheingans R. Using structural-nested models to estimate the effect of cluster-level adherence on individual-level outcomes with a three-armed cluster-randomized trial. Stat Med 2013; 33:1490-502. [PMID: 24288357 DOI: 10.1002/sim.6049] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2013] [Revised: 10/10/2013] [Accepted: 10/29/2013] [Indexed: 11/07/2022]
Abstract
Much attention has been paid to estimating the causal effect of adherence to a randomized protocol using instrumental variables to adjust for unmeasured confounding. Researchers tend to use the instrumental variable within one of the three main frameworks: regression with an endogenous variable, principal stratification, or structural-nested modeling. We found in our literature review that even in simple settings, causal interpretations of analyses with endogenous regressors can be ambiguous or rely on a strong assumption that can be difficult to interpret. Principal stratification and structural-nested modeling are alternative frameworks that render unambiguous causal interpretations based on assumptions that are, arguably, easier to interpret. Our interest stems from a wish to estimate the effect of cluster-level adherence on individual-level binary outcomes with a three-armed cluster-randomized trial and polytomous adherence. Principal stratification approaches to this problem are quite challenging because of the sheer number of principal strata involved. Therefore, we developed a structural-nested modeling approach and, in the process, extended the methodology to accommodate cluster-randomized trials with unequal probability of selecting individuals. Furthermore, we developed a method to implement the approach with relatively simple programming. The approach works quite well, but when the structural-nested model does not fit the data, there is no solution to the estimating equation. We investigate the performance of the approach using simulated data, and we also use the approach to estimate the effect on pupil absence of school-level adherence to a randomized water, sanitation, and hygiene intervention in western Kenya.
Collapse
Affiliation(s)
- Babette A Brumback
- Department of Biostatistics, University of Florida, Gainesville, FL 32611, U.S.A
| | | | | | | | | |
Collapse
|
6
|
|
7
|
Vansteelandt S, Bowden J, Babanezhad M, Goetghebeur E. On Instrumental Variables Estimation of Causal Odds Ratios. Stat Sci 2011. [DOI: 10.1214/11-sts360] [Citation(s) in RCA: 80] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
8
|
Abstract
Background Departures from randomized treatments
complicate the analysis of many randomized controlled trials. Intention-to-treat
analysis estimates the effect of being allocated to treatment.
It is now possible to estimate the effect of receiving
treatment without assuming comparability of groups defined by actual treatment.
However, the methodology is largely confined to trials where the only treatment
changes were switches to other trial treatments. Purpose To propose a method for comparing the
effects of receiving trial treatments in an active-controlled clinical trial
where some participants received nontrial treatments and others received no
treatment at all, and to illustrate the method in the PENTA 5 trial in
HIV-infected children. Methods We combine the instrumental variables
approach, which forms unbiased estimating equations based on the randomization
but does not fully identify the contrasts of trial treatment effects, with prior
information about the distribution of possible effects of nontrial treatments
and of one trial treatment; we do not need to use prior information about the
comparisons of trial treatments. Prior information in PENTA 5 was elicited from
the investigators. Results In PENTA 5, the prior information suggested
that all treatments were beneficial, but there was uncertainty about the degree
of benefit. Allowing for this prior information changed point estimates and
increased standard errors compared with ignoring nontrial treatments. Limitations The method depends on the correct
specification of the causal effect of treatment: in PENTA 5, this assumed a
linear effect of dose and no interactions between treatments. This specification
is hard to check from the data but can be explored in sensitivity analyses.
Prior information would be better derived from the literature whenever
possible. Conclusions The use of partial prior information
gives a way to adjust for complex patterns of departures from randomized
treatments. It should be useful in all trials where nontrial treatments are used
and in active-controlled trials where trial treatments are not universally
used.
Collapse
Affiliation(s)
- Simon J Bond
- MRC Biostatistics Unit, Institute of Public Health, Cambridge, UK
| | | |
Collapse
|
9
|
Goetghebeur E. Commentary: To cause or not to cause confusion vs transparency with Mendelian Randomization. Int J Epidemiol 2010; 39:918-20. [DOI: 10.1093/ije/dyq100] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
10
|
Comté L, Vansteelandt S, Tousset E, Baxter G, Vrijens B. Linear and loglinear structural mean models to evaluate the benefits of an on-demand dosing regimen. Clin Trials 2009; 6:403-15. [DOI: 10.1177/1740774509344777] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
Background Structural mean models can be used to estimate treatment efficacy when drug exposure varies. We applied stuctural mean model to evaluate the clinical benefits of a proton pump inhibitor prescribed to be taken as needed to alleviate epigastric pain. We also investigated a new diagnostic approach to evaluate model assumptions. Methods All patients were suffering from nonerosive reflux disease or functional ulcer-like dyspepsia and were prescribed a proton pump inhibitor to be taken as needed for relief of epigastric pain. The primary endpoint was a score variable that expresses the magnitude of gastro-intestinal symptoms at 8 weeks after randomization. We developed linear and loglinear versions of the structural mean models to derive an unbiased estimator of the reduction in symptom score as a function of exposure to the test drug. Semi-parametric models based on splines and corresponding simultaneous confidence bands identified the presence of potential interactions between drug exposure and baseline covariates. Results The on-demand dosing regimen generated a wide range of drug exposure. Application of SMM showed that the potential treatment-induced reduction in symptom score was much greater than the average treatment reduction observed in this population of patients. Our diagnostic tool was useful for detecting the interaction between drug exposure and baseline covariates. Limitations Analysis could only be performed over the first 2 months after randomization because, afterwards, many patients dropped out from the placebo group. Conclusions The structural mean model approach allows one to estimate treatment efficacy in the presence of variable drug exposure. Similar results were obtained using linear and loglinear structural mean model. Clinical Trials 2009; 6: 403—415. http://ctj.sagepub.com
Collapse
Affiliation(s)
- Laetitia Comté
- Department of Mathematics, University of Liège, Grande Traverse 12, B37, Liège 4000, Belgium
| | - Stijn Vansteelandt
- Department of Applied Mathematics and Computer Science, Ghent University, Ghent, Belgium
| | | | | | | |
Collapse
|
11
|
Sjölander A, Humphreys K, Vansteelandt S, Bellocco R, Palmgren J. Sensitivity analysis for principal stratum direct effects, with an application to a study of physical activity and coronary heart disease. Biometrics 2009; 65:514-20. [PMID: 18759834 DOI: 10.1111/j.1541-0420.2008.01108.x] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
SUMMARY In many studies, the aim is to learn about the direct exposure effect, that is, the effect not mediated through an intermediate variable. For example, in circulation disease studies it may be of interest to assess whether a suitable level of physical activity can prevent disease, even if it fails to prevent obesity. It is well known that stratification on the intermediate may introduce a so-called posttreatment selection bias. To handle this problem, we use the framework of principal stratification (Frangakis and Rubin, 2002, Biometrics 58, 21-29) to define a causally relevant estimand--the principal stratum direct effect (PSDE). The PSDE is not identified in our setting. We propose a method of sensitivity analysis that yields a range of plausible values for the causal estimand. We compare our work to similar methods proposed in the literature for handling the related problem of "truncation by death."
Collapse
Affiliation(s)
- Arvid Sjölander
- Department of Medical Epidemiology and Biostatistics, Karolinska Institutet, Stockholm, Sweden.
| | | | | | | | | |
Collapse
|
12
|
Cain LE, Cole SR, Greenland S, Brown TT, Chmiel JS, Kingsley L, Detels R. Effect of highly active antiretroviral therapy on incident AIDS using calendar period as an instrumental variable. Am J Epidemiol 2009; 169:1124-32. [PMID: 19318615 DOI: 10.1093/aje/kwp002] [Citation(s) in RCA: 27] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
Human immunodeficiency virus (HIV) researchers often use calendar periods as an imperfect proxy for highly active antiretroviral therapy (HAART) when estimating the effect of HAART on HIV disease progression. The authors report on 614 HIV-positive homosexual men followed from 1984 to 2007 in 4 US cities. During 5,321 person-years, 268 of 614 men incurred acquired immunodeficiency syndrome, 49 died, and 90 were lost to follow-up. Comparing the pre-HAART calendar period (<1996) with the HAART calendar period (>or=1996) resulted in a naive rate ratio of 3.62 (95% confidence limits: 2.67, 4.92). However, this estimate is likely biased because of misclassification of HAART use by calendar period. Simple calendar period approaches may circumvent confounding by indication at the cost of inducing exposure misclassification. To correct this misclassification, the authors propose an instrumental-variable estimator analogous to ones previously used for noncompliance corrections in randomized clinical trials. When the pre-HAART calendar period was compared with the HAART calendar period, the instrumental-variable rate ratio was 5.02 (95% confidence limits: 3.45, 7.31), 39% higher than the naive result. Weighting by the inverse probability of calendar period given age at seroconversion, race/ethnicity, and time since seroconversion did not appreciably alter the results. These methods may help resolve discrepancies between observational and randomized evidence.
Collapse
Affiliation(s)
- Lauren E Cain
- Department of Epidemiology, Bloomberg School of Public Health, Johns Hopkins University, Baltimore, Maryland, USA.
| | | | | | | | | | | | | |
Collapse
|
13
|
Cain LE, Cole SR. Inverse probability-of-censoring weights for the correction of time-varying noncompliance in the effect of randomized highly active antiretroviral therapy on incident AIDS or death. Stat Med 2009; 28:1725-38. [DOI: 10.1002/sim.3585] [Citation(s) in RCA: 68] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
|
14
|
Chen H, Geng Z, Zhou XH. Identifiability and estimation of causal effects in randomized trials with noncompliance and completely nonignorable missing data. Biometrics 2008; 65:675-82. [PMID: 18759847 DOI: 10.1111/j.1541-0420.2008.01120.x] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
In this article, we first study parameter identifiability in randomized clinical trials with noncompliance and missing outcomes. We show that under certain conditions the parameters of interest are identifiable even under different types of completely nonignorable missing data: that is, the missing mechanism depends on the outcome. We then derive their maximum likelihood and moment estimators and evaluate their finite-sample properties in simulation studies in terms of bias, efficiency, and robustness. Our sensitivity analysis shows that the assumed nonignorable missing-data model has an important impact on the estimated complier average causal effect (CACE) parameter. Our new method provides some new and useful alternative nonignorable missing-data models over the existing latent ignorable model, which guarantees parameter identifiability, for estimating the CACE in a randomized clinical trial with noncompliance and missing data.
Collapse
Affiliation(s)
- Hua Chen
- School of Mathematical Sciences, Peking University, Beijing 100871, China
| | | | | |
Collapse
|
15
|
Abstract
PURPOSE We review and compare two causal modeling approaches that correspond to two major and distinct classes of inference - efficacy and intervention-based inference - in the context of randomized trials with subject noncompliance. METHODS We review the definitions of efficacy and intervention-based effects in the clinical trials literature and relate these to two separate and distinct causal modeling approaches: the structural mean modeling (SMM) approach and the principal stratification, instrumental variable approach. RESULTS The SMM-based efficacy approach focuses on the effect of actually receiving treatment. In contrast, the principal stratification method addresses the effect of treatment assignment within partially unobserved latent subgroups defined by compliance behavior. While these approaches differ in terms of philosophy, model definitions, and estimation, they estimate the same causal effect under certain assumptions, but estimate very different causal effects when those assumptions are relaxed. We illustrate these results using a randomized psychiatry trial where the focus is physician compliance to the designated protocol and the other examines patient compliance to the designated protocol, both from the same trial. LIMITATIONS The validity of the models under the instrumental variable, SMM and principal stratification approaches depends on modeling assumptions, some of which may not be verifiable from the observed data and potentially less realistic than the no-confounding assumption made by non-causal approaches. CONCLUSIONS This comparison in terms of efficacy versus intervention-based effects in causal modeling parallels the explanatory versus pragmatic approaches in clinical trials research; therefore researchers should weigh carefully when choosing causal modeling methodology based on whether efficacy or intervention-based effects are of interest.
Collapse
Affiliation(s)
- Scarlett L Bellamy
- Department of Biostatistics and Epidemiology, University of Pennsylvania, School of Medicine, Philadelphia, PA 19104-6021, USA.
| | | | | |
Collapse
|
16
|
Goetghebeur E, Stijn V. Structural mean models for compliance analysis in randomized clinical trials and the impact of errors on measures of exposure. Stat Methods Med Res 2006; 14:397-415. [PMID: 16178139 DOI: 10.1191/0962280205sm407oa] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Abstract
Partial compliance with assigned treatment regimes is common in drug trials and calls for a causal analysis of the effect of treatment actually received. As such observed exposure is no longer randomized, selection bias must be carefully accounted for. The framework of potential outcomes allows this by defining a subject-specific treatment-free reference outcome, which may be latent and is modelled in relation to the observed (treated) data. Causal parameters enter these structural models explicitly. In this paper we review recent progress in randomization-based inference for structural mean modelling, from the additive linear model to the structural generalized linear models. An arsenal of tools currently available for standard association regression has steadily been developed in the structural setting, providing many parallel features to help randomization-based inference. We argue that measurement error on exposure is an important practical complication that has, however, not yet been addressed. We show how standard additive linear structural mean models are robust against unbiased measurement error and how efficient, asymptotically unbiased inference can be drawn when the degree of measurement error bias is known. The impact of measurement error is illustrated in a blood pressure example and finite sample properties are verified by simulation. We end with a plea for more and careful use of this methodology and point to directions for further development.
Collapse
|