Systolic Blood Pressure
m1sBP = lmer(adjSBP~Age * LifeIPV * Sex * Race + BMIcat + PovStat + (Age|HNDid) + (1|subclass),LifeIPVbp)
(st = step(m1sBP))
Final Model by backwards elimination
m1sBP = lmer(adjSBP ~Age + LifeIPV + Sex + Race + BMIcat + Age:LifeIPV + Age:Sex + LifeIPV:Sex + Age:LifeIPV:Sex + (Age|HNDid) + (1|subclass),data = LifeIPVbp)
summary(m1sBP)
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: adjSBP ~ Age + LifeIPV + Sex + Race + BMIcat + Age:LifeIPV + Age:Sex + LifeIPV:Sex + Age:LifeIPV:Sex + (Age | HNDid) + (1 | subclass)
## Data: LifeIPVbp
##
## REML criterion at convergence: 5533
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 165.906 12.880
## Age 0.031 0.176 1.00
## subclass (Intercept) 0.000 0.000
## Residual 133.399 11.550
## Number of obs: 666, groups: HNDid, 334; subclass, 222
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 113.413 2.272 386.000 49.93 < 2e-16
## Age 0.795 0.149 477.000 5.34 1.4e-07
## LifeIPV1 -0.860 2.217 280.000 -0.39 0.6984
## SexMen -2.414 2.095 278.000 -1.15 0.2502
## RaceAfrAm 4.542 1.793 314.000 2.53 0.0118
## BMIcatObese 8.975 1.663 604.000 5.40 9.8e-08
## Age:LifeIPV1 -0.211 0.224 477.000 -0.94 0.3472
## Age:SexMen -0.389 0.220 489.000 -1.77 0.0772
## LifeIPV1:SexMen 15.249 4.823 288.000 3.16 0.0017
## Age:LifeIPV1:SexMen 1.370 0.513 502.000 2.67 0.0078
##
## Correlation of Fixed Effects:
## (Intr) Age LfIPV1 SexMen RcAfrA BMIctO Ag:LIPV1 Ag:SxM
## Age 0.201
## LifeIPV1 -0.499 -0.212
## SexMen -0.451 -0.220 0.453
## RaceAfrAm -0.519 -0.010 0.194 -0.051
## BMIcatObese -0.511 0.015 -0.071 0.048 -0.089
## Age:LifIPV1 -0.122 -0.664 0.289 0.143 0.032 -0.056
## Age:SexMen -0.115 -0.676 0.135 0.287 -0.034 -0.007 0.448
## LfIPV1:SxMn 0.163 0.097 -0.447 -0.429 -0.010 0.071 -0.133 -0.124
## Ag:LIPV1:SM 0.001 0.290 -0.118 -0.122 0.044 0.060 -0.437 -0.430
## LIPV1:
## Age
## LifeIPV1
## SexMen
## RaceAfrAm
## BMIcatObese
## Age:LifIPV1
## Age:SexMen
## LfIPV1:SxMn
## Ag:LIPV1:SM 0.414
Age * LifeIPV * Sex Interaction Plot (sBP)
Diastolic Blood Pressure
m1dBP = lmer(adjDBP~Age * LifeIPV * Sex * Race + PovStat + BMIcat + (Age|HNDid) + (1|subclass),LifeIPVbp)
(st = step(m1dBP))
Final Model by backwards elimination
m1dBP = lmer(adjDBP~Age + LifeIPV + Sex + Race + PovStat + BMIcat + Age:LifeIPV + Age:Sex + LifeIPV:Sex + Age:LifeIPV:Sex + (Age|HNDid) + (1|subclass), LifeIPVbp)
summary(m1dBP)
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: adjDBP ~ Age + LifeIPV + Sex + Race + PovStat + BMIcat + Age:LifeIPV + Age:Sex + LifeIPV:Sex + Age:LifeIPV:Sex + (Age | HNDid) + (1 | subclass)
## Data: LifeIPVbp
##
## REML criterion at convergence: 4887
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 3.38e+01 5.8133
## Age 7.28e-04 0.0270 -1.00
## subclass (Intercept) 5.72e-04 0.0239
## Residual 6.38e+01 7.9866
## Number of obs: 666, groups: HNDid, 334; subclass, 222
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 64.0079 1.2794 359.0000 50.03 < 2e-16
## Age -0.0262 0.0835 416.0000 -0.31 0.7535
## LifeIPV1 1.1155 1.1695 306.0000 0.95 0.3409
## SexMen 2.3086 1.1020 304.0000 2.09 0.0370
## RaceAfrAm 3.2575 0.9821 321.0000 3.32 0.0010
## PovStatBelow 2.8070 0.9415 320.0000 2.98 0.0031
## BMIcatObese 5.3263 0.9808 536.0000 5.43 8.5e-08
## Age:LifeIPV1 -0.0370 0.1254 413.0000 -0.30 0.7680
## Age:SexMen 0.1214 0.1235 421.0000 0.98 0.3261
## LifeIPV1:SexMen 4.4797 2.5417 302.0000 1.76 0.0790
## Age:LifeIPV1:SexMen 0.6703 0.2987 428.0000 2.24 0.0253
##
## Correlation of Fixed Effects:
## (Intr) Age LfIPV1 SexMen RcAfrA PvSttB BMIctO Ag:LIPV1
## Age 0.036
## LifeIPV1 -0.435 -0.094
## SexMen -0.438 -0.081 0.439
## RaceAfrAm -0.476 -0.001 0.204 -0.057
## PovStatBelw -0.251 0.100 -0.092 0.068 -0.042
## BMIcatObese -0.531 0.025 -0.083 0.055 -0.109 0.022
## Age:LifIPV1 -0.009 -0.667 0.108 0.050 0.024 -0.070 -0.062
## Age:SexMen -0.011 -0.679 0.062 0.097 -0.016 -0.089 -0.014 0.452
## LfIPV1:SxMn 0.158 0.033 -0.435 -0.430 -0.001 -0.070 0.074 -0.043
## Ag:LIPV1:SM -0.036 0.277 -0.032 -0.043 0.053 -0.009 0.052 -0.418
## Ag:SxM LIPV1:
## Age
## LifeIPV1
## SexMen
## RaceAfrAm
## PovStatBelw
## BMIcatObese
## Age:LifIPV1
## Age:SexMen
## LfIPV1:SxMn -0.039
## Ag:LIPV1:SM -0.411 0.282
Age * LifeIPV * Sex Interaction Plot (dBP)