Model 1 - Age * Incident IPV * Sex Interaction Plot (MAP)
Mean Arterial Pressure - Model 1
m1MAP = lmer(MAP~Age * IncidIPV * Sex * Race + BMIcat + PovStat + (1|subclass) + (Age|HNDid),IncidIPVbp)
(st = step(m1MAP))
Mean Arterial Pressure - Final Model 1 (after backwards elimination)
m1MAP = lmer(MAP~Age + IncidIPV + Sex + BMIcat + PovStat + Age:IncidIPV + Age:Sex + IncidIPV:Sex + Age:IncidIPV:Sex + (1|subclass) + (Age|HNDid),IncidIPVbp)
summary(m1MAP)
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: MAP ~ Age + IncidIPV + Sex + BMIcat + PovStat + Age:IncidIPV + Age:Sex + IncidIPV:Sex + Age:IncidIPV:Sex + (1 | subclass) + (Age | HNDid)
## Data: IncidIPVbp
##
## REML criterion at convergence: 1002
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 45.2012 6.723
## Age 0.0264 0.163 0.12
## subclass (Intercept) 3.2316 1.798
## Residual 63.2176 7.951
## Number of obs: 138, groups: HNDid, 71; subclass, 46
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 87.533 3.032 24.700 28.87 <2e-16
## Age 0.923 0.273 20.200 3.38 0.0029
## IncidIPV1 -2.710 4.825 29.900 -0.56 0.5785
## SexMen -4.769 3.685 28.400 -1.29 0.2060
## BMIcatObese 5.232 2.229 101.800 2.35 0.0208
## PovStatBelow 4.770 2.242 59.700 2.13 0.0376
## Age:IncidIPV1 -0.805 0.505 32.700 -1.59 0.1207
## Age:SexMen -0.722 0.386 28.500 -1.87 0.0715
## IncidIPV1:SexMen 13.923 6.880 31.200 2.02 0.0516
## Age:IncidIPV1:SexMen 1.781 0.702 30.200 2.54 0.0166
##
## Correlation of Fixed Effects:
## (Intr) Age InIPV1 SexMen BMIctO PvSttB Ag:IIPV1 Ag:SxM
## Age 0.488
## IncidIPV1 -0.386 -0.340
## SexMen -0.547 -0.435 0.390
## BMIcatObese -0.436 0.082 -0.187 -0.113
## PovStatBelw -0.325 0.016 0.120 0.103 -0.030
## Ag:IncdIPV1 -0.234 -0.537 0.761 0.253 -0.132 0.037
## Age:SexMen -0.337 -0.698 0.256 0.697 -0.093 0.033 0.392
## IncdIPV1:SM 0.272 0.244 -0.736 -0.572 0.214 -0.188 -0.553 -0.396
## Ag:IIPV1:SM 0.182 0.390 -0.576 -0.412 0.150 -0.137 -0.737 -0.569
## IIPV1:
## Age
## IncidIPV1
## SexMen
## BMIcatObese
## PovStatBelw
## Ag:IncdIPV1
## Age:SexMen
## IncdIPV1:SM
## Ag:IIPV1:SM 0.725
Model 1 - Age * Incident IPV * Sex Interaction Plot (MAP)
Mean Arterial Pressure - Model 2
sBP and dBP adjustments for MAP
m2MAP = lmer(adjMAP~Age * IncidIPV * Sex * Race + BMIcat + PovStat + (Age|HNDid) + (1|subclass),IncidIPVbp)
(st = step(m2MAP))
Mean Arterial Pressure - Final Model 2
m2MAP = lmer(adjMAP~Age + IncidIPV + Sex + BMIcat + Age:IncidIPV + Age:Sex + IncidIPV:Sex + Age:IncidIPV:Sex + (Age|HNDid) + (1|subclass),IncidIPVbp)
summary(m2MAP)
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: adjMAP ~ Age + IncidIPV + Sex + BMIcat + Age:IncidIPV + Age:Sex + IncidIPV:Sex + Age:IncidIPV:Sex + (Age | HNDid) + (1 | subclass)
## Data: IncidIPVbp
##
## REML criterion at convergence: 1016
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 50.4183 7.101
## Age 0.0863 0.294 -0.05
## subclass (Intercept) 8.4806 2.912
## Residual 58.1547 7.626
## Number of obs: 138, groups: HNDid, 71; subclass, 46
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 90.896 2.978 24.300 30.52 <2e-16
## Age 0.979 0.291 24.100 3.37 0.0025
## IncidIPV1 -3.004 4.898 31.100 -0.61 0.5441
## SexMen -7.187 3.754 30.700 -1.91 0.0649
## BMIcatObese 6.223 2.324 103.300 2.68 0.0086
## Age:IncidIPV1 -0.841 0.526 35.700 -1.60 0.1191
## Age:SexMen -0.928 0.406 31.400 -2.28 0.0294
## IncidIPV1:SexMen 18.538 6.967 30.500 2.66 0.0123
## Age:IncidIPV1:SexMen 2.270 0.737 31.700 3.08 0.0042
##
## Correlation of Fixed Effects:
## (Intr) Age InIPV1 SexMen BMIctO Ag:IIPV1 Ag:SxM IIPV1:
## Age 0.487
## IncidIPV1 -0.360 -0.320
## SexMen -0.535 -0.407 0.392
## BMIcatObese -0.473 0.083 -0.190 -0.119
## Ag:IncdIPV1 -0.212 -0.537 0.734 0.244 -0.133
## Age:SexMen -0.310 -0.693 0.245 0.661 -0.099 0.403
## IncdIPV1:SM 0.223 0.233 -0.736 -0.571 0.213 -0.538 -0.380
## Ag:IIPV1:SM 0.131 0.388 -0.545 -0.389 0.140 -0.736 -0.575 0.684
Model 2 - Age * Incident IPV * Sex Interaction Plot (MAP)
Mean Arterial Pressure - Model 3
sBP adjustments for MAP
m3MAP = lmer(MAP_adjsBP~Age * IncidIPV * Sex * Race + BMIcat + PovStat + (1|subclass) + (Age|HNDid),IncidIPVbp)
(st = step(m3MAP))
Mean Arterial Pressure - Final Model 3
m3MAP = lmer(MAP_adjsBP~Age + IncidIPV + Sex + BMIcat + PovStat + Age:IncidIPV + Age:Sex + IncidIPV:Sex + Age:IncidIPV:Sex + (Age|HNDid) + (1|subclass),IncidIPVbp)
summary(m3MAP)
## [1] "Asymptotic covariance matrix A is not positive!"
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: MAP_adjsBP ~ Age + IncidIPV + Sex + BMIcat + PovStat + Age:IncidIPV + Age:Sex + IncidIPV:Sex + Age:IncidIPV:Sex + (Age | HNDid) + (1 | subclass)
## Data: IncidIPVbp
##
## REML criterion at convergence: 1003
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 39.2701 6.267
## Age 0.0106 0.103 -1.00
## subclass (Intercept) 5.4041 2.325
## Residual 60.7393 7.794
## Number of obs: 138, groups: HNDid, 71; subclass, 46
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 88.073 2.923 31.400 30.14 < 2e-16
## Age 0.953 0.270 84.800 3.53 0.00067
## IncidIPV1 -2.130 4.622 31.700 -0.46 0.64812
## SexMen -5.712 3.528 31.000 -1.62 0.11556
## BMIcatObese 5.581 2.236 99.000 2.50 0.01420
## PovStatBelow 5.013 2.285 60.400 2.19 0.03209
## Age:IncidIPV1 -0.796 0.496 81.900 -1.60 0.11258
## Age:SexMen -0.835 0.383 101.100 -2.18 0.03145
## IncidIPV1:SexMen 15.032 6.634 34.300 2.27 0.02988
## Age:IncidIPV1:SexMen 1.950 0.699 102.300 2.79 0.00629
##
## Correlation of Fixed Effects:
## (Intr) Age InIPV1 SexMen BMIctO PvSttB Ag:IIPV1 Ag:SxM
## Age 0.426
## IncidIPV1 -0.372 -0.304
## SexMen -0.524 -0.387 0.390
## BMIcatObese -0.437 0.099 -0.203 -0.135
## PovStatBelw -0.351 -0.006 0.146 0.119 -0.043
## Ag:IncdIPV1 -0.201 -0.534 0.724 0.233 -0.141 0.060
## Age:SexMen -0.285 -0.690 0.234 0.650 -0.112 0.048 0.394
## IncdIPV1:SM 0.269 0.221 -0.742 -0.578 0.226 -0.219 -0.530 -0.375
## Ag:IIPV1:SM 0.165 0.386 -0.548 -0.392 0.152 -0.163 -0.734 -0.573
## IIPV1:
## Age
## IncidIPV1
## SexMen
## BMIcatObese
## PovStatBelw
## Ag:IncdIPV1
## Age:SexMen
## IncdIPV1:SM
## Ag:IIPV1:SM 0.684
Model 3 - Age * Incident IPV * Sex Interaction Plot (MAP)
Mean Arterial Pressure - Model 4
dBP adjustments for MAP
m4MAP = lmer(MAP_adjdBP~Age * IncidIPV * Sex * Race + BMIcat + PovStat + (Age|HNDid) + (1|subclass),IncidIPVbp)
(st = step(m4MAP))
Mean Arterial Pressure - Final Model 4
m4MAP = lmer(MAP_adjdBP~Age + IncidIPV + Sex + BMIcat + PovStat + Age:IncidIPV + Age:Sex + IncidIPV:Sex + Age:IncidIPV:Sex + (Age|HNDid) + (1|subclass),IncidIPVbp)
summary(m4MAP)
## [1] "Asymptotic covariance matrix A is not positive!"
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: MAP_adjdBP ~ Age + IncidIPV + Sex + BMIcat + PovStat + Age:IncidIPV + Age:Sex + IncidIPV:Sex + Age:IncidIPV:Sex + (Age | HNDid) + (1 | subclass)
## Data: IncidIPVbp
##
## REML criterion at convergence: 1003
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 39.2701 6.267
## Age 0.0106 0.103 -1.00
## subclass (Intercept) 5.4041 2.325
## Residual 60.7393 7.794
## Number of obs: 138, groups: HNDid, 71; subclass, 46
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 88.073 2.923 31.400 30.14 < 2e-16
## Age 0.953 0.270 84.800 3.53 0.00067
## IncidIPV1 -2.130 4.622 31.700 -0.46 0.64812
## SexMen -5.712 3.528 31.000 -1.62 0.11556
## BMIcatObese 5.581 2.236 99.000 2.50 0.01420
## PovStatBelow 5.013 2.285 60.400 2.19 0.03209
## Age:IncidIPV1 -0.796 0.496 81.900 -1.60 0.11258
## Age:SexMen -0.835 0.383 101.100 -2.18 0.03145
## IncidIPV1:SexMen 15.032 6.634 34.300 2.27 0.02988
## Age:IncidIPV1:SexMen 1.950 0.699 102.300 2.79 0.00629
##
## Correlation of Fixed Effects:
## (Intr) Age InIPV1 SexMen BMIctO PvSttB Ag:IIPV1 Ag:SxM
## Age 0.426
## IncidIPV1 -0.372 -0.304
## SexMen -0.524 -0.387 0.390
## BMIcatObese -0.437 0.099 -0.203 -0.135
## PovStatBelw -0.351 -0.006 0.146 0.119 -0.043
## Ag:IncdIPV1 -0.201 -0.534 0.724 0.233 -0.141 0.060
## Age:SexMen -0.285 -0.690 0.234 0.650 -0.112 0.048 0.394
## IncdIPV1:SM 0.269 0.221 -0.742 -0.578 0.226 -0.219 -0.530 -0.375
## Ag:IIPV1:SM 0.165 0.386 -0.548 -0.392 0.152 -0.163 -0.734 -0.573
## IIPV1:
## Age
## IncidIPV1
## SexMen
## BMIcatObese
## PovStatBelw
## Ag:IncdIPV1
## Age:SexMen
## IncdIPV1:SM
## Ag:IIPV1:SM 0.684
Model 3 - Age * Incident IPV * Sex Interaction Plot (MAP)
