Standard statistical methods were used for descriptive analyses, taking into account cohort sampling based on twin pairs. As twins within a pair are not statistically independent observations, robust estimates of standard errors were obtained for individual-based analyses (Williams, 2000).