the cardiovascular health study ( chs ) is a national institutes of health national heart , lung , and blood institute  funded , population - based longitudinal study of adults aged 65 years ( 7 ) . ( a full list of principal chs investigators and institutions can be found at http://www.chs-nhlbi.org/pi.htm . ) the main objective of the study was to study the onset and course of coronary heart disease and stroke . communities : forsyth county , north carolina ; sacramento county , california ; washington county , maryland ; and pittsburgh , pennsylvania . subsequent biannual in - person or telephone contacts were used to ascertain and verify the incidence of outcome events . for our analysis , we excluded 78 women with missing information on parity , 59 with missing information on baseline diabetes status , and 45 with missing information on other covariates , leaving a final eligible sample of 3,211 women . our main exposure of interest was parity , assessed in women at the baseline interview , with the question ,  how many live births have you had ? the response was modeled as a categorical variable : nulliparous ( 0 births ) , 12 livebirths , 34 livebirths , and grand multiparity ( 5 livebirths ) . information on medication use was ascertained annually with detailed medication inventories ( 8) . fasting serum blood glucose and insulin were measured at the central laboratory at the university of vermont in 19891990 , 19921993 , and 19961997 ( 9 ) . in this cohort of older women , we studied both the cross - sectional relationship between parity and prevalent diabetes at baseline and the association of parity with incident diabetes only among those who were nondiabetic at baseline ( n = 2,761 ) . baseline diabetes was assessed using data from the baseline clinic examination and defined with american diabetes association criteria as the use of hypoglycemic medications or a fasting blood glucose of 126 mg / dl . incident diabetes was defined as the first use of diabetes medications or by a criterion of fasting blood glucose of 126 mg / dl during the follow - up years for which fasting blood glucose was available . insulin resistance was estimated at baseline using homa - ir , calculated by the following formula : homa - ir = fasting serum insulin ( microunits per milliliter )  fasting plasma glucose ( millimoles per liter)/22.5 . we used descriptive statistics to characterize participant 's demographic characteristics at baseline . to examine the relationship between parity and prevalent diabetes , we estimated the prevalence ratios using generalized linear models with a log - link and poisson distribution ( 10 ) . models using a binomial link yielded similar point estimates and confidence intervals but did not converge in all cases . to examine the association between parity and incident diabetes in those without diabetes at baseline , we used cox proportional hazard models to estimate hazard risk ratios related to parity . for both the incident and prevalent diabetes outcomes , we examined the relationship between parity and outcomes in multivariable sequentially adjusted models . the next models controlled for several possible confounders including age , race , income ( three categories ) , education ( more than high school , high school , or less than high school ) , marital status ( married , single , widowed , or divorced ) , height ( centimeters ) smoking ( current , never , or former ) , study site , and alcohol intake ( user / abstainer ) . our final models contained all of the previously listed potential confounders and in addition controlled for body anthropometrics as potential mediators ( baseline bmi [ calculated from measured weight and height ] , height and measured waist circumference , and self - report of bmi at age 50 years ; because of missing data , age 50 years bmi was imputed from age , race , and baseline bmi in 128 women ) . in these final models , we examined the degree of mediation using the sas mediate macro ( 11 ) . we also conducted multivariable linear regression analyses to examine whether parity was associated with baseline fasting serum glucose , insulin levels , and homa - ir only among those without diabetes at baseline . because insulin and homa - ir were , as expected , right - skewed , both were log - transformed to improve normality . our main exposure of interest was parity , assessed in women at the baseline interview , with the question ,  how many live births have you had ? the response was modeled as a categorical variable : nulliparous ( 0 births ) , 12 livebirths , 34 livebirths , and grand multiparity ( 5 livebirths ) . information on medication use was ascertained annually with detailed medication inventories ( 8) . fasting serum blood glucose and insulin were measured at the central laboratory at the university of vermont in 19891990 , 19921993 , and 19961997 ( 9 ) . in this cohort of older women , we studied both the cross - sectional relationship between parity and prevalent diabetes at baseline and the association of parity with incident diabetes only among those who were nondiabetic at baseline ( n = 2,761 ) . baseline diabetes was assessed using data from the baseline clinic examination and defined with american diabetes association criteria as the use of hypoglycemic medications or a fasting blood glucose of 126 mg / dl . incident diabetes was defined as the first use of diabetes medications or by a criterion of fasting blood glucose of 126 mg / dl during the follow - up years for which fasting blood glucose was available . insulin resistance was estimated at baseline using homa - ir , calculated by the following formula : homa - ir = fasting serum insulin ( microunits per milliliter )  fasting plasma glucose ( millimoles per liter)/22.5 . we used descriptive statistics to characterize participant 's demographic characteristics at baseline . to examine the relationship between parity and prevalent diabetes , we estimated the prevalence ratios using generalized linear models with a log - link and poisson distribution ( 10 ) . models using a binomial link yielded similar point estimates and confidence intervals but did not converge in all cases . to examine the association between parity and incident diabetes in those without diabetes at baseline , we used cox proportional hazard models to estimate hazard risk ratios related to parity . for both the incident and prevalent diabetes outcomes , we examined the relationship between parity and outcomes in multivariable sequentially adjusted models . the next models controlled for several possible confounders including age , race , income ( three categories ) , education ( more than high school , high school , or less than high school ) , marital status ( married , single , widowed , or divorced ) , height ( centimeters ) smoking ( current , never , or former ) , study site , and alcohol intake ( user / abstainer ) . our final models contained all of the previously listed potential confounders and in addition controlled for body anthropometrics as potential mediators ( baseline bmi [ calculated from measured weight and height ] , height and measured waist circumference , and self - report of bmi at age 50 years ; because of missing data , age 50 years bmi was imputed from age , race , and baseline bmi in 128 women ) . in these final models , we examined the degree of mediation using the sas mediate macro ( 11 ) . we also conducted multivariable linear regression analyses to examine whether parity was associated with baseline fasting serum glucose , insulin levels , and homa - ir only among those without diabetes at baseline . because insulin and homa - ir were , as expected , right - skewed , both were log - transformed to improve normality . those with grand multiparity ( 5 livebirths ) were more likely to be black , have larger bmi and waist circumferences at baseline , more likely to abstain from alcohol , and less likely to have graduated from high school . in unadjusted analyses , diabetes was considerably more prevalent in women with grand multiparity ( 25% ) compared with those with fewer livebirths ( 12% ) and with those who were nulliparous ( 15% ; p < 0.001 ) . baseline characteristics of women enrolled in chs data are means  sd or n ( % ) . in regression models controlling for age and race , we observed a nearly 60% increase in prevalent diabetes associated with grand multiparity compared with women who were nulliparous ( table 2 ) . after addition of demographic and clinical factors ( including measures of socioeconomic status ) to the model , the observed increase in prevalent diabetes associated with grand multiparity was reduced to 33% but remained statistically significant . in a final model that controlled for body anthropometrics , the prevalence of diabetes associated with grand multiparity was 21% higher than that for the referent nulliparous group and no longer significant . the degree of attenuation was somewhat less if we adjusted only for bmi recalled from age 50 years rather than measured bmi also ( prevalence ratio 1.27 [ 95% ci , 0.961.68 ] ) , consistent with a greater degree of measurement error in recalled bmi . comparing generalized estimating equation models with and without adjustment for anthropometrics , these factors explained 44% ( [ 286% ] ; p = 0.04 ) of the adjusted association of grand multiparity with prevalent diabetes . multivariable models of parity and prevalent diabetes data are prevalence ratio ( 95% ci ) . * model 2 : age , race , marital status ( 4 ) , income ( 3 ) , education ( 3 ) , height ( in centimeters ) , alcohol yes / no , clinic ( 4 ) , and smoking ( 3 ) . model 3 : also for bmi , bmi at age 50 years , and waist circumference . we observed similar associations between parity and prevalent diabetes in blacks and whites ( table 2 ) . although the association of grand multiparity with prevalent diabetes was numerically stronger in blacks , interaction terms of race with grand multiparity were not significant , even in models that included anthropometrics ( p > 0.16 ) . we next examined the relationship of parity with incident diabetes . despite the strong association of grand multiparity with diabetes at baseline , it was not associated with the incident development of diabetes in older women ( table 3 ) . multivariable models of parity and incident diabetes data are n or hrs ( 95% ci ) . * model 2 : age , race , marital status ( 4 ) , income ( 3 ) , education ( 3 ) , height ( in centimeters ) , alcohol yes / no , clinic ( 4 ) , and smoking ( 3 ) . unadjusted analysis suggested a modest trend toward increased fasting insulin and homa - ir associated with grandmultiparity , whereas fasting glucose was not significantly associated with parity ( table 4 ) . in multivariable models , controlling for age , race , marital status , income , education , height , alcohol use , clinic site , and smoking , we again observed no association with fasting glucose but a trend toward increased insulin and homa - ir with higher categories of parity ( ptrend = 0.01 for insulin and 0.02 for homa - ir ) . parity and metabolic markers among those without diabetes data are means  sd or  coefficient ( sem ) . * multivariable linear regression : models include age , race , marital status ( 4 ) , income ( 3 ) , education ( 3 ) , height ( in centimeters ) , alcohol yes / no , clinic ( 4 ) , and smoking ( 3 ) . in analyses stratified by race , there was no interaction between parity and race and fasting glucose ( p = 0.3 ) . intermediate degrees of parity appeared to have a greater impact on insulin and homa - ir in black women than in white women , whereas the reverse was true for grandmultiparity ( table 4 ) . interaction terms of parity ( defined categorically ) and race were statistically significant in the adjusted models with log(insulin ) and log(homa - ir ) as the outcome ( p = 0.02 for both ) . when we tested the interaction specifically between race and grandmultiparity , the p values for effect modification by race were 0.90 , 0.06 , and 0.09 for glucose , log(insulin ) , and log(homa - ir ) , respectively . the apparent differences by race in the relationship of parity with insulin resistance were not markedly different in analyses that further adjusted for anthropometrics ( data not shown ) . in regression models controlling for age and race , we observed a nearly 60% increase in prevalent diabetes associated with grand multiparity compared with women who were nulliparous ( table 2 ) . after addition of demographic and clinical factors ( including measures of socioeconomic status ) to the model , the observed increase in prevalent diabetes associated with grand multiparity was reduced to 33% but remained statistically significant . in a final model that controlled for body anthropometrics , the prevalence of diabetes associated with grand multiparity was 21% higher than that for the referent nulliparous group and no longer significant . the degree of attenuation was somewhat less if we adjusted only for bmi recalled from age 50 years rather than measured bmi also ( prevalence ratio 1.27 [ 95% ci , 0.961.68 ] ) , consistent with a greater degree of measurement error in recalled bmi . comparing generalized estimating equation models with and without adjustment for anthropometrics , these factors explained 44% ( [ 286% ] ; p = 0.04 ) of the adjusted association of grand multiparity with prevalent diabetes . multivariable models of parity and prevalent diabetes data are prevalence ratio ( 95% ci ) . * model 2 : age , race , marital status ( 4 ) , income ( 3 ) , education ( 3 ) , height ( in centimeters ) , alcohol yes / no , clinic ( 4 ) , and smoking ( 3 ) . model 3 : also for bmi , bmi at age 50 years , and waist circumference . we observed similar associations between parity and prevalent diabetes in blacks and whites ( table 2 ) . although the association of grand multiparity with prevalent diabetes was numerically stronger in blacks , interaction terms of race with grand multiparity were not significant , even in models that included anthropometrics ( p > 0.16 ) . we next examined the relationship of parity with incident diabetes . despite the strong association of grand multiparity with diabetes at baseline , it was not associated with the incident development of diabetes in older women ( table 3 ) . multivariable models of parity and incident diabetes data are n or hrs ( 95% ci ) . * model 2 : age , race , marital status ( 4 ) , income ( 3 ) , education ( 3 ) , height ( in centimeters ) , alcohol yes / no , clinic ( 4 ) , and smoking ( 3 ) . unadjusted analysis suggested a modest trend toward increased fasting insulin and homa - ir associated with grandmultiparity , whereas fasting glucose was not significantly associated with parity ( table 4 ) . in multivariable models , controlling for age , race , marital status , income , education , height , alcohol use , clinic site , and smoking , we again observed no association with fasting glucose but a trend toward increased insulin and homa - ir with higher categories of parity ( ptrend = 0.01 for insulin and 0.02 for homa - ir ) . parity and metabolic markers among those without diabetes data are means  sd or  coefficient ( sem ) . * multivariable linear regression : models include age , race , marital status ( 4 ) , income ( 3 ) , education ( 3 ) , height ( in centimeters ) , alcohol yes / no , clinic ( 4 ) , and smoking ( 3 ) . in analyses stratified by race , there was no interaction between parity and race and fasting glucose ( p = 0.3 ) . intermediate degrees of parity appeared to have a greater impact on insulin and homa - ir in black women than in white women , whereas the reverse was true for grandmultiparity ( table 4 ) . interaction terms of parity ( defined categorically ) and race were statistically significant in the adjusted models with log(insulin ) and log(homa - ir ) as the outcome ( p = 0.02 for both ) . when we tested the interaction specifically between race and grandmultiparity , the p values for effect modification by race were 0.90 , 0.06 , and 0.09 for glucose , log(insulin ) , and log(homa - ir ) , respectively . the apparent differences by race in the relationship of parity with insulin resistance were not markedly different in analyses that further adjusted for anthropometrics ( data not shown ) . our study , among the first in a population - based cohort of older women , helps to elucidate the relationship of parity with diabetes among older women . we found that grand multiparity ( 5 live births ) was associated with higher risk of prevalent diabetes in analyses controlling for age and race . this association was attenuated but continued to be statistically significant after adjustment for demographic and clinical factors . when anthropometric measures were added to the models , the magnitude of association was further attenuated and became statistically insignificant . thus , the association of grand multiparity with increased prevalence of diabetes seems to be confounded or mediated , in large part , by variation in sociodemographic factors and higher body weight associated with grand multiparity . further , we did not observe any significant variation in the association of parity with diabetes between blacks and whites . on the other hand , among women who remained free of diabetes into older age , parity was not associated with incident diabetes nor with fasting glucose among those free of diabetes . however , there was a statistically significant , albeit weak in magnitude , association of parity and fasting insulin and homa - ir at baseline in women free of diabetes . our data are generally consistent with the hypothesis that parity appears to play an ongoing role in the development of diabetes through middle age , but its effect wanes among older women who have remained nondiabetic . previous studies examining the relationship between parity and diabetes have come to discordant conclusions . some have found that parity , particularly higher levels of parity , is associated with increased risk of diabetes . for example , nicholson et al . ( 5 ) analyzed a cohort of middle - aged black and white women and found that those who had borne 5 children had twice the risk of incident diabetes in unadjusted analyses ( hazard ratio [ hr ] 2.10 [ 95% ci 1.732.53 ] ) . in their study , the increased risk of diabetes was attenuated but persisted after adjustment for factors including sociodemographic factors and anthropometrics ( adjusted hr 1.27 [ 1.021.57 ] ) . consistent with our findings , an earlier study examining the relationship between parity and glucose homeostasis in middle - aged to older women found that after adjustment for covariates , each pregnancy was associated with increased fasting insulin and decreased insulin sensitivity that was not explained by obesity and body composition measures ( 12 ) . this study suggested that changes in insulin sensitivity related to parity persist many years after child - bearing . ( 6 ) examined a cohort of > 120,000 registered nurses , aged 3055 years at baseline . similar to our findings , their unadjusted analyses suggested that women with high parity ( defined in their study as 6 live births ) had a 50% higher risk of incident diabetes over 12 years of follow - up . however , after adjustment for age and bmi , they observed no significant relationship between parity and incident diabetes . their results are concordant with ours in connecting high levels of parity with diabetes in middle age and in the finding that the effects of high parity are largely related to the impact on midlife weight . together , ours and the previous study emphasize the importance of weight management targeted to multiparous women . we also had standardized measures of several potential confounding and mediating variables , including directly measured weight and waist circumference , allowing us to examine the mechanism of association of our relationship of interest . first , some information that may have been helpful in the interpretation of the results , such as history of gestational diabetes mellitus , was not available . next , our results may not be generalizable to other races and ethnicities , such as native americans and asians ( 13,14 ) , and we may have had limited power to detect differences between blacks and whites . in addition , we had a limited number of cases of incident diabetes , particularly among grand multiparous women , 25% of whom already had diabetes at baseline , limiting our power to detect a small effect of grand multiparity on incident diabetes . finally , because our analysis examines the relationship between a relatively early life exposure with outcomes much later in life , our results may be influenced by survivor bias . accordingly , selective survival of women into older age may result in an attenuated measured association of parity with later diabetes outcomes . in summary , grand multiparity , but not lesser degrees of parity , much of the higher prevalence of diabetes associated with past child - bearing seems to be mediated ( or confounded ) by the heavier bmi associated with grand multiparity . this finding presents an opportunity for education and intervention related to weight control among grand multiparous women to reduce diabetes prevalence . in nondiabetic women who have reached older age , higher parity is not associated with increased risk of developing diabetes but tended to be associated with small increases in fasting insulin and homa - ir . additional longitudinal studies designed specifically to study the long - term medical effects of child - bearing should be conducted . future studies should include multiple longitudinal measures of diabetes related factors over time to better understand the interplay of risk factors over time .