“…We fit four quantile g-computation models with SRS total raw score as the outcome and the 6 phthalate metabolites (i.e., ΣDEHP, MBP, MBzP, MCPP, MEP, and MiBP) as the independent variables. 39,40 Model 1 was a joint marginal structural model given by falsedouble-struckE(YXq|Z,ψ,η)=g(ψ0+ψ1Sq+ηZ) where falseY was SRS total raw score, falseXq was the phthalates decile, falseboldZ were the covariates, falseg(⋅) was a link function in a generalized linear model (e.g., the identity function in the case of a linear regression model for the expected value of falseY), falseψ0 was the model intercept, falseψ1 was the expected change in the outcome, given a one quantile increase in all exposures simultaneously, falseSq is an “index” that represented a joint value of the phthalates, and falseη was a set of model coefficients for the covariates …”