Statistical approaches for estimating sex-specific effects in endocrine disruptors research

Abstract

BACKGROUND: When a biologic mechanism of interest is anticipated to operate differentially according to sex, as is often the case in endocrine disruptors research, investigators routinely estimate sex-specific associations. Less attention has been given to potential sexual heterogeneity of confounder associations with outcomes. When relationships of covariates with outcomes differ according to sex, commonly applied statistical approaches for estimating sex-specific endocrine disruptor effects may produce divergent estimates. OBJECTIVES: We discuss underlying assumptions and evaluate the performance of two traditional approaches for estimating sex-specific effects, stratification and product terms, and introduce a simple modeling alternative: an augmented product term approach. METHODS: We describe the impact of assumptions regarding sexual heterogeneity of confounder relationships on estimates of sex-specific effects of the exposure of interest for three approaches: stratification, traditional product terms, and augmented product terms. Using simulated and applied examples, we demonstrate properties of each approach under a range of scenarios. RESULTS: In simulations, sex-specific exposure effects estimated using the traditional product term approach were biased when confounders had sex-dependent associations with the outcome. Sex-specific estimates from stratification and the augmented product term approach were unbiased but less precise. In the applied example, the three approaches yielded similar estimates, but resulted in some meaningful differences in conclusions based on statistical significance. CONCLUSIONS: Investigators should consider sexual heterogeneity of confounder associations when choosing an analytic approach to estimate sex-specific effects of endocrine disruptors on health. In the presence of sex-dependent confounding, our augmented product term approach may be advantageous over stratification when there is prior knowledge available to fit reduced models or when investigators seek an automated test for effect measure modification.