Analyses on complete sample regardless of 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 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)

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.