Model 1 - Age * LifeIPV * Sex Interaction Plot (MAP)
Mean Arterial Pressure - Model 1
m1MAP = lmer(MAP~Age * LifeIPV * Sex * Race + BMIcat + PovStat + (1|subclass) + (Age|HNDid),LifeIPVbp)
(st = step(m1MAP))
Mean Arterial Pressure - Final Model 1 (after backwards elimination)
m1MAP = lmer(MAP~Age + LifeIPV + Sex + Race + BMIcat + PovStat + Age:LifeIPV + Age:Sex + LifeIPV:Sex + Age:LifeIPV:Sex + (1|subclass) + (Age|HNDid),LifeIPVbp)
summary(m1MAP)
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: MAP ~ Age + LifeIPV + Sex + Race + BMIcat + PovStat + Age:LifeIPV + Age:Sex + LifeIPV:Sex + Age:LifeIPV:Sex + (1 | subclass) + (Age | HNDid)
## Data: LifeIPVbp
##
## REML criterion at convergence: 4983
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 4.24e+01 6.5114
## Age 1.07e-04 0.0104 -1.00
## subclass (Intercept) 0.00e+00 0.0000
## Residual 7.20e+01 8.4839
## Number of obs: 666, groups: HNDid, 334; subclass, 222
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 78.6060 1.3953 355.0000 56.34 < 2e-16
## Age 0.1371 0.0909 415.0000 1.51 0.13215
## LifeIPV1 0.6015 1.2788 296.0000 0.47 0.63844
## SexMen 1.4964 1.2051 294.0000 1.24 0.21534
## RaceAfrAm 3.5934 1.0705 322.0000 3.36 0.00088
## BMIcatObese 5.4789 1.0616 546.0000 5.16 3.4e-07
## PovStatBelow 3.1802 1.0266 320.0000 3.10 0.00212
## Age:LifeIPV1 -0.0819 0.1365 412.0000 -0.60 0.54885
## Age:SexMen 0.0058 0.1343 420.0000 0.04 0.96558
## LifeIPV1:SexMen 6.1225 2.7805 292.0000 2.20 0.02845
## Age:LifeIPV1:SexMen 0.7231 0.3232 441.0000 2.24 0.02576
##
## Correlation of Fixed Effects:
## (Intr) Age LfIPV1 SexMen RcAfrA BMIctO PvSttB Ag:LIPV1
## Age 0.052
## LifeIPV1 -0.437 -0.110
## SexMen -0.439 -0.098 0.440
## RaceAfrAm -0.479 -0.003 0.203 -0.056
## BMIcatObese -0.528 0.024 -0.082 0.054 -0.107
## PovStatBelw -0.252 0.102 -0.091 0.068 -0.039 0.022
## Age:LifIPV1 -0.020 -0.667 0.132 0.062 0.025 -0.061 -0.072
## Age:SexMen -0.021 -0.679 0.072 0.121 -0.018 -0.013 -0.090 0.452
## LfIPV1:SxMn 0.159 0.040 -0.435 -0.430 -0.002 0.074 -0.070 -0.054
## Ag:LIPV1:SM -0.032 0.279 -0.042 -0.053 0.052 0.053 -0.008 -0.420
## Ag:SxM LIPV1:
## Age
## LifeIPV1
## SexMen
## RaceAfrAm
## BMIcatObese
## PovStatBelw
## Age:LifIPV1
## Age:SexMen
## LfIPV1:SxMn -0.049
## Ag:LIPV1:SM -0.413 0.299
Model 1 - Age * LifeIPV * Sex Interaction Plot (MAP)
Mean Arterial Pressure - Model 2
sBP and dBP adjustments for MAP
m2MAP = lmer(adjMAP~Age * LifeIPV * Sex * Race + BMIcat + PovStat + (Age|HNDid) + (1|subclass),LifeIPVbp)
(st = step(m2MAP))
Mean Arterial Pressure - Final Model 2
m2MAP = lmer(adjMAP~Age + LifeIPV + Sex + Race + BMIcat + Age:LifeIPV + Age:Sex + LifeIPV:Sex + Age:LifeIPV:Sex + (Age|HNDid) + (1|subclass),LifeIPVbp)
summary(m2MAP)
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: adjMAP ~ Age + LifeIPV + Sex + Race + BMIcat + Age:LifeIPV + Age:Sex + LifeIPV:Sex + Age:LifeIPV:Sex + (Age | HNDid) + (1 | subclass)
## Data: LifeIPVbp
##
## REML criterion at convergence: 5049
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 5.72e+01 7.565638
## Age 2.18e-05 0.004671 1.00
## subclass (Intercept) 1.06e-07 0.000326
## Residual 7.33e+01 8.561082
## Number of obs: 666, groups: HNDid, 334; subclass, 222
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 80.9383 1.4673 368.0000 55.16 < 2e-16
## Age 0.2312 0.0976 435.0000 2.37 0.01833
## LifeIPV1 0.7436 1.3946 295.0000 0.53 0.59430
## SexMen 0.6252 1.3167 294.0000 0.47 0.63525
## RaceAfrAm 3.9346 1.1681 322.0000 3.37 0.00085
## BMIcatObese 6.5423 1.1319 574.0000 5.78 1.2e-08
## Age:LifeIPV1 -0.0810 0.1471 433.0000 -0.55 0.58219
## Age:SexMen -0.0229 0.1444 442.0000 -0.16 0.87430
## LifeIPV1:SexMen 8.2907 3.0349 295.0000 2.73 0.00668
## Age:LifeIPV1:SexMen 0.8997 0.3467 462.0000 2.60 0.00976
##
## Correlation of Fixed Effects:
## (Intr) Age LfIPV1 SexMen RcAfrA BMIctO Ag:LIPV1 Ag:SxM
## Age 0.095
## LifeIPV1 -0.483 -0.114
## SexMen -0.439 -0.118 0.449
## RaceAfrAm -0.511 0.000 0.200 -0.053
## BMIcatObese -0.531 0.020 -0.078 0.052 -0.102
## Age:LifIPV1 -0.049 -0.665 0.145 0.075 0.023 -0.060
## Age:SexMen -0.053 -0.676 0.072 0.145 -0.023 -0.010 0.449
## LfIPV1:SxMn 0.150 0.053 -0.445 -0.428 -0.006 0.074 -0.067 -0.063
## Ag:LIPV1:SM -0.031 0.283 -0.052 -0.060 0.050 0.054 -0.425 -0.418
## LIPV1:
## Age
## LifeIPV1
## SexMen
## RaceAfrAm
## BMIcatObese
## Age:LifIPV1
## Age:SexMen
## LfIPV1:SxMn
## Ag:LIPV1:SM 0.310
Model 2 - Age * LifeIPV * Sex Interaction Plot (MAP)
Mean Arterial Pressure - Model 3
sBP adjustments for MAP
m3MAP = lmer(MAP_adjsBP~Age * LifeIPV * Sex * Race + BMIcat + PovStat + (1|subclass) + (Age|HNDid),LifeIPVbp)
(st = step(m3MAP))
Mean Arterial Pressure - Final Model 3
m3MAP = lmer(MAP_adjsBP~Age + LifeIPV + Sex + Race + BMIcat + PovStat + Age:LifeIPV + Age:Sex + LifeIPV:Sex + Age:LifeIPV:Sex + (Age|HNDid) + (1|subclass),LifeIPVbp)
summary(m3MAP)
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: MAP_adjsBP ~ Age + LifeIPV + Sex + Race + BMIcat + PovStat + Age:LifeIPV + Age:Sex + LifeIPV:Sex + Age:LifeIPV:Sex + (Age | HNDid) + (1 | subclass)
## Data: LifeIPVbp
##
## REML criterion at convergence: 5002
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 4.76e+01 6.90e+00
## Age 6.43e-05 8.02e-03 -1.00
## subclass (Intercept) 1.91e-10 1.38e-05
## Residual 7.18e+01 8.47e+00
## Number of obs: 666, groups: HNDid, 334; subclass, 222
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 79.2004 1.4364 357.0000 55.14 < 2e-16
## Age 0.1971 0.0933 422.0000 2.11 0.03525
## LifeIPV1 0.5071 1.3199 296.0000 0.38 0.70113
## SexMen 1.1950 1.2439 294.0000 0.96 0.33750
## RaceAfrAm 3.7121 1.1045 321.0000 3.36 0.00087
## BMIcatObese 6.0758 1.0849 557.0000 5.60 3.4e-08
## PovStatBelow 3.1524 1.0593 320.0000 2.98 0.00314
## Age:LifeIPV1 -0.0930 0.1402 419.0000 -0.66 0.50721
## Age:SexMen -0.0243 0.1379 428.0000 -0.18 0.86044
## LifeIPV1:SexMen 6.8847 2.8687 294.0000 2.40 0.01702
## Age:LifeIPV1:SexMen 0.7996 0.3312 449.0000 2.41 0.01618
##
## Correlation of Fixed Effects:
## (Intr) Age LfIPV1 SexMen RcAfrA BMIctO PvSttB Ag:LIPV1
## Age 0.054
## LifeIPV1 -0.439 -0.111
## SexMen -0.440 -0.099 0.440
## RaceAfrAm -0.481 -0.003 0.203 -0.056
## BMIcatObese -0.524 0.023 -0.081 0.054 -0.106
## PovStatBelw -0.253 0.102 -0.091 0.068 -0.039 0.021
## Age:LifIPV1 -0.021 -0.667 0.135 0.063 0.025 -0.062 -0.072
## Age:SexMen -0.022 -0.679 0.073 0.123 -0.018 -0.013 -0.090 0.452
## LfIPV1:SxMn 0.160 0.041 -0.435 -0.430 -0.002 0.073 -0.070 -0.055
## Ag:LIPV1:SM -0.031 0.279 -0.044 -0.054 0.052 0.053 -0.007 -0.421
## Ag:SxM LIPV1:
## Age
## LifeIPV1
## SexMen
## RaceAfrAm
## BMIcatObese
## PovStatBelw
## Age:LifIPV1
## Age:SexMen
## LfIPV1:SxMn -0.050
## Ag:LIPV1:SM -0.414 0.300
Model 3 - Age * LifeIPV * Sex Interaction Plot (MAP)
Mean Arterial Pressure - Model 4
dBP adjustments for MAP
m4MAP = lmer(MAP_adjdBP~Age * LifeIPV * Sex * Race + BMIcat + PovStat + (Age|HNDid) + (1|subclass),LifeIPVbp)
(st = step(m4MAP))
Mean Arterial Pressure - Final Model 4
m4MAP = lmer(MAP_adjdBP~Age + LifeIPV + Sex + Race + BMIcat + PovStat + Age:LifeIPV + Age:Sex + LifeIPV:Sex + Age:LifeIPV:Sex + (Age|HNDid) + (1|subclass),LifeIPVbp)
summary(m4MAP)
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: MAP_adjdBP ~ Age + LifeIPV + Sex + Race + BMIcat + PovStat + Age:LifeIPV + Age:Sex + LifeIPV:Sex + Age:LifeIPV:Sex + (Age | HNDid) + (1 | subclass)
## Data: LifeIPVbp
##
## REML criterion at convergence: 5002
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 4.76e+01 6.90e+00
## Age 6.43e-05 8.02e-03 -1.00
## subclass (Intercept) 1.91e-10 1.38e-05
## Residual 7.18e+01 8.47e+00
## Number of obs: 666, groups: HNDid, 334; subclass, 222
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 79.2004 1.4364 357.0000 55.14 < 2e-16
## Age 0.1971 0.0933 422.0000 2.11 0.03525
## LifeIPV1 0.5071 1.3199 296.0000 0.38 0.70113
## SexMen 1.1950 1.2439 294.0000 0.96 0.33750
## RaceAfrAm 3.7121 1.1045 321.0000 3.36 0.00087
## BMIcatObese 6.0758 1.0849 557.0000 5.60 3.4e-08
## PovStatBelow 3.1524 1.0593 320.0000 2.98 0.00314
## Age:LifeIPV1 -0.0930 0.1402 419.0000 -0.66 0.50721
## Age:SexMen -0.0243 0.1379 428.0000 -0.18 0.86044
## LifeIPV1:SexMen 6.8847 2.8687 294.0000 2.40 0.01702
## Age:LifeIPV1:SexMen 0.7996 0.3312 449.0000 2.41 0.01618
##
## Correlation of Fixed Effects:
## (Intr) Age LfIPV1 SexMen RcAfrA BMIctO PvSttB Ag:LIPV1
## Age 0.054
## LifeIPV1 -0.439 -0.111
## SexMen -0.440 -0.099 0.440
## RaceAfrAm -0.481 -0.003 0.203 -0.056
## BMIcatObese -0.524 0.023 -0.081 0.054 -0.106
## PovStatBelw -0.253 0.102 -0.091 0.068 -0.039 0.021
## Age:LifIPV1 -0.021 -0.667 0.135 0.063 0.025 -0.062 -0.072
## Age:SexMen -0.022 -0.679 0.073 0.123 -0.018 -0.013 -0.090 0.452
## LfIPV1:SxMn 0.160 0.041 -0.435 -0.430 -0.002 0.073 -0.070 -0.055
## Ag:LIPV1:SM -0.031 0.279 -0.044 -0.054 0.052 0.053 -0.007 -0.421
## Ag:SxM LIPV1:
## Age
## LifeIPV1
## SexMen
## RaceAfrAm
## BMIcatObese
## PovStatBelw
## Age:LifIPV1
## Age:SexMen
## LfIPV1:SxMn -0.050
## Ag:LIPV1:SM -0.414 0.300
Model 3 - Age * LifeIPV * Sex Interaction Plot (MAP)
