Clock Total Regression Model 1
## Loading required package: Matrix
## 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
## Warning: number of observations <= rank(Z); variance-covariance matrix
## will be unidentifiable
## Linear mixed model fit by maximum likelihood ['merModLmerTest']
## Formula: ClockTotal ~ (Age + IPVstatus + Sex + Race)^4 + (Age | HNDid) + (1 | subclass)
## Data: IPVandCognitionDataSet2
## AIC BIC logLik deviance
## 409.4 469.0 -183.7 367.4
## Random effects:
## Groups Name Std.Dev. Corr
## HNDid (Intercept) 0.1878
## Age 0.0171 -1.00
## subclass (Intercept) 0.3120
## Residual 0.9555
## Number of obs: 126, groups: HNDid, 63; subclass, 21
## Fixed Effects:
## (Intercept) Age
## 9.8576 0.0473
## IPVstatus1 SexMen
## -1.0731 -0.3281
## RaceAfrAm Age:IPVstatus1
## -1.7616 -0.0891
## Age:SexMen Age:RaceAfrAm
## -0.0289 -0.0412
## IPVstatus1:SexMen IPVstatus1:RaceAfrAm
## 1.3641 1.6606
## SexMen:RaceAfrAm Age:IPVstatus1:SexMen
## 1.6400 0.1263
## Age:IPVstatus1:RaceAfrAm Age:SexMen:RaceAfrAm
## 0.0322 0.0326
## IPVstatus1:SexMen:RaceAfrAm Age:IPVstatus1:SexMen:RaceAfrAm
## -2.6329 -0.0784
## Warning: number of observations <= rank(Z); variance-covariance matrix will be unidentifiable
## Warning:
## model has been refitted with REML=TRUE
##
## Warning: number of observations <= rank(Z); variance-covariance matrix will be unidentifiable
## Warning: number of observations <= rank(Z); variance-covariance matrix will be unidentifiable
## Warning: number of observations <= rank(Z); variance-covariance matrix will be unidentifiable
## 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 | HNDid) 0.05 1 1 0.8182
## (1 | HNDid) 0.08 1 2 0.7831
## (Age + 0 | HNDid) 0.69 1 3 0.4071
## (1 | subclass) 3.76 1 kept 0.0524
##
## Fixed effects:
## Sum Sq Mean Sq NumDF DenDF F.value elim.num Pr(>F)
## Age:IPVstatus:Sex:Race 0.2073 0.2073 1 105.70 0.1884 1 0.6651
## Age:IPVstatus:Race 0.0009 0.0009 1 108.97 0.0045 2 0.9468
## Age:Sex:Race 0.0182 0.0182 1 108.80 0.0223 3 0.8816
## Age:Race 1.4114 1.4114 1 99.69 0.8625 4 0.3553
## Age:IPVstatus:Sex 0.3715 0.3715 1 113.03 1.0625 5 0.3049
## Age:IPVstatus 1.3709 1.3709 1 100.08 0.5637 6 0.4545
## Age:Sex 0.9769 0.9769 1 105.92 0.3437 7 0.5590
## Age 0.1959 0.1959 1 32.06 0.0612 8 0.8062
## IPVstatus 0.0433 0.0433 1 104.11 0.0001 kept 0.9920
## Sex 7.1067 7.1067 1 117.92 1.8623 kept 0.1750
## Race 5.6901 5.6901 1 109.06 6.4736 kept 0.0123
## IPVstatus:Sex 3.3284 3.3284 1 105.57 1.2272 kept 0.2705
## IPVstatus:Race 0.4005 0.4005 1 116.65 0.4375 kept 0.5096
## Sex:Race 2.0465 2.0465 1 111.62 0.8826 kept 0.3495
## IPVstatus:Sex:Race 5.8356 5.8356 1 116.61 5.7222 kept 0.0183
##
## Least squares means:
## IPVstatus Sex Race Estimate
## IPVstatus 0 1.0 NA NA 9.089
## IPVstatus 1 2.0 NA NA 9.087
## Sex Women NA 2.0 NA 8.936
## Sex Men NA 1.0 NA 9.240
## Race White NA NA 2.0 9.382
## Race AfrAm NA NA 1.0 8.794
## IPVstatus:Sex 0 Women 1.0 2.0 NA 8.809
## IPVstatus:Sex 1 Women 2.0 2.0 NA 9.064
## IPVstatus:Sex 0 Men 1.0 1.0 NA 9.370
## IPVstatus:Sex 1 Men 2.0 1.0 NA 9.111
## IPVstatus:Race 0 White 1.0 NA 2.0 9.458
## IPVstatus:Race 1 White 2.0 NA 2.0 9.307
## IPVstatus:Race 0 AfrAm 1.0 NA 1.0 8.721
## IPVstatus:Race 1 AfrAm 2.0 NA 1.0 8.867
## Sex:Race Women White NA 2.0 2.0 9.338
## Sex:Race Men White NA 1.0 2.0 9.426
## Sex:Race Women AfrAm NA 2.0 1.0 8.534
## Sex:Race Men AfrAm NA 1.0 1.0 9.054
## IPVstatus:Sex:Race 0 Women White 1.0 2.0 2.0 9.555
## IPVstatus:Sex:Race 1 Women White 2.0 2.0 2.0 9.121
## IPVstatus:Sex:Race 0 Men White 1.0 1.0 2.0 9.360
## IPVstatus:Sex:Race 1 Men White 2.0 1.0 2.0 9.493
## IPVstatus:Sex:Race 0 Women AfrAm 1.0 2.0 1.0 8.062
## IPVstatus:Sex:Race 1 Women AfrAm 2.0 2.0 1.0 9.007
## IPVstatus:Sex:Race 0 Men AfrAm 1.0 1.0 1.0 9.379
## IPVstatus:Sex:Race 1 Men AfrAm 2.0 1.0 1.0 8.729
## Standard Error DF t-value Lower CI
## IPVstatus 0 0.152 33.8 59.6 8.78
## IPVstatus 1 0.191 68.1 47.5 8.71
## Sex Women 0.174 45.1 51.2 8.59
## Sex Men 0.179 49.3 51.5 8.88
## Race White 0.200 53.9 47.0 8.98
## Race AfrAm 0.157 32.2 56.0 8.47
## IPVstatus:Sex 0 Women 0.220 74.3 40.0 8.37
## IPVstatus:Sex 1 Women 0.241 85.0 37.7 8.59
## IPVstatus:Sex 0 Men 0.183 46.7 51.2 9.00
## IPVstatus:Sex 1 Men 0.288 94.8 31.6 8.54
## IPVstatus:Race 0 White 0.232 69.9 40.8 8.99
## IPVstatus:Race 1 White 0.297 98.7 31.4 8.72
## IPVstatus:Race 0 AfrAm 0.172 42.2 50.8 8.37
## IPVstatus:Race 1 AfrAm 0.229 83.4 38.8 8.41
## Sex:Race Women White 0.267 77.0 35.0 8.81
## Sex:Race Men White 0.269 94.1 35.0 8.89
## Sex:Race Women AfrAm 0.205 54.3 41.6 8.12
## Sex:Race Men AfrAm 0.213 56.9 42.5 8.63
## IPVstatus:Sex:Race 0 Women White 0.361 105.2 26.4 8.84
## IPVstatus:Sex:Race 1 Women White 0.362 104.0 25.2 8.40
## IPVstatus:Sex:Race 0 Men White 0.263 81.7 35.5 8.84
## IPVstatus:Sex:Race 1 Men White 0.465 109.1 20.4 8.57
## IPVstatus:Sex:Race 0 Women AfrAm 0.244 55.8 33.0 7.57
## IPVstatus:Sex:Race 1 Women AfrAm 0.307 99.8 29.3 8.40
## IPVstatus:Sex:Race 0 Men AfrAm 0.225 71.9 41.6 8.93
## IPVstatus:Sex:Race 1 Men AfrAm 0.331 101.6 26.4 8.07
## Upper CI p-value
## IPVstatus 0 9.40 <2e-16
## IPVstatus 1 9.47 <2e-16
## Sex Women 9.29 <2e-16
## Sex Men 9.60 <2e-16
## Race White 9.78 <2e-16
## Race AfrAm 9.11 <2e-16
## IPVstatus:Sex 0 Women 9.25 <2e-16
## IPVstatus:Sex 1 Women 9.54 <2e-16
## IPVstatus:Sex 0 Men 9.74 <2e-16
## IPVstatus:Sex 1 Men 9.68 <2e-16
## IPVstatus:Race 0 White 9.92 <2e-16
## IPVstatus:Race 1 White 9.90 <2e-16
## IPVstatus:Race 0 AfrAm 9.07 <2e-16
## IPVstatus:Race 1 AfrAm 9.32 <2e-16
## Sex:Race Women White 9.87 <2e-16
## Sex:Race Men White 9.96 <2e-16
## Sex:Race Women AfrAm 8.95 <2e-16
## Sex:Race Men AfrAm 9.48 <2e-16
## IPVstatus:Sex:Race 0 Women White 10.27 <2e-16
## IPVstatus:Sex:Race 1 Women White 9.84 <2e-16
## IPVstatus:Sex:Race 0 Men White 9.88 <2e-16
## IPVstatus:Sex:Race 1 Men White 10.41 <2e-16
## IPVstatus:Sex:Race 0 Women AfrAm 8.55 <2e-16
## IPVstatus:Sex:Race 1 Women AfrAm 9.62 <2e-16
## IPVstatus:Sex:Race 0 Men AfrAm 9.83 <2e-16
## IPVstatus:Sex:Race 1 Men AfrAm 9.39 <2e-16
##
## Differences of LSMEANS:
## Estimate Standard Error DF t-value
## IPVstatus 0-1 0.0 0.2098 104.1 0.01
## Sex Women-Men -0.3 0.2226 117.9 -1.36
## Race White-AfrAm 0.6 0.2311 109.1 2.54
## IPVstatus:Sex 0 Women- 1 Women -0.3 0.3022 117.1 -0.84
## IPVstatus:Sex 0 Women- 0 Men -0.6 0.2665 115.5 -2.10
## IPVstatus:Sex 0 Women- 1 Men -0.3 0.3371 116.0 -0.90
## IPVstatus:Sex 1 Women- 0 Men -0.3 0.2712 113.4 -1.13
## IPVstatus:Sex 1 Women- 1 Men 0.0 0.3685 112.8 -0.13
## IPVstatus:Sex 0 Men- 1 Men 0.3 0.3231 116.2 0.80
## IPVstatus:Race 0 White- 1 White 0.2 0.3524 117.6 0.43
## IPVstatus:Race 0 White- 0 AfrAm 0.7 0.2710 108.8 2.72
## IPVstatus:Race 0 White- 1 AfrAm 0.6 0.3047 117.9 1.94
## IPVstatus:Race 1 White- 0 AfrAm 0.6 0.3195 117.7 1.83
## IPVstatus:Race 1 White- 1 AfrAm 0.4 0.3668 115.4 1.20
## IPVstatus:Race 0 AfrAm- 1 AfrAm -0.1 0.2550 112.4 -0.58
## Sex:Race Women White- Men White -0.1 0.3575 117.9 -0.25
## Sex:Race Women White- Women AfrAm 0.8 0.3245 102.7 2.48
## Sex:Race Women White- Men AfrAm 0.3 0.3251 108.9 0.87
## Sex:Race Men White- Women AfrAm 0.9 0.3167 117.9 2.82
## Sex:Race Men White- Men AfrAm 0.4 0.3268 115.7 1.14
## Sex:Race Women AfrAm- Men AfrAm -0.5 0.2766 106.9 -1.88
## Lower CI Upper CI p-value
## IPVstatus 0-1 -0.4140 0.4182 0.992
## Sex Women-Men -0.7447 0.1370 0.175
## Race White-AfrAm 0.1300 1.0462 0.012
## IPVstatus:Sex 0 Women- 1 Women -0.8534 0.3435 0.401
## IPVstatus:Sex 0 Women- 0 Men -1.0887 -0.0330 0.037
## IPVstatus:Sex 0 Women- 1 Men -0.9693 0.3659 0.373
## IPVstatus:Sex 1 Women- 0 Men -0.8432 0.2314 0.262
## IPVstatus:Sex 1 Women- 1 Men -0.7768 0.6832 0.899
## IPVstatus:Sex 0 Men- 1 Men -0.3809 0.8992 0.424
## IPVstatus:Race 0 White- 1 White -0.5470 0.8487 0.669
## IPVstatus:Race 0 White- 0 AfrAm 0.1998 1.2739 0.008
## IPVstatus:Race 0 White- 1 AfrAm -0.0132 1.1936 0.055
## IPVstatus:Race 1 White- 0 AfrAm -0.0467 1.2187 0.069
## IPVstatus:Race 1 White- 1 AfrAm -0.2873 1.1660 0.234
## IPVstatus:Race 0 AfrAm- 1 AfrAm -0.6519 0.3586 0.566
## Sex:Race Women White- Men White -0.7963 0.6197 0.805
## Sex:Race Women White- Women AfrAm 0.1601 1.4471 0.015
## Sex:Race Women White- Men AfrAm -0.3600 0.9285 0.384
## Sex:Race Men White- Women AfrAm 0.2647 1.5191 0.006
## Sex:Race Men White- Men AfrAm -0.2748 1.0199 0.257
## Sex:Race Women AfrAm- Men AfrAm -1.0677 0.0291 0.063
##
## Final model:
## lme4::lmer(formula = ClockTotal ~ IPVstatus + Sex + Race + (1 |
## subclass) + IPVstatus:Sex + IPVstatus:Race + Sex:Race + IPVstatus:Sex:Race,
## data = IPVandCognitionDataSet2, REML = reml, contrasts = l)
Re-run suggested final Model 1
(mm1 = lmer(ClockTotal ~ IPVstatus + Sex + Race + (1 | subclass) + IPVstatus:Sex +
IPVstatus:Race + Sex:Race + IPVstatus:Sex:Race, data = IPVandCognitionDataSet2,
REML = F))
## Linear mixed model fit by maximum likelihood ['merModLmerTest']
## Formula: ClockTotal ~ IPVstatus + Sex + Race + (1 | subclass) + IPVstatus:Sex + IPVstatus:Race + Sex:Race + IPVstatus:Sex:Race
## Data: IPVandCognitionDataSet2
## AIC BIC logLik deviance
## 391.6 419.9 -185.8 371.6
## Random effects:
## Groups Name Std.Dev.
## subclass (Intercept) 0.367
## Residual 1.007
## Number of obs: 126, groups: subclass, 21
## Fixed Effects:
## (Intercept) IPVstatus1
## 9.559 -0.430
## SexMen RaceAfrAm
## -0.199 -1.491
## IPVstatus1:SexMen IPVstatus1:RaceAfrAm
## 0.549 1.375
## SexMen:RaceAfrAm IPVstatus1:SexMen:RaceAfrAm
## 1.505 -2.145
summary(mm1)
## Linear mixed model fit by maximum likelihood ['merModLmerTest']
## Formula: ClockTotal ~ IPVstatus + Sex + Race + (1 | subclass) + IPVstatus:Sex + IPVstatus:Race + Sex:Race + IPVstatus:Sex:Race
## Data: IPVandCognitionDataSet2
##
## AIC BIC logLik deviance
## 391.6 419.9 -185.8 371.6
##
## Random effects:
## Groups Name Variance Std.Dev.
## subclass (Intercept) 0.135 0.367
## Residual 1.013 1.007
## Number of obs: 126, groups: subclass, 21
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 9.559 0.348 113.200 27.45 < 2e-16
## IPVstatus1 -0.430 0.473 124.900 -0.91 0.36429
## SexMen -0.199 0.416 125.900 -0.48 0.63337
## RaceAfrAm -1.491 0.415 110.300 -3.59 0.00049
## IPVstatus1:SexMen 0.549 0.711 119.100 0.77 0.44163
## IPVstatus1:RaceAfrAm 1.375 0.602 124.800 2.28 0.02415
## SexMen:RaceAfrAm 1.505 0.521 119.100 2.89 0.00458
## IPVstatus1:SexMen:RaceAfrAm -2.145 0.873 124.300 -2.46 0.01541
##
## Correlation of Fixed Effects:
## (Intr) IPVst1 SexMen RcAfrA IPVs1:SM IPV1:R SM:RAA
## IPVstatus1 -0.677
## SexMen -0.795 0.586
## RaceAfrAm -0.828 0.583 0.671
## IPVstts1:SM 0.458 -0.694 -0.608 -0.368
## IPVstt1:RAA 0.546 -0.802 -0.445 -0.662 0.522
## SxMn:RcAfrA 0.661 -0.486 -0.806 -0.795 0.479 0.533
## IPV1:SM:RAA -0.387 0.574 0.482 0.456 -0.799 -0.693 -0.585
plot(st)
plot(mm1)
Clock Total Regression Model 2 (with CES)
load("~/Desktop/Megan/Research/IPV and Cognition Paper/IPV R Output/IPVandCognitionDataSet2.rda")
library(lme4)
library(lmerTest)
(mm2 = lmer(ClockTotal ~ (Age + IPVstatus + Sex + Race + CES1)^5 + (Age | HNDid) +
(1 | subclass), data = IPVandCognitionDataSet2, REML = F))
## Warning: number of observations <= rank(Z); variance-covariance matrix
## will be unidentifiable
## Linear mixed model fit by maximum likelihood ['merModLmerTest']
## Formula: ClockTotal ~ (Age + IPVstatus + Sex + Race + CES1)^5 + (Age | HNDid) + (1 | subclass)
## Data: IPVandCognitionDataSet2
## AIC BIC logLik deviance
## 424.8 529.7 -175.4 350.8
## Random effects:
## Groups Name Std.Dev. Corr
## HNDid (Intercept) 0.0000
## Age 0.0343 NaN
## subclass (Intercept) 0.0000
## Residual 0.9271
## Number of obs: 126, groups: HNDid, 63; subclass, 21
## Fixed Effects:
## (Intercept)
## 9.83453
## Age
## 0.08222
## IPVstatus1
## -5.83453
## SexMen
## 0.13461
## RaceAfrAm
## -1.48418
## CES11
## 0.16547
## Age:IPVstatus1
## -0.58222
## Age:SexMen
## -0.05494
## Age:RaceAfrAm
## -0.02764
## Age:CES11
## -0.08222
## IPVstatus1:SexMen
## 5.49586
## IPVstatus1:RaceAfrAm
## 5.65603
## IPVstatus1:CES11
## 5.33420
## SexMen:RaceAfrAm
## 0.82196
## SexMen:CES11
## -1.76667
## RaceAfrAm:CES11
## -0.49528
## Age:IPVstatus1:SexMen
## 0.62312
## Age:IPVstatus1:RaceAfrAm
## 0.40785
## Age:IPVstatus1:CES11
## 0.59925
## Age:SexMen:RaceAfrAm
## 0.02402
## Age:SexMen:CES11
## 0.04313
## Age:RaceAfrAm:CES11
## 0.01727
## IPVstatus1:SexMen:RaceAfrAm
## -5.62428
## IPVstatus1:SexMen:CES11
## -1.03014
## IPVstatus1:RaceAfrAm:CES11
## -4.34699
## SexMen:RaceAfrAm:CES11
## 2.32420
## Age:IPVstatus1:SexMen:RaceAfrAm
## -0.47241
## Age:IPVstatus1:SexMen:CES11
## -0.29501
## Age:IPVstatus1:RaceAfrAm:CES11
## -0.40121
## Age:SexMen:RaceAfrAm:CES11
## 0.00676
## IPVstatus1:SexMen:RaceAfrAm:CES11
## -0.46247
## Age:IPVstatus1:SexMen:RaceAfrAm:CES11
## 0.14321
(st = step(mm2))
## Warning: number of observations <= rank(Z); variance-covariance matrix will be unidentifiable
## Warning:
## model has been refitted with REML=TRUE
##
## Warning: number of observations <= rank(Z); variance-covariance matrix will be unidentifiable
## Warning: number of observations <= rank(Z); variance-covariance matrix will be unidentifiable
## 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 | HNDid) 0.0 1 1 1.000
## (1 | subclass) 0.0 1 2 1.000
## (1 | HNDid) 0.0 1 3 1.000
## (Age + 0 | HNDid) 4.6 1 kept 0.032
##
## Fixed effects:
## Sum Sq Mean Sq NumDF DenDF F.value elim.num
## Age:IPVstatus:Sex:Race:CES1 0.0385 0.0385 1 85.94 0.0364 1
## Age:Sex:Race:CES1 0.0093 0.0093 1 94.65 0.0024 2
## IPVstatus:Sex:Race:CES1 0.2395 0.2395 1 80.30 0.2359 3
## Age:IPVstatus:Sex:CES1 0.1316 0.1316 1 89.13 0.2429 4
## Age:Sex:CES1 0.0920 0.0920 1 75.13 0.1548 5
## IPVstatus:Sex:CES1 0.4234 0.4234 1 43.61 0.5353 6
## Age:IPVstatus:Race:CES1 1.3819 1.3819 1 99.99 1.6120 7
## Age:Race:CES1 0.5974 0.5974 1 93.62 0.3477 8
## Age:IPVstatus:Sex:Race 3.0818 3.0818 1 71.09 1.8027 9
## Age:Sex:Race 0.0680 0.0680 1 73.66 0.3903 10
## Age:IPVstatus:Race 0.1522 0.1522 1 71.83 0.3939 11
## Age:Race 0.4992 0.4992 1 81.81 0.4027 12
## Age:IPVstatus:Sex 0.6449 0.6449 1 79.67 1.7548 13
## Age:Sex 0.4163 0.4163 1 76.63 0.1170 14
## Age:IPVstatus:CES1 1.3252 1.3252 1 69.27 1.3834 15
## Age:IPVstatus 0.7330 0.7330 1 68.29 0.1435 16
## Age:CES1 0.0037 0.0037 1 71.87 0.4228 17
## Age 0.7665 0.7665 1 80.24 1.3275 18
## IPVstatus:Race:CES1 3.2657 3.2657 1 62.65 1.8196 19
## IPVstatus:CES1 0.0670 0.0670 1 63.00 0.0375 20
## IPVstatus 0.0437 0.0437 1 76.56 0.4414 kept
## Sex 4.6437 4.6437 1 92.21 0.2917 kept
## Race 8.2411 8.2411 1 92.21 5.9343 kept
## CES1 1.1001 1.1001 1 74.23 0.0492 kept
## IPVstatus:Sex 3.2449 3.2449 1 76.56 0.2264 kept
## IPVstatus:Race 0.3025 0.3025 1 76.56 0.7677 kept
## Sex:Race 2.1462 2.1462 1 92.21 1.6725 kept
## Sex:CES1 0.5738 0.5738 1 74.23 1.8208 kept
## Race:CES1 0.5823 0.5823 1 74.23 0.0428 kept
## IPVstatus:Sex:Race 4.8728 4.8728 1 76.56 7.0240 kept
## Sex:Race:CES1 2.1571 2.1571 1 74.23 4.2264 kept
## Pr(>F)
## Age:IPVstatus:Sex:Race:CES1 0.8491
## Age:Sex:Race:CES1 0.9609
## IPVstatus:Sex:Race:CES1 0.6285
## Age:IPVstatus:Sex:CES1 0.6233
## Age:Sex:CES1 0.6951
## IPVstatus:Sex:CES1 0.4683
## Age:IPVstatus:Race:CES1 0.2072
## Age:Race:CES1 0.5568
## Age:IPVstatus:Sex:Race 0.1837
## Age:Sex:Race 0.5341
## Age:IPVstatus:Race 0.5322
## Age:Race 0.5275
## Age:IPVstatus:Sex 0.1891
## Age:Sex 0.7333
## Age:IPVstatus:CES1 0.2436
## Age:IPVstatus 0.7060
## Age:CES1 0.5176
## Age 0.2527
## IPVstatus:Race:CES1 0.1822
## IPVstatus:CES1 0.8471
## IPVstatus 0.5085
## Sex 0.5905
## Race 0.0168
## CES1 0.8251
## IPVstatus:Sex 0.6356
## IPVstatus:Race 0.3837
## Sex:Race 0.1992
## Sex:CES1 0.1813
## Race:CES1 0.8367
## IPVstatus:Sex:Race 0.0098
## Sex:Race:CES1 0.0433
##
## Least squares means:
## IPVstatus Sex Race CES1 Estimate
## IPVstatus 0 1.0 NA NA NA 9.13
## IPVstatus 1 2.0 NA NA NA 8.96
## Sex Women NA 2.0 NA NA 8.98
## Sex Men NA 1.0 NA NA 9.11
## Race White NA NA 2.0 NA 9.34
## Race AfrAm NA NA 1.0 NA 8.76
## CES1 0 NA NA NA 1.0 9.08
## CES1 1 NA NA NA 2.0 9.02
## IPVstatus:Sex 0 Women 1.0 2.0 NA NA 9.01
## IPVstatus:Sex 1 Women 2.0 2.0 NA NA 8.96
## IPVstatus:Sex 0 Men 1.0 1.0 NA NA 9.26
## IPVstatus:Sex 1 Men 2.0 1.0 NA NA 8.96
## IPVstatus:Race 0 White 1.0 NA 2.0 NA 9.54
## IPVstatus:Race 1 White 2.0 NA 2.0 NA 9.13
## IPVstatus:Race 0 AfrAm 1.0 NA 1.0 NA 8.73
## IPVstatus:Race 1 AfrAm 2.0 NA 1.0 NA 8.78
## Sex:Race Women White NA 2.0 2.0 NA 9.43
## Sex:Race Men White NA 1.0 2.0 NA 9.25
## Sex:Race Women AfrAm NA 2.0 1.0 NA 8.54
## Sex:Race Men AfrAm NA 1.0 1.0 NA 8.98
## Sex:CES1 Women 0 NA 2.0 NA 1.0 8.84
## Sex:CES1 Men 0 NA 1.0 NA 1.0 9.31
## Sex:CES1 Women 1 NA 2.0 NA 2.0 9.12
## Sex:CES1 Men 1 NA 1.0 NA 2.0 8.91
## Race:CES1 White 0 NA NA 2.0 1.0 9.39
## Race:CES1 AfrAm 0 NA NA 1.0 1.0 8.76
## Race:CES1 White 1 NA NA 2.0 2.0 9.28
## Race:CES1 AfrAm 1 NA NA 1.0 2.0 8.76
## IPVstatus:Sex:Race 0 Women White 1.0 2.0 2.0 NA 9.92
## IPVstatus:Sex:Race 1 Women White 2.0 2.0 2.0 NA 8.94
## IPVstatus:Sex:Race 0 Men White 1.0 1.0 2.0 NA 9.16
## IPVstatus:Sex:Race 1 Men White 2.0 1.0 2.0 NA 9.33
## IPVstatus:Sex:Race 0 Women AfrAm 1.0 2.0 1.0 NA 8.10
## IPVstatus:Sex:Race 1 Women AfrAm 2.0 2.0 1.0 NA 8.98
## IPVstatus:Sex:Race 0 Men AfrAm 1.0 1.0 1.0 NA 9.36
## IPVstatus:Sex:Race 1 Men AfrAm 2.0 1.0 1.0 NA 8.59
## Sex:Race:CES1 Women White 0 NA 2.0 2.0 1.0 9.05
## Sex:Race:CES1 Men White 0 NA 1.0 2.0 1.0 9.73
## Sex:Race:CES1 Women AfrAm 0 NA 2.0 1.0 1.0 8.63
## Sex:Race:CES1 Men AfrAm 0 NA 1.0 1.0 1.0 8.89
## Sex:Race:CES1 Women White 1 NA 2.0 2.0 2.0 9.80
## Sex:Race:CES1 Men White 1 NA 1.0 2.0 2.0 8.77
## Sex:Race:CES1 Women AfrAm 1 NA 2.0 1.0 2.0 8.45
## Sex:Race:CES1 Men AfrAm 1 NA 1.0 1.0 2.0 9.06
## Standard Error DF t-value Lower CI
## IPVstatus 0 0.151 80.6 60.6 8.83
## IPVstatus 1 0.201 85.3 44.6 8.56
## Sex Women 0.161 94.2 55.7 8.66
## Sex Men 0.176 90.6 51.9 8.76
## Race White 0.194 86.1 48.2 8.95
## Race AfrAm 0.139 103.2 63.1 8.48
## CES1 0 0.166 88.5 54.8 8.75
## CES1 1 0.181 77.8 49.9 8.66
## IPVstatus:Sex 0 Women 0.244 73.1 36.9 8.52
## IPVstatus:Sex 1 Women 0.261 81.4 34.3 8.44
## IPVstatus:Sex 0 Men 0.177 95.2 52.4 8.91
## IPVstatus:Sex 1 Men 0.305 88.1 29.4 8.35
## IPVstatus:Race 0 White 0.257 72.3 37.2 9.03
## IPVstatus:Race 1 White 0.324 77.7 28.2 8.49
## IPVstatus:Race 0 AfrAm 0.159 102.2 55.0 8.41
## IPVstatus:Race 1 AfrAm 0.238 99.1 37.0 8.31
## Sex:Race Women White 0.258 94.7 36.6 8.92
## Sex:Race Men White 0.290 79.2 31.9 8.67
## Sex:Race Women AfrAm 0.194 93.1 44.0 8.15
## Sex:Race Men AfrAm 0.198 110.5 45.2 8.58
## Sex:CES1 Women 0 0.243 86.6 36.4 8.36
## Sex:CES1 Men 0 0.225 90.8 41.4 8.86
## Sex:CES1 Women 1 0.261 67.7 34.9 8.60
## Sex:CES1 Men 1 0.250 89.4 35.6 8.42
## Race:CES1 White 0 0.262 84.0 35.9 8.87
## Race:CES1 AfrAm 0 0.203 95.9 43.2 8.36
## Race:CES1 White 1 0.310 67.6 30.0 8.67
## Race:CES1 AfrAm 1 0.187 105.6 46.8 8.38
## IPVstatus:Sex:Race 0 Women White 0.428 69.9 23.2 9.06
## IPVstatus:Sex:Race 1 Women White 0.423 73.7 21.1 8.09
## IPVstatus:Sex:Race 0 Men White 0.283 77.9 32.4 8.60
## IPVstatus:Sex:Race 1 Men White 0.490 80.7 19.0 8.36
## IPVstatus:Sex:Race 0 Women AfrAm 0.236 84.1 34.3 7.63
## IPVstatus:Sex:Race 1 Women AfrAm 0.306 96.4 29.3 8.37
## IPVstatus:Sex:Race 0 Men AfrAm 0.212 113.9 44.1 8.94
## IPVstatus:Sex:Race 1 Men AfrAm 0.363 100.9 23.6 7.87
## Sex:Race:CES1 Women White 0 0.414 85.3 21.9 8.23
## Sex:Race:CES1 Men White 0 0.320 81.7 30.4 9.09
## Sex:Race:CES1 Women AfrAm 0 0.254 89.8 34.0 8.13
## Sex:Race:CES1 Men AfrAm 0 0.316 99.7 28.1 8.26
## Sex:Race:CES1 Women White 1 0.438 59.1 22.4 8.93
## Sex:Race:CES1 Men White 1 0.438 77.1 20.0 7.89
## Sex:Race:CES1 Women AfrAm 1 0.285 91.0 29.6 7.88
## Sex:Race:CES1 Men AfrAm 1 0.242 113.9 37.5 8.58
## Upper CI p-value
## IPVstatus 0 9.43 <2e-16
## IPVstatus 1 9.36 <2e-16
## Sex Women 9.30 <2e-16
## Sex Men 9.46 <2e-16
## Race White 9.72 <2e-16
## Race AfrAm 9.03 <2e-16
## CES1 0 9.40 <2e-16
## CES1 1 9.38 <2e-16
## IPVstatus:Sex 0 Women 9.49 <2e-16
## IPVstatus:Sex 1 Women 9.48 <2e-16
## IPVstatus:Sex 0 Men 9.61 <2e-16
## IPVstatus:Sex 1 Men 9.57 <2e-16
## IPVstatus:Race 0 White 10.05 <2e-16
## IPVstatus:Race 1 White 9.78 <2e-16
## IPVstatus:Race 0 AfrAm 9.04 <2e-16
## IPVstatus:Race 1 AfrAm 9.26 <2e-16
## Sex:Race Women White 9.94 <2e-16
## Sex:Race Men White 9.82 <2e-16
## Sex:Race Women AfrAm 8.92 <2e-16
## Sex:Race Men AfrAm 9.37 <2e-16
## Sex:CES1 Women 0 9.32 <2e-16
## Sex:CES1 Men 0 9.76 <2e-16
## Sex:CES1 Women 1 9.65 <2e-16
## Sex:CES1 Men 1 9.41 <2e-16
## Race:CES1 White 0 9.91 <2e-16
## Race:CES1 AfrAm 0 9.16 <2e-16
## Race:CES1 White 1 9.90 <2e-16
## Race:CES1 AfrAm 1 9.13 <2e-16
## IPVstatus:Sex:Race 0 Women White 10.77 <2e-16
## IPVstatus:Sex:Race 1 Women White 9.78 <2e-16
## IPVstatus:Sex:Race 0 Men White 9.73 <2e-16
## IPVstatus:Sex:Race 1 Men White 10.31 <2e-16
## IPVstatus:Sex:Race 0 Women AfrAm 8.57 <2e-16
## IPVstatus:Sex:Race 1 Women AfrAm 9.59 <2e-16
## IPVstatus:Sex:Race 0 Men AfrAm 9.78 <2e-16
## IPVstatus:Sex:Race 1 Men AfrAm 9.31 <2e-16
## Sex:Race:CES1 Women White 0 9.88 <2e-16
## Sex:Race:CES1 Men White 0 10.37 <2e-16
## Sex:Race:CES1 Women AfrAm 0 9.13 <2e-16
## Sex:Race:CES1 Men AfrAm 0 9.52 <2e-16
## Sex:Race:CES1 Women White 1 10.68 <2e-16
## Sex:Race:CES1 Men White 1 9.64 <2e-16
## Sex:Race:CES1 Women AfrAm 1 9.01 <2e-16
## Sex:Race:CES1 Men AfrAm 1 9.54 <2e-16
##
## Differences of LSMEANS:
## Estimate Standard Error DF t-value
## IPVstatus 0-1 0.2 0.2635 76.6 0.66
## Sex Women-Men -0.1 0.2384 92.2 -0.54
## Race White-AfrAm 0.6 0.2384 92.2 2.44
## CES1 0-1 0.1 0.2518 74.2 0.22
## IPVstatus:Sex 0 Women- 1 Women 0.0 0.3900 66.7 0.13
## IPVstatus:Sex 0 Women- 0 Men -0.3 0.3017 80.6 -0.84
## IPVstatus:Sex 0 Women- 1 Men 0.0 0.3910 82.1 0.12
## IPVstatus:Sex 1 Women- 0 Men -0.3 0.3155 85.8 -0.96
## IPVstatus:Sex 1 Women- 1 Men 0.0 0.4018 85.3 -0.01
## IPVstatus:Sex 0 Men- 1 Men 0.3 0.3543 89.2 0.85
## IPVstatus:Race 0 White- 1 White 0.4 0.4373 67.7 0.93
## IPVstatus:Race 0 White- 0 AfrAm 0.8 0.3017 80.6 2.69
## IPVstatus:Race 0 White- 1 AfrAm 0.8 0.3497 84.7 2.16
## IPVstatus:Race 1 White- 0 AfrAm 0.4 0.3608 82.5 1.12
## IPVstatus:Race 1 White- 1 AfrAm 0.3 0.4018 85.3 0.87
## IPVstatus:Race 0 AfrAm- 1 AfrAm -0.1 0.2940 97.0 -0.19
## Sex:Race Women White- Men White 0.2 0.3878 86.1 0.46
## Sex:Race Women White- Women AfrAm 0.9 0.3223 94.2 2.76
## Sex:Race Women White- Men AfrAm 0.5 0.3250 101.6 1.39
## Sex:Race Men White- Women AfrAm 0.7 0.3488 83.5 2.03
## Sex:Race Men White- Men AfrAm 0.3 0.3513 90.6 0.78
## Sex:Race Women AfrAm- Men AfrAm -0.4 0.2774 103.2 -1.58
## Sex:CES1 Women 0- Men 0 -0.5 0.3310 88.5 -1.42
## Sex:CES1 Women 0- Women 1 -0.3 0.3880 64.7 -0.73
## Sex:CES1 Women 0- Men 1 -0.1 0.3485 88.0 -0.21
## Sex:CES1 Men 0- Women 1 0.2 0.3449 77.3 0.54
## Sex:CES1 Men 0- Men 1 0.4 0.3209 89.4 1.23
## Sex:CES1 Women 1- Men 1 0.2 0.3617 77.8 0.58
## Race:CES1 White 0- AfrAm 0 0.6 0.3310 88.5 1.91
## Race:CES1 White 0- White 1 0.1 0.4223 65.0 0.26
## Race:CES1 White 0- AfrAm 1 0.6 0.3215 91.8 1.98
## Race:CES1 AfrAm 0- White 1 -0.5 0.3702 75.8 -1.42
## Race:CES1 AfrAm 0- AfrAm 1 0.0 0.2742 97.8 0.01
## Race:CES1 White 1- AfrAm 1 0.5 0.3617 77.8 1.46
## Lower CI Upper CI p-value
## IPVstatus 0-1 -0.3496 0.700 0.508
## Sex Women-Men -0.6022 0.345 0.591
## Race White-AfrAm 0.1073 1.054 0.017
## CES1 0-1 -0.4458 0.557 0.825
## IPVstatus:Sex 0 Women- 1 Women -0.7289 0.828 0.899
## IPVstatus:Sex 0 Women- 0 Men -0.8543 0.346 0.402
## IPVstatus:Sex 0 Women- 1 Men -0.7316 0.824 0.906
## IPVstatus:Sex 1 Women- 0 Men -0.9311 0.324 0.338
## IPVstatus:Sex 1 Women- 1 Men -0.8023 0.795 0.993
## IPVstatus:Sex 0 Men- 1 Men -0.4036 1.004 0.399
## IPVstatus:Race 0 White- 1 White -0.4667 1.278 0.357
## IPVstatus:Race 0 White- 0 AfrAm 0.2113 1.412 0.009
## IPVstatus:Race 0 White- 1 AfrAm 0.0603 1.451 0.034
## IPVstatus:Race 1 White- 0 AfrAm -0.3120 1.123 0.264
## IPVstatus:Race 1 White- 1 AfrAm -0.4491 1.149 0.386
## IPVstatus:Race 0 AfrAm- 1 AfrAm -0.6393 0.528 0.850
## Sex:Race Women White- Men White -0.5913 0.950 0.644
## Sex:Race Women White- Women AfrAm 0.2490 1.529 0.007
## Sex:Race Women White- Men AfrAm -0.1928 1.097 0.167
## Sex:Race Men White- Women AfrAm 0.0158 1.403 0.045
## Sex:Race Men White- Men AfrAm -0.4254 0.970 0.440
## Sex:Race Women AfrAm- Men AfrAm -0.9871 0.113 0.118
## Sex:CES1 Women 0- Men 0 -1.1262 0.189 0.161
## Sex:CES1 Women 0- Women 1 -1.0588 0.491 0.467
## Sex:CES1 Women 0- Men 1 -0.7655 0.620 0.835
## Sex:CES1 Men 0- Women 1 -0.5021 0.871 0.594
## Sex:CES1 Men 0- Men 1 -0.2421 1.033 0.221
## Sex:CES1 Women 1- Men 1 -0.5092 0.931 0.561
## Race:CES1 White 0- AfrAm 0 -0.0250 1.290 0.059
## Race:CES1 White 0- White 1 -0.7355 0.951 0.799
## Race:CES1 White 0- AfrAm 1 -0.0021 1.275 0.051
## Race:CES1 AfrAm 0- White 1 -1.2622 0.212 0.160
## Race:CES1 AfrAm 0- AfrAm 1 -0.5404 0.548 0.989
## Race:CES1 White 1- AfrAm 1 -0.1915 1.249 0.148
##
## Final model:
## lme4::lmer(formula = ClockTotal ~ IPVstatus + Sex + Race + CES1 +
## (Age + 0 | HNDid) + IPVstatus:Sex + IPVstatus:Race + Sex:Race +
## Sex:CES1 + Race:CES1 + IPVstatus:Sex:Race + Sex:Race:CES1,
## data = IPVandCognitionDataSet2, REML = reml, contrasts = l)
Re-run suggested final Model 2
(mm2 = lmer(ClockTotal ~ IPVstatus + Sex + Race + CES1 + (Age + 0 | HNDid) +
IPVstatus:Sex + IPVstatus:Race + Sex:Race + Sex:CES1 + Race:CES1 + IPVstatus:Sex:Race +
Sex:Race:CES1, data = IPVandCognitionDataSet2, REML = F))
## Linear mixed model fit by maximum likelihood ['merModLmerTest']
## Formula: ClockTotal ~ IPVstatus + Sex + Race + CES1 + (Age + 0 | HNDid) + IPVstatus:Sex + IPVstatus:Race + Sex:Race + Sex:CES1 + Race:CES1 + IPVstatus:Sex:Race + Sex:Race:CES1
## Data: IPVandCognitionDataSet2
## AIC BIC logLik deviance
## 394.7 434.4 -183.3 366.7
## Random effects:
## Groups Name Std.Dev.
## HNDid Age 0.0407
## Residual 0.9768
## Number of obs: 126, groups: HNDid, 63
## Fixed Effects:
## (Intercept) IPVstatus1
## 9.527 -0.942
## SexMen RaceAfrAm
## 0.114 -1.344
## CES11 IPVstatus1:SexMen
## 0.731 1.086
## IPVstatus1:RaceAfrAm SexMen:RaceAfrAm
## 1.839 0.975
## SexMen:CES11 RaceAfrAm:CES11
## -1.678 -0.899
## IPVstatus1:SexMen:RaceAfrAm SexMen:RaceAfrAm:CES11
## -2.744 2.016
summary(mm2)
## Linear mixed model fit by maximum likelihood ['merModLmerTest']
## Formula: ClockTotal ~ IPVstatus + Sex + Race + CES1 + (Age + 0 | HNDid) + IPVstatus:Sex + IPVstatus:Race + Sex:Race + Sex:CES1 + Race:CES1 + IPVstatus:Sex:Race + Sex:Race:CES1
## Data: IPVandCognitionDataSet2
##
## AIC BIC logLik deviance
## 394.7 434.4 -183.3 366.7
##
## Random effects:
## Groups Name Variance Std.Dev.
## HNDid Age 0.00166 0.0407
## Residual 0.95415 0.9768
## Number of obs: 126, groups: HNDid, 63
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 9.527 0.353 124.700 27.02 <2e-16
## IPVstatus1 -0.942 0.631 69.600 -1.49 0.1399
## SexMen 0.114 0.454 113.700 0.25 0.8028
## RaceAfrAm -1.344 0.441 121.900 -3.05 0.0028
## CES11 0.731 0.632 68.200 1.16 0.2515
## IPVstatus1:SexMen 1.086 0.816 77.500 1.33 0.1874
## IPVstatus1:RaceAfrAm 1.839 0.728 77.000 2.53 0.0136
## SexMen:RaceAfrAm 0.975 0.575 117.600 1.69 0.0928
## SexMen:CES11 -1.678 0.788 75.100 -2.13 0.0364
## RaceAfrAm:CES11 -0.899 0.724 74.300 -1.24 0.2182
## IPVstatus1:SexMen:RaceAfrAm -2.744 0.986 86.700 -2.78 0.0066
## SexMen:RaceAfrAm:CES11 2.016 0.942 84.800 2.14 0.0352
##
## Correlation of Fixed Effects:
## (Intr) IPVst1 SexMen RcAfrA CES11 IPVs1:SM IPV1:R SxM:RAA
## IPVstatus1 -0.327
## SexMen -0.777 0.254
## RaceAfrAm -0.800 0.262 0.622
## CES11 -0.288 -0.641 0.224 0.231
## IPVstts1:SM 0.253 -0.773 -0.354 -0.202 0.496
## IPVstt1:RAA 0.284 -0.867 -0.221 -0.370 0.556 0.670
## SxMn:RcAfrA 0.613 -0.201 -0.788 -0.766 -0.177 0.279 0.284
## SexMn:CES11 0.231 0.515 -0.366 -0.185 -0.802 -0.413 -0.446 0.289
## RcAfA:CES11 0.252 0.560 -0.196 -0.365 -0.873 -0.433 -0.495 0.280
## IPV1:SM:RAA -0.209 0.640 0.293 0.273 -0.410 -0.827 -0.738 -0.330
## SM:RAA:CES1 -0.193 -0.430 0.306 0.280 0.671 0.346 0.380 -0.430
## SM:CES RAA:CE IPV1:SM:
## IPVstatus1
## SexMen
## RaceAfrAm
## CES11
## IPVstts1:SM
## IPVstt1:RAA
## SxMn:RcAfrA
## SexMn:CES11
## RcAfA:CES11 0.701
## IPV1:SM:RAA 0.342 0.365
## SM:RAA:CES1 -0.836 -0.768 -0.376
plot(st)
plot(mm2)
