Trails B-A 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: IPVandCognitionDataSet2$"TrailsB-A" ~ (Age + IPVstatus + Sex + Race)^4 + (Age | HNDid) + (1 | subclass)
## Data: IPVandCognitionDataSet2
## AIC BIC logLik deviance
## 1587.0 1646.5 -772.5 1545.0
## Random effects:
## Groups Name Std.Dev. Corr
## HNDid (Intercept) 1.36e+02
## Age 5.65e+00 1.00
## subclass (Intercept) 2.89e-03
## Residual 7.83e+01
## Number of obs: 126, groups: HNDid, 63; subclass, 21
## Fixed Effects:
## (Intercept) Age
## 31.25 -1.21
## IPVstatus1 SexMen
## 11.83 9.49
## RaceAfrAm Age:IPVstatus1
## 130.31 1.69
## Age:SexMen Age:RaceAfrAm
## 2.31 2.95
## IPVstatus1:SexMen IPVstatus1:RaceAfrAm
## 251.18 -135.67
## SexMen:RaceAfrAm Age:IPVstatus1:SexMen
## -30.80 14.53
## Age:IPVstatus1:RaceAfrAm Age:SexMen:RaceAfrAm
## -3.14 3.67
## IPVstatus1:SexMen:RaceAfrAm Age:IPVstatus1:SexMen:RaceAfrAm
## -40.40 -28.58
## 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
## 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 effects:
## Chi.sq Chi.DF elim.num p.value
## (1 | subclass) 0.05 1 1 0.8249
## (1 | HNDid) 22.85 1 kept 0
##
## Fixed effects:
## Sum Sq Mean Sq NumDF DenDF F.value elim.num
## Age:IPVstatus:Sex:Race 18083.3 18083.3 1 105.45 2.3725 1
## Age:IPVstatus:Sex 137.0 137.0 1 107.43 0.0017 2
## Age:Sex:Race 4236.6 4236.6 1 111.20 0.4524 3
## IPVstatus:Sex:Race 3899.5 3899.5 1 50.88 0.5066 4
## Sex:Race 179.9 179.9 1 53.60 0.0000 5
## Age:Sex 9915.8 9915.8 1 112.99 1.3441 6
## Age:IPVstatus:Race 29626.6 29626.6 1 113.01 3.3392 7
## Age:IPVstatus 1746.3 1746.3 1 113.57 0.0385 8
## Age:Race 7060.9 7060.9 1 117.61 0.2334 9
## IPVstatus:Race 4731.1 4731.1 1 55.35 0.7611 10
## Age 20013.6 20013.6 1 118.12 2.1321 11
## Race 35928.9 35928.9 1 58.00 3.5779 12
## IPVstatus 9966.6 9966.6 1 59.00 1.9029 kept
## Sex 17661.2 17661.2 1 59.00 6.4676 kept
## IPVstatus:Sex 79169.6 79169.6 1 59.00 10.8175 kept
## Pr(>F)
## Age:IPVstatus:Sex:Race 0.1265
## Age:IPVstatus:Sex 0.9669
## Age:Sex:Race 0.5026
## IPVstatus:Sex:Race 0.4799
## Sex:Race 0.9997
## Age:Sex 0.2488
## Age:IPVstatus:Race 0.0703
## Age:IPVstatus 0.8448
## Age:Race 0.6299
## IPVstatus:Race 0.3867
## Age 0.1469
## Race 0.0635
## IPVstatus 0.1730
## Sex 0.0136
## IPVstatus:Sex 0.0017
##
## Least squares means:
## IPVstatus Sex Estimate Standard Error DF t-value
## IPVstatus 0 1 NA 92.6 20.6 59 4.49
## IPVstatus 1 2 NA 141.8 29.1 59 4.87
## Sex Women NA 2 71.8 24.6 59 2.92
## Sex Men NA 1 162.6 25.8 59 6.29
## IPVstatus:Sex 0 Women 1 2 105.9 31.2 59 3.40
## IPVstatus:Sex 1 Women 2 2 37.8 38.2 59 0.99
## IPVstatus:Sex 0 Men 1 1 79.3 27.0 59 2.94
## IPVstatus:Sex 1 Men 2 1 245.9 44.1 59 5.58
## Lower CI Upper CI p-value
## IPVstatus 0 51.4 134 <2e-16
## IPVstatus 1 83.5 200 <2e-16
## Sex Women 22.5 121 0.0050
## Sex Men 110.9 214 <2e-16
## IPVstatus:Sex 0 Women 43.6 168 0.0012
## IPVstatus:Sex 1 Women -38.6 114 0.3266
## IPVstatus:Sex 0 Men 25.3 133 0.0047
## IPVstatus:Sex 1 Men 157.8 334 <2e-16
##
## Differences of LSMEANS:
## Estimate Standard Error DF t-value
## IPVstatus 0-1 -49.2 35.7 59.0 -1.38
## Sex Women-Men -90.8 35.7 59.0 -2.54
## IPVstatus:Sex 0 Women- 1 Women 68.2 49.3 59.0 1.38
## IPVstatus:Sex 0 Women- 0 Men 26.6 41.2 59.0 0.65
## IPVstatus:Sex 0 Women- 1 Men -140.0 54.0 59.0 -2.59
## IPVstatus:Sex 1 Women- 0 Men -41.5 46.7 59.0 -0.89
## IPVstatus:Sex 1 Women- 1 Men -208.2 58.3 59.0 -3.57
## IPVstatus:Sex 0 Men- 1 Men -166.7 51.7 59.0 -3.23
## Lower CI Upper CI p-value
## IPVstatus 0-1 -120.7 22.2 0.173
## Sex Women-Men -162.2 -19.4 0.014
## IPVstatus:Sex 0 Women- 1 Women -30.4 166.8 0.172
## IPVstatus:Sex 0 Women- 0 Men -55.9 109.1 0.521
## IPVstatus:Sex 0 Women- 1 Men -248.0 -32.0 0.012
## IPVstatus:Sex 1 Women- 0 Men -135.1 52.0 0.378
## IPVstatus:Sex 1 Women- 1 Men -324.8 -91.5 7e-04
## IPVstatus:Sex 0 Men- 1 Men -270.0 -63.3 0.002
##
## Final model:
## lme4::lmer(formula = IPVandCognitionDataSet2$"TrailsB-A" ~ IPVstatus +
## Sex + (1 | HNDid) + IPVstatus:Sex, data = IPVandCognitionDataSet2,
## REML = reml, contrasts = l)
Re-run the suggested final Model 1
(mm1 = lmer(IPVandCognitionDataSet2$"TrailsB-A" ~ IPVstatus + Sex + (1 | HNDid) +
IPVstatus:Sex, data = IPVandCognitionDataSet2, REML = F))
## Linear mixed model fit by maximum likelihood ['merModLmerTest']
## Formula: IPVandCognitionDataSet2$"TrailsB-A" ~ IPVstatus + Sex + (1 | HNDid) + IPVstatus:Sex
## Data: IPVandCognitionDataSet2
## AIC BIC logLik deviance
## 1582.3 1599.3 -785.1 1570.3
## Random effects:
## Groups Name Std.Dev.
## HNDid (Intercept) 113.4
## Residual 83.7
## Number of obs: 126, groups: HNDid, 63
## Fixed Effects:
## (Intercept) IPVstatus1 SexMen
## 105.9 -68.2 -26.6
## IPVstatus1:SexMen
## 234.8
summary(mm1)
## Linear mixed model fit by maximum likelihood ['merModLmerTest']
## Formula: IPVandCognitionDataSet2$"TrailsB-A" ~ IPVstatus + Sex + (1 | HNDid) + IPVstatus:Sex
## Data: IPVandCognitionDataSet2
##
## AIC BIC logLik deviance
## 1582.3 1599.3 -785.1 1570.3
##
## Random effects:
## Groups Name Variance Std.Dev.
## HNDid (Intercept) 12868 113.4
## Residual 6997 83.7
## Number of obs: 126, groups: HNDid, 63
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 105.9 30.2 63.0 3.51 0.00083
## IPVstatus1 -68.2 47.7 63.0 -1.43 0.15773
## SexMen -26.6 39.9 63.0 -0.67 0.50692
## IPVstatus1:SexMen 234.8 69.1 63.0 3.40 0.00118
##
## Correlation of Fixed Effects:
## (Intr) IPVst1 SexMen
## IPVstatus1 -0.632
## SexMen -0.756 0.478
## IPVstts1:SM 0.436 -0.690 -0.577
plot(st)
plot(mm1)
Trails B-A Regression Model 2 (with CES)
load("~/Desktop/Megan/Research/IPV and Cognition Paper/IPV R Output/IPVandCognitionDataSet2.rda")
library(lme4)
library(lmerTest)
(mm2 = lmer(IPVandCognitionDataSet2$"TrailsB-A" ~ (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: IPVandCognitionDataSet2$"TrailsB-A" ~ (Age + IPVstatus + Sex + Race + CES1)^5 + (Age | HNDid) + (1 | subclass)
## Data: IPVandCognitionDataSet2
## AIC BIC logLik deviance
## 1599.2 1704.2 -762.6 1525.2
## Random effects:
## Groups Name Std.Dev. Corr
## HNDid (Intercept) 1.50e+02
## Age 7.79e+00 1.00
## subclass (Intercept) 1.73e-03
## Residual 6.97e+01
## Number of obs: 126, groups: HNDid, 63; subclass, 21
## Fixed Effects:
## (Intercept)
## 31.54
## Age
## -1.87
## IPVstatus1
## 29.46
## SexMen
## 9.18
## RaceAfrAm
## 239.70
## CES11
## 5.46
## Age:IPVstatus1
## 4.20
## Age:SexMen
## 3.37
## Age:RaceAfrAm
## 8.13
## Age:CES11
## 1.87
## IPVstatus1:SexMen
## 314.97
## IPVstatus1:RaceAfrAm
## -279.15
## IPVstatus1:CES11
## -27.90
## SexMen:RaceAfrAm
## -160.70
## SexMen:CES11
## -12.71
## RaceAfrAm:CES11
## -282.77
## Age:IPVstatus1:SexMen
## 14.76
## Age:IPVstatus1:RaceAfrAm
## -11.76
## Age:IPVstatus1:CES11
## -4.27
## Age:SexMen:RaceAfrAm
## -5.17
## Age:SexMen:CES11
## -3.89
## Age:RaceAfrAm:CES11
## -15.07
## IPVstatus1:SexMen:RaceAfrAm
## 290.00
## IPVstatus1:SexMen:CES11
## -260.00
## IPVstatus1:RaceAfrAm:CES11
## 330.67
## SexMen:RaceAfrAm:CES11
## 358.45
## Age:IPVstatus1:SexMen:RaceAfrAm
## -11.67
## Age:IPVstatus1:SexMen:CES11
## -8.17
## Age:IPVstatus1:RaceAfrAm:CES11
## 21.23
## Age:SexMen:RaceAfrAm:CES11
## 27.66
## IPVstatus1:SexMen:RaceAfrAm:CES11
## -317.21
## Age:IPVstatus1:SexMen:RaceAfrAm:CES11
## -10.38
(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
## 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 effects:
## Chi.sq Chi.DF elim.num p.value
## (1 | subclass) 1.06 1 1 0.3034
## (1 | HNDid) 22.29 1 kept 0
##
## Fixed effects:
## Sum Sq Mean Sq NumDF DenDF F.value
## Age:IPVstatus:Sex:Race:CES1 6.453e+01 6.453e+01 1 56.67 0.0086
## IPVstatus:Sex:Race:CES1 5.221e+02 5.221e+02 1 66.60 0.0689
## Age:IPVstatus:Sex:CES1 1.372e+03 1.372e+03 1 77.51 0.1140
## Age:IPVstatus:Race:CES1 3.388e+02 3.388e+02 1 68.66 0.2573
## IPVstatus:Race:CES1 4.700e-03 4.700e-03 1 45.91 0.0479
## Age:IPVstatus:CES1 3.357e+02 3.357e+02 1 97.84 0.1720
## Age:Sex:Race:CES1 4.455e+03 4.455e+03 1 88.91 0.5259
## Sex:Race:CES1 2.289e+02 2.289e+02 1 43.28 0.0013
## Age:Race:CES1 5.904e+03 5.904e+03 1 94.09 0.1283
## Race:CES1 2.853e+03 2.853e+03 1 43.07 0.1403
## Age:Sex:CES1 1.303e+04 1.303e+04 1 101.45 0.7982
## Age:CES1 1.309e+03 1.309e+03 1 99.68 0.0118
## Age:IPVstatus:Sex:Race 5.760e+03 5.760e+03 1 101.74 1.0488
## Age:IPVstatus:Sex 7.245e+02 7.245e+02 1 103.96 0.0610
## Age:Sex:Race 3.866e+03 3.866e+03 1 106.05 0.1338
## Age:IPVstatus:Race 1.453e+04 1.453e+04 1 104.52 2.0156
## Age:Race 2.788e+03 2.788e+03 1 108.79 0.0015
## Age:IPVstatus 2.322e+03 2.322e+03 1 107.68 0.5199
## IPVstatus:Sex:Race 6.492e+03 6.492e+03 1 49.38 2.8579
## IPVstatus:Race 4.625e+03 4.625e+03 1 50.09 0.0014
## Sex:Race 3.201e+02 3.201e+02 1 51.23 0.0251
## Age:Sex 8.294e+03 8.294e+03 1 110.65 2.7394
## Age 2.085e+04 2.085e+04 1 114.10 2.0428
## IPVstatus 8.191e+03 8.191e+03 1 54.00 3.9893
## Sex 1.423e+04 1.423e+04 1 54.00 10.6464
## Race 2.880e+04 2.880e+04 1 54.00 6.1147
## CES1 5.494e+03 5.494e+03 1 54.00 2.8659
## IPVstatus:Sex 6.789e+04 6.789e+04 1 54.00 11.8153
## IPVstatus:CES1 4.938e+03 4.938e+03 1 54.00 0.7132
## Sex:CES1 6.993e+02 6.993e+02 1 54.00 1.3620
## IPVstatus:Sex:CES1 2.632e+04 2.632e+04 1 54.00 4.7126
## elim.num Pr(>F)
## Age:IPVstatus:Sex:Race:CES1 1 0.9265
## IPVstatus:Sex:Race:CES1 2 0.7937
## Age:IPVstatus:Sex:CES1 3 0.7366
## Age:IPVstatus:Race:CES1 4 0.6136
## IPVstatus:Race:CES1 5 0.8277
## Age:IPVstatus:CES1 6 0.6792
## Age:Sex:Race:CES1 7 0.4703
## Sex:Race:CES1 8 0.9715
## Age:Race:CES1 9 0.7210
## Race:CES1 10 0.7099
## Age:Sex:CES1 11 0.3737
## Age:CES1 12 0.9136
## Age:IPVstatus:Sex:Race 13 0.3082
## Age:IPVstatus:Sex 14 0.8054
## Age:Sex:Race 15 0.7153
## Age:IPVstatus:Race 16 0.1587
## Age:Race 17 0.9693
## Age:IPVstatus 18 0.4724
## IPVstatus:Sex:Race 19 0.0972
## IPVstatus:Race 20 0.9704
## Sex:Race 21 0.8747
## Age:Sex 22 0.1007
## Age 23 0.1557
## IPVstatus kept 0.0508
## Sex kept 0.0019
## Race kept 0.0166
## CES1 kept 0.0962
## IPVstatus:Sex kept 0.0011
## IPVstatus:CES1 kept 0.4021
## Sex:CES1 kept 0.2483
## IPVstatus:Sex:CES1 kept 0.0344
##
## Least squares means:
## IPVstatus Sex Race CES1 Estimate
## IPVstatus 0 1 NA NA NA 76.3
## IPVstatus 1 2 NA NA NA 147.5
## Sex Women NA 2 NA NA 53.6
## Sex Men NA 1 NA NA 170.2
## Race White NA NA 2 NA 68.9
## Race AfrAm NA NA 1 NA 154.9
## CES1 0 NA NA NA 1 142.1
## CES1 1 NA NA NA 2 81.6
## IPVstatus:Sex 0 Women 1 2 NA NA 79.2
## IPVstatus:Sex 1 Women 2 2 NA NA 28.0
## IPVstatus:Sex 0 Men 1 1 NA NA 73.3
## IPVstatus:Sex 1 Men 2 1 NA NA 267.0
## IPVstatus:CES1 0 0 1 NA NA 1 91.5
## IPVstatus:CES1 1 0 2 NA NA 1 192.8
## IPVstatus:CES1 0 1 1 NA NA 2 61.0
## IPVstatus:CES1 1 1 2 NA NA 2 102.2
## Sex:CES1 Women 0 NA 2 NA 1 62.8
## Sex:CES1 Men 0 NA 1 NA 1 221.5
## Sex:CES1 Women 1 NA 2 NA 2 44.4
## Sex:CES1 Men 1 NA 1 NA 2 118.8
## IPVstatus:Sex:CES1 0 Women 0 1 2 NA 1 113.2
## IPVstatus:Sex:CES1 1 Women 0 2 2 NA 1 12.4
## IPVstatus:Sex:CES1 0 Men 0 1 1 NA 1 69.8
## IPVstatus:Sex:CES1 1 Men 0 2 1 NA 1 373.2
## IPVstatus:Sex:CES1 0 Women 1 1 2 NA 2 45.3
## IPVstatus:Sex:CES1 1 Women 1 2 2 NA 2 43.6
## IPVstatus:Sex:CES1 0 Men 1 1 1 NA 2 76.8
## IPVstatus:Sex:CES1 1 Men 1 2 1 NA 2 160.8
## Standard Error DF t-value Lower CI Upper CI
## IPVstatus 0 21.4 54 3.57 33.42 119
## IPVstatus 1 29.2 54 5.04 88.85 206
## Sex Women 24.8 54 2.16 3.83 103
## Sex Men 26.4 54 6.44 117.19 223
## Race White 28.3 54 2.43 12.03 126
## Race AfrAm 21.9 54 7.08 111.06 199
## CES1 0 26.4 54 5.39 89.26 195
## CES1 1 24.9 54 3.28 31.71 132
## IPVstatus:Sex 0 Women 31.8 54 2.49 15.54 143
## IPVstatus:Sex 1 Women 37.3 54 0.75 -46.72 103
## IPVstatus:Sex 0 Men 27.6 54 2.66 17.98 129
## IPVstatus:Sex 1 Men 44.8 54 5.96 177.11 357
## IPVstatus:CES1 0 0 25.0 54 3.66 41.33 142
## IPVstatus:CES1 1 0 46.3 54 4.17 100.02 286
## IPVstatus:CES1 0 1 33.8 54 1.81 -6.73 129
## IPVstatus:CES1 1 1 35.5 54 2.88 31.03 173
## Sex:CES1 Women 0 35.0 54 1.80 -7.29 133
## Sex:CES1 Men 0 39.9 54 5.56 141.58 301
## Sex:CES1 Women 1 34.2 54 1.30 -24.11 113
## Sex:CES1 Men 1 35.1 54 3.39 48.46 189
## IPVstatus:Sex:CES1 0 Women 0 38.4 54 2.94 36.14 190
## IPVstatus:Sex:CES1 1 Women 0 57.5 54 0.22 -102.99 128
## IPVstatus:Sex:CES1 0 Men 0 31.7 54 2.20 6.20 133
## IPVstatus:Sex:CES1 1 Men 0 73.3 54 5.09 226.23 520
## IPVstatus:Sex:CES1 0 Women 1 49.4 54 0.92 -53.80 144
## IPVstatus:Sex:CES1 1 Women 1 47.9 54 0.91 -52.45 140
## IPVstatus:Sex:CES1 0 Men 1 45.0 54 1.71 -13.32 167
## IPVstatus:Sex:CES1 1 Men 1 52.9 54 3.04 54.67 267
## p-value
## IPVstatus 0 0.0008
## IPVstatus 1 <2e-16
## Sex Women 0.0353
## Sex Men <2e-16
## Race White 0.0185
## Race AfrAm <2e-16
## CES1 0 <2e-16
## CES1 1 0.0018
## IPVstatus:Sex 0 Women 0.0157
## IPVstatus:Sex 1 Women 0.4559
## IPVstatus:Sex 0 Men 0.0104
## IPVstatus:Sex 1 Men <2e-16
## IPVstatus:CES1 0 0 0.0006
## IPVstatus:CES1 1 0 0.0001
## IPVstatus:CES1 0 1 0.0765
## IPVstatus:CES1 1 1 0.0057
## Sex:CES1 Women 0 0.0780
## Sex:CES1 Men 0 <2e-16
## Sex:CES1 Women 1 0.1993
## Sex:CES1 Men 1 0.0013
## IPVstatus:Sex:CES1 0 Women 0 0.0048
## IPVstatus:Sex:CES1 1 Women 0 0.8304
## IPVstatus:Sex:CES1 0 Men 0 0.0321
## IPVstatus:Sex:CES1 1 Men 0 <2e-16
## IPVstatus:Sex:CES1 0 Women 1 0.3637
## IPVstatus:Sex:CES1 1 Women 1 0.3670
## IPVstatus:Sex:CES1 0 Men 1 0.0933
## IPVstatus:Sex:CES1 1 Men 1 0.0037
##
## Differences of LSMEANS:
## Estimate Standard Error DF t-value
## IPVstatus 0-1 -71.2 35.66 54.0 -2.00
## Sex Women-Men -116.5 35.71 54.0 -3.26
## Race White-AfrAm -86.0 34.79 54.0 -2.47
## CES1 0-1 60.5 35.75 54.0 1.69
## IPVstatus:Sex 0 Women- 1 Women 51.3 48.27 54.0 1.06
## IPVstatus:Sex 0 Women- 0 Men 6.0 41.41 54.0 0.14
## IPVstatus:Sex 0 Women- 1 Men -187.7 54.50 54.0 -3.45
## IPVstatus:Sex 1 Women- 0 Men -45.3 46.09 54.0 -0.98
## IPVstatus:Sex 1 Women- 1 Men -239.0 58.10 54.0 -4.11
## IPVstatus:Sex 0 Men- 1 Men -193.7 52.47 54.0 -3.69
## IPVstatus:CES1 0 0- 1 0 -101.3 52.45 54.0 -1.93
## IPVstatus:CES1 0 0- 0 1 30.5 41.36 54.0 0.74
## IPVstatus:CES1 0 0- 1 1 -10.7 43.06 54.0 -0.25
## IPVstatus:CES1 1 0- 0 1 131.7 56.96 54.0 2.31
## IPVstatus:CES1 1 0- 1 1 90.6 58.13 54.0 1.56
## IPVstatus:CES1 0 1- 1 1 -41.2 48.23 54.0 -0.85
## Sex:CES1 Women 0- Men 0 -158.7 53.27 54.0 -2.98
## Sex:CES1 Women 0- Women 1 18.4 48.12 54.0 0.38
## Sex:CES1 Women 0- Men 1 -56.0 48.30 54.0 -1.16
## Sex:CES1 Men 0- Women 1 177.1 52.67 54.0 3.36
## Sex:CES1 Men 0- Men 1 102.7 53.37 54.0 1.92
## Sex:CES1 Women 1- Men 1 -74.4 48.17 54.0 -1.54
## Lower CI Upper CI p-value
## IPVstatus 0-1 -142.71 0.27 0.051
## Sex Women-Men -188.13 -44.93 0.002
## Race White-AfrAm -155.79 -16.28 0.017
## CES1 0-1 -11.15 132.20 0.096
## IPVstatus:Sex 0 Women- 1 Women -45.51 148.04 0.293
## IPVstatus:Sex 0 Women- 0 Men -77.07 88.98 0.886
## IPVstatus:Sex 0 Women- 1 Men -297.01 -78.49 0.001
## IPVstatus:Sex 1 Women- 0 Men -137.71 47.08 0.330
## IPVstatus:Sex 1 Women- 1 Men -355.50 -122.53 1e-04
## IPVstatus:Sex 0 Men- 1 Men -298.89 -88.52 5e-04
## IPVstatus:CES1 0 0- 1 0 -206.43 3.87 0.059
## IPVstatus:CES1 0 0- 0 1 -52.47 113.39 0.465
## IPVstatus:CES1 0 0- 1 1 -97.04 75.64 0.805
## IPVstatus:CES1 1 0- 0 1 17.54 245.94 0.025
## IPVstatus:CES1 1 0- 1 1 -25.97 207.13 0.125
## IPVstatus:CES1 0 1- 1 1 -137.86 55.54 0.397
## Sex:CES1 Women 0- Men 0 -265.47 -51.87 0.004
## Sex:CES1 Women 0- Women 1 -78.09 114.85 0.704
## Sex:CES1 Women 0- Men 1 -152.84 40.82 0.251
## Sex:CES1 Men 0- Women 1 71.45 282.66 0.001
## Sex:CES1 Men 0- Men 1 -4.34 209.66 0.060
## Sex:CES1 Women 1- Men 1 -170.97 22.19 0.128
##
## Final model:
## lme4::lmer(formula = IPVandCognitionDataSet2$"TrailsB-A" ~ IPVstatus +
## Sex + Race + CES1 + (1 | HNDid) + IPVstatus:Sex + IPVstatus:CES1 +
## Sex:CES1 + IPVstatus:Sex:CES1, data = IPVandCognitionDataSet2,
## REML = reml, contrasts = l)
Re-run the suggested final Model 2
(mm2 = lmer(IPVandCognitionDataSet2$"TrailsB-A" ~ IPVstatus + Sex + Race + CES1 +
(1 | HNDid) + IPVstatus:Sex + IPVstatus:CES1 + Sex:CES1 + IPVstatus:Sex:CES1,
data = IPVandCognitionDataSet2, REML = F))
## Linear mixed model fit by maximum likelihood ['merModLmerTest']
## Formula: IPVandCognitionDataSet2$"TrailsB-A" ~ IPVstatus + Sex + Race + CES1 + (1 | HNDid) + IPVstatus:Sex + IPVstatus:CES1 + Sex:CES1 + IPVstatus:Sex:CES1
## Data: IPVandCognitionDataSet2
## AIC BIC logLik deviance
## 1581.2 1612.4 -779.6 1559.2
## Random effects:
## Groups Name Std.Dev.
## HNDid (Intercept) 101.1
## Residual 83.6
## Number of obs: 126, groups: HNDid, 63
## Fixed Effects:
## (Intercept) IPVstatus1 SexMen
## 70.2 -100.8 -43.4
## RaceAfrAm CES11 IPVstatus1:SexMen
## 86.0 -67.9 404.2
## IPVstatus1:CES11 SexMen:CES11 IPVstatus1:SexMen:CES11
## 99.1 75.0 -318.5
summary(mm2)
## Linear mixed model fit by maximum likelihood ['merModLmerTest']
## Formula: IPVandCognitionDataSet2$"TrailsB-A" ~ IPVstatus + Sex + Race + CES1 + (1 | HNDid) + IPVstatus:Sex + IPVstatus:CES1 + Sex:CES1 + IPVstatus:Sex:CES1
## Data: IPVandCognitionDataSet2
##
## AIC BIC logLik deviance
## 1581.2 1612.4 -779.6 1559.2
##
## Random effects:
## Groups Name Variance Std.Dev.
## HNDid (Intercept) 10228 101.1
## Residual 6997 83.6
## Number of obs: 126, groups: HNDid, 63
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 70.2 40.8 63.0 1.72 0.09053
## IPVstatus1 -100.8 63.4 63.0 -1.59 0.11681
## SexMen -43.4 46.0 63.0 -0.95 0.34805
## RaceAfrAm 86.0 32.2 63.0 2.67 0.00961
## CES11 -67.9 57.1 63.0 -1.19 0.23844
## IPVstatus1:SexMen 404.2 97.9 63.0 4.13 0.00011
## IPVstatus1:CES11 99.1 91.0 63.0 1.09 0.28008
## SexMen:CES11 75.0 76.2 63.0 0.98 0.32898
## IPVstatus1:SexMen:CES11 -318.5 135.8 63.0 -2.34 0.02220
##
## Correlation of Fixed Effects:
## (Intr) IPVst1 SexMen RcAfrA CES11 IPVs1:SM IPV1:C SM:CES
## IPVstatus1 -0.440
## SexMen -0.691 0.424
## RaceAfrAm -0.502 -0.083 0.052
## CES11 -0.473 0.355 0.469 -0.125
## IPVstts1:SM 0.247 -0.654 -0.461 0.129 -0.239
## IPVs1:CES11 0.231 -0.710 -0.288 0.210 -0.644 0.475
## SexMn:CES11 0.367 -0.264 -0.598 0.067 -0.746 0.291 0.477
## IPV1:SM:CES -0.102 0.484 0.325 -0.244 0.444 -0.740 -0.692 -0.575
plot(st)
plot(mm2)
