Statistical associations between the continuous and categorical predictors and continuous and binary outcome variables were estimated using Generalized Linear Mixed Effects Models (GLMM) in the PROC GLIMMIX procedure of SAS (SAS Institute, Inc., Cary, NC). This procedure allowed us to account for the dependent nature of the twin observations (all twins were included in the analyses) similarly to what has been shown previously (49).
