Trails 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: TrailsAtestSec ~ (Age + IPVstatus + Sex + PovStat)^4 + (Age | HNDid) + (1 | subclass)
## Data: IPVandCognitionDataSet2
## AIC BIC logLik deviance
## 960.9 1020.5 -459.5 918.9
## Random effects:
## Groups Name Std.Dev. Corr
## HNDid (Intercept) 6.623
## Age 0.187 1.00
## subclass (Intercept) 4.917
## Residual 7.190
## Number of obs: 126, groups: HNDid, 63; subclass, 21
## Fixed Effects:
## (Intercept) Age
## 36.42042 0.71803
## IPVstatus1 SexMen
## -1.94609 -7.47931
## PovStatBelow Age:IPVstatus1
## 3.48355 -0.04630
## Age:SexMen Age:PovStatBelow
## -0.57003 -0.20116
## IPVstatus1:SexMen IPVstatus1:PovStatBelow
## 3.02145 23.49130
## SexMen:PovStatBelow Age:IPVstatus1:SexMen
## 9.52528 0.00969
## Age:IPVstatus1:PovStatBelow Age:SexMen:PovStatBelow
## 2.81165 0.96526
## IPVstatus1:SexMen:PovStatBelow Age:IPVstatus1:SexMen:PovStatBelow
## -10.30524 -2.19099
## 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) 4.34 1 kept 0.0372
## (1 | HNDid) 5.59 1 kept 0.0181
##
## Fixed effects:
## Sum Sq Mean Sq NumDF DenDF F.value elim.num
## Age:IPVstatus:Sex:PovStat 125.498 125.498 1 97.94 2.3080 1
## IPVstatus:Sex:PovStat 4.208 4.208 1 55.51 0.0792 2
## Age:Sex:PovStat 9.024 9.024 1 80.78 0.1809 3
## Age:IPVstatus:Sex 2.843 2.843 1 96.76 0.7380 4
## IPVstatus:Sex 88.424 88.424 1 44.37 0.5778 5
## Age:Sex 1.252 1.252 1 77.70 1.1182 6
## Sex:PovStat 60.154 60.154 1 45.56 1.0679 7
## Sex 26.385 26.385 1 56.47 0.0143 8
## Age:IPVstatus:PovStat 206.510 206.510 1 101.76 2.9369 9
## IPVstatus:PovStat 19.548 19.548 1 57.87 0.3642 10
## Age:PovStat 182.417 182.417 1 100.63 3.3887 11
## Age:IPVstatus 262.071 262.071 1 82.95 3.8159 12
## IPVstatus 52.281 52.281 1 39.69 0.6913 13
## Age 1357.222 1357.222 1 64.17 19.6560 kept
## PovStat 635.262 635.262 1 61.17 8.8820 kept
## Pr(>F)
## Age:IPVstatus:Sex:PovStat 0.1319
## IPVstatus:Sex:PovStat 0.7795
## Age:Sex:PovStat 0.6718
## Age:IPVstatus:Sex 0.3924
## IPVstatus:Sex 0.4512
## Age:Sex 0.2936
## Sex:PovStat 0.3069
## Sex 0.9054
## Age:IPVstatus:PovStat 0.0896
## IPVstatus:PovStat 0.5485
## Age:PovStat 0.0686
## Age:IPVstatus 0.0541
## IPVstatus 0.4107
## Age 0
## PovStat 0.0041
##
## Least squares means:
## PovStat Estimate Standard Error DF t-value Lower CI
## PovStat Above 1.0 30.25 1.57 29.9 19.25 27.0
## PovStat Below 2.0 38.21 2.25 52.3 16.97 33.7
## Upper CI p-value
## PovStat Above 33.5 <2e-16
## PovStat Below 42.7 <2e-16
##
## Differences of LSMEANS:
## Estimate Standard Error DF t-value Lower CI Upper CI
## PovStat Above-Below -8.0 2.67 61.2 -2.98 -13.3 -2.62
## p-value
## PovStat Above-Below 0.004
##
## Final model:
## lme4::lmer(formula = TrailsAtestSec ~ Age + PovStat + (1 | subclass) +
## (1 | HNDid), data = IPVandCognitionDataSet2, REML = reml,
## contrasts = l)
Re-run final Model 1
(mm1 = lmer(TrailsAtestSec ~ Age + PovStat + (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: TrailsAtestSec ~ Age + PovStat + (Age | HNDid) + (1 | subclass)
## Data: IPVandCognitionDataSet2
## AIC BIC logLik deviance
## 950.4 973.1 -467.2 934.4
## Random effects:
## Groups Name Std.Dev. Corr
## HNDid (Intercept) 9.433
## Age 0.338 1.00
## subclass (Intercept) 3.356
## Residual 7.414
## Number of obs: 126, groups: HNDid, 63; subclass, 21
## Fixed Effects:
## (Intercept) Age PovStatBelow
## 35.431 0.678 6.717
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 maximum likelihood ['merModLmerTest']
## Formula: TrailsAtestSec ~ Age + PovStat + (Age | HNDid) + (1 | subclass)
## Data: IPVandCognitionDataSet2
##
## AIC BIC logLik deviance
## 950.4 973.1 -467.2 934.4
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 88.989 9.433
## Age 0.114 0.338 1.00
## subclass (Intercept) 11.265 3.356
## Residual 54.970 7.414
## Number of obs: 126, groups: HNDid, 63; subclass, 21
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 35.431 2.306 40.100 15.36 < 2e-16
## Age 0.678 0.174 60.800 3.90 0.00024
## PovStatBelow 6.717 2.547 64.300 2.64 0.01048
##
## Correlation of Fixed Effects:
## (Intr) Age
## Age 0.770
## PovStatBelw -0.443 -0.183
plot(st)
plot(mm1)
Trails 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(TrailsAtestSec ~ (Age + IPVstatus + PovStat + Sex + 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: TrailsAtestSec ~ (Age + IPVstatus + PovStat + Sex + CES1)^5 + (Age | HNDid) + (1 | subclass)
## Data: IPVandCognitionDataSet2
## REML criterion at convergence: 790.7
## Random effects:
## Groups Name Std.Dev. Corr
## HNDid (Intercept) 8.064
## Age 0.395 1.00
## subclass (Intercept) 7.397
## Residual 6.980
## Number of obs: 126, groups: HNDid, 63; subclass, 21
## Fixed Effects:
## (Intercept)
## 30.365
## Age
## 0.487
## IPVstatus1
## 14.351
## PovStatBelow
## 12.095
## SexMen
## 2.020
## CES11
## 13.566
## Age:IPVstatus1
## 1.446
## Age:PovStatBelow
## -0.054
## Age:SexMen
## -0.202
## Age:CES11
## 0.542
## IPVstatus1:PovStatBelow
## 44.256
## IPVstatus1:SexMen
## -11.188
## IPVstatus1:CES11
## -31.685
## PovStatBelow:SexMen
## -2.804
## PovStatBelow:CES11
## -21.172
## SexMen:CES11
## -24.971
## Age:IPVstatus1:PovStatBelow
## 4.121
## Age:IPVstatus1:SexMen
## -1.731
## Age:IPVstatus1:CES11
## -3.004
## Age:PovStatBelow:SexMen
## 0.177
## Age:PovStatBelow:CES11
## -0.636
## Age:SexMen:CES11
## -1.076
## IPVstatus1:PovStatBelow:SexMen
## -40.719
## IPVstatus1:PovStatBelow:CES11
## -5.895
## IPVstatus1:SexMen:CES11
## 44.076
## PovStatBelow:SexMen:CES11
## 27.872
## Age:IPVstatus1:PovStatBelow:SexMen
## -3.250
## Age:IPVstatus1:PovStatBelow:CES11
## 0.257
## Age:IPVstatus1:SexMen:CES11
## 5.012
## Age:PovStatBelow:SexMen:CES11
## 2.264
## IPVstatus1:PovStatBelow:SexMen:CES11
## 12.776
## Age:IPVstatus1:PovStatBelow:SexMen:CES11
## -1.919
(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 effects:
## Chi.sq Chi.DF elim.num p.value
## (1 | HNDid) 0.00 1 1 1
## (1 | subclass) 21.12 1 kept 0
##
## Fixed effects:
## Sum Sq Mean Sq NumDF DenDF F.value
## Age:IPVstatus:PovStat:Sex:CES1 42.8865 42.8865 1 72.59 0.6859
## Age:IPVstatus:PovStat:CES1 174.3430 174.3430 1 74.57 0.0016
## Age:PovStat:Sex:CES1 13.7889 13.7889 1 79.65 0.7095
## Age:PovStat:CES1 37.4892 37.4892 1 81.03 0.4787
## IPVstatus:PovStat:Sex:CES1 62.4128 62.4128 1 84.69 0.9357
## IPVstatus:PovStat:CES1 133.2820 133.2820 1 90.14 0.0034
## Age 2103.1625 2103.1625 1 96.80 27.5107
## IPVstatus 126.8242 126.8242 1 86.48 16.4815
## PovStat 1759.5480 1759.5480 1 82.47 33.9089
## Sex 92.2956 92.2956 1 90.37 4.8045
## CES1 2.9722 2.9722 1 89.00 2.2859
## Age:IPVstatus 658.2019 658.2019 1 86.56 13.7845
## Age:PovStat 379.5660 379.5660 1 88.18 11.1758
## Age:Sex 14.7744 14.7744 1 81.87 3.7534
## Age:CES1 0.9615 0.9615 1 88.48 0.0750
## IPVstatus:PovStat 118.8368 118.8368 1 97.07 11.9257
## IPVstatus:Sex 264.5152 264.5152 1 84.80 0.5870
## IPVstatus:CES1 21.4711 21.4711 1 95.97 1.2055
## PovStat:Sex 130.9547 130.9547 1 86.86 0.2393
## PovStat:CES1 191.5774 191.5774 1 97.52 2.3473
## Sex:CES1 249.4829 249.4829 1 97.03 3.3226
## Age:IPVstatus:PovStat 653.1834 653.1834 1 97.84 10.0966
## Age:IPVstatus:Sex 0.0969 0.0969 1 84.29 1.8685
## Age:IPVstatus:CES1 82.1379 82.1379 1 92.73 0.5804
## Age:PovStat:Sex 92.9537 92.9537 1 90.51 0.7839
## Age:Sex:CES1 519.6226 519.6226 1 94.44 5.1850
## IPVstatus:PovStat:Sex 0.5051 0.5051 1 90.42 4.2837
## IPVstatus:Sex:CES1 258.2273 258.2273 1 92.67 11.9337
## PovStat:Sex:CES1 126.6146 126.6146 1 95.21 3.9738
## Age:IPVstatus:PovStat:Sex 364.9849 364.9849 1 86.40 6.4192
## Age:IPVstatus:Sex:CES1 345.1971 345.1971 1 90.64 9.0169
## elim.num Pr(>F)
## Age:IPVstatus:PovStat:Sex:CES1 1 0.4103
## Age:IPVstatus:PovStat:CES1 2 0.9679
## Age:PovStat:Sex:CES1 3 0.4021
## Age:PovStat:CES1 4 0.4910
## IPVstatus:PovStat:Sex:CES1 5 0.3361
## IPVstatus:PovStat:CES1 6 0.9535
## Age kept 0
## IPVstatus kept 0.0001
## PovStat kept 0
## Sex kept 0.0310
## CES1 kept 0.1341
## Age:IPVstatus kept 0.0004
## Age:PovStat kept 0.0012
## Age:Sex kept 0.0561
## Age:CES1 kept 0.7849
## IPVstatus:PovStat kept 0.0008
## IPVstatus:Sex kept 0.4457
## IPVstatus:CES1 kept 0.2750
## PovStat:Sex kept 0.6260
## PovStat:CES1 kept 0.1287
## Sex:CES1 kept 0.0714
## Age:IPVstatus:PovStat kept 0.0020
## Age:IPVstatus:Sex kept 0.1753
## Age:IPVstatus:CES1 kept 0.4481
## Age:PovStat:Sex kept 0.3783
## Age:Sex:CES1 kept 0.0250
## IPVstatus:PovStat:Sex kept 0.0413
## IPVstatus:Sex:CES1 kept 0.0008
## PovStat:Sex:CES1 kept 0.0491
## Age:IPVstatus:PovStat:Sex kept 0.0131
## Age:IPVstatus:Sex:CES1 kept 0.0035
##
## Least squares means:
## IPVstatus PovStat Sex CES1 Estimate
## IPVstatus 0 1.0 NA NA NA 31.9
## IPVstatus 1 2.0 NA NA NA 37.0
## PovStat Above NA 1.0 NA NA 29.1
## PovStat Below NA 2.0 NA NA 39.8
## Sex Women NA NA 2.0 NA 36.0
## Sex Men NA NA 1.0 NA 33.0
## CES1 0 NA NA NA 1.0 36.8
## CES1 1 NA NA NA 2.0 32.1
## IPVstatus:PovStat 0 Above 1.0 1.0 NA NA 29.1
## IPVstatus:PovStat 1 Above 2.0 1.0 NA NA 29.1
## IPVstatus:PovStat 0 Below 1.0 2.0 NA NA 34.7
## IPVstatus:PovStat 1 Below 2.0 2.0 NA NA 44.9
## IPVstatus:Sex 0 Women 1.0 NA 2.0 NA 34.0
## IPVstatus:Sex 1 Women 2.0 NA 2.0 NA 37.9
## IPVstatus:Sex 0 Men 1.0 NA 1.0 NA 29.8
## IPVstatus:Sex 1 Men 2.0 NA 1.0 NA 36.1
## IPVstatus:CES1 0 0 1.0 NA NA 1.0 33.3
## IPVstatus:CES1 1 0 2.0 NA NA 1.0 40.4
## IPVstatus:CES1 0 1 1.0 NA NA 2.0 30.5
## IPVstatus:CES1 1 1 2.0 NA NA 2.0 33.7
## PovStat:Sex Above Women NA 1.0 2.0 NA 31.0
## PovStat:Sex Below Women NA 2.0 2.0 NA 40.9
## PovStat:Sex Above Men NA 1.0 1.0 NA 27.2
## PovStat:Sex Below Men NA 2.0 1.0 NA 38.7
## PovStat:CES1 Above 0 NA 1.0 NA 1.0 29.3
## PovStat:CES1 Below 0 NA 2.0 NA 1.0 44.4
## PovStat:CES1 Above 1 NA 1.0 NA 2.0 28.9
## PovStat:CES1 Below 1 NA 2.0 NA 2.0 35.3
## Sex:CES1 Women 0 NA NA 2.0 1.0 38.4
## Sex:CES1 Men 0 NA NA 1.0 1.0 35.3
## Sex:CES1 Women 1 NA NA 2.0 2.0 33.5
## Sex:CES1 Men 1 NA NA 1.0 2.0 30.6
## IPVstatus:PovStat:Sex 0 Above Women 1.0 1.0 2.0 NA 32.2
## IPVstatus:PovStat:Sex 1 Above Women 2.0 1.0 2.0 NA 29.8
## IPVstatus:PovStat:Sex 0 Below Women 1.0 2.0 2.0 NA 35.8
## IPVstatus:PovStat:Sex 1 Below Women 2.0 2.0 2.0 NA 46.1
## IPVstatus:PovStat:Sex 0 Above Men 1.0 1.0 1.0 NA 25.9
## IPVstatus:PovStat:Sex 1 Above Men 2.0 1.0 1.0 NA 28.5
## IPVstatus:PovStat:Sex 0 Below Men 1.0 2.0 1.0 NA 33.6
## IPVstatus:PovStat:Sex 1 Below Men 2.0 2.0 1.0 NA 43.7
## IPVstatus:Sex:CES1 0 Women 0 1.0 NA 2.0 1.0 33.4
## IPVstatus:Sex:CES1 1 Women 0 2.0 NA 2.0 1.0 43.4
## IPVstatus:Sex:CES1 0 Men 0 1.0 NA 1.0 1.0 33.2
## IPVstatus:Sex:CES1 1 Men 0 2.0 NA 1.0 1.0 37.4
## IPVstatus:Sex:CES1 0 Women 1 1.0 NA 2.0 2.0 34.6
## IPVstatus:Sex:CES1 1 Women 1 2.0 NA 2.0 2.0 32.5
## IPVstatus:Sex:CES1 0 Men 1 1.0 NA 1.0 2.0 26.4
## IPVstatus:Sex:CES1 1 Men 1 2.0 NA 1.0 2.0 34.9
## PovStat:Sex:CES1 Above Women 0 NA 1.0 2.0 1.0 28.6
## PovStat:Sex:CES1 Below Women 0 NA 2.0 2.0 1.0 48.2
## PovStat:Sex:CES1 Above Men 0 NA 1.0 1.0 1.0 30.1
## PovStat:Sex:CES1 Below Men 0 NA 2.0 1.0 1.0 40.5
## PovStat:Sex:CES1 Above Women 1 NA 1.0 2.0 2.0 33.4
## PovStat:Sex:CES1 Below Women 1 NA 2.0 2.0 2.0 33.7
## PovStat:Sex:CES1 Above Men 1 NA 1.0 1.0 2.0 24.3
## PovStat:Sex:CES1 Below Men 1 NA 2.0 1.0 2.0 36.9
## Standard Error DF t-value Lower CI
## IPVstatus 0 2.28 21.2 14.01 27.2
## IPVstatus 1 2.74 40.5 13.52 31.5
## PovStat Above 2.26 21.0 12.86 24.4
## PovStat Below 2.76 39.2 14.40 34.2
## Sex Women 2.35 23.4 15.30 31.1
## Sex Men 2.88 43.9 11.42 27.1
## CES1 0 2.56 31.3 14.40 31.6
## CES1 1 2.45 27.8 13.06 27.0
## IPVstatus:PovStat 0 Above 2.29 21.8 12.70 24.3
## IPVstatus:PovStat 1 Above 3.06 56.8 9.51 23.0
## IPVstatus:PovStat 0 Below 3.04 54.1 11.44 28.6
## IPVstatus:PovStat 1 Below 3.69 73.4 12.18 37.6
## IPVstatus:Sex 0 Women 2.40 25.1 14.15 29.0
## IPVstatus:Sex 1 Women 3.05 53.6 12.43 31.8
## IPVstatus:Sex 0 Men 3.07 50.2 9.69 23.6
## IPVstatus:Sex 1 Men 3.89 83.4 9.27 28.4
## IPVstatus:CES1 0 0 2.47 27.6 13.51 28.3
## IPVstatus:CES1 1 0 3.74 77.4 10.78 32.9
## IPVstatus:CES1 0 1 2.72 39.4 11.22 25.0
## IPVstatus:CES1 1 1 2.99 50.4 11.27 27.7
## PovStat:Sex Above Women 2.35 23.4 13.20 26.1
## PovStat:Sex Below Women 3.06 53.1 13.37 34.8
## PovStat:Sex Above Men 3.01 50.0 9.04 21.2
## PovStat:Sex Below Men 3.83 78.7 10.10 31.1
## PovStat:CES1 Above 0 3.09 53.1 9.49 23.1
## PovStat:CES1 Below 0 3.14 53.3 14.11 38.0
## PovStat:CES1 Above 1 2.52 28.9 11.45 23.7
## PovStat:CES1 Below 1 3.60 72.5 9.81 28.1
## Sex:CES1 Women 0 2.93 47.6 13.09 32.5
## Sex:CES1 Men 0 3.36 65.2 10.50 28.6
## Sex:CES1 Women 1 2.77 35.6 12.11 27.9
## Sex:CES1 Men 1 3.34 64.2 9.18 23.9
## IPVstatus:PovStat:Sex 0 Above Women 2.59 31.8 12.44 26.9
## IPVstatus:PovStat:Sex 1 Above Women 2.87 47.3 10.39 24.0
## IPVstatus:PovStat:Sex 0 Below Women 3.38 58.9 10.58 29.0
## IPVstatus:PovStat:Sex 1 Below Women 4.68 94.1 9.84 36.8
## IPVstatus:PovStat:Sex 0 Above Men 2.92 45.3 8.89 20.1
## IPVstatus:PovStat:Sex 1 Above Men 4.84 98.3 5.89 18.9
## IPVstatus:PovStat:Sex 0 Below Men 4.65 96.0 7.23 24.4
## IPVstatus:PovStat:Sex 1 Below Men 5.22 95.1 8.38 33.4
## IPVstatus:Sex:CES1 0 Women 0 2.71 39.9 12.33 27.9
## IPVstatus:Sex:CES1 1 Women 0 4.35 89.7 9.97 34.7
## IPVstatus:Sex:CES1 0 Men 0 3.24 57.4 10.24 26.7
## IPVstatus:Sex:CES1 1 Men 0 5.64 99.5 6.62 26.2
## IPVstatus:Sex:CES1 0 Women 1 3.09 49.3 11.18 28.4
## IPVstatus:Sex:CES1 1 Women 1 3.33 61.3 9.74 25.8
## IPVstatus:Sex:CES1 0 Men 1 3.99 82.8 6.61 18.4
## IPVstatus:Sex:CES1 1 Men 1 4.40 90.7 7.92 26.1
## PovStat:Sex:CES1 Above Women 0 2.79 41.6 10.25 23.0
## PovStat:Sex:CES1 Below Women 0 4.18 88.9 11.53 39.9
## PovStat:Sex:CES1 Above Men 0 4.75 95.9 6.33 20.7
## PovStat:Sex:CES1 Below Men 0 3.90 78.5 10.38 32.7
## PovStat:Sex:CES1 Above Women 1 2.80 39.1 11.93 27.7
## PovStat:Sex:CES1 Below Women 1 3.63 70.8 9.27 26.4
## PovStat:Sex:CES1 Above Men 1 3.57 59.6 6.82 17.2
## PovStat:Sex:CES1 Below Men 1 5.71 99.3 6.46 25.5
## Upper CI p-value
## IPVstatus 0 36.6 <2e-16
## IPVstatus 1 42.5 <2e-16
## PovStat Above 33.8 <2e-16
## PovStat Below 45.4 <2e-16
## Sex Women 40.8 <2e-16
## Sex Men 38.8 <2e-16
## CES1 0 42.1 <2e-16
## CES1 1 37.1 <2e-16
## IPVstatus:PovStat 0 Above 33.8 <2e-16
## IPVstatus:PovStat 1 Above 35.3 <2e-16
## IPVstatus:PovStat 0 Below 40.8 <2e-16
## IPVstatus:PovStat 1 Below 52.2 <2e-16
## IPVstatus:Sex 0 Women 38.9 <2e-16
## IPVstatus:Sex 1 Women 44.1 <2e-16
## IPVstatus:Sex 0 Men 36.0 <2e-16
## IPVstatus:Sex 1 Men 43.8 <2e-16
## IPVstatus:CES1 0 0 38.4 <2e-16
## IPVstatus:CES1 1 0 47.8 <2e-16
## IPVstatus:CES1 0 1 36.0 <2e-16
## IPVstatus:CES1 1 1 39.7 <2e-16
## PovStat:Sex Above Women 35.8 <2e-16
## PovStat:Sex Below Women 47.1 <2e-16
## PovStat:Sex Above Men 33.3 <2e-16
## PovStat:Sex Below Men 46.3 <2e-16
## PovStat:CES1 Above 0 35.5 <2e-16
## PovStat:CES1 Below 0 50.7 <2e-16
## PovStat:CES1 Above 1 34.0 <2e-16
## PovStat:CES1 Below 1 42.4 <2e-16
## Sex:CES1 Women 0 44.3 <2e-16
## Sex:CES1 Men 0 42.0 <2e-16
## Sex:CES1 Women 1 39.1 <2e-16
## Sex:CES1 Men 1 37.3 <2e-16
## IPVstatus:PovStat:Sex 0 Above Women 37.4 <2e-16
## IPVstatus:PovStat:Sex 1 Above Women 35.6 <2e-16
## IPVstatus:PovStat:Sex 0 Below Women 42.6 <2e-16
## IPVstatus:PovStat:Sex 1 Below Women 55.4 <2e-16
## IPVstatus:PovStat:Sex 0 Above Men 31.8 <2e-16
## IPVstatus:PovStat:Sex 1 Above Men 38.1 <2e-16
## IPVstatus:PovStat:Sex 0 Below Men 42.9 <2e-16
## IPVstatus:PovStat:Sex 1 Below Men 54.1 <2e-16
## IPVstatus:Sex:CES1 0 Women 0 38.9 <2e-16
## IPVstatus:Sex:CES1 1 Women 0 52.0 <2e-16
## IPVstatus:Sex:CES1 0 Men 0 39.7 <2e-16
## IPVstatus:Sex:CES1 1 Men 0 48.5 <2e-16
## IPVstatus:Sex:CES1 0 Women 1 40.8 <2e-16
## IPVstatus:Sex:CES1 1 Women 1 39.1 <2e-16
## IPVstatus:Sex:CES1 0 Men 1 34.3 <2e-16
## IPVstatus:Sex:CES1 1 Men 1 43.6 <2e-16
## PovStat:Sex:CES1 Above Women 0 34.2 <2e-16
## PovStat:Sex:CES1 Below Women 0 56.5 <2e-16
## PovStat:Sex:CES1 Above Men 0 39.5 <2e-16
## PovStat:Sex:CES1 Below Men 0 48.3 <2e-16
## PovStat:Sex:CES1 Above Women 1 39.0 <2e-16
## PovStat:Sex:CES1 Below Women 1 40.9 <2e-16
## PovStat:Sex:CES1 Above Men 1 31.5 <2e-16
## PovStat:Sex:CES1 Below Men 1 48.2 <2e-16
##
## Differences of LSMEANS:
## Estimate Standard Error DF
## IPVstatus 0-1 -5.1 2.3586 86.7
## PovStat Above-Below -10.7 2.3906 94.1
## Sex Women-Men 3.0 2.8084 96.3
## CES1 0-1 4.8 2.3122 92.7
## IPVstatus:PovStat 0 Above- 1 Above -0.1 2.9624 86.1
## IPVstatus:PovStat 0 Above- 0 Below -5.7 2.8579 87.1
## IPVstatus:PovStat 0 Above- 1 Below -15.8 3.5316 92.6
## IPVstatus:PovStat 1 Above- 0 Below -5.6 3.1755 88.1
## IPVstatus:PovStat 1 Above- 1 Below -15.8 3.9922 95.0
## IPVstatus:PovStat 0 Below- 1 Below -10.2 3.8762 90.9
## IPVstatus:Sex 0 Women- 1 Women -3.9 2.8424 89.8
## IPVstatus:Sex 0 Women- 0 Men 4.2 3.1166 98.5
## IPVstatus:Sex 0 Women- 1 Men -2.1 3.7006 89.3
## IPVstatus:Sex 1 Women- 0 Men 8.1 3.6341 96.0
## IPVstatus:Sex 1 Women- 1 Men 1.8 4.3522 92.2
## IPVstatus:Sex 0 Men- 1 Men -6.3 3.9908 90.2
## IPVstatus:CES1 0 0- 1 0 -7.1 3.7437 92.6
## IPVstatus:CES1 0 0- 0 1 2.8 2.4860 90.9
## IPVstatus:CES1 0 0- 1 1 -0.3 2.9332 90.2
## IPVstatus:CES1 1 0- 0 1 9.9 3.6353 89.5
## IPVstatus:CES1 1 0- 1 1 6.7 3.9863 95.4
## IPVstatus:CES1 0 1- 1 1 -3.2 2.9134 90.0
## PovStat:Sex Above Women- Below Women -10.0 2.7679 90.7
## PovStat:Sex Above Women- Above Men 3.8 2.9469 96.3
## PovStat:Sex Above Women- Below Men -7.7 3.7790 95.5
## PovStat:Sex Below Women- Above Men 13.7 3.5951 95.4
## PovStat:Sex Below Women- Below Men 2.3 4.1857 91.9
## PovStat:Sex Above Men- Below Men -11.5 3.7681 91.9
## PovStat:CES1 Above 0- Below 0 -15.0 3.5602 96.6
## PovStat:CES1 Above 0- Above 1 0.5 3.3694 97.4
## PovStat:CES1 Above 0- Below 1 -5.9 3.5987 94.1
## PovStat:CES1 Below 0- Above 1 15.5 3.0286 92.4
## PovStat:CES1 Below 0- Below 1 9.1 3.8816 94.4
## PovStat:CES1 Above 1- Below 1 -6.4 3.8056 95.9
## Sex:CES1 Women 0- Men 0 3.1 3.6859 91.0
## Sex:CES1 Women 0- Women 1 4.9 3.2297 98.5
## Sex:CES1 Women 0- Men 1 7.8 3.6260 90.6
## Sex:CES1 Men 0- Women 1 1.8 3.6497 98.1
## Sex:CES1 Men 0- Men 1 4.7 3.3982 86.9
## Sex:CES1 Women 1- Men 1 2.9 3.6708 98.6
## t-value Lower CI Upper CI p-value
## IPVstatus 0-1 -2.17 -9.813 -0.4365 0.033
## PovStat Above-Below -4.48 -15.462 -5.9686 <2e-16
## Sex Women-Men 1.07 -2.564 8.5847 0.286
## CES1 0-1 2.07 0.184 9.3675 0.042
## IPVstatus:PovStat 0 Above- 1 Above -0.03 -5.968 5.8104 0.979
## IPVstatus:PovStat 0 Above- 0 Below -1.98 -11.349 0.0112 0.050
## IPVstatus:PovStat 0 Above- 1 Below -4.49 -22.853 -8.8265 <2e-16
## IPVstatus:PovStat 1 Above- 0 Below -1.76 -11.901 0.7200 0.082
## IPVstatus:PovStat 1 Above- 1 Below -3.95 -23.687 -7.8358 2e-04
## IPVstatus:PovStat 0 Below- 1 Below -2.62 -17.870 -2.4710 0.010
## IPVstatus:Sex 0 Women- 1 Women -1.39 -9.594 1.6999 0.168
## IPVstatus:Sex 0 Women- 0 Men 1.34 -1.996 10.3721 0.182
## IPVstatus:Sex 0 Women- 1 Men -0.57 -9.467 5.2381 0.569
## IPVstatus:Sex 1 Women- 0 Men 2.24 0.921 15.3484 0.028
## IPVstatus:Sex 1 Women- 1 Men 0.42 -6.811 10.4762 0.675
## IPVstatus:Sex 0 Men- 1 Men -1.58 -14.230 1.6260 0.118
## IPVstatus:CES1 0 0- 1 0 -1.88 -14.487 0.3821 0.063
## IPVstatus:CES1 0 0- 0 1 1.15 -2.091 7.7857 0.255
## IPVstatus:CES1 0 0- 1 1 -0.12 -6.176 5.4780 0.905
## IPVstatus:CES1 1 0- 0 1 2.72 2.678 17.1228 0.008
## IPVstatus:CES1 1 0- 1 1 1.68 -1.210 14.6170 0.096
## IPVstatus:CES1 0 1- 1 1 -1.10 -8.985 2.5914 0.276
## PovStat:Sex Above Women- Below Women -3.60 -15.462 -4.4648 5e-04
## PovStat:Sex Above Women- Above Men 1.28 -2.087 9.6115 0.205
## PovStat:Sex Above Women- Below Men -2.04 -15.207 -0.2033 0.044
## PovStat:Sex Below Women- Above Men 3.82 6.589 20.8622 2e-04
## PovStat:Sex Below Women- Below Men 0.54 -6.055 10.5715 0.591
## PovStat:Sex Above Men- Below Men -3.04 -18.951 -3.9833 0.003
## PovStat:CES1 Above 0- Below 0 -4.22 -22.078 -7.9451 1e-04
## PovStat:CES1 Above 0- Above 1 0.14 -6.208 7.1662 0.887
## PovStat:CES1 Above 0- Below 1 -1.65 -13.085 1.2056 0.102
## PovStat:CES1 Below 0- Above 1 5.11 9.476 21.5055 <2e-16
## PovStat:CES1 Below 0- Below 1 2.34 1.365 16.7785 0.021
## PovStat:CES1 Above 1- Below 1 -1.69 -13.973 1.1352 0.095
## Sex:CES1 Women 0- Men 0 0.84 -4.217 10.4261 0.402
## Sex:CES1 Women 0- Women 1 1.51 -1.539 11.2788 0.135
## Sex:CES1 Women 0- Men 1 2.15 0.583 14.9887 0.034
## Sex:CES1 Men 0- Women 1 0.48 -5.477 9.0079 0.630
## Sex:CES1 Men 0- Men 1 1.38 -2.073 11.4356 0.172
## Sex:CES1 Women 1- Men 1 0.79 -4.368 10.1998 0.429
##
## Final model:
## lme4::lmer(formula = TrailsAtestSec ~ Age + IPVstatus + PovStat +
## Sex + CES1 + (1 | subclass) + Age:IPVstatus + Age:PovStat +
## Age:Sex + Age:CES1 + IPVstatus:PovStat + IPVstatus:Sex +
## IPVstatus:CES1 + PovStat:Sex + PovStat:CES1 + Sex:CES1 +
## Age:IPVstatus:PovStat + Age:IPVstatus:Sex + Age:IPVstatus:CES1 +
## Age:PovStat:Sex + Age:Sex:CES1 + IPVstatus:PovStat:Sex +
## IPVstatus:Sex:CES1 + PovStat:Sex:CES1 + Age:IPVstatus:PovStat:Sex +
## Age:IPVstatus:Sex:CES1, data = IPVandCognitionDataSet2, REML = reml,
## contrasts = l)
Re-run suggested final Model 2
(mm2 = lmer(TrailsAtestSec ~ Age + IPVstatus + PovStat + Sex + CES1 + (Age |
HNDid) + (1 | subclass) + Age:IPVstatus + Age:PovStat + Age:Sex + Age:CES1 +
IPVstatus:PovStat + IPVstatus:Sex + IPVstatus:CES1 + PovStat:Sex + PovStat:CES1 +
Sex:CES1 + Age:IPVstatus:PovStat + Age:IPVstatus:Sex + Age:IPVstatus:CES1 +
Age:PovStat:Sex + Age:Sex:CES1 + IPVstatus:PovStat:Sex + IPVstatus:Sex:CES1 +
PovStat:Sex:CES1 + Age:IPVstatus:PovStat:Sex + Age:IPVstatus:Sex:CES1, data = IPVandCognitionDataSet2))
## Warning: number of observations <= rank(Z); variance-covariance matrix
## will be unidentifiable
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: TrailsAtestSec ~ Age + IPVstatus + PovStat + Sex + CES1 + (Age | HNDid) + (1 | subclass) + Age:IPVstatus + Age:PovStat + Age:Sex + Age:CES1 + IPVstatus:PovStat + IPVstatus:Sex + IPVstatus:CES1 + PovStat:Sex + PovStat:CES1 + Sex:CES1 + Age:IPVstatus:PovStat + Age:IPVstatus:Sex + Age:IPVstatus:CES1 + Age:PovStat:Sex + Age:Sex:CES1 + IPVstatus:PovStat:Sex + IPVstatus:Sex:CES1 + PovStat:Sex:CES1 + Age:IPVstatus:PovStat:Sex + Age:IPVstatus:Sex:CES1
## Data: IPVandCognitionDataSet2
## REML criterion at convergence: 820.9
## Random effects:
## Groups Name Std.Dev. Corr
## HNDid (Intercept) 7.545
## Age 0.364 1.00
## subclass (Intercept) 7.557
## Residual 6.825
## Number of obs: 126, groups: HNDid, 63; subclass, 21
## Fixed Effects:
## (Intercept) Age
## 30.723 0.586
## IPVstatus1 PovStatBelow
## 15.276 9.076
## SexMen CES11
## 1.239 12.807
## Age:IPVstatus1 Age:PovStatBelow
## 1.431 -0.483
## Age:SexMen Age:CES11
## -0.352 0.325
## IPVstatus1:PovStatBelow IPVstatus1:SexMen
## 39.325 -22.233
## IPVstatus1:CES11 PovStatBelow:SexMen
## -32.968 0.560
## PovStatBelow:CES11 SexMen:CES11
## -16.996 -21.380
## Age:IPVstatus1:PovStatBelow Age:IPVstatus1:SexMen
## 4.099 -2.155
## Age:IPVstatus1:CES11 Age:PovStatBelow:SexMen
## -2.901 0.871
## Age:SexMen:CES11 IPVstatus1:PovStatBelow:SexMen
## -0.454 -24.420
## IPVstatus1:SexMen:CES11 PovStatBelow:SexMen:CES11
## 52.708 18.418
## Age:IPVstatus1:PovStatBelow:SexMen Age:IPVstatus1:SexMen:CES11
## -2.997 4.912
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: TrailsAtestSec ~ Age + IPVstatus + PovStat + Sex + CES1 + (Age | HNDid) + (1 | subclass) + Age:IPVstatus + Age:PovStat + Age:Sex + Age:CES1 + IPVstatus:PovStat + IPVstatus:Sex + IPVstatus:CES1 + PovStat:Sex + PovStat:CES1 + Sex:CES1 + Age:IPVstatus:PovStat + Age:IPVstatus:Sex + Age:IPVstatus:CES1 + Age:PovStat:Sex + Age:Sex:CES1 + IPVstatus:PovStat:Sex + IPVstatus:Sex:CES1 + PovStat:Sex:CES1 + Age:IPVstatus:PovStat:Sex + Age:IPVstatus:Sex:CES1
## Data: IPVandCognitionDataSet2
##
## REML criterion at convergence: 820.9
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 56.924 7.545
## Age 0.132 0.364 1.00
## subclass (Intercept) 57.113 7.557
## Residual 46.576 6.825
## Number of obs: 126, groups: HNDid, 63; subclass, 21
##
## Fixed effects:
## Estimate Std. Error df t value
## (Intercept) 30.723 5.077 25.300 6.05
## Age 0.586 0.361 38.900 1.62
## IPVstatus1 15.276 9.216 45.600 1.66
## PovStatBelow 9.076 7.454 31.800 1.22
## SexMen 1.239 6.706 19.300 0.18
## CES11 12.807 6.478 27.800 1.98
## Age:IPVstatus1 1.431 0.801 78.400 1.79
## Age:PovStatBelow -0.483 0.550 33.500 -0.88
## Age:SexMen -0.352 0.515 29.900 -0.68
## Age:CES11 0.325 0.518 63.800 0.63
## IPVstatus1:PovStatBelow 39.325 14.214 47.000 2.77
## IPVstatus1:SexMen -22.233 18.424 34.500 -1.21
## IPVstatus1:CES11 -32.968 11.556 41.800 -2.85
## PovStatBelow:SexMen 0.560 9.896 24.300 0.06
## PovStatBelow:CES11 -16.996 5.358 25.900 -3.17
## SexMen:CES11 -21.380 11.678 50.300 -1.83
## Age:IPVstatus1:PovStatBelow 4.099 1.285 73.300 3.19
## Age:IPVstatus1:SexMen -2.155 1.302 55.600 -1.66
## Age:IPVstatus1:CES11 -2.901 1.011 71.300 -2.87
## Age:PovStatBelow:SexMen 0.871 0.907 32.700 0.96
## Age:SexMen:CES11 -0.454 1.042 80.400 -0.44
## IPVstatus1:PovStatBelow:SexMen -24.420 19.490 37.700 -1.25
## IPVstatus1:SexMen:CES11 52.708 17.911 25.600 2.94
## PovStatBelow:SexMen:CES11 18.418 11.846 57.600 1.55
## Age:IPVstatus1:PovStatBelow:SexMen -2.997 1.731 59.500 -1.73
## Age:IPVstatus1:SexMen:CES11 4.912 1.558 61.800 3.15
## Pr(>|t|)
## (Intercept) 2.4e-06
## Age 0.1132
## IPVstatus1 0.1043
## PovStatBelow 0.2323
## SexMen 0.8553
## CES11 0.0580
## Age:IPVstatus1 0.0777
## Age:PovStatBelow 0.3859
## Age:SexMen 0.4999
## Age:CES11 0.5329
## IPVstatus1:PovStatBelow 0.0081
## IPVstatus1:SexMen 0.2357
## IPVstatus1:CES11 0.0067
## PovStatBelow:SexMen 0.9553
## PovStatBelow:CES11 0.0039
## SexMen:CES11 0.0731
## Age:IPVstatus1:PovStatBelow 0.0021
## Age:IPVstatus1:SexMen 0.1035
## Age:IPVstatus1:CES11 0.0054
## Age:PovStatBelow:SexMen 0.3438
## Age:SexMen:CES11 0.6640
## IPVstatus1:PovStatBelow:SexMen 0.2179
## IPVstatus1:SexMen:CES11 0.0068
## PovStatBelow:SexMen:CES11 0.1255
## Age:IPVstatus1:PovStatBelow:SexMen 0.0885
## Age:IPVstatus1:SexMen:CES11 0.0025
##
## Correlation matrix not shown by default, as p = 26 > 20.
## Use print(x, correlation=TRUE) or
## vcov(x) if you need it
plot(st)
plot(mm2)
