Analyses excluding participants with osteophytes

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 1442 African Americans and whites with mean age of 48.3 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.

125 participants with osteophytes were excluded from these analyses.

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          1442 8162651617.72 43065677.53 8031117901.00 8224521902.00
HNDwave        1442          1.00        0.00          1.00          1.00
HNDvisit       1442          1.00        0.00          1.00          1.00
Age            1442         48.29        8.94         30.02         66.22
DXAlumbTBMD    1442          1.16        0.19          0.60          2.15
DXAhipTotalBMD 1442          1.09        0.17          0.57          2.13
PhysBMI        1440         29.54        7.00         14.36         56.08
ageDecade      1442         -0.17        0.89         -2.00          1.62
osteo             0           NaN          NA           Inf          -Inf
  PovStat Age.n Age.mean Age.sd Age.min Age.max
1   Above   868    48.69  9.036   30.02   66.22
2   Below   574    47.68  8.768   30.10   64.85

    Welch Two Sample t-test

data:  Age by PovStat 
t = 2.107, df = 1252, p-value = 0.03536
alternative hypothesis: true difference in means is not equal to 0 
95 percent confidence interval:
 0.06905 1.94261 
sample estimates:
mean in group Above mean in group Below 
              48.69               47.68 
       Sex
PovStat Women  Men   
  Above "0.58" "0.42"
  Below "0.61" "0.39"
Call: xtabs(formula = ~PovStat + Sex, data = dxa[dxa$HNDwave == 1, 
    ])
Number of cases in table: 1442 
Number of factors: 2 
Test for independence of all factors:
    Chisq = 1, df = 1, p-value = 0.3
       Race
PovStat White  AfrAm 
  Above "0.51" "0.49"
  Below "0.30" "0.70"
Call: xtabs(formula = ~PovStat + Race, data = dxa[dxa$HNDwave == 1, 
    ])
Number of cases in table: 1442 
Number of factors: 2 
Test for independence of all factors:
    Chisq = 63, df = 1, p-value = 2e-15
  PovStat PhysBMI.n PhysBMI.mean PhysBMI.sd PhysBMI.min PhysBMI.max
1   Above       867        29.93      6.858       14.36       56.08
2   Below       573        28.96      7.184       14.70       54.65

    Welch Two Sample t-test

data:  PhysBMI by PovStat 
t = 2.551, df = 1185, p-value = 0.01087
alternative hypothesis: true difference in means is not equal to 0 
95 percent confidence interval:
 0.2238 1.7145 
sample estimates:
mean in group Above mean in group Below 
              29.93               28.96 
       MedHxAlcCurr
PovStat No     Yes   
  Above "0.76" "0.24"
  Below "0.75" "0.25"
Call: xtabs(formula = ~PovStat + MedHxAlcCurr, data = dxa[dxa$HNDwave == 
    1, ])
Number of cases in table: 1323 
Number of factors: 2 
Test for independence of all factors:
    Chisq = 0.023, df = 1, p-value = 0.9
       MedHxCigaretteCurr
PovStat No     Yes   
  Above "0.62" "0.38"
  Below "0.41" "0.59"
Call: xtabs(formula = ~PovStat + MedHxCigaretteCurr, data = dxa[dxa$HNDwave == 
    1, ])
Number of cases in table: 1370 
Number of factors: 2 
Test for independence of all factors:
    Chisq = 59, df = 1, p-value = 1e-14
   Race Age.n Age.mean Age.sd Age.min Age.max
1 White   621    48.50  9.012   30.10   66.22
2 AfrAm   821    48.12  8.889   30.02   65.98

    Welch Two Sample t-test

data:  Age by Race 
t = 0.7957, df = 1326, p-value = 0.4263
alternative hypothesis: true difference in means is not equal to 0 
95 percent confidence interval:
 -0.5556  1.3138 
sample estimates:
mean in group White mean in group AfrAm 
              48.50               48.12 
       Sex
Race    Women  Men   
  White "0.60" "0.40"
  AfrAm "0.59" "0.41"
Call: xtabs(formula = ~Race + Sex, data = dxa[dxa$HNDwave == 1, ])
Number of cases in table: 1442 
Number of factors: 2 
Test for independence of all factors:
    Chisq = 0.12, df = 1, p-value = 0.7
       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: 1442 
Number of factors: 2 
Test for independence of all factors:
    Chisq = 63, df = 1, p-value = 2e-15
   Race PhysBMI.n PhysBMI.mean PhysBMI.sd PhysBMI.min PhysBMI.max
1 White       620        29.43      6.979       14.36       56.08
2 AfrAm       820        29.62      7.025       15.21       55.36

    Welch Two Sample t-test

data:  PhysBMI by Race 
t = -0.5058, df = 1338, p-value = 0.6131
alternative hypothesis: true difference in means is not equal to 0 
95 percent confidence interval:
 -0.9191  0.5423 
sample estimates:
mean in group White mean in group AfrAm 
              29.43               29.62 
       MedHxAlcCurr
Race    No     Yes   
  White "0.76" "0.24"
  AfrAm "0.75" "0.25"
Call: xtabs(formula = ~Race + MedHxAlcCurr, data = dxa[dxa$HNDwave == 
    1, ])
Number of cases in table: 1323 
Number of factors: 2 
Test for independence of all factors:
    Chisq = 0.002, df = 1, p-value = 1
       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: 1370 
Number of factors: 2 
Test for independence of all factors:
    Chisq = 9, df = 1, p-value = 0.002

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 were significant interactions between age by sex and age by 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.0% 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.7% 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.273082   0.017273  -15.81  < 2e-16
ageDecade == 0             -0.047604   0.005188   -9.18  < 2e-16
SexMen == 0                 0.063841   0.006819    9.36  < 2e-16
RaceAfrAm == 0              0.055725   0.006764    8.24  2.2e-16
PovStatBelow == 0          -0.003705   0.006966   -0.53   0.5948
PhysBMI == 0                0.010098   0.000502   20.11  < 2e-16
MedHxAlcCurrYes == 0       -0.020366   0.007694   -2.65   0.0081
MedHxCigaretteCurrYes == 0 -0.010702   0.006930   -1.54   0.1225
ageDecade:SexMen == 0       0.014287   0.006072    2.35   0.0186
ageDecade:RaceAfrAm == 0    0.017123   0.006046    2.83   0.0046
(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.166092   0.021632   -7.68  1.6e-14
ageDecade == 0             -0.079888   0.006803  -11.74  < 2e-16
SexMen == 0                 0.039378   0.008608    4.57  4.8e-06
RaceAfrAm == 0              0.057883   0.008542    6.78  1.2e-11
PovStatBelow == 0           0.018890   0.008761    2.16  0.03107
PhysBMI == 0                0.007225   0.000629   11.49  < 2e-16
MedHxAlcCurrYes == 0        0.007831   0.009789    0.80  0.42370
MedHxCigaretteCurrYes == 0 -0.012515   0.008671   -1.44  0.14893
ageDecade:SexMen == 0       0.028170   0.007948    3.54  0.00039
ageDecade:RaceAfrAm == 0   -0.002683   0.007920   -0.34  0.73483
(Univariate p values reported)

Figure 1. Rate of change in log transformed hip bone mineral density for African Americans and whites above and below 125% of the poverty level
Rate of change in log transformed hip bone mineral density for African Americans and whites above and below 125% of the poverty level

Figure 2. Rate of change in log transformed lumbar bone mineral density for African Americans and whites above and below 125% of the poverty level.
Rate of change in log transformed lumbar bone mineral density for African Americans and whites above and below 125% of the poverty level.