Methods
Participants were drawn from the Healthy Aging in Neighborhoods of Diversity across the Life Span (HANDLS), a longitudinal population-based study of health disparities based on an area probability sample of Baltimore City initiated in 2004. HANDLS was designed to disentangle the effects of race and SES on risk factors for morbidity and mortality, to examine the incidence and progression of pre-clinical disease, and to follow-up the development and persistence of health disparities, longitudinal health status, and health risks.{Evans, 2004} Participants were examined initially in 2004-2009 and first follow-up examinations began in 2009, approximately 4-5 years after the initial examination.
In the present study we examined 1546 African Americans and whites with mean age of 48.6 years of age at initial examination who have thus far completed two DXA assessments.
Densitometry measurements. Describe DXA.
The first examination measured hip and lumber BMD by DXA (Lunar DPX-iq). The second examination measured BMD by DXA (Hologic QDR Discovery-A). Bland-Altman statistics for cross-calibration between the Lunar and Hologic machines showed a strong correspondence between the two devices.
Bone density measurements were obtained for the total body, hip and lumbar spine.
Covariates.
Smoking and daily alcohol consumption were obtained through self-report.
Race/ethnicity was determined by self-identification. Two racial/ethnic groups were included in this study African American and white.
Analyses. We log-transformed hip and lumber BMD to normalize their distributions and to facilitate interpretation of percent change over time. We performed separate analyses of change in hip and lumber BMD using mixed-model regressions because this technique accounts appropriately for intra-individual correlations over time.{Singer, 2003} We evaluated the effects of race (African-American and white) and socio-economic status (below and above poverty status) on BMD, and adjusted for the effects of age, sex, alcohol consumption, present cigarette smoking, and body mass index. Coefficients from analyses of outcomes transformed by the natural log are interpretable as the percent change in average value of the outcome per unit change in the predictor according to the formula 100(exp(b)-1), where b is a fixed-effect coefficient.{Singer, 2003; Vittinghoff, 2005} We performed analyses with R 2.15.2{R Development Core Team, 2012} and we considered p<.05 for significant differences.
Results
In this subsample of HANDLS, participants who had family incomes above the Federal poverty limit were a year older (p<.05) than participants below the poverty limit (Table 1). There was significantly greater body mass index and fewer current cigarette smokers among those who were above the Federal poverty limit compared with those below the poverty limit. There were no age differences or differences in sex distributions between African Americans and whites, but there were significantly greater proportion of African Americans who were below the poverty limit and were current cigarette smokers. There were no differences in the distribution of sex by poverty status or race. There were also no differences in body mass index associated with race, and neither poverty status nor race was associated with current alcohol.
N Mean SD Min Max
HNDid 1546 8162090533.81 44184614.34 8031107801.00 8224521902.00
HNDwave 1546 1.00 0.00 1.00 1.00
HNDvisit 1546 1.00 0.00 1.00 1.00
Age 1546 48.57 8.92 30.02 66.22
DXAlumbTBMD 1546 1.16 0.19 0.60 2.15
DXAhipTotalBMD 1546 1.09 0.17 0.57 2.13
PhysBMI 1544 29.55 7.01 14.36 56.08
ageDecade 1546 -0.14 0.89 -2.00 1.62
PovStat Age.n Age.mean Age.sd Age.min Age.max
1 Above 926 48.97 9.037 30.02 66.22
2 Below 620 47.97 8.723 30.10 64.85
Welch Two Sample t-test
data: Age by PovStat
t = 2.184, df = 1359, p-value = 0.02913
alternative hypothesis: true difference in means is not equal to 0
95 percent confidence interval:
0.1021 1.9039
sample estimates:
mean in group Above mean in group Below
48.97 47.97
Sex
PovStat Women Men
Above "0.57" "0.43"
Below "0.61" "0.39"
Call: xtabs(formula = ~PovStat + Sex, data = dxa[dxa$HNDwave == 1,
])
Number of cases in table: 1546
Number of factors: 2
Test for independence of all factors:
Chisq = 2.2, df = 1, p-value = 0.1
Race
PovStat White AfrAm
Above "0.52" "0.48"
Below "0.30" "0.70"
Call: xtabs(formula = ~PovStat + Race, data = dxa[dxa$HNDwave == 1,
])
Number of cases in table: 1546
Number of factors: 2
Test for independence of all factors:
Chisq = 74, df = 1, p-value = 6e-18
PovStat PhysBMI.n PhysBMI.mean PhysBMI.sd PhysBMI.min PhysBMI.max
1 Above 925 29.94 6.840 14.36 56.08
2 Below 619 28.97 7.223 14.70 54.65
Welch Two Sample t-test
data: PhysBMI by PovStat
t = 2.647, df = 1275, p-value = 0.008214
alternative hypothesis: true difference in means is not equal to 0
95 percent confidence interval:
0.2517 1.6926
sample estimates:
mean in group Above mean in group Below
29.94 28.97
MedHxAlcCurr
PovStat No Yes
Above "0.75" "0.25"
Below "0.75" "0.25"
Call: xtabs(formula = ~PovStat + MedHxAlcCurr, data = dxa[dxa$HNDwave ==
1, ])
Number of cases in table: 1409
Number of factors: 2
Test for independence of all factors:
Chisq = 0.031, df = 1, p-value = 0.9
MedHxCigaretteCurr
PovStat No Yes
Above "0.63" "0.37"
Below "0.41" "0.59"
Call: xtabs(formula = ~PovStat + MedHxCigaretteCurr, data = dxa[dxa$HNDwave ==
1, ])
Number of cases in table: 1469
Number of factors: 2
Test for independence of all factors:
Chisq = 66, df = 1, p-value = 4e-16
Race Age.n Age.mean Age.sd Age.min Age.max
1 White 664 48.86 9.026 30.10 66.22
2 AfrAm 882 48.35 8.843 30.02 65.98
Welch Two Sample t-test
data: Age by Race
t = 1.118, df = 1412, p-value = 0.2636
alternative hypothesis: true difference in means is not equal to 0
95 percent confidence interval:
-0.3877 1.4160
sample estimates:
mean in group White mean in group AfrAm
48.86 48.35
Sex
Race Women Men
White "0.59" "0.41"
AfrAm "0.59" "0.41"
Call: xtabs(formula = ~Race + Sex, data = dxa[dxa$HNDwave == 1, ])
Number of cases in table: 1546
Number of factors: 2
Test for independence of all factors:
Chisq = 0.05, df = 1, p-value = 0.8
PovStat
Race Above Below
White "0.72" "0.28"
AfrAm "0.51" "0.49"
Call: xtabs(formula = ~Race + PovStat, data = dxa[dxa$HNDwave == 1,
])
Number of cases in table: 1546
Number of factors: 2
Test for independence of all factors:
Chisq = 74, df = 1, p-value = 6e-18
Race PhysBMI.n PhysBMI.mean PhysBMI.sd PhysBMI.min PhysBMI.max
1 White 663 29.52 6.933 14.36 56.08
2 AfrAm 881 29.58 7.071 15.21 55.36
Welch Two Sample t-test
data: PhysBMI by Race
t = -0.1672, df = 1440, p-value = 0.8672
alternative hypothesis: true difference in means is not equal to 0
95 percent confidence interval:
-0.7654 0.6451
sample estimates:
mean in group White mean in group AfrAm
29.52 29.58
MedHxAlcCurr
Race No Yes
White "0.75" "0.25"
AfrAm "0.74" "0.26"
Call: xtabs(formula = ~Race + MedHxAlcCurr, data = dxa[dxa$HNDwave ==
1, ])
Number of cases in table: 1409
Number of factors: 2
Test for independence of all factors:
Chisq = 0.18, df = 1, p-value = 0.7
MedHxCigaretteCurr
Race No Yes
White "0.59" "0.41"
AfrAm "0.50" "0.50"
Call: xtabs(formula = ~Race + MedHxCigaretteCurr, data = dxa[dxa$HNDwave ==
1, ])
Number of cases in table: 1469
Number of factors: 2
Test for independence of all factors:
Chisq = 13, df = 1, p-value = 0.0004
In repeated measures analyses age, sex, race, body mass index, and current alcohol use were associated with bone mineral density in the hip (Table 2), and there was a significant interaction between age and race. Hip BMD was not associated with poverty status or current cigarette smoking. Overall, there was a significant decline in hip BMD over time and whites and African Americans declined at significantly different rates as indicated by a significant interaction between age and race. There were no differences in rates of decline by poverty status. Whites declined by 4.6% per decade; African Americans declined by 3.2% per decade (Figure 1).
Age, sex, race, and body mass index were associated with bone mineral density in the lumbar spine (Table 2), and there was a significant interaction between age and sex. Lumbar BMD was not associated with poverty status, current alcohol consumption, or current smoking. Overall, there was a significant decline in lumbar BMD over time and men and women declined at significantly different rates as indicated by a significant interaction between age and sex. There were no differences in the rates of decline in lumbar spine poverty status. There was no difference in the rate of decline by race: Whites declined by 7.3% per decade; African Americans declined by 7.9% per decade (Figure 2).
Simultaneous Tests for General Linear Hypotheses
Fit: lmer(formula = logHip ~ ageDecade * (Sex + Race) + PovStat +
PhysBMI + MedHxAlcCurr + MedHxCigaretteCurr + (ageDecade |
HNDid), data = dxa, na.action = na.omit)
Linear Hypotheses:
Estimate Std. Error z value Pr(>|z|)
(Intercept) == 0 -0.283211 0.016940 -16.72 < 2e-16
ageDecade == 0 -0.046979 0.005064 -9.28 < 2e-16
SexMen == 0 0.063768 0.006678 9.55 < 2e-16
RaceAfrAm == 0 0.052136 0.006647 7.84 4.4e-15
PovStatBelow == 0 -0.006024 0.006849 -0.88 0.379
PhysBMI == 0 0.010500 0.000492 21.35 < 2e-16
MedHxAlcCurrYes == 0 -0.023253 0.007507 -3.10 0.002
MedHxCigaretteCurrYes == 0 -0.011129 0.006798 -1.64 0.102
ageDecade:SexMen == 0 0.011409 0.005906 1.93 0.053
ageDecade:RaceAfrAm == 0 0.014431 0.005888 2.45 0.014
(Univariate p values reported)
Simultaneous Tests for General Linear Hypotheses
Fit: lmer(formula = logLum ~ ageDecade * (Sex + Race) + PovStat +
PhysBMI + MedHxAlcCurr + MedHxCigaretteCurr + (ageDecade |
HNDid), data = dxa, na.action = na.omit)
Linear Hypotheses:
Estimate Std. Error z value Pr(>|z|)
(Intercept) == 0 -0.160506 0.021034 -7.63 2.3e-14
ageDecade == 0 -0.075408 0.006730 -11.20 < 2e-16
SexMen == 0 0.038667 0.008352 4.63 3.7e-06
RaceAfrAm == 0 0.053153 0.008319 6.39 1.7e-10
PovStatBelow == 0 0.016175 0.008526 1.90 0.058
PhysBMI == 0 0.007262 0.000611 11.89 < 2e-16
MedHxAlcCurrYes == 0 0.004626 0.009437 0.49 0.624
MedHxCigaretteCurrYes == 0 -0.015175 0.008429 -1.80 0.072
ageDecade:SexMen == 0 0.031260 0.007848 3.98 6.8e-05
ageDecade:RaceAfrAm == 0 -0.007018 0.007828 -0.90 0.370
(Univariate p values reported)