Systolic Blood Pressure Final Model (BMI as categorical)load("~/Desktop/IPV and Hypertension/Data Analysis/IncidIPVbp.rda")
library(lme4)
## Loading required package: Matrix
library(lmerTest)
## KernSmooth 2.23 loaded
## Copyright M. P. Wand 1997-2009
##
## Attaching package: 'lmerTest'
##
## The following object is masked from 'package:lme4':
##
## lmer
##
## The following object is masked from 'package:stats':
##
## step
library("lme4")
library("multcomp")
## Loading required package: mvtnorm
## Loading required package: TH.data
library("zStat")
Systolic Blood Pressure (BMI as categorical)
load("~/Desktop/IPV and Hypertension/Data Analysis/IncidIPVbp.rda")
library(lme4)
library(lmerTest)
(mm2 = lmer(adjSBP ~(Age + IncidIPV + Sex + PovStat + BMIcat2)^5 + (Age|HNDid) + (1|subclass) ,IncidIPVbp))
## Warning: number of observations <= rank(Z); variance-covariance matrix
## will be unidentifiable
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: adjSBP ~ (Age + IncidIPV + Sex + PovStat + BMIcat2)^5 + (Age | HNDid) + (1 | subclass)
## Data: IncidIPVbp
## REML criterion at convergence: 977.3
## Random effects:
## Groups Name Std.Dev. Corr
## HNDid (Intercept) 1.14e+01
## Age 1.49e-02 -1.00
## subclass (Intercept) 4.09e-05
## Residual 1.12e+01
## Number of obs: 138, groups: HNDid, 70; subclass, 46
## Fixed Effects:
## (Intercept)
## 117.383
## Age
## 1.558
## IncidIPV1
## 10.962
## SexMen
## -7.372
## PovStatBelow
## 7.246
## BMIcat2Obese
## 8.815
## Age:IncidIPV1
## 0.539
## Age:SexMen
## -2.445
## Age:PovStatBelow
## -0.220
## Age:BMIcat2Obese
## -1.119
## IncidIPV1:SexMen
## -4.982
## IncidIPV1:PovStatBelow
## -3.191
## IncidIPV1:BMIcat2Obese
## -2.130
## SexMen:PovStatBelow
## 10.717
## SexMen:BMIcat2Obese
## -2.307
## PovStatBelow:BMIcat2Obese
## 8.991
## Age:IncidIPV1:SexMen
## -0.101
## Age:IncidIPV1:PovStatBelow
## -0.077
## Age:IncidIPV1:BMIcat2Obese
## 0.454
## Age:SexMen:PovStatBelow
## 2.336
## Age:SexMen:BMIcat2Obese
## 1.376
## Age:PovStatBelow:BMIcat2Obese
## 2.425
## IncidIPV1:SexMen:PovStatBelow
## 9.408
## IncidIPV1:SexMen:BMIcat2Obese
## 16.589
## IncidIPV1:PovStatBelow:BMIcat2Obese
## -26.600
## SexMen:PovStatBelow:BMIcat2Obese
## -14.137
## Age:IncidIPV1:SexMen:PovStatBelow
## 5.068
## Age:IncidIPV1:SexMen:BMIcat2Obese
## 1.420
## Age:IncidIPV1:PovStatBelow:BMIcat2Obese
## -6.941
## Age:SexMen:PovStatBelow:BMIcat2Obese
## -2.512
## IncidIPV1:SexMen:PovStatBelow:BMIcat2Obese
## 11.947
## Age:IncidIPV1:SexMen:PovStatBelow:BMIcat2Obese
## -0.435
(st = step(mm2))
## Warning: number of observations <= rank(Z); variance-covariance matrix will be unidentifiable
## Warning: number of observations <= rank(Z); variance-covariance matrix will be unidentifiable
## Random term (Age | HNDid) was eliminated because of having correlation +-1 or NaN
## Warning: number of observations <= rank(Z); variance-covariance matrix
## will be unidentifiable
## Random term (Age + 0 | HNDid) was eliminated because of standard deviation being equal to 0
##
## Random term (1 | subclass) was eliminated because of standard deviation being equal to 0
##
## Random effects:
## Chi.sq Chi.DF elim.num p.value
## (1 | HNDid) 14.19 1 kept 2e-04
##
## Fixed effects:
## Sum Sq Mean Sq NumDF DenDF F.value
## Age:IncidIPV:Sex:PovStat:BMIcat2 0.5741 0.5741 1 99.11 0.0046
## IncidIPV:Sex:PovStat:BMIcat2 20.6063 20.6063 1 91.13 0.1681
## Age:IncidIPV:Sex:BMIcat2 3.9527 3.9527 1 107.96 0.2342
## IncidIPV:Sex:BMIcat2 0.2107 0.2107 1 83.66 0.2476
## Age:Sex:PovStat:BMIcat2 107.5097 107.5097 1 108.79 0.6983
## Sex:PovStat:BMIcat2 3.2473 3.2473 1 93.30 0.5626
## Age:Sex:BMIcat2 12.3043 12.3043 1 111.26 0.5800
## Sex:BMIcat2 441.8659 441.8659 1 89.23 1.6245
## Age:IncidIPV:Sex:PovStat 497.5082 497.5082 1 111.77 2.0858
## IncidIPV:Sex:PovStat 17.0158 17.0158 1 55.80 1.7853
## Age 1919.9581 1919.9581 1 107.99 16.3612
## IncidIPV 600.6325 600.6325 1 85.18 0.9962
## Sex 3.7749 3.7749 1 93.60 0.0509
## PovStat 662.7933 662.7933 1 86.04 4.8118
## BMIcat2 1233.7989 1233.7989 1 83.07 3.3153
## Age:IncidIPV 191.2687 191.2687 1 112.23 0.1640
## Age:Sex 412.1801 412.1801 1 110.57 0.1885
## Age:PovStat 835.1758 835.1758 1 111.86 3.1034
## Age:BMIcat2 18.0002 18.0002 1 107.40 0.7618
## IncidIPV:Sex 46.3274 46.3274 1 91.59 0.9851
## IncidIPV:PovStat 411.1877 411.1877 1 89.47 0.3146
## IncidIPV:BMIcat2 21.3280 21.3280 1 84.49 0.0658
## Sex:PovStat 56.5951 56.5951 1 84.53 1.0014
## PovStat:BMIcat2 2.8045 2.8045 1 82.67 0.3669
## Age:IncidIPV:Sex 205.9596 205.9596 1 110.65 4.5914
## Age:IncidIPV:PovStat 183.6904 183.6904 1 113.07 1.1025
## Age:IncidIPV:BMIcat2 78.0728 78.0728 1 105.41 3.4977
## Age:Sex:PovStat 675.4175 675.4175 1 108.99 5.3102
## Age:PovStat:BMIcat2 6.8082 6.8082 1 104.06 2.6315
## IncidIPV:PovStat:BMIcat2 202.4580 202.4580 1 82.21 1.4673
## Age:IncidIPV:PovStat:BMIcat2 875.6610 875.6610 1 106.03 7.8463
## elim.num Pr(>F)
## Age:IncidIPV:Sex:PovStat:BMIcat2 1 0.9460
## IncidIPV:Sex:PovStat:BMIcat2 2 0.6828
## Age:IncidIPV:Sex:BMIcat2 3 0.6294
## IncidIPV:Sex:BMIcat2 4 0.6201
## Age:Sex:PovStat:BMIcat2 5 0.4052
## Sex:PovStat:BMIcat2 6 0.4551
## Age:Sex:BMIcat2 7 0.4479
## Sex:BMIcat2 8 0.2058
## Age:IncidIPV:Sex:PovStat 9 0.1515
## IncidIPV:Sex:PovStat 10 0.1869
## Age kept 1e-04
## IncidIPV kept 0.3211
## Sex kept 0.8221
## PovStat kept 0.0310
## BMIcat2 kept 0.0722
## Age:IncidIPV kept 0.6863
## Age:Sex kept 0.6650
## Age:PovStat kept 0.0809
## Age:BMIcat2 kept 0.3847
## IncidIPV:Sex kept 0.3236
## IncidIPV:PovStat kept 0.5763
## IncidIPV:BMIcat2 kept 0.7982
## Sex:PovStat kept 0.3198
## PovStat:BMIcat2 kept 0.5464
## Age:IncidIPV:Sex kept 0.0343
## Age:IncidIPV:PovStat kept 0.2959
## Age:IncidIPV:BMIcat2 kept 0.0642
## Age:Sex:PovStat kept 0.0231
## Age:PovStat:BMIcat2 kept 0.1078
## IncidIPV:PovStat:BMIcat2 kept 0.2292
## Age:IncidIPV:PovStat:BMIcat2 kept 0.0061
##
## Least squares means:
## IncidIPV Sex PovStat BMIcat2
## IncidIPV 0 1.0 NA NA NA
## IncidIPV 1 2.0 NA NA NA
## Sex Women NA 2.0 NA NA
## Sex Men NA 1.0 NA NA
## PovStat Above NA NA 1.0 NA
## PovStat Below NA NA 2.0 NA
## BMIcat2 nonObese NA NA NA 1.0
## BMIcat2 Obese NA NA NA 2.0
## IncidIPV:Sex 0 Women 1.0 2.0 NA NA
## IncidIPV:Sex 1 Women 2.0 2.0 NA NA
## IncidIPV:Sex 0 Men 1.0 1.0 NA NA
## IncidIPV:Sex 1 Men 2.0 1.0 NA NA
## IncidIPV:PovStat 0 Above 1.0 NA 1.0 NA
## IncidIPV:PovStat 1 Above 2.0 NA 1.0 NA
## IncidIPV:PovStat 0 Below 1.0 NA 2.0 NA
## IncidIPV:PovStat 1 Below 2.0 NA 2.0 NA
## IncidIPV:BMIcat2 0 nonObese 1.0 NA NA 1.0
## IncidIPV:BMIcat2 1 nonObese 2.0 NA NA 1.0
## IncidIPV:BMIcat2 0 Obese 1.0 NA NA 2.0
## IncidIPV:BMIcat2 1 Obese 2.0 NA NA 2.0
## Sex:PovStat Women Above NA 2.0 1.0 NA
## Sex:PovStat Men Above NA 1.0 1.0 NA
## Sex:PovStat Women Below NA 2.0 2.0 NA
## Sex:PovStat Men Below NA 1.0 2.0 NA
## PovStat:BMIcat2 Above nonObese NA NA 1.0 1.0
## PovStat:BMIcat2 Below nonObese NA NA 2.0 1.0
## PovStat:BMIcat2 Above Obese NA NA 1.0 2.0
## PovStat:BMIcat2 Below Obese NA NA 2.0 2.0
## IncidIPV:PovStat:BMIcat2 0 Above nonObese 1.0 NA 1.0 1.0
## IncidIPV:PovStat:BMIcat2 1 Above nonObese 2.0 NA 1.0 1.0
## IncidIPV:PovStat:BMIcat2 0 Below nonObese 1.0 NA 2.0 1.0
## IncidIPV:PovStat:BMIcat2 1 Below nonObese 2.0 NA 2.0 1.0
## IncidIPV:PovStat:BMIcat2 0 Above Obese 1.0 NA 1.0 2.0
## IncidIPV:PovStat:BMIcat2 1 Above Obese 2.0 NA 1.0 2.0
## IncidIPV:PovStat:BMIcat2 0 Below Obese 1.0 NA 2.0 2.0
## IncidIPV:PovStat:BMIcat2 1 Below Obese 2.0 NA 2.0 2.0
## Estimate Standard Error DF
## IncidIPV 0 119 2.19 51.4
## IncidIPV 1 123 3.45 56.5
## Sex Women 121 2.77 54.4
## Sex Men 121 3.09 59.0
## PovStat Above 118 2.36 50.1
## PovStat Below 124 3.35 57.2
## BMIcat2 nonObese 114 3.06 66.1
## BMIcat2 Obese 128 2.60 58.7
## IncidIPV:Sex 0 Women 118 3.13 49.9
## IncidIPV:Sex 1 Women 125 4.54 56.2
## IncidIPV:Sex 0 Men 120 3.12 53.9
## IncidIPV:Sex 1 Men 122 5.33 60.9
## IncidIPV:PovStat 0 Above 117 2.73 49.5
## IncidIPV:PovStat 1 Above 120 3.88 50.3
## IncidIPV:PovStat 0 Below 121 3.42 52.5
## IncidIPV:PovStat 1 Below 126 5.80 60.1
## IncidIPV:BMIcat2 0 nonObese 115 2.72 61.0
## IncidIPV:BMIcat2 1 nonObese 113 5.45 66.1
## IncidIPV:BMIcat2 0 Obese 123 3.30 65.6
## IncidIPV:BMIcat2 1 Obese 133 4.10 57.6
## Sex:PovStat Women Above 117 3.17 51.3
## Sex:PovStat Men Above 120 3.69 54.4
## Sex:PovStat Women Below 125 4.42 53.6
## Sex:PovStat Men Below 122 4.79 60.8
## PovStat:BMIcat2 Above nonObese 113 3.16 72.2
## PovStat:BMIcat2 Below nonObese 115 5.16 61.0
## PovStat:BMIcat2 Above Obese 124 3.14 71.6
## PovStat:BMIcat2 Below Obese 132 4.22 55.3
## IncidIPV:PovStat:BMIcat2 0 Above nonObese 112 3.71 70.4
## IncidIPV:PovStat:BMIcat2 1 Above nonObese 114 5.13 73.3
## IncidIPV:PovStat:BMIcat2 0 Below nonObese 118 3.95 53.8
## IncidIPV:PovStat:BMIcat2 1 Below nonObese 113 9.52 62.3
## IncidIPV:PovStat:BMIcat2 0 Above Obese 121 3.56 61.2
## IncidIPV:PovStat:BMIcat2 1 Above Obese 126 5.24 77.1
## IncidIPV:PovStat:BMIcat2 0 Below Obese 124 5.59 67.5
## IncidIPV:PovStat:BMIcat2 1 Below Obese 140 6.29 48.2
## t-value Lower CI Upper CI
## IncidIPV 0 54.5 114.7 124
## IncidIPV 1 35.6 116.2 130
## Sex Women 43.8 115.7 127
## Sex Men 39.2 114.8 127
## PovStat Above 50.2 113.6 123
## PovStat Below 37.0 117.2 131
## BMIcat2 nonObese 37.3 108.1 120
## BMIcat2 Obese 49.3 122.9 133
## IncidIPV:Sex 0 Women 37.7 111.7 124
## IncidIPV:Sex 1 Women 27.5 115.5 134
## IncidIPV:Sex 0 Men 38.5 114.0 126
## IncidIPV:Sex 1 Men 22.8 111.0 132
## IncidIPV:PovStat 0 Above 42.8 111.4 122
## IncidIPV:PovStat 1 Above 30.9 112.0 128
## IncidIPV:PovStat 0 Below 35.4 114.5 128
## IncidIPV:PovStat 1 Below 21.8 114.9 138
## IncidIPV:BMIcat2 0 nonObese 42.4 109.8 121
## IncidIPV:BMIcat2 1 nonObese 20.8 102.3 124
## IncidIPV:BMIcat2 0 Obese 37.2 116.4 130
## IncidIPV:BMIcat2 1 Obese 32.5 124.9 141
## Sex:PovStat Women Above 37.0 110.8 124
## Sex:PovStat Men Above 32.4 112.2 127
## Sex:PovStat Women Below 28.4 116.6 134
## Sex:PovStat Men Below 25.5 112.8 132
## PovStat:BMIcat2 Above nonObese 35.8 106.7 119
## PovStat:BMIcat2 Below nonObese 22.4 105.1 126
## PovStat:BMIcat2 Above Obese 39.4 117.5 130
## PovStat:BMIcat2 Below Obese 31.4 123.9 141
## IncidIPV:PovStat:BMIcat2 0 Above nonObese 30.3 104.9 120
## IncidIPV:PovStat:BMIcat2 1 Above nonObese 22.1 103.5 124
## IncidIPV:PovStat:BMIcat2 0 Below nonObese 29.9 110.3 126
## IncidIPV:PovStat:BMIcat2 1 Below nonObese 11.8 93.6 132
## IncidIPV:PovStat:BMIcat2 0 Above Obese 34.1 114.3 129
## IncidIPV:PovStat:BMIcat2 1 Above Obese 24.0 115.6 136
## IncidIPV:PovStat:BMIcat2 0 Below Obese 22.3 113.3 136
## IncidIPV:PovStat:BMIcat2 1 Below Obese 22.3 127.6 153
## p-value
## IncidIPV 0 <2e-16
## IncidIPV 1 <2e-16
## Sex Women <2e-16
## Sex Men <2e-16
## PovStat Above <2e-16
## PovStat Below <2e-16
## BMIcat2 nonObese <2e-16
## BMIcat2 Obese <2e-16
## IncidIPV:Sex 0 Women <2e-16
## IncidIPV:Sex 1 Women <2e-16
## IncidIPV:Sex 0 Men <2e-16
## IncidIPV:Sex 1 Men <2e-16
## IncidIPV:PovStat 0 Above <2e-16
## IncidIPV:PovStat 1 Above <2e-16
## IncidIPV:PovStat 0 Below <2e-16
## IncidIPV:PovStat 1 Below <2e-16
## IncidIPV:BMIcat2 0 nonObese <2e-16
## IncidIPV:BMIcat2 1 nonObese <2e-16
## IncidIPV:BMIcat2 0 Obese <2e-16
## IncidIPV:BMIcat2 1 Obese <2e-16
## Sex:PovStat Women Above <2e-16
## Sex:PovStat Men Above <2e-16
## Sex:PovStat Women Below <2e-16
## Sex:PovStat Men Below <2e-16
## PovStat:BMIcat2 Above nonObese <2e-16
## PovStat:BMIcat2 Below nonObese <2e-16
## PovStat:BMIcat2 Above Obese <2e-16
## PovStat:BMIcat2 Below Obese <2e-16
## IncidIPV:PovStat:BMIcat2 0 Above nonObese <2e-16
## IncidIPV:PovStat:BMIcat2 1 Above nonObese <2e-16
## IncidIPV:PovStat:BMIcat2 0 Below nonObese <2e-16
## IncidIPV:PovStat:BMIcat2 1 Below nonObese <2e-16
## IncidIPV:PovStat:BMIcat2 0 Above Obese <2e-16
## IncidIPV:PovStat:BMIcat2 1 Above Obese <2e-16
## IncidIPV:PovStat:BMIcat2 0 Below Obese <2e-16
## IncidIPV:PovStat:BMIcat2 1 Below Obese <2e-16
##
## Differences of LSMEANS:
## Estimate Standard Error
## IncidIPV 0-1 -4.0 4.108
## Sex Women-Men 0.3 4.228
## PovStat Above-Below -5.5 4.126
## BMIcat2 nonObese-Obese -13.8 3.956
## IncidIPV:Sex 0 Women- 1 Women -6.6 5.484
## IncidIPV:Sex 0 Women- 0 Men -2.2 4.472
## IncidIPV:Sex 0 Women- 1 Men -3.7 6.160
## IncidIPV:Sex 1 Women- 0 Men 4.4 5.617
## IncidIPV:Sex 1 Women- 1 Men 2.9 7.094
## IncidIPV:Sex 0 Men- 1 Men -1.5 6.185
## IncidIPV:PovStat 0 Above- 1 Above -3.0 4.767
## IncidIPV:PovStat 0 Above- 0 Below -4.5 4.385
## IncidIPV:PovStat 0 Above- 1 Below -9.6 6.368
## IncidIPV:PovStat 1 Above- 0 Below -1.5 5.221
## IncidIPV:PovStat 1 Above- 1 Below -6.6 7.046
## IncidIPV:PovStat 0 Below- 1 Below -5.1 6.771
## IncidIPV:BMIcat2 0 nonObese- 1 nonObese 2.1 6.065
## IncidIPV:BMIcat2 0 nonObese- 0 Obese -7.7 4.180
## IncidIPV:BMIcat2 0 nonObese- 1 Obese -17.9 4.892
## IncidIPV:BMIcat2 1 nonObese- 0 Obese -9.8 6.413
## IncidIPV:BMIcat2 1 nonObese- 1 Obese -20.0 6.736
## IncidIPV:BMIcat2 0 Obese- 1 Obese -10.2 5.339
## Sex:PovStat Women Above- Men Above -2.4 5.006
## Sex:PovStat Women Above- Women Below -8.3 5.336
## Sex:PovStat Women Above- Men Below -5.2 5.721
## Sex:PovStat Men Above- Women Below -5.9 6.088
## Sex:PovStat Men Above- Men Below -2.8 5.917
## Sex:PovStat Women Below- Men Below 3.1 6.333
## PovStat:BMIcat2 Above nonObese- Below nonObese -2.4 5.974
## PovStat:BMIcat2 Above nonObese- Above Obese -10.7 4.181
## PovStat:BMIcat2 Above nonObese- Below Obese -19.4 5.240
## PovStat:BMIcat2 Below nonObese- Above Obese -8.3 6.156
## PovStat:BMIcat2 Below nonObese- Below Obese -17.0 6.624
## PovStat:BMIcat2 Above Obese- Below Obese -8.7 5.330
## DF t-value Lower CI
## IncidIPV 0-1 55.1 -0.98 -12.27
## Sex Women-Men 58.8 0.08 -8.12
## PovStat Above-Below 54.6 -1.34 -13.82
## BMIcat2 nonObese-Obese 72.4 -3.50 -21.73
## IncidIPV:Sex 0 Women- 1 Women 53.8 -1.20 -17.60
## IncidIPV:Sex 0 Women- 0 Men 52.3 -0.50 -11.20
## IncidIPV:Sex 0 Women- 1 Men 57.8 -0.60 -16.03
## IncidIPV:Sex 1 Women- 0 Men 56.1 0.78 -6.88
## IncidIPV:Sex 1 Women- 1 Men 61.2 0.41 -11.28
## IncidIPV:Sex 0 Men- 1 Men 59.0 -0.24 -13.84
## IncidIPV:PovStat 0 Above- 1 Above 49.9 -0.62 -12.54
## IncidIPV:PovStat 0 Above- 0 Below 51.2 -1.02 -13.28
## IncidIPV:PovStat 0 Above- 1 Below 57.0 -1.50 -22.33
## IncidIPV:PovStat 1 Above- 0 Below 51.7 -0.29 -11.99
## IncidIPV:PovStat 1 Above- 1 Below 57.2 -0.94 -20.73
## IncidIPV:PovStat 0 Below- 1 Below 58.9 -0.75 -18.65
## IncidIPV:BMIcat2 0 nonObese- 1 nonObese 64.1 0.34 -10.03
## IncidIPV:BMIcat2 0 nonObese- 0 Obese 80.0 -1.85 -16.04
## IncidIPV:BMIcat2 0 nonObese- 1 Obese 57.5 -3.65 -27.67
## IncidIPV:BMIcat2 1 nonObese- 0 Obese 66.3 -1.53 -22.61
## IncidIPV:BMIcat2 1 nonObese- 1 Obese 70.4 -2.96 -33.40
## IncidIPV:BMIcat2 0 Obese- 1 Obese 62.6 -1.90 -20.83
## Sex:PovStat Women Above- Men Above 55.8 -0.48 -12.44
## Sex:PovStat Women Above- Women Below 51.1 -1.56 -19.01
## Sex:PovStat Women Above- Men Below 57.6 -0.91 -16.66
## Sex:PovStat Men Above- Women Below 56.0 -0.97 -18.08
## Sex:PovStat Men Above- Men Below 57.6 -0.47 -14.64
## Sex:PovStat Women Below- Men Below 57.6 0.49 -9.59
## PovStat:BMIcat2 Above nonObese- Below nonObese 61.6 -0.41 -14.38
## PovStat:BMIcat2 Above nonObese- Above Obese 105.2 -2.57 -19.02
## PovStat:BMIcat2 Above nonObese- Below Obese 59.8 -3.70 -29.87
## PovStat:BMIcat2 Below nonObese- Above Obese 64.6 -1.35 -20.59
## PovStat:BMIcat2 Below nonObese- Below Obese 60.2 -2.56 -30.20
## PovStat:BMIcat2 Above Obese- Below Obese 62.7 -1.62 -19.31
## Upper CI p-value
## IncidIPV 0-1 4.198 0.330
## Sex Women-Men 8.799 0.936
## PovStat Above-Below 2.722 0.184
## BMIcat2 nonObese-Obese -5.954 8e-04
## IncidIPV:Sex 0 Women- 1 Women 4.392 0.234
## IncidIPV:Sex 0 Women- 0 Men 6.743 0.620
## IncidIPV:Sex 0 Women- 1 Men 8.637 0.551
## IncidIPV:Sex 1 Women- 0 Men 15.625 0.440
## IncidIPV:Sex 1 Women- 1 Men 17.092 0.683
## IncidIPV:Sex 0 Men- 1 Men 10.911 0.814
## IncidIPV:PovStat 0 Above- 1 Above 6.612 0.537
## IncidIPV:PovStat 0 Above- 0 Below 4.324 0.312
## IncidIPV:PovStat 0 Above- 1 Below 3.168 0.138
## IncidIPV:PovStat 1 Above- 0 Below 8.963 0.773
## IncidIPV:PovStat 1 Above- 1 Below 7.488 0.351
## IncidIPV:PovStat 0 Below- 1 Below 8.444 0.454
## IncidIPV:BMIcat2 0 nonObese- 1 nonObese 14.203 0.732
## IncidIPV:BMIcat2 0 nonObese- 0 Obese 0.598 0.069
## IncidIPV:BMIcat2 0 nonObese- 1 Obese -8.081 6e-04
## IncidIPV:BMIcat2 1 nonObese- 0 Obese 2.997 0.131
## IncidIPV:BMIcat2 1 nonObese- 1 Obese -6.528 0.004
## IncidIPV:BMIcat2 0 Obese- 1 Obese 0.515 0.062
## Sex:PovStat Women Above- Men Above 7.617 0.632
## Sex:PovStat Women Above- Women Below 2.411 0.126
## Sex:PovStat Women Above- Men Below 6.245 0.366
## Sex:PovStat Men Above- Women Below 6.308 0.338
## Sex:PovStat Men Above- Men Below 9.048 0.638
## Sex:PovStat Women Below- Men Below 15.770 0.627
## PovStat:BMIcat2 Above nonObese- Below nonObese 9.507 0.685
## PovStat:BMIcat2 Above nonObese- Above Obese -2.439 0.012
## PovStat:BMIcat2 Above nonObese- Below Obese -8.906 5e-04
## PovStat:BMIcat2 Below nonObese- Above Obese 4.003 0.183
## PovStat:BMIcat2 Below nonObese- Below Obese -3.703 0.013
## PovStat:BMIcat2 Above Obese- Below Obese 1.992 0.109
##
## Final model:
## lme4::lmer(formula = adjSBP ~ Age + IncidIPV + Sex + PovStat +
## BMIcat2 + (1 | HNDid) + Age:IncidIPV + Age:Sex + Age:PovStat +
## Age:BMIcat2 + IncidIPV:Sex + IncidIPV:PovStat + IncidIPV:BMIcat2 +
## Sex:PovStat + PovStat:BMIcat2 + Age:IncidIPV:Sex + Age:IncidIPV:PovStat +
## Age:IncidIPV:BMIcat2 + Age:Sex:PovStat + Age:PovStat:BMIcat2 +
## IncidIPV:PovStat:BMIcat2 + Age:IncidIPV:PovStat:BMIcat2,
## data = IncidIPVbp, REML = reml, contrasts = l)
plot(st)
Systolic Blood Pressure Final Model (BMI as categorical)
(mm2 = lmer(adjSBP ~Age + IncidIPV + Sex + PovStat + BMIcat2 + Age:IncidIPV + Age:Sex + Age:PovStat + Age:BMIcat2 + IncidIPV:Sex + IncidIPV:PovStat + IncidIPV:BMIcat2 + Sex:PovStat + PovStat:BMIcat2 + Age:IncidIPV:Sex + Age:IncidIPV:PovStat + Age:IncidIPV:BMIcat2 + Age:Sex:PovStat + Age:PovStat:BMIcat2 + IncidIPV:PovStat:BMIcat2 + Age:IncidIPV:PovStat:BMIcat2 + (Age|HNDid) + (1|subclass) ,data = IncidIPVbp))
## Warning: number of observations <= rank(Z); variance-covariance matrix
## will be unidentifiable
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: adjSBP ~ Age + IncidIPV + Sex + PovStat + BMIcat2 + Age:IncidIPV + Age:Sex + Age:PovStat + Age:BMIcat2 + IncidIPV:Sex + IncidIPV:PovStat + IncidIPV:BMIcat2 + Sex:PovStat + PovStat:BMIcat2 + Age:IncidIPV:Sex + Age:IncidIPV:PovStat + Age:IncidIPV:BMIcat2 + Age:Sex:PovStat + Age:PovStat:BMIcat2 + IncidIPV:PovStat:BMIcat2 + Age:IncidIPV:PovStat:BMIcat2 + (Age | HNDid) + (1 | subclass)
## Data: IncidIPVbp
## REML criterion at convergence: 1041
## Random effects:
## Groups Name Std.Dev. Corr
## HNDid (Intercept) 10.74305
## Age 0.09616 -1.00
## subclass (Intercept) 0.00019
## Residual 11.01309
## Number of obs: 138, groups: HNDid, 70; subclass, 46
## Fixed Effects:
## (Intercept)
## 119.128
## Age
## 1.359
## IncidIPV1
## -2.318
## SexMen
## -10.555
## PovStatBelow
## 7.644
## BMIcat2Obese
## 8.562
## Age:IncidIPV1
## -0.936
## Age:SexMen
## -2.284
## Age:PovStatBelow
## -0.139
## Age:BMIcat2Obese
## -0.068
## IncidIPV1:SexMen
## 13.106
## IncidIPV1:PovStatBelow
## 5.707
## IncidIPV1:BMIcat2Obese
## 9.662
## SexMen:PovStatBelow
## 11.266
## PovStatBelow:BMIcat2Obese
## 6.856
## Age:IncidIPV1:SexMen
## 2.705
## Age:IncidIPV1:PovStatBelow
## 1.908
## Age:IncidIPV1:BMIcat2Obese
## 0.962
## Age:SexMen:PovStatBelow
## 2.434
## Age:PovStatBelow:BMIcat2Obese
## 1.411
## IncidIPV1:PovStatBelow:BMIcat2Obese
## -25.797
## Age:IncidIPV1:PovStatBelow:BMIcat2Obese
## -6.525
summary(mm2)
## Warning: number of observations <= rank(Z); variance-covariance matrix will be unidentifiable
## Warning: number of observations <= rank(Z); variance-covariance matrix will be unidentifiable
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: adjSBP ~ Age + IncidIPV + Sex + PovStat + BMIcat2 + Age:IncidIPV + Age:Sex + Age:PovStat + Age:BMIcat2 + IncidIPV:Sex + IncidIPV:PovStat + IncidIPV:BMIcat2 + Sex:PovStat + PovStat:BMIcat2 + Age:IncidIPV:Sex + Age:IncidIPV:PovStat + Age:IncidIPV:BMIcat2 + Age:Sex:PovStat + Age:PovStat:BMIcat2 + IncidIPV:PovStat:BMIcat2 + Age:IncidIPV:PovStat:BMIcat2 + (Age | HNDid) + (1 | subclass)
## Data: IncidIPVbp
##
## REML criterion at convergence: 1041
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 1.15e+02 10.74305
## Age 9.25e-03 0.09616 -1.00
## subclass (Intercept) 3.60e-08 0.00019
## Residual 1.21e+02 11.01309
## Number of obs: 138, groups: HNDid, 70; subclass, 46
##
## Fixed effects:
## Estimate Std. Error df
## (Intercept) 119.128 5.541 16.700
## Age 1.359 0.655 62.700
## IncidIPV1 -2.318 13.545 59.300
## SexMen -10.555 6.834 25.100
## PovStatBelow 7.644 9.104 17.700
## BMIcat2Obese 8.562 7.010 26.000
## Age:IncidIPV1 -0.936 1.450 104.400
## Age:SexMen -2.284 0.754 78.500
## Age:PovStatBelow -0.139 0.895 62.700
## Age:BMIcat2Obese -0.068 0.755 57.100
## IncidIPV1:SexMen 13.106 12.974 57.100
## IncidIPV1:PovStatBelow 5.707 15.197 31.300
## IncidIPV1:BMIcat2Obese 9.662 15.537 64.100
## SexMen:PovStatBelow 11.266 10.819 30.000
## PovStatBelow:BMIcat2Obese 6.856 11.870 20.200
## Age:IncidIPV1:SexMen 2.705 1.253 109.100
## Age:IncidIPV1:PovStatBelow 1.908 1.752 107.200
## Age:IncidIPV1:BMIcat2Obese 0.962 1.594 104.700
## Age:SexMen:PovStatBelow 2.434 1.055 97.900
## Age:PovStatBelow:BMIcat2Obese 1.411 1.199 53.300
## IncidIPV1:PovStatBelow:BMIcat2Obese -25.797 20.512 29.900
## Age:IncidIPV1:PovStatBelow:BMIcat2Obese -6.525 2.282 96.000
## t value Pr(>|t|)
## (Intercept) 21.50 1.3e-13
## Age 2.08 0.0420
## IncidIPV1 -0.17 0.8647
## SexMen -1.54 0.1350
## PovStatBelow 0.84 0.4123
## BMIcat2Obese 1.22 0.2329
## Age:IncidIPV1 -0.65 0.5198
## Age:SexMen -3.03 0.0033
## Age:PovStatBelow -0.15 0.8774
## Age:BMIcat2Obese -0.09 0.9285
## IncidIPV1:SexMen 1.01 0.3167
## IncidIPV1:PovStatBelow 0.38 0.7098
## IncidIPV1:BMIcat2Obese 0.62 0.5363
## SexMen:PovStatBelow 1.04 0.3060
## PovStatBelow:BMIcat2Obese 0.58 0.5699
## Age:IncidIPV1:SexMen 2.16 0.0331
## Age:IncidIPV1:PovStatBelow 1.09 0.2786
## Age:IncidIPV1:BMIcat2Obese 0.60 0.5473
## Age:SexMen:PovStatBelow 2.31 0.0232
## Age:PovStatBelow:BMIcat2Obese 1.18 0.2447
## IncidIPV1:PovStatBelow:BMIcat2Obese -1.26 0.2183
## Age:IncidIPV1:PovStatBelow:BMIcat2Obese -2.86 0.0052
##
## Correlation matrix not shown by default, as p = 22 > 20.
## Use print(x, correlation=TRUE) or
## vcov(x) if you need it
SBP and IPV: Age x IncidIPV x Sex
## Age Sex IncidIPV adjSBP PovStat BMIcat2 hat
## 1 30 Women 0 0 0.3913 0.4928 174.0
## 2 31 Women 0 0 0.3913 0.4928 175.5
## 3 32 Women 0 0 0.3913 0.4928 177.1
## 4 33 Women 0 0 0.3913 0.4928 178.6
## 5 34 Women 0 0 0.3913 0.4928 180.1
## 6 35 Women 0 0 0.3913 0.4928 181.7
SBP and IPV: Age x IncidIPV x PovStat x BMI
## Age IncidIPV PovStat BMIcat2 adjSBP Sex hat
## 1 30 0 Above nonObese 0 0.4928 120.9
## 2 31 0 Above nonObese 0 0.4928 121.2
## 3 32 0 Above nonObese 0 0.4928 121.4
## 4 33 0 Above nonObese 0 0.4928 121.6
## 5 34 0 Above nonObese 0 0.4928 121.9
## 6 35 0 Above nonObese 0 0.4928 122.1
SBP and IPV: Age x IncidIPV x PovStat x BMI (Plots broken up by poverty status)
## Age IncidIPV PovStat BMIcat2 adjSBP Sex hat
## 1 30 0 Above nonObese 0 0.4928 120.9
## 2 31 0 Above nonObese 0 0.4928 121.2
## 3 32 0 Above nonObese 0 0.4928 121.4
## 4 33 0 Above nonObese 0 0.4928 121.6
## 5 34 0 Above nonObese 0 0.4928 121.9
## 6 35 0 Above nonObese 0 0.4928 122.1
Diastolic Blood Pressure (BMI as categorical)
(mm4 = lmer(adjDBP ~(Age + IncidIPV + Sex + PovStat + BMIcat2)^5 + (Age|HNDid) + (1|subclass) ,IncidIPVbp))
## Warning: number of observations <= rank(Z); variance-covariance matrix
## will be unidentifiable
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: adjDBP ~ (Age + IncidIPV + Sex + PovStat + BMIcat2)^5 + (Age | HNDid) + (1 | subclass)
## Data: IncidIPVbp
## REML criterion at convergence: 859.3
## Random effects:
## Groups Name Std.Dev. Corr
## HNDid (Intercept) 5.636
## Age 0.121 -1.00
## subclass (Intercept) 1.179
## Residual 6.333
## Number of obs: 138, groups: HNDid, 70; subclass, 46
## Fixed Effects:
## (Intercept)
## 70.8643
## Age
## 0.5661
## IncidIPV1
## 12.1346
## SexMen
## 1.0032
## PovStatBelow
## 4.9500
## BMIcat2Obese
## 5.5728
## Age:IncidIPV1
## 1.3441
## Age:SexMen
## -0.7740
## Age:PovStatBelow
## 0.0218
## Age:BMIcat2Obese
## -0.5869
## IncidIPV1:SexMen
## -18.3981
## IncidIPV1:PovStatBelow
## 11.2074
## IncidIPV1:BMIcat2Obese
## -14.6964
## SexMen:PovStatBelow
## 0.4771
## SexMen:BMIcat2Obese
## -4.6951
## PovStatBelow:BMIcat2Obese
## -2.2832
## Age:IncidIPV1:SexMen
## -1.5549
## Age:IncidIPV1:PovStatBelow
## 0.2967
## Age:IncidIPV1:BMIcat2Obese
## -1.3722
## Age:SexMen:PovStatBelow
## 0.7028
## Age:SexMen:BMIcat2Obese
## 0.6419
## Age:PovStatBelow:BMIcat2Obese
## 0.7831
## IncidIPV1:SexMen:PovStatBelow
## 11.6144
## IncidIPV1:SexMen:BMIcat2Obese
## 33.4345
## IncidIPV1:PovStatBelow:BMIcat2Obese
## -24.9780
## SexMen:PovStatBelow:BMIcat2Obese
## -12.9211
## Age:IncidIPV1:SexMen:PovStatBelow
## 2.9659
## Age:IncidIPV1:SexMen:BMIcat2Obese
## 3.0754
## Age:IncidIPV1:PovStatBelow:BMIcat2Obese
## -4.2474
## Age:SexMen:PovStatBelow:BMIcat2Obese
## -1.6385
## IncidIPV1:SexMen:PovStatBelow:BMIcat2Obese
## 12.9543
## Age:IncidIPV1:SexMen:PovStatBelow:BMIcat2Obese
## 0.5618
(st = step(mm4))
## Warning: number of observations <= rank(Z); variance-covariance matrix will be unidentifiable
## Warning: number of observations <= rank(Z); variance-covariance matrix will be unidentifiable
## Random term (Age | HNDid) was eliminated because of having correlation +-1 or NaN
## Warning: number of observations <= rank(Z); variance-covariance matrix will be unidentifiable
## Warning: number of observations <= rank(Z); variance-covariance matrix will be unidentifiable
##
## Random effects:
## Chi.sq Chi.DF elim.num p.value
## (Age + 0 | HNDid) 0.00 1 1 1.0000
## (1 | subclass) 0.19 1 2 0.6606
## (1 | HNDid) 11.42 1 kept 0.0007
##
## Fixed effects:
## Sum Sq Mean Sq NumDF DenDF F.value
## Age:IncidIPV:Sex:PovStat:BMIcat2 3.0422 3.0422 1 98.84 0.0727
## IncidIPV:Sex:PovStat:BMIcat2 5.4796 5.4796 1 90.37 0.1317
## Age:Sex:PovStat:BMIcat2 50.3662 50.3662 1 107.10 1.2299
## Sex:PovStat:BMIcat2 0.4912 0.4912 1 91.13 0.1937
## Age:IncidIPV:Sex:PovStat 252.8939 252.8939 1 108.10 2.6539
## IncidIPV:Sex:PovStat 4.2955 4.2955 1 56.43 0.6814
## Age:Sex:PovStat 128.3729 128.3729 1 97.14 2.3401
## Sex:PovStat 0.8682 0.8682 1 51.46 1.7706
## Age 375.0425 375.0425 1 97.48 12.6682
## IncidIPV 113.4509 113.4509 1 68.68 7.5730
## Sex 6.0040 6.0040 1 90.78 0.0072
## PovStat 267.2578 267.2578 1 82.01 5.4394
## BMIcat2 204.7792 204.7792 1 74.36 0.7200
## Age:IncidIPV 50.3611 50.3611 1 97.48 2.6657
## Age:Sex 6.8708 6.8708 1 109.93 0.0033
## Age:PovStat 194.8877 194.8877 1 104.17 1.3372
## Age:BMIcat2 83.7275 83.7275 1 101.39 5.7278
## IncidIPV:Sex 8.6957 8.6957 1 90.78 0.2237
## IncidIPV:PovStat 443.6150 443.6150 1 82.01 2.3879
## IncidIPV:BMIcat2 15.7783 15.7783 1 74.36 2.1464
## Sex:BMIcat2 207.0763 207.0763 1 94.82 1.9573
## PovStat:BMIcat2 79.6665 79.6665 1 86.02 6.5079
## Age:IncidIPV:Sex 154.0713 154.0713 1 109.93 0.7616
## Age:IncidIPV:PovStat 12.4582 12.4582 1 104.17 0.1425
## Age:IncidIPV:BMIcat2 100.3472 100.3472 1 101.39 6.0432
## Age:Sex:BMIcat2 47.0881 47.0881 1 110.91 5.3944
## Age:PovStat:BMIcat2 7.6851 7.6851 1 104.90 5.3781
## IncidIPV:Sex:BMIcat2 6.2230 6.2230 1 94.82 5.2930
## IncidIPV:PovStat:BMIcat2 22.2187 22.2187 1 86.02 3.6983
## Age:IncidIPV:Sex:BMIcat2 118.1944 118.1944 1 110.91 5.6254
## Age:IncidIPV:PovStat:BMIcat2 273.1737 273.1737 1 104.90 7.0224
## elim.num Pr(>F)
## Age:IncidIPV:Sex:PovStat:BMIcat2 1 0.7881
## IncidIPV:Sex:PovStat:BMIcat2 2 0.7175
## Age:Sex:PovStat:BMIcat2 3 0.2699
## Sex:PovStat:BMIcat2 4 0.6609
## Age:IncidIPV:Sex:PovStat 5 0.1062
## IncidIPV:Sex:PovStat 6 0.4126
## Age:Sex:PovStat 7 0.1293
## Sex:PovStat 8 0.1892
## Age kept 0.0006
## IncidIPV kept 0.0076
## Sex kept 0.9327
## PovStat kept 0.0221
## BMIcat2 kept 0.3989
## Age:IncidIPV kept 0.1058
## Age:Sex kept 0.9541
## Age:PovStat kept 0.2502
## Age:BMIcat2 kept 0.0185
## IncidIPV:Sex kept 0.6374
## IncidIPV:PovStat kept 0.1261
## IncidIPV:BMIcat2 kept 0.1471
## Sex:BMIcat2 kept 0.1651
## PovStat:BMIcat2 kept 0.0125
## Age:IncidIPV:Sex kept 0.3847
## Age:IncidIPV:PovStat kept 0.7065
## Age:IncidIPV:BMIcat2 kept 0.0157
## Age:Sex:BMIcat2 kept 0.0220
## Age:PovStat:BMIcat2 kept 0.0223
## IncidIPV:Sex:BMIcat2 kept 0.0236
## IncidIPV:PovStat:BMIcat2 kept 0.0578
## Age:IncidIPV:Sex:BMIcat2 kept 0.0194
## Age:IncidIPV:PovStat:BMIcat2 kept 0.0093
##
## Least squares means:
## IncidIPV Sex PovStat BMIcat2
## IncidIPV 0 1.0 NA NA NA
## IncidIPV 1 2.0 NA NA NA
## Sex Women NA 2.0 NA NA
## Sex Men NA 1.0 NA NA
## PovStat Above NA NA 1.0 NA
## PovStat Below NA NA 2.0 NA
## BMIcat2 nonObese NA NA NA 1.0
## BMIcat2 Obese NA NA NA 2.0
## IncidIPV:Sex 0 Women 1.0 2.0 NA NA
## IncidIPV:Sex 1 Women 2.0 2.0 NA NA
## IncidIPV:Sex 0 Men 1.0 1.0 NA NA
## IncidIPV:Sex 1 Men 2.0 1.0 NA NA
## IncidIPV:PovStat 0 Above 1.0 NA 1.0 NA
## IncidIPV:PovStat 1 Above 2.0 NA 1.0 NA
## IncidIPV:PovStat 0 Below 1.0 NA 2.0 NA
## IncidIPV:PovStat 1 Below 2.0 NA 2.0 NA
## IncidIPV:BMIcat2 0 nonObese 1.0 NA NA 1.0
## IncidIPV:BMIcat2 1 nonObese 2.0 NA NA 1.0
## IncidIPV:BMIcat2 0 Obese 1.0 NA NA 2.0
## IncidIPV:BMIcat2 1 Obese 2.0 NA NA 2.0
## Sex:BMIcat2 Women nonObese NA 2.0 NA 1.0
## Sex:BMIcat2 Men nonObese NA 1.0 NA 1.0
## Sex:BMIcat2 Women Obese NA 2.0 NA 2.0
## Sex:BMIcat2 Men Obese NA 1.0 NA 2.0
## PovStat:BMIcat2 Above nonObese NA NA 1.0 1.0
## PovStat:BMIcat2 Below nonObese NA NA 2.0 1.0
## PovStat:BMIcat2 Above Obese NA NA 1.0 2.0
## PovStat:BMIcat2 Below Obese NA NA 2.0 2.0
## IncidIPV:Sex:BMIcat2 0 Women nonObese 1.0 2.0 NA 1.0
## IncidIPV:Sex:BMIcat2 1 Women nonObese 2.0 2.0 NA 1.0
## IncidIPV:Sex:BMIcat2 0 Men nonObese 1.0 1.0 NA 1.0
## IncidIPV:Sex:BMIcat2 1 Men nonObese 2.0 1.0 NA 1.0
## IncidIPV:Sex:BMIcat2 0 Women Obese 1.0 2.0 NA 2.0
## IncidIPV:Sex:BMIcat2 1 Women Obese 2.0 2.0 NA 2.0
## IncidIPV:Sex:BMIcat2 0 Men Obese 1.0 1.0 NA 2.0
## IncidIPV:Sex:BMIcat2 1 Men Obese 2.0 1.0 NA 2.0
## IncidIPV:PovStat:BMIcat2 0 Above nonObese 1.0 NA 1.0 1.0
## IncidIPV:PovStat:BMIcat2 1 Above nonObese 2.0 NA 1.0 1.0
## IncidIPV:PovStat:BMIcat2 0 Below nonObese 1.0 NA 2.0 1.0
## IncidIPV:PovStat:BMIcat2 1 Below nonObese 2.0 NA 2.0 1.0
## IncidIPV:PovStat:BMIcat2 0 Above Obese 1.0 NA 1.0 2.0
## IncidIPV:PovStat:BMIcat2 1 Above Obese 2.0 NA 1.0 2.0
## IncidIPV:PovStat:BMIcat2 0 Below Obese 1.0 NA 2.0 2.0
## IncidIPV:PovStat:BMIcat2 1 Below Obese 2.0 NA 2.0 2.0
## Estimate Standard Error DF
## IncidIPV 0 72.6 1.22 44.5
## IncidIPV 1 78.1 2.10 52.8
## Sex Women 75.6 1.59 48.5
## Sex Men 75.1 1.77 52.5
## PovStat Above 72.7 1.38 44.8
## PovStat Below 78.0 1.95 52.3
## BMIcat2 nonObese 73.9 1.87 59.6
## BMIcat2 Obese 76.8 1.52 47.6
## IncidIPV:Sex 0 Women 72.6 1.74 41.6
## IncidIPV:Sex 1 Women 78.6 2.66 51.8
## IncidIPV:Sex 0 Men 72.6 1.73 48.9
## IncidIPV:Sex 1 Men 77.6 3.09 53.7
## IncidIPV:PovStat 0 Above 73.1 1.53 41.8
## IncidIPV:PovStat 1 Above 72.3 2.31 46.2
## IncidIPV:PovStat 0 Below 72.0 1.99 48.8
## IncidIPV:PovStat 1 Below 83.9 3.35 53.6
## IncidIPV:BMIcat2 0 nonObese 71.6 1.54 48.6
## IncidIPV:BMIcat2 1 nonObese 76.3 3.40 62.1
## IncidIPV:BMIcat2 0 Obese 73.5 1.82 54.1
## IncidIPV:BMIcat2 1 Obese 80.0 2.43 44.4
## Sex:BMIcat2 Women nonObese 73.8 2.54 61.1
## Sex:BMIcat2 Men nonObese 74.0 2.37 58.3
## Sex:BMIcat2 Women Obese 77.4 1.79 47.2
## Sex:BMIcat2 Men Obese 76.1 2.58 53.5
## PovStat:BMIcat2 Above nonObese 69.7 1.84 56.1
## PovStat:BMIcat2 Below nonObese 78.1 3.06 57.5
## PovStat:BMIcat2 Above Obese 75.7 1.94 55.0
## PovStat:BMIcat2 Below Obese 77.8 2.37 48.9
## IncidIPV:Sex:BMIcat2 0 Women nonObese 69.9 2.40 47.9
## IncidIPV:Sex:BMIcat2 1 Women nonObese 77.7 4.48 65.6
## IncidIPV:Sex:BMIcat2 0 Men nonObese 73.3 1.92 49.7
## IncidIPV:Sex:BMIcat2 1 Men nonObese 74.8 4.34 60.2
## IncidIPV:Sex:BMIcat2 0 Women Obese 75.2 2.41 47.3
## IncidIPV:Sex:BMIcat2 1 Women Obese 79.5 2.63 47.1
## IncidIPV:Sex:BMIcat2 0 Men Obese 71.8 2.73 70.4
## IncidIPV:Sex:BMIcat2 1 Men Obese 80.5 4.39 48.1
## IncidIPV:PovStat:BMIcat2 0 Above nonObese 70.6 2.10 57.0
## IncidIPV:PovStat:BMIcat2 1 Above nonObese 68.8 3.02 55.7
## IncidIPV:PovStat:BMIcat2 0 Below nonObese 72.6 2.20 42.9
## IncidIPV:PovStat:BMIcat2 1 Below nonObese 83.7 5.72 60.1
## IncidIPV:PovStat:BMIcat2 0 Above Obese 75.5 2.03 54.1
## IncidIPV:PovStat:BMIcat2 1 Above Obese 75.8 3.30 55.4
## IncidIPV:PovStat:BMIcat2 0 Below Obese 71.5 3.23 62.6
## IncidIPV:PovStat:BMIcat2 1 Below Obese 84.1 3.48 39.8
## t-value Lower CI Upper CI
## IncidIPV 0 59.6 70.1 75.0
## IncidIPV 1 37.1 73.9 82.3
## Sex Women 47.6 72.4 78.8
## Sex Men 42.4 71.5 78.6
## PovStat Above 52.5 69.9 75.5
## PovStat Below 40.0 74.1 81.9
## BMIcat2 nonObese 39.6 70.2 77.7
## BMIcat2 Obese 50.5 73.7 79.8
## IncidIPV:Sex 0 Women 41.8 69.1 76.1
## IncidIPV:Sex 1 Women 29.6 73.3 83.9
## IncidIPV:Sex 0 Men 41.9 69.1 76.0
## IncidIPV:Sex 1 Men 25.1 71.4 83.8
## IncidIPV:PovStat 0 Above 47.8 70.0 76.2
## IncidIPV:PovStat 1 Above 31.3 67.7 77.0
## IncidIPV:PovStat 0 Below 36.2 68.1 76.0
## IncidIPV:PovStat 1 Below 25.0 77.2 90.6
## IncidIPV:BMIcat2 0 nonObese 46.5 68.5 74.7
## IncidIPV:BMIcat2 1 nonObese 22.4 69.5 83.1
## IncidIPV:BMIcat2 0 Obese 40.3 69.9 77.2
## IncidIPV:BMIcat2 1 Obese 32.9 75.1 84.9
## Sex:BMIcat2 Women nonObese 29.1 68.7 78.9
## Sex:BMIcat2 Men nonObese 31.2 69.3 78.8
## Sex:BMIcat2 Women Obese 43.3 73.8 80.9
## Sex:BMIcat2 Men Obese 29.5 71.0 81.3
## PovStat:BMIcat2 Above nonObese 37.9 66.0 73.4
## PovStat:BMIcat2 Below nonObese 25.5 72.0 84.3
## PovStat:BMIcat2 Above Obese 39.0 71.8 79.6
## PovStat:BMIcat2 Below Obese 32.8 73.0 82.6
## IncidIPV:Sex:BMIcat2 0 Women nonObese 29.1 65.1 74.7
## IncidIPV:Sex:BMIcat2 1 Women nonObese 17.4 68.8 86.7
## IncidIPV:Sex:BMIcat2 0 Men nonObese 38.2 69.4 77.1
## IncidIPV:Sex:BMIcat2 1 Men nonObese 17.2 66.1 83.5
## IncidIPV:Sex:BMIcat2 0 Women Obese 31.2 70.4 80.1
## IncidIPV:Sex:BMIcat2 1 Women Obese 30.2 74.2 84.8
## IncidIPV:Sex:BMIcat2 0 Men Obese 26.3 66.4 77.3
## IncidIPV:Sex:BMIcat2 1 Men Obese 18.3 71.6 89.3
## IncidIPV:PovStat:BMIcat2 0 Above nonObese 33.6 66.4 74.8
## IncidIPV:PovStat:BMIcat2 1 Above nonObese 22.8 62.8 74.9
## IncidIPV:PovStat:BMIcat2 0 Below nonObese 32.9 68.1 77.0
## IncidIPV:PovStat:BMIcat2 1 Below nonObese 14.6 72.3 95.1
## IncidIPV:PovStat:BMIcat2 0 Above Obese 37.2 71.5 79.6
## IncidIPV:PovStat:BMIcat2 1 Above Obese 23.0 69.2 82.4
## IncidIPV:PovStat:BMIcat2 0 Below Obese 22.2 65.1 78.0
## IncidIPV:PovStat:BMIcat2 1 Below Obese 24.1 77.1 91.1
## p-value
## IncidIPV 0 <2e-16
## IncidIPV 1 <2e-16
## Sex Women <2e-16
## Sex Men <2e-16
## PovStat Above <2e-16
## PovStat Below <2e-16
## BMIcat2 nonObese <2e-16
## BMIcat2 Obese <2e-16
## IncidIPV:Sex 0 Women <2e-16
## IncidIPV:Sex 1 Women <2e-16
## IncidIPV:Sex 0 Men <2e-16
## IncidIPV:Sex 1 Men <2e-16
## IncidIPV:PovStat 0 Above <2e-16
## IncidIPV:PovStat 1 Above <2e-16
## IncidIPV:PovStat 0 Below <2e-16
## IncidIPV:PovStat 1 Below <2e-16
## IncidIPV:BMIcat2 0 nonObese <2e-16
## IncidIPV:BMIcat2 1 nonObese <2e-16
## IncidIPV:BMIcat2 0 Obese <2e-16
## IncidIPV:BMIcat2 1 Obese <2e-16
## Sex:BMIcat2 Women nonObese <2e-16
## Sex:BMIcat2 Men nonObese <2e-16
## Sex:BMIcat2 Women Obese <2e-16
## Sex:BMIcat2 Men Obese <2e-16
## PovStat:BMIcat2 Above nonObese <2e-16
## PovStat:BMIcat2 Below nonObese <2e-16
## PovStat:BMIcat2 Above Obese <2e-16
## PovStat:BMIcat2 Below Obese <2e-16
## IncidIPV:Sex:BMIcat2 0 Women nonObese <2e-16
## IncidIPV:Sex:BMIcat2 1 Women nonObese <2e-16
## IncidIPV:Sex:BMIcat2 0 Men nonObese <2e-16
## IncidIPV:Sex:BMIcat2 1 Men nonObese <2e-16
## IncidIPV:Sex:BMIcat2 0 Women Obese <2e-16
## IncidIPV:Sex:BMIcat2 1 Women Obese <2e-16
## IncidIPV:Sex:BMIcat2 0 Men Obese <2e-16
## IncidIPV:Sex:BMIcat2 1 Men Obese <2e-16
## IncidIPV:PovStat:BMIcat2 0 Above nonObese <2e-16
## IncidIPV:PovStat:BMIcat2 1 Above nonObese <2e-16
## IncidIPV:PovStat:BMIcat2 0 Below nonObese <2e-16
## IncidIPV:PovStat:BMIcat2 1 Below nonObese <2e-16
## IncidIPV:PovStat:BMIcat2 0 Above Obese <2e-16
## IncidIPV:PovStat:BMIcat2 1 Above Obese <2e-16
## IncidIPV:PovStat:BMIcat2 0 Below Obese <2e-16
## IncidIPV:PovStat:BMIcat2 1 Below Obese <2e-16
##
## Differences of LSMEANS:
## Estimate Standard Error
## IncidIPV 0-1 -5.5 2.4315
## Sex Women-Men 0.5 2.3220
## PovStat Above-Below -5.3 2.3470
## BMIcat2 nonObese-Obese -2.8 2.3821
## IncidIPV:Sex 0 Women- 1 Women -6.0 3.1756
## IncidIPV:Sex 0 Women- 0 Men 0.0 2.4672
## IncidIPV:Sex 0 Women- 1 Men -5.1 3.5420
## IncidIPV:Sex 1 Women- 0 Men 6.0 3.1720
## IncidIPV:Sex 1 Women- 1 Men 1.0 3.9344
## IncidIPV:Sex 0 Men- 1 Men -5.1 3.5388
## IncidIPV:PovStat 0 Above- 1 Above 0.8 2.7687
## IncidIPV:PovStat 0 Above- 0 Below 1.0 2.5764
## IncidIPV:PovStat 0 Above- 1 Below -10.8 3.6820
## IncidIPV:PovStat 1 Above- 0 Below 0.3 3.0470
## IncidIPV:PovStat 1 Above- 1 Below -11.6 3.9238
## IncidIPV:PovStat 0 Below- 1 Below -11.8 3.8956
## IncidIPV:BMIcat2 0 nonObese- 1 nonObese -4.7 3.7325
## IncidIPV:BMIcat2 0 nonObese- 0 Obese -1.9 2.3380
## IncidIPV:BMIcat2 0 nonObese- 1 Obese -8.4 2.8782
## IncidIPV:BMIcat2 1 nonObese- 0 Obese 2.7 3.8585
## IncidIPV:BMIcat2 1 nonObese- 1 Obese -3.7 4.1509
## IncidIPV:BMIcat2 0 Obese- 1 Obese -6.4 3.0399
## Sex:BMIcat2 Women nonObese- Men nonObese -0.2 3.1997
## Sex:BMIcat2 Women nonObese- Women Obese -3.5 3.0305
## Sex:BMIcat2 Women nonObese- Men Obese -2.3 3.6306
## Sex:BMIcat2 Men nonObese- Women Obese -3.3 2.9917
## Sex:BMIcat2 Men nonObese- Men Obese -2.1 3.4799
## Sex:BMIcat2 Women Obese- Men Obese 1.2 3.2389
## PovStat:BMIcat2 Above nonObese- Below nonObese -8.4 3.4065
## PovStat:BMIcat2 Above nonObese- Above Obese -6.0 2.5701
## PovStat:BMIcat2 Above nonObese- Below Obese -8.1 3.0243
## PovStat:BMIcat2 Below nonObese- Above Obese 2.5 3.6358
## PovStat:BMIcat2 Below nonObese- Below Obese 0.3 3.8568
## PovStat:BMIcat2 Above Obese- Below Obese -2.1 3.0895
## DF t-value Lower CI
## IncidIPV 0-1 50.6 -2.28 -10.43
## Sex Women-Men 50.8 0.21 -4.17
## PovStat Above-Below 48.7 -2.25 -9.99
## BMIcat2 nonObese-Obese 58.9 -1.18 -7.59
## IncidIPV:Sex 0 Women- 1 Women 48.5 -1.90 -12.41
## IncidIPV:Sex 0 Women- 0 Men 45.6 0.00 -4.96
## IncidIPV:Sex 0 Women- 1 Men 50.4 -1.43 -12.17
## IncidIPV:Sex 1 Women- 0 Men 50.9 1.90 -0.33
## IncidIPV:Sex 1 Women- 1 Men 53.0 0.25 -6.92
## IncidIPV:Sex 0 Men- 1 Men 52.5 -1.43 -12.17
## IncidIPV:PovStat 0 Above- 1 Above 44.8 0.27 -4.83
## IncidIPV:PovStat 0 Above- 0 Below 47.5 0.40 -4.16
## IncidIPV:PovStat 0 Above- 1 Below 51.3 -2.94 -18.22
## IncidIPV:PovStat 1 Above- 0 Below 47.3 0.09 -5.86
## IncidIPV:PovStat 1 Above- 1 Below 49.2 -2.95 -19.46
## IncidIPV:PovStat 0 Below- 1 Below 52.3 -3.04 -19.67
## IncidIPV:BMIcat2 0 nonObese- 1 nonObese 59.6 -1.25 -12.14
## IncidIPV:BMIcat2 0 nonObese- 0 Obese 60.9 -0.83 -6.62
## IncidIPV:BMIcat2 0 nonObese- 1 Obese 45.5 -2.91 -14.16
## IncidIPV:BMIcat2 1 nonObese- 0 Obese 60.2 0.71 -4.99
## IncidIPV:BMIcat2 1 nonObese- 1 Obese 58.3 -0.89 -12.00
## IncidIPV:BMIcat2 0 Obese- 1 Obese 47.6 -2.11 -12.54
## Sex:BMIcat2 Women nonObese- Men nonObese 60.1 -0.07 -6.64
## Sex:BMIcat2 Women nonObese- Women Obese 65.9 -1.17 -9.60
## Sex:BMIcat2 Women nonObese- Men Obese 56.4 -0.64 -9.60
## Sex:BMIcat2 Men nonObese- Women Obese 52.4 -1.11 -9.31
## Sex:BMIcat2 Men nonObese- Men Obese 59.1 -0.60 -9.06
## Sex:BMIcat2 Women Obese- Men Obese 54.9 0.38 -5.28
## PovStat:BMIcat2 Above nonObese- Below nonObese 54.3 -2.47 -15.24
## PovStat:BMIcat2 Above nonObese- Above Obese 71.2 -2.32 -11.08
## PovStat:BMIcat2 Above nonObese- Below Obese 49.7 -2.68 -14.17
## PovStat:BMIcat2 Below nonObese- Above Obese 56.5 0.68 -4.83
## PovStat:BMIcat2 Below nonObese- Below Obese 56.0 0.08 -7.41
## PovStat:BMIcat2 Above Obese- Below Obese 55.0 -0.69 -8.33
## Upper CI p-value
## IncidIPV 0-1 -0.667 0.027
## Sex Women-Men 5.152 0.834
## PovStat Above-Below -0.559 0.029
## BMIcat2 nonObese-Obese 1.947 0.241
## IncidIPV:Sex 0 Women- 1 Women 0.355 0.064
## IncidIPV:Sex 0 Women- 0 Men 4.978 0.997
## IncidIPV:Sex 0 Women- 1 Men 2.053 0.159
## IncidIPV:Sex 1 Women- 0 Men 12.407 0.063
## IncidIPV:Sex 1 Women- 1 Men 8.860 0.806
## IncidIPV:Sex 0 Men- 1 Men 2.030 0.158
## IncidIPV:PovStat 0 Above- 1 Above 6.328 0.787
## IncidIPV:PovStat 0 Above- 0 Below 6.205 0.693
## IncidIPV:PovStat 0 Above- 1 Below -3.434 0.005
## IncidIPV:PovStat 1 Above- 0 Below 6.402 0.929
## IncidIPV:PovStat 1 Above- 1 Below -3.692 0.005
## IncidIPV:PovStat 0 Below- 1 Below -4.033 0.004
## IncidIPV:BMIcat2 0 nonObese- 1 nonObese 2.793 0.215
## IncidIPV:BMIcat2 0 nonObese- 0 Obese 2.731 0.409
## IncidIPV:BMIcat2 0 nonObese- 1 Obese -2.574 0.006
## IncidIPV:BMIcat2 1 nonObese- 0 Obese 10.447 0.482
## IncidIPV:BMIcat2 1 nonObese- 1 Obese 4.614 0.377
## IncidIPV:BMIcat2 0 Obese- 1 Obese -0.310 0.040
## Sex:BMIcat2 Women nonObese- Men nonObese 6.163 0.941
## Sex:BMIcat2 Women nonObese- Women Obese 2.505 0.246
## Sex:BMIcat2 Women nonObese- Men Obese 4.942 0.524
## Sex:BMIcat2 Men nonObese- Women Obese 2.693 0.274
## Sex:BMIcat2 Men nonObese- Men Obese 4.870 0.550
## Sex:BMIcat2 Women Obese- Men Obese 7.707 0.709
## PovStat:BMIcat2 Above nonObese- Below nonObese -1.586 0.017
## PovStat:BMIcat2 Above nonObese- Above Obese -0.834 0.023
## PovStat:BMIcat2 Above nonObese- Below Obese -2.020 0.010
## PovStat:BMIcat2 Below nonObese- Above Obese 9.739 0.502
## PovStat:BMIcat2 Below nonObese- Below Obese 8.045 0.934
## PovStat:BMIcat2 Above Obese- Below Obese 4.054 0.492
##
## Final model:
## lme4::lmer(formula = adjDBP ~ Age + IncidIPV + Sex + PovStat +
## BMIcat2 + (1 | HNDid) + Age:IncidIPV + Age:Sex + Age:PovStat +
## Age:BMIcat2 + IncidIPV:Sex + IncidIPV:PovStat + IncidIPV:BMIcat2 +
## Sex:BMIcat2 + PovStat:BMIcat2 + Age:IncidIPV:Sex + Age:IncidIPV:PovStat +
## Age:IncidIPV:BMIcat2 + Age:Sex:BMIcat2 + Age:PovStat:BMIcat2 +
## IncidIPV:Sex:BMIcat2 + IncidIPV:PovStat:BMIcat2 + Age:IncidIPV:Sex:BMIcat2 +
## Age:IncidIPV:PovStat:BMIcat2, data = IncidIPVbp, REML = reml,
## contrasts = l)
summary(mm4)
## Warning: number of observations <= rank(Z); variance-covariance matrix will be unidentifiable
## Warning: number of observations <= rank(Z); variance-covariance matrix will be unidentifiable
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: adjDBP ~ (Age + IncidIPV + Sex + PovStat + BMIcat2)^5 + (Age | HNDid) + (1 | subclass)
## Data: IncidIPVbp
##
## REML criterion at convergence: 859.3
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 31.7614 5.636
## Age 0.0145 0.121 -1.00
## subclass (Intercept) 1.3901 1.179
## Residual 40.1047 6.333
## Number of obs: 138, groups: HNDid, 70; subclass, 46
##
## Fixed effects:
## Estimate Std. Error
## (Intercept) 70.8643 3.3520
## Age 0.5661 0.4521
## IncidIPV1 12.1346 10.4527
## SexMen 1.0032 4.8841
## PovStatBelow 4.9500 5.5676
## BMIcat2Obese 5.5728 6.6320
## Age:IncidIPV1 1.3441 1.1564
## Age:SexMen -0.7740 0.5749
## Age:PovStatBelow 0.0218 0.5953
## Age:BMIcat2Obese -0.5869 0.8268
## IncidIPV1:SexMen -18.3981 12.8348
## IncidIPV1:PovStatBelow 11.2074 23.3907
## IncidIPV1:BMIcat2Obese -14.6964 12.3643
## SexMen:PovStatBelow 0.4771 7.6918
## SexMen:BMIcat2Obese -4.6951 8.2296
## PovStatBelow:BMIcat2Obese -2.2832 8.7571
## Age:IncidIPV1:SexMen -1.5549 1.4455
## Age:IncidIPV1:PovStatBelow 0.2967 2.1925
## Age:IncidIPV1:BMIcat2Obese -1.3722 1.3885
## Age:SexMen:PovStatBelow 0.7028 0.8005
## Age:SexMen:BMIcat2Obese 0.6419 0.9686
## Age:PovStatBelow:BMIcat2Obese 0.7831 1.0282
## IncidIPV1:SexMen:PovStatBelow 11.6144 25.3854
## IncidIPV1:SexMen:BMIcat2Obese 33.4345 21.0401
## IncidIPV1:PovStatBelow:BMIcat2Obese -24.9780 26.1762
## SexMen:PovStatBelow:BMIcat2Obese -12.9211 19.2301
## Age:IncidIPV1:SexMen:PovStatBelow 2.9659 2.6722
## Age:IncidIPV1:SexMen:BMIcat2Obese 3.0754 2.0170
## Age:IncidIPV1:PovStatBelow:BMIcat2Obese -4.2474 2.5788
## Age:SexMen:PovStatBelow:BMIcat2Obese -1.6385 1.8893
## IncidIPV1:SexMen:PovStatBelow:BMIcat2Obese 12.9543 36.1037
## Age:IncidIPV1:SexMen:PovStatBelow:BMIcat2Obese 0.5618 3.5704
## df t value Pr(>|t|)
## (Intercept) 12.3000 21.14 4.6e-11
## Age 58.4000 1.25 0.22
## IncidIPV1 56.2000 1.16 0.25
## SexMen 15.3000 0.21 0.84
## PovStatBelow 11.3000 0.89 0.39
## BMIcat2Obese 34.8000 0.84 0.41
## Age:IncidIPV1 76.5000 1.16 0.25
## Age:SexMen 57.1000 -1.35 0.18
## Age:PovStatBelow 54.1000 0.04 0.97
## Age:BMIcat2Obese 94.5000 -0.71 0.48
## IncidIPV1:SexMen 46.1000 -1.43 0.16
## IncidIPV1:PovStatBelow 74.0000 0.48 0.63
## IncidIPV1:BMIcat2Obese 57.7000 -1.19 0.24
## SexMen:PovStatBelow 15.8000 0.06 0.95
## SexMen:BMIcat2Obese 26.5000 -0.57 0.57
## PovStatBelow:BMIcat2Obese 20.4000 -0.26 0.80
## Age:IncidIPV1:SexMen 90.8000 -1.08 0.28
## Age:IncidIPV1:PovStatBelow 83.0000 0.14 0.89
## Age:IncidIPV1:BMIcat2Obese 90.8000 -0.99 0.33
## Age:SexMen:PovStatBelow 72.0000 0.88 0.38
## Age:SexMen:BMIcat2Obese 74.3000 0.66 0.51
## Age:PovStatBelow:BMIcat2Obese 78.1000 0.76 0.45
## IncidIPV1:SexMen:PovStatBelow 68.7000 0.46 0.65
## IncidIPV1:SexMen:BMIcat2Obese 77.2000 1.59 0.12
## IncidIPV1:PovStatBelow:BMIcat2Obese 71.2000 -0.95 0.34
## SexMen:PovStatBelow:BMIcat2Obese 80.8000 -0.67 0.50
## Age:IncidIPV1:SexMen:PovStatBelow 85.0000 1.11 0.27
## Age:IncidIPV1:SexMen:BMIcat2Obese 89.3000 1.52 0.13
## Age:IncidIPV1:PovStatBelow:BMIcat2Obese 91.1000 -1.65 0.10
## Age:SexMen:PovStatBelow:BMIcat2Obese 94.3000 -0.87 0.39
## IncidIPV1:SexMen:PovStatBelow:BMIcat2Obese 88.6000 0.36 0.72
## Age:IncidIPV1:SexMen:PovStatBelow:BMIcat2Obese 95.8000 0.16 0.88
##
## Correlation matrix not shown by default, as p = 32 > 20.
## Use print(x, correlation=TRUE) or
## vcov(x) if you need it
plot(st)
Diastolic Blood Pressure Final Model (BMI as categorical)
(mm4 = lmer(adjDBP ~Age + IncidIPV + Sex + PovStat + BMIcat2 + Age:IncidIPV + Age:Sex + Age:PovStat + Age:BMIcat2 + IncidIPV:Sex + IncidIPV:PovStat + IncidIPV:BMIcat2 + Sex:BMIcat2 + PovStat:BMIcat2 + Age:IncidIPV:Sex + Age:IncidIPV:PovStat + Age:IncidIPV:BMIcat2 + Age:Sex:BMIcat2 + Age:PovStat:BMIcat2 + IncidIPV:Sex:BMIcat2 + IncidIPV:PovStat:BMIcat2 + Age:IncidIPV:Sex:BMIcat2 + Age:IncidIPV:PovStat:BMIcat2 + (Age|HNDid) + (1|subclass) ,data = IncidIPVbp))
## Warning: number of observations <= rank(Z); variance-covariance matrix
## will be unidentifiable
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: adjDBP ~ Age + IncidIPV + Sex + PovStat + BMIcat2 + Age:IncidIPV + Age:Sex + Age:PovStat + Age:BMIcat2 + IncidIPV:Sex + IncidIPV:PovStat + IncidIPV:BMIcat2 + Sex:BMIcat2 + PovStat:BMIcat2 + Age:IncidIPV:Sex + Age:IncidIPV:PovStat + Age:IncidIPV:BMIcat2 + Age:Sex:BMIcat2 + Age:PovStat:BMIcat2 + IncidIPV:Sex:BMIcat2 + IncidIPV:PovStat:BMIcat2 + Age:IncidIPV:Sex:BMIcat2 + Age:IncidIPV:PovStat:BMIcat2 + (Age | HNDid) + (1 | subclass)
## Data: IncidIPVbp
## REML criterion at convergence: 903
## Random effects:
## Groups Name Std.Dev. Corr
## HNDid (Intercept) 4.706559
## Age 0.250493 -1.00
## subclass (Intercept) 0.000219
## Residual 6.400323
## Number of obs: 138, groups: HNDid, 70; subclass, 46
## Fixed Effects:
## (Intercept)
## 70.910
## Age
## 0.342
## IncidIPV1
## 8.121
## SexMen
## 1.373
## PovStatBelow
## 4.680
## BMIcat2Obese
## 7.139
## Age:IncidIPV1
## 0.945
## Age:SexMen
## -0.428
## Age:PovStatBelow
## 0.386
## Age:BMIcat2Obese
## -0.249
## IncidIPV1:SexMen
## -14.094
## IncidIPV1:PovStatBelow
## 25.134
## IncidIPV1:BMIcat2Obese
## -12.632
## SexMen:BMIcat2Obese
## -7.326
## PovStatBelow:BMIcat2Obese
## -4.483
## Age:IncidIPV1:SexMen
## -0.937
## Age:IncidIPV1:PovStatBelow
## 1.918
## Age:IncidIPV1:BMIcat2Obese
## -1.327
## Age:SexMen:BMIcat2Obese
## 0.145
## Age:PovStatBelow:BMIcat2Obese
## 0.216
## IncidIPV1:SexMen:BMIcat2Obese
## 39.427
## IncidIPV1:PovStatBelow:BMIcat2Obese
## -30.278
## Age:IncidIPV1:SexMen:BMIcat2Obese
## 3.794
## Age:IncidIPV1:PovStatBelow:BMIcat2Obese
## -4.438
summary(mm4)
## Warning: number of observations <= rank(Z); variance-covariance matrix will be unidentifiable
## Warning: number of observations <= rank(Z); variance-covariance matrix will be unidentifiable
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: adjDBP ~ Age + IncidIPV + Sex + PovStat + BMIcat2 + Age:IncidIPV + Age:Sex + Age:PovStat + Age:BMIcat2 + IncidIPV:Sex + IncidIPV:PovStat + IncidIPV:BMIcat2 + Sex:BMIcat2 + PovStat:BMIcat2 + Age:IncidIPV:Sex + Age:IncidIPV:PovStat + Age:IncidIPV:BMIcat2 + Age:Sex:BMIcat2 + Age:PovStat:BMIcat2 + IncidIPV:Sex:BMIcat2 + IncidIPV:PovStat:BMIcat2 + Age:IncidIPV:Sex:BMIcat2 + Age:IncidIPV:PovStat:BMIcat2 + (Age | HNDid) + (1 | subclass)
## Data: IncidIPVbp
##
## REML criterion at convergence: 903
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 2.22e+01 4.706559
## Age 6.27e-02 0.250493 -1.00
## subclass (Intercept) 4.78e-08 0.000219
## Residual 4.10e+01 6.400323
## Number of obs: 138, groups: HNDid, 70; subclass, 46
##
## Fixed effects:
## Estimate Std. Error df
## (Intercept) 70.910 2.684 18.800
## Age 0.342 0.354 90.500
## IncidIPV1 8.121 9.140 77.200
## SexMen 1.373 3.406 24.300
## PovStatBelow 4.680 3.450 22.700
## BMIcat2Obese 7.139 5.827 52.200
## Age:IncidIPV1 0.945 1.024 100.200
## Age:SexMen -0.428 0.389 98.600
## Age:PovStatBelow 0.386 0.386 97.800
## Age:BMIcat2Obese -0.249 0.712 111.800
## IncidIPV1:SexMen -14.094 10.461 71.300
## IncidIPV1:PovStatBelow 25.134 7.846 32.500
## IncidIPV1:BMIcat2Obese -12.632 10.927 74.200
## SexMen:BMIcat2Obese -7.326 6.744 47.500
## PovStatBelow:BMIcat2Obese -4.483 6.873 43.400
## Age:IncidIPV1:SexMen -0.937 1.195 113.100
## Age:IncidIPV1:PovStatBelow 1.918 0.982 107.300
## Age:IncidIPV1:BMIcat2Obese -1.327 1.233 108.700
## Age:SexMen:BMIcat2Obese 0.145 0.772 110.400
## Age:PovStatBelow:BMIcat2Obese 0.216 0.765 106.500
## IncidIPV1:SexMen:BMIcat2Obese 39.427 14.609 71.800
## IncidIPV1:PovStatBelow:BMIcat2Obese -30.278 12.555 53.000
## Age:IncidIPV1:SexMen:BMIcat2Obese 3.794 1.563 112.500
## Age:IncidIPV1:PovStatBelow:BMIcat2Obese -4.438 1.432 109.100
## t value Pr(>|t|)
## (Intercept) 26.42 2.2e-16
## Age 0.96 0.3375
## IncidIPV1 0.89 0.3770
## SexMen 0.40 0.6905
## PovStatBelow 1.36 0.1882
## BMIcat2Obese 1.23 0.2260
## Age:IncidIPV1 0.92 0.3581
## Age:SexMen -1.10 0.2743
## Age:PovStatBelow 1.00 0.3188
## Age:BMIcat2Obese -0.35 0.7276
## IncidIPV1:SexMen -1.35 0.1821
## IncidIPV1:PovStatBelow 3.20 0.0030
## IncidIPV1:BMIcat2Obese -1.16 0.2514
## SexMen:BMIcat2Obese -1.09 0.2828
## PovStatBelow:BMIcat2Obese -0.65 0.5176
## Age:IncidIPV1:SexMen -0.78 0.4345
## Age:IncidIPV1:PovStatBelow 1.95 0.0534
## Age:IncidIPV1:BMIcat2Obese -1.08 0.2844
## Age:SexMen:BMIcat2Obese 0.19 0.8516
## Age:PovStatBelow:BMIcat2Obese 0.28 0.7782
## IncidIPV1:SexMen:BMIcat2Obese 2.70 0.0087
## IncidIPV1:PovStatBelow:BMIcat2Obese -2.41 0.0194
## Age:IncidIPV1:SexMen:BMIcat2Obese 2.43 0.0168
## Age:IncidIPV1:PovStatBelow:BMIcat2Obese -3.10 0.0025
##
## Correlation matrix not shown by default, as p = 24 > 20.
## Use print(x, correlation=TRUE) or
## vcov(x) if you need it
DBP and IPV: Age x IncidIPV x PovStat x BMI
## Age IncidIPV PovStat BMIcat2 adjDBP Sex hat
## 1 30 0 Above nonObese 0 0.4928 75.51
## 2 31 0 Above nonObese 0 0.4928 75.64
## 3 32 0 Above nonObese 0 0.4928 75.77
## 4 33 0 Above nonObese 0 0.4928 75.90
## 5 34 0 Above nonObese 0 0.4928 76.04
## 6 35 0 Above nonObese 0 0.4928 76.17
DBP and IPV: Age x IncidIPV x PovStat x BMI (Plots broken up by poverty status)
## Age IncidIPV PovStat BMIcat2 adjDBP Sex hat
## 1 30 0 Above nonObese 0 0.4928 75.51
## 2 31 0 Above nonObese 0 0.4928 75.64
## 3 32 0 Above nonObese 0 0.4928 75.77
## 4 33 0 Above nonObese 0 0.4928 75.90
## 5 34 0 Above nonObese 0 0.4928 76.04
## 6 35 0 Above nonObese 0 0.4928 76.17
DBP and IPV: Age x IncidIPV x Sex x BMI
## Age IncidIPV Sex BMIcat2 adjDBP PovStat hat
## 1 30 0 Women nonObese 0 0.3913 87.53
## 2 31 0 Women nonObese 0 0.3913 88.02
## 3 32 0 Women nonObese 0 0.3913 88.51
## 4 33 0 Women nonObese 0 0.3913 89.01
## 5 34 0 Women nonObese 0 0.3913 89.50
## 6 35 0 Women nonObese 0 0.3913 89.99
DBP and IPV: Age x IncidIPV x Sex x BMI (broken up by gender)
## Age IncidIPV Sex BMIcat2 adjDBP PovStat hat
## 1 30 0 Women nonObese 0 0.3913 87.53
## 2 31 0 Women nonObese 0 0.3913 88.02
## 3 32 0 Women nonObese 0 0.3913 88.51
## 4 33 0 Women nonObese 0 0.3913 89.01
## 5 34 0 Women nonObese 0 0.3913 89.50
## 6 35 0 Women nonObese 0 0.3913 89.99