Trails B 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 REML ['merModLmerTest']
## Formula: TrailsBtestSec ~ (Age + IPVstatus + Sex + PovStat)^4 + (Age | HNDid) + (1 | subclass)
## Data: IPVandCognitionDataSet2
## REML criterion at convergence: 1411
## Random effects:
## Groups Name Std.Dev. Corr
## HNDid (Intercept) 1.74e+02
## Age 5.95e+00 1.00
## subclass (Intercept) 1.25e-04
## Residual 6.60e+01
## Number of obs: 126, groups: HNDid, 63; subclass, 21
## Fixed Effects:
## (Intercept) Age
## 127.854 0.952
## IPVstatus1 SexMen
## -51.092 -54.043
## PovStatBelow Age:IPVstatus1
## 124.524 -0.129
## Age:SexMen Age:PovStatBelow
## -0.535 11.093
## IPVstatus1:SexMen IPVstatus1:PovStatBelow
## 164.986 -101.976
## SexMen:PovStatBelow Age:IPVstatus1:SexMen
## 71.367 -11.199
## Age:IPVstatus1:PovStatBelow Age:SexMen:PovStatBelow
## -7.372 4.763
## IPVstatus1:SexMen:PovStatBelow Age:IPVstatus1:SexMen:PovStatBelow
## 62.968 21.256
## 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 (1 | subclass) was eliminated because of standard deviation being equal to 0
## 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 | HNDid) 52.6 1 kept < 1e-07
##
## Fixed effects:
## Sum Sq Mean Sq NumDF DenDF F.value elim.num
## Age:IPVstatus:Sex:PovStat 2974.3 2974.3 1 102.77 0.6242 1
## Age:IPVstatus:Sex 430.0 430.0 1 106.66 0.1167 2
## Age:IPVstatus:PovStat 2480.0 2480.0 1 108.43 0.1613 3
## Age:IPVstatus 346.7 346.7 1 110.58 0.0162 4
## IPVstatus:Sex:PovStat 2396.8 2396.8 1 51.85 0.5536 5
## IPVstatus:PovStat 4826.9 4826.9 1 52.83 0.9661 6
## Age:Sex:PovStat 10702.8 10702.8 1 115.74 2.6506 7
## Sex:PovStat 678.6 678.6 1 60.98 0.1041 8
## Age:Sex 6648.0 6648.0 1 117.55 1.0530 9
## Age 42780.8 42780.8 1 119.00 11.0562 kept
## IPVstatus 3978.6 3978.6 1 55.63 1.8726 kept
## Sex 14437.0 14437.0 1 56.78 4.0286 kept
## PovStat 11330.8 11330.8 1 77.79 8.0366 kept
## Age:PovStat 55425.7 55425.7 1 118.97 12.1417 kept
## IPVstatus:Sex 34937.8 34937.8 1 55.63 7.2835 kept
## Pr(>F)
## Age:IPVstatus:Sex:PovStat 0.4313
## Age:IPVstatus:Sex 0.7333
## Age:IPVstatus:PovStat 0.6887
## Age:IPVstatus 0.8990
## IPVstatus:Sex:PovStat 0.4602
## IPVstatus:PovStat 0.3301
## Age:Sex:PovStat 0.1062
## Sex:PovStat 0.7481
## Age:Sex 0.3069
## Age 0.0012
## IPVstatus 0.1767
## Sex 0.0495
## PovStat 0.0058
## Age:PovStat 0.0007
## IPVstatus:Sex 0.0092
##
## Least squares means:
## IPVstatus Sex PovStat Estimate Standard Error DF
## IPVstatus 0 1.0 NA NA 118.2 23.2 57.9
## IPVstatus 1 2.0 NA NA 169.6 31.5 56.4
## Sex Women NA 2.0 NA 105.7 26.4 55.7
## Sex Men NA 1.0 NA 182.1 29.3 59.0
## PovStat Above NA NA 1.0 127.6 22.7 56.6
## PovStat Below NA NA 2.0 160.2 32.7 58.8
## IPVstatus:Sex 0 Women 1.0 2.0 NA 130.9 30.6 55.9
## IPVstatus:Sex 1 Women 2.0 2.0 NA 80.4 41.4 55.6
## IPVstatus:Sex 0 Men 1.0 1.0 NA 105.5 32.9 58.3
## IPVstatus:Sex 1 Men 2.0 1.0 NA 258.8 47.0 56.8
## t-value Lower CI Upper CI p-value
## IPVstatus 0 5.09 71.73 165 <2e-16
## IPVstatus 1 5.38 106.45 233 <2e-16
## Sex Women 4.01 52.84 158 0.0002
## Sex Men 6.22 123.54 241 <2e-16
## PovStat Above 5.62 82.12 173 <2e-16
## PovStat Below 4.90 94.80 226 <2e-16
## IPVstatus:Sex 0 Women 4.28 69.59 192 0.0001
## IPVstatus:Sex 1 Women 1.94 -2.48 163 0.0570
## IPVstatus:Sex 0 Men 3.21 39.60 171 0.0022
## IPVstatus:Sex 1 Men 5.50 164.61 353 <2e-16
##
## Differences of LSMEANS:
## Estimate Standard Error DF t-value
## IPVstatus 0-1 -51.4 37.6 55.6 -1.37
## Sex Women-Men -76.5 38.1 56.8 -2.01
## PovStat Above-Below -32.6 38.9 58.1 -0.84
## IPVstatus:Sex 0 Women- 1 Women 50.5 50.2 55.8 1.01
## IPVstatus:Sex 0 Women- 0 Men 25.5 43.4 56.4 0.59
## IPVstatus:Sex 0 Women- 1 Men -127.9 55.6 56.0 -2.30
## IPVstatus:Sex 1 Women- 0 Men -25.0 51.3 56.5 -0.49
## IPVstatus:Sex 1 Women- 1 Men -178.4 62.2 56.1 -2.87
## IPVstatus:Sex 0 Men- 1 Men -153.3 56.2 55.5 -2.73
## Lower CI Upper CI p-value
## IPVstatus 0-1 -126.7 23.864 0.177
## Sex Women-Men -152.8 -0.172 0.050
## PovStat Above-Below -110.4 45.303 0.406
## IPVstatus:Sex 0 Women- 1 Women -50.0 150.996 0.319
## IPVstatus:Sex 0 Women- 0 Men -61.5 112.403 0.560
## IPVstatus:Sex 0 Women- 1 Men -239.3 -16.519 0.025
## IPVstatus:Sex 1 Women- 0 Men -127.9 77.789 0.628
## IPVstatus:Sex 1 Women- 1 Men -303.0 -53.762 0.006
## IPVstatus:Sex 0 Men- 1 Men -266.0 -40.713 0.009
##
## Final model:
## lme4::lmer(formula = TrailsBtestSec ~ Age + IPVstatus + Sex +
## PovStat + (1 | HNDid) + Age:PovStat + IPVstatus:Sex, data = IPVandCognitionDataSet2,
## REML = reml, contrasts = l)
Re-run the suggested final Model 1
(mm1 = lmer(TrailsBtestSec ~ Age + IPVstatus + Sex + PovStat + (Age | HNDid) +
(1 | subclass) + Age:PovStat + IPVstatus:Sex, data = IPVandCognitionDataSet2))
## Warning: number of observations <= rank(Z); variance-covariance matrix
## will be unidentifiable
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: TrailsBtestSec ~ Age + IPVstatus + Sex + PovStat + (Age | HNDid) + (1 | subclass) + Age:PovStat + IPVstatus:Sex
## Data: IPVandCognitionDataSet2
## REML criterion at convergence: 1487
## Random effects:
## Groups Name Std.Dev. Corr
## HNDid (Intercept) 163.1942
## Age 4.8649 0.93
## subclass (Intercept) 0.0252
## Residual 63.8126
## Number of obs: 126, groups: HNDid, 63; subclass, 21
## Fixed Effects:
## (Intercept) Age IPVstatus1
## 119.945 -0.423 -57.966
## SexMen PovStatBelow Age:PovStatBelow
## -44.287 141.345 14.785
## IPVstatus1:SexMen
## 232.299
summary(mm1)
## 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: TrailsBtestSec ~ Age + IPVstatus + Sex + PovStat + (Age | HNDid) + (1 | subclass) + Age:PovStat + IPVstatus:Sex
## Data: IPVandCognitionDataSet2
##
## REML criterion at convergence: 1487
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 2.66e+04 163.1942
## Age 2.37e+01 4.8649 0.93
## subclass (Intercept) 6.33e-04 0.0252
## Residual 4.07e+03 63.8126
## Number of obs: 126, groups: HNDid, 63; subclass, 21
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 119.945 41.258 68.400 2.91 0.00491
## Age -0.423 2.436 58.200 -0.17 0.86285
## IPVstatus1 -57.966 46.782 55.700 -1.24 0.22052
## SexMen -44.287 40.603 55.000 -1.09 0.28014
## PovStatBelow 141.345 54.448 35.700 2.60 0.01361
## Age:PovStatBelow 14.785 4.076 43.300 3.63 0.00075
## IPVstatus1:SexMen 232.299 71.520 56.300 3.25 0.00196
##
## Correlation of Fixed Effects:
## (Intr) Age IPVst1 SexMen PvSttB Ag:PSB
## Age 0.709
## IPVstatus1 -0.447 -0.091
## SexMen -0.479 -0.030 0.387
## PovStatBelw -0.511 -0.514 0.052 -0.011
## Ag:PvSttBlw -0.363 -0.595 0.017 -0.109 0.732
## IPVstts1:SM 0.319 0.062 -0.659 -0.570 -0.075 0.011
plot(st)
plot(mm1)
Trails B Regression Model 2 (with CES)
load("~/Desktop/Megan/Research/IPV and Cognition Paper/IPV R Output/IPVandCognitionDataSet2.rda")
library(lme4)
library(lmerTest)
(mm2 = lmer(TrailsBtestSec ~ (Age + IPVstatus + Sex + PovStat + CES1)^5 + (Age |
HNDid) + (1 | subclass), data = IPVandCognitionDataSet2))
## Warning: number of observations <= rank(Z); variance-covariance matrix
## will be unidentifiable
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: TrailsBtestSec ~ (Age + IPVstatus + Sex + PovStat + CES1)^5 + (Age | HNDid) + (1 | subclass)
## Data: IPVandCognitionDataSet2
## REML criterion at convergence: 1243
## Random effects:
## Groups Name Std.Dev. Corr
## HNDid (Intercept) 1.90e+02
## Age 7.31e+00 1.00
## subclass (Intercept) 1.89e-04
## Residual 6.16e+01
## Number of obs: 126, groups: HNDid, 63; subclass, 21
## Fixed Effects:
## (Intercept)
## 62.863
## Age
## -1.119
## IPVstatus1
## 18.196
## SexMen
## -19.207
## PovStatBelow
## 387.716
## CES11
## 132.144
## Age:IPVstatus1
## 2.308
## Age:SexMen
## -0.482
## Age:PovStatBelow
## 27.329
## Age:CES11
## 3.974
## IPVstatus1:SexMen
## 451.148
## IPVstatus1:PovStatBelow
## -388.775
## IPVstatus1:CES11
## -140.541
## SexMen:PovStatBelow
## -181.767
## SexMen:CES11
## 10.535
## PovStatBelow:CES11
## -485.694
## Age:IPVstatus1:SexMen
## -0.707
## Age:IPVstatus1:PovStatBelow
## -26.518
## Age:IPVstatus1:CES11
## -4.763
## Age:SexMen:PovStatBelow
## -15.055
## Age:SexMen:CES11
## 6.058
## Age:PovStatBelow:CES11
## -27.443
## IPVstatus1:SexMen:PovStatBelow
## 52.129
## IPVstatus1:SexMen:CES11
## -290.507
## IPVstatus1:PovStatBelow:CES11
## 523.740
## SexMen:PovStatBelow:CES11
## 426.090
## Age:IPVstatus1:SexMen:PovStatBelow
## 30.018
## Age:IPVstatus1:SexMen:CES11
## -3.703
## Age:IPVstatus1:PovStatBelow:CES11
## 32.779
## Age:SexMen:PovStatBelow:CES11
## 37.822
## IPVstatus1:SexMen:PovStatBelow:CES11
## -319.157
## Age:IPVstatus1:SexMen:PovStatBelow:CES11
## -17.745
(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 (1 | subclass) was eliminated because of standard deviation being equal to 0
## 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 | HNDid) 51.67 1 kept < 1e-07
##
## Fixed effects:
## Sum Sq Mean Sq NumDF DenDF F.value
## Age:IPVstatus:Sex:PovStat:CES1 2.440e+02 2.440e+02 1 63.08 0.0569
## IPVstatus:Sex:PovStat:CES1 1.047e+03 1.047e+03 1 77.22 0.2504
## Age:IPVstatus:Sex:PovStat 5.931e+03 5.931e+03 1 82.77 0.6064
## Age:IPVstatus:Sex:CES1 2.370e+02 2.370e+02 1 77.92 0.8054
## Age:IPVstatus:Sex 4.542e+01 4.542e+01 1 70.96 0.1278
## Age:IPVstatus:PovStat:CES1 6.921e+03 6.921e+03 1 85.41 0.3899
## Age:IPVstatus:CES1 3.239e+03 3.239e+03 1 93.52 0.0005
## Age:IPVstatus:PovStat 2.512e+03 2.512e+03 1 92.40 0.0189
## Age:IPVstatus 2.118e+02 2.118e+02 1 93.68 0.0002
## IPVstatus:PovStat:CES1 8.882e+02 8.882e+02 1 47.24 0.4988
## IPVstatus:Sex:PovStat 2.793e+03 2.793e+03 1 48.62 2.2838
## IPVstatus:PovStat 4.086e+03 4.086e+03 1 49.06 1.4550
## IPVstatus:Sex:CES1 7.732e+03 7.732e+03 1 50.17 1.6633
## IPVstatus:CES1 9.016e+03 9.016e+03 1 50.75 2.5003
## Age:Sex:PovStat:CES1 9.370e+03 9.370e+03 1 88.53 2.9682
## Sex:PovStat:CES1 1.374e+03 1.374e+03 1 65.23 0.0246
## Age:PovStat:CES1 1.968e+03 1.968e+03 1 101.24 0.3115
## Age:Sex:PovStat 8.945e+03 8.945e+03 1 108.62 1.6973
## Sex:PovStat 3.107e+01 3.107e+01 1 60.29 0.0233
## Age 4.103e+04 4.103e+04 1 107.55 17.0879
## IPVstatus 3.270e+03 3.270e+03 1 53.71 1.3545
## Sex 1.183e+04 1.183e+04 1 81.95 4.1830
## PovStat 9.115e+03 9.115e+03 1 80.06 7.3063
## CES1 2.349e+02 2.349e+02 1 77.23 0.3959
## Age:Sex 5.742e+03 5.742e+03 1 106.60 1.9424
## Age:PovStat 5.378e+04 5.378e+04 1 111.09 15.2043
## Age:CES1 9.716e-01 9.716e-01 1 104.07 0.4527
## IPVstatus:Sex 2.923e+04 2.923e+04 1 54.13 4.8678
## Sex:CES1 1.935e+03 1.935e+03 1 85.59 1.2459
## PovStat:CES1 7.072e+03 7.072e+03 1 59.47 4.1587
## Age:Sex:CES1 2.204e+04 2.204e+04 1 103.56 5.3869
## elim.num Pr(>F)
## Age:IPVstatus:Sex:PovStat:CES1 1 0.8123
## IPVstatus:Sex:PovStat:CES1 2 0.6182
## Age:IPVstatus:Sex:PovStat 3 0.4384
## Age:IPVstatus:Sex:CES1 4 0.3723
## Age:IPVstatus:Sex 5 0.7218
## Age:IPVstatus:PovStat:CES1 6 0.5340
## Age:IPVstatus:CES1 7 0.9821
## Age:IPVstatus:PovStat 8 0.8910
## Age:IPVstatus 9 0.9876
## IPVstatus:PovStat:CES1 10 0.4835
## IPVstatus:Sex:PovStat 11 0.1372
## IPVstatus:PovStat 12 0.2335
## IPVstatus:Sex:CES1 13 0.2031
## IPVstatus:CES1 14 0.1200
## Age:Sex:PovStat:CES1 15 0.0884
## Sex:PovStat:CES1 16 0.8760
## Age:PovStat:CES1 17 0.5780
## Age:Sex:PovStat 18 0.1954
## Sex:PovStat 19 0.8791
## Age kept 1e-04
## IPVstatus kept 0.2496
## Sex kept 0.0440
## PovStat kept 0.0084
## CES1 kept 0.5311
## Age:Sex kept 0.1663
## Age:PovStat kept 0.0002
## Age:CES1 kept 0.5025
## IPVstatus:Sex kept 0.0316
## Sex:CES1 kept 0.2675
## PovStat:CES1 kept 0.0459
## Age:Sex:CES1 kept 0.0222
##
## Least squares means:
## IPVstatus Sex PovStat CES1 Estimate Standard Error
## IPVstatus 0 1.0 NA NA NA 111.1 24.8
## IPVstatus 1 2.0 NA NA NA 157.9 32.7
## Sex Women NA 2.0 NA NA 105.8 26.9
## Sex Men NA 1.0 NA NA 163.1 30.2
## PovStat Above NA NA 1.0 NA 125.8 23.0
## PovStat Below NA NA 2.0 NA 143.1 33.8
## CES1 0 NA NA NA 1.0 159.1 28.0
## CES1 1 NA NA NA 2.0 109.8 30.3
## IPVstatus:Sex 0 Women 1.0 2.0 NA NA 127.3 31.1
## IPVstatus:Sex 1 Women 2.0 2.0 NA NA 84.3 42.1
## IPVstatus:Sex 0 Men 1.0 1.0 NA NA 94.8 36.5
## IPVstatus:Sex 1 Men 2.0 1.0 NA NA 231.4 49.7
## Sex:CES1 Women 0 NA 2.0 NA 1.0 124.1 38.1
## Sex:CES1 Men 0 NA 1.0 NA 1.0 194.2 39.9
## Sex:CES1 Women 1 NA 2.0 NA 2.0 87.5 36.6
## Sex:CES1 Men 1 NA 1.0 NA 2.0 132.1 45.1
## PovStat:CES1 Above 0 NA NA 1.0 1.0 109.3 32.3
## PovStat:CES1 Below 0 NA NA 2.0 1.0 208.9 44.4
## PovStat:CES1 Above 1 NA NA 1.0 2.0 142.3 33.1
## PovStat:CES1 Below 1 NA NA 2.0 2.0 77.3 50.2
## DF t-value Lower CI Upper CI p-value
## IPVstatus 0 56.0 4.49 61.475 161 <2e-16
## IPVstatus 1 56.2 4.83 92.326 223 <2e-16
## Sex Women 54.4 3.94 51.957 160 0.0002
## Sex Men 60.5 5.39 102.653 224 <2e-16
## PovStat Above 54.3 5.48 79.747 172 <2e-16
## PovStat Below 60.1 4.23 75.468 211 0.0001
## CES1 0 54.4 5.68 102.930 215 <2e-16
## CES1 1 59.9 3.62 49.156 170 0.0006
## IPVstatus:Sex 0 Women 53.9 4.09 64.957 190 0.0001
## IPVstatus:Sex 1 Women 53.8 2.00 -0.193 169 0.0505
## IPVstatus:Sex 0 Men 55.9 2.60 21.626 168 0.0121
## IPVstatus:Sex 1 Men 58.0 4.66 131.990 331 <2e-16
## Sex:CES1 Women 0 55.7 3.25 47.651 200 0.0019
## Sex:CES1 Men 0 54.8 4.87 114.258 274 <2e-16
## Sex:CES1 Women 1 54.1 2.39 14.070 161 0.0204
## Sex:CES1 Men 1 60.2 2.93 41.889 222 0.0048
## PovStat:CES1 Above 0 55.3 3.38 44.552 174 0.0013
## PovStat:CES1 Below 0 56.0 4.70 119.908 298 <2e-16
## PovStat:CES1 Above 1 53.0 4.29 75.808 209 0.0001
## PovStat:CES1 Below 1 63.1 1.54 -22.878 178 0.1281
##
## Differences of LSMEANS:
## Estimate Standard Error DF t-value
## IPVstatus 0-1 -46.8 40.21 53.7 -1.16
## Sex Women-Men -57.4 39.01 56.9 -1.47
## PovStat Above-Below -17.3 39.93 57.8 -0.43
## CES1 0-1 49.3 40.75 56.1 1.21
## IPVstatus:Sex 0 Women- 1 Women 43.0 50.98 53.2 0.84
## IPVstatus:Sex 0 Women- 0 Men 32.5 46.40 53.9 0.70
## IPVstatus:Sex 0 Women- 1 Men -104.2 57.82 55.8 -1.80
## IPVstatus:Sex 1 Women- 0 Men -10.6 54.17 54.5 -0.19
## IPVstatus:Sex 1 Women- 1 Men -147.2 64.85 56.2 -2.27
## IPVstatus:Sex 0 Men- 1 Men -136.6 62.84 54.4 -2.17
## Sex:CES1 Women 0- Men 0 -70.1 54.28 56.2 -1.29
## Sex:CES1 Women 0- Women 1 36.6 52.04 55.5 0.70
## Sex:CES1 Women 0- Men 1 -8.1 59.24 59.2 -0.14
## Sex:CES1 Men 0- Women 1 106.7 53.44 53.2 2.00
## Sex:CES1 Men 0- Men 1 62.1 59.92 55.2 1.04
## Sex:CES1 Women 1- Men 1 -44.6 55.46 55.1 -0.80
## PovStat:CES1 Above 0- Below 0 -99.6 53.81 57.3 -1.85
## PovStat:CES1 Above 0- Above 1 -33.0 46.63 54.0 -0.71
## PovStat:CES1 Above 0- Below 1 32.0 58.80 58.9 0.54
## PovStat:CES1 Below 0- Above 1 66.7 55.25 54.8 1.21
## PovStat:CES1 Below 0- Below 1 131.6 66.35 59.7 1.98
## PovStat:CES1 Above 1- Below 1 64.9 59.58 59.8 1.09
## Lower CI Upper CI p-value
## IPVstatus 0-1 -127.429 33.83 0.25
## Sex Women-Men -135.488 20.77 0.15
## PovStat Above-Below -97.279 62.59 0.67
## CES1 0-1 -32.322 130.93 0.23
## IPVstatus:Sex 0 Women- 1 Women -59.224 145.24 0.40
## IPVstatus:Sex 0 Women- 0 Men -60.579 125.48 0.49
## IPVstatus:Sex 0 Women- 1 Men -219.997 11.68 0.08
## IPVstatus:Sex 1 Women- 0 Men -119.149 98.03 0.85
## IPVstatus:Sex 1 Women- 1 Men -277.064 -17.27 0.03
## IPVstatus:Sex 0 Men- 1 Men -262.582 -10.64 0.03
## Sex:CES1 Women 0- Men 0 -178.833 38.62 0.20
## Sex:CES1 Women 0- Women 1 -67.714 140.84 0.49
## Sex:CES1 Women 0- Men 1 -126.589 110.48 0.89
## Sex:CES1 Men 0- Women 1 -0.504 213.84 0.05
## Sex:CES1 Men 0- Men 1 -58.023 182.12 0.30
## Sex:CES1 Women 1- Men 1 -155.748 66.52 0.42
## PovStat:CES1 Above 0- Below 0 -207.369 8.11 0.07
## PovStat:CES1 Above 0- Above 1 -126.474 60.52 0.48
## PovStat:CES1 Above 0- Below 1 -85.703 149.63 0.59
## PovStat:CES1 Below 0- Above 1 -44.079 177.38 0.23
## PovStat:CES1 Below 0- Below 1 -1.146 264.33 0.05
## PovStat:CES1 Above 1- Below 1 -54.244 184.12 0.28
##
## Final model:
## lme4::lmer(formula = TrailsBtestSec ~ Age + IPVstatus + Sex +
## PovStat + CES1 + (1 | HNDid) + Age:Sex + Age:PovStat + Age:CES1 +
## IPVstatus:Sex + Sex:CES1 + PovStat:CES1 + Age:Sex:CES1, data = IPVandCognitionDataSet2,
## REML = reml, contrasts = l)
Re-run the suggested final Model 2
(mm2 = lmer(TrailsBtestSec ~ Age + IPVstatus + Sex + PovStat + CES1 + (Age |
HNDid) + (1 | subclass) + Age:Sex + Age:PovStat + Age:CES1 + IPVstatus:Sex +
Sex:CES1 + PovStat:CES1 + Age:Sex:CES1, data = IPVandCognitionDataSet2))
## Warning: number of observations <= rank(Z); variance-covariance matrix
## will be unidentifiable
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: TrailsBtestSec ~ Age + IPVstatus + Sex + PovStat + CES1 + (Age | HNDid) + (1 | subclass) + Age:Sex + Age:PovStat + Age:CES1 + IPVstatus:Sex + Sex:CES1 + PovStat:CES1 + Age:Sex:CES1
## Data: IPVandCognitionDataSet2
## REML criterion at convergence: 1434
## Random effects:
## Groups Name Std.Dev. Corr
## HNDid (Intercept) 169.93
## Age 6.79 0.81
## subclass (Intercept) 0.00
## Residual 59.51
## Number of obs: 126, groups: HNDid, 63; subclass, 21
## Fixed Effects:
## (Intercept) Age IPVstatus1
## 115.48 1.42 -54.49
## SexMen PovStatBelow CES11
## -70.96 204.51 2.33
## Age:SexMen Age:PovStatBelow Age:CES11
## -4.85 15.54 -5.26
## IPVstatus1:SexMen SexMen:CES11 PovStatBelow:CES11
## 216.59 98.21 -139.52
## Age:SexMen:CES11
## 17.63
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: TrailsBtestSec ~ Age + IPVstatus + Sex + PovStat + CES1 + (Age | HNDid) + (1 | subclass) + Age:Sex + Age:PovStat + Age:CES1 + IPVstatus:Sex + Sex:CES1 + PovStat:CES1 + Age:Sex:CES1
## Data: IPVandCognitionDataSet2
##
## REML criterion at convergence: 1434
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 28876.8 169.93
## Age 46.1 6.79 0.81
## subclass (Intercept) 0.0 0.00
## Residual 3541.0 59.51
## Number of obs: 126, groups: HNDid, 63; subclass, 21
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 115.48 61.03 48.80 1.89 0.06440
## Age 1.42 4.03 46.60 0.35 0.72612
## IPVstatus1 -54.49 48.22 51.00 -1.13 0.26372
## SexMen -70.96 77.17 46.30 -0.92 0.36253
## PovStatBelow 204.51 67.66 55.80 3.02 0.00378
## CES11 2.33 79.53 47.30 0.03 0.97671
## Age:SexMen -4.85 5.33 48.20 -0.91 0.36721
## Age:PovStatBelow 15.54 4.22 44.60 3.68 0.00062
## Age:CES11 -5.26 5.48 46.90 -0.96 0.34233
## IPVstatus1:SexMen 216.59 77.64 50.70 2.79 0.00741
## SexMen:CES11 98.21 113.06 44.10 0.87 0.38972
## PovStatBelow:CES11 -139.52 78.17 57.90 -1.78 0.07953
## Age:SexMen:CES11 17.63 8.41 55.40 2.10 0.04075
##
## Correlation of Fixed Effects:
## (Intr) Age IPVst1 SexMen PvSttB CES11 Ag:SxM Ag:PSB A:CES1
## Age 0.782
## IPVstatus1 -0.284 -0.047
## SexMen -0.673 -0.542 0.206
## PovStatBelw -0.329 -0.179 0.058 -0.107
## CES11 -0.653 -0.539 -0.072 0.483 0.118
## Age:SexMen -0.499 -0.673 0.023 0.780 -0.166 0.396
## Ag:PvSttBlw -0.184 -0.265 0.011 -0.139 0.659 -0.043 -0.111
## Age:CES11 -0.502 -0.658 -0.048 0.422 -0.058 0.771 0.525 -0.087
## IPVstts1:SM 0.216 0.027 -0.631 -0.191 -0.141 0.008 0.061 0.006 0.029
## SexMn:CES11 0.395 0.375 0.068 -0.609 0.093 -0.614 -0.552 0.042 -0.545
## PvStB:CES11 0.198 0.025 -0.052 0.020 -0.547 -0.278 0.092 -0.079 0.011
## Ag:SM:CES11 0.303 0.429 0.038 -0.477 0.101 -0.460 -0.634 0.051 -0.651
## IPV1:S SM:CES PSB:CE
## Age
## IPVstatus1
## SexMen
## PovStatBelw
## CES11
## Age:SexMen
## Ag:PvSttBlw
## Age:CES11
## IPVstts1:SM
## SexMn:CES11 -0.236
## PvStB:CES11 0.155 -0.123
## Ag:SM:CES11 -0.103 0.761 -0.158
plot(st)
plot(mm2)
