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 REML ['merModLmerTest']
## Formula: logTrailsA ~ (Age + IPVstatus + Sex + PovStat)^4 + (Age | HNDid) + (1 | subclass)
## Data: IPVandCognitionDataSet2
## REML criterion at convergence: 105.1
## Random effects:
## Groups Name Std.Dev. Corr
## HNDid (Intercept) 0.24226
## Age 0.00363 1.00
## subclass (Intercept) 0.16528
## Residual 0.20530
## Number of obs: 126, groups: HNDid, 63; subclass, 21
## Fixed Effects:
## (Intercept) Age
## 3.49590 0.01943
## IPVstatus1 SexMen
## 0.01224 -0.12694
## PovStatBelow Age:IPVstatus1
## 0.13647 0.00401
## Age:SexMen Age:PovStatBelow
## -0.00852 -0.00467
## IPVstatus1:SexMen IPVstatus1:PovStatBelow
## 0.03726 0.59477
## SexMen:PovStatBelow Age:IPVstatus1:SexMen
## 0.17288 -0.00661
## Age:IPVstatus1:PovStatBelow Age:SexMen:PovStatBelow
## 0.07783 0.00854
## IPVstatus1:SexMen:PovStatBelow Age:IPVstatus1:SexMen:PovStatBelow
## -0.28503 -0.05589
## 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.39 1 kept 0.0362
## (1 | HNDid) 9.34 1 kept 0.0022
##
## Fixed effects:
## Sum Sq Mean Sq NumDF DenDF F.value elim.num
## Age:IPVstatus:Sex:PovStat 0.0718 0.0718 1 101.05 1.6908 1
## Age:Sex:PovStat 0.0014 0.0014 1 89.96 0.0222 2
## IPVstatus:Sex:PovStat 0.0016 0.0016 1 57.48 0.0486 3
## Age:IPVstatus:Sex 0.0088 0.0088 1 99.84 0.9785 4
## IPVstatus:Sex 0.0561 0.0561 1 44.26 0.5451 5
## Age:Sex 0.0081 0.0081 1 81.40 1.4800 6
## Sex:PovStat 0.0450 0.0450 1 45.43 0.9178 7
## Sex 0.0182 0.0182 1 56.60 0.0813 8
## Age:IPVstatus:PovStat 0.1344 0.1344 1 106.00 2.0650 9
## IPVstatus:PovStat 0.0029 0.0029 1 57.59 0.1146 10
## Age:PovStat 0.0578 0.0578 1 105.57 1.4302 11
## Age:IPVstatus 0.1567 0.1567 1 85.63 3.0904 12
## IPVstatus 0.0432 0.0432 1 39.50 0.8161 13
## Age 0.9870 0.9870 1 72.43 17.7015 kept
## PovStat 0.4068 0.4068 1 60.46 8.1005 kept
## Pr(>F)
## Age:IPVstatus:Sex:PovStat 0.1964
## Age:Sex:PovStat 0.8819
## IPVstatus:Sex:PovStat 0.8264
## Age:IPVstatus:Sex 0.3250
## IPVstatus:Sex 0.4642
## Age:Sex 0.2273
## Sex:PovStat 0.3431
## Sex 0.7766
## Age:IPVstatus:PovStat 0.1537
## IPVstatus:PovStat 0.7362
## Age:PovStat 0.2344
## Age:IPVstatus 0.0823
## IPVstatus 0.3718
## Age 1e-04
## PovStat 0.0060
##
## Least squares means:
## PovStat Estimate Standard Error DF t-value Lower CI
## PovStat Above 1.0 3.3456 0.0482 29.6 69.3800 3.25
## PovStat Below 2.0 3.5723 0.0681 52.5 52.4400 3.44
## Upper CI p-value
## PovStat Above 3.44 <2e-16
## PovStat Below 3.71 <2e-16
##
## Differences of LSMEANS:
## Estimate Standard Error DF t-value Lower CI Upper CI
## PovStat Above-Below -0.2 0.0796 60.5 -2.85 -0.386 -0.0674
## p-value
## PovStat Above-Below 0.006
##
## Final model:
## lme4::lmer(formula = logTrailsA ~ Age + PovStat + (1 | subclass) +
## (1 | HNDid), data = IPVandCognitionDataSet2, REML = reml,
## contrasts = l)
Re-run final Model 1
(mm1 = lmer(logTrailsA ~ Age + PovStat + (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: logTrailsA ~ Age + PovStat + (Age | HNDid) + (1 | subclass)
## Data: IPVandCognitionDataSet2
## REML criterion at convergence: 60.73
## Random effects:
## Groups Name Std.Dev. Corr
## HNDid (Intercept) 0.25622
## Age 0.00428 1.00
## subclass (Intercept) 0.11868
## Residual 0.20765
## Number of obs: 126, groups: HNDid, 63; subclass, 21
## Fixed Effects:
## (Intercept) Age PovStatBelow
## 3.493 0.020 0.217
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: logTrailsA ~ Age + PovStat + (Age | HNDid) + (1 | subclass)
## Data: IPVandCognitionDataSet2
##
## REML criterion at convergence: 60.73
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 6.56e-02 0.25622
## Age 1.83e-05 0.00428 1.00
## subclass (Intercept) 1.41e-02 0.11868
## Residual 4.31e-02 0.20765
## Number of obs: 126, groups: HNDid, 63; subclass, 21
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 3.49339 0.06713 29.40000 52.04 < 2e-16
## Age 0.02004 0.00502 55.50000 4.00 0.00019
## PovStatBelow 0.21736 0.07990 60.90000 2.72 0.00849
##
## Correlation of Fixed Effects:
## (Intr) Age
## Age 0.690
## PovStatBelw -0.479 -0.195
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(logTrailsA ~ (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: logTrailsA ~ (Age + IPVstatus + PovStat + Sex + CES1)^5 + (Age | HNDid) + (1 | subclass)
## Data: IPVandCognitionDataSet2
## REML criterion at convergence: 126.3
## Random effects:
## Groups Name Std.Dev. Corr
## HNDid (Intercept) 0.2484
## Age 0.0108 1.00
## subclass (Intercept) 0.2400
## Residual 0.1948
## Number of obs: 126, groups: HNDid, 63; subclass, 21
## Fixed Effects:
## (Intercept)
## 3.27571
## Age
## 0.01132
## IPVstatus1
## 0.56229
## PovStatBelow
## 0.51101
## SexMen
## 0.18535
## CES11
## 0.53124
## Age:IPVstatus1
## 0.05488
## Age:PovStatBelow
## 0.00817
## Age:SexMen
## 0.00300
## Age:CES11
## 0.02339
## IPVstatus1:PovStatBelow
## 0.79088
## IPVstatus1:SexMen
## -0.48505
## IPVstatus1:CES11
## -1.11796
## PovStatBelow:SexMen
## -0.27201
## PovStatBelow:CES11
## -0.84280
## SexMen:CES11
## -0.88872
## Age:IPVstatus1:PovStatBelow
## 0.08332
## Age:IPVstatus1:SexMen
## -0.06921
## Age:IPVstatus1:CES11
## -0.10798
## Age:PovStatBelow:SexMen
## -0.01737
## Age:PovStatBelow:CES11
## -0.03368
## Age:SexMen:CES11
## -0.04307
## IPVstatus1:PovStatBelow:SexMen
## -0.68797
## IPVstatus1:PovStatBelow:CES11
## 0.36562
## IPVstatus1:SexMen:CES11
## 1.50619
## PovStatBelow:SexMen:CES11
## 0.94029
## Age:IPVstatus1:PovStatBelow:SexMen
## -0.04672
## Age:IPVstatus1:PovStatBelow:CES11
## 0.05031
## Age:IPVstatus1:SexMen:CES11
## 0.17332
## Age:PovStatBelow:SexMen:CES11
## 0.06501
## IPVstatus1:PovStatBelow:SexMen:CES11
## -0.12875
## Age:IPVstatus1:PovStatBelow:SexMen:CES11
## -0.10717
(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.08 1 1 0.7715
## (1 | subclass) 27.02 1 kept 0
##
## Fixed effects:
## Sum Sq Mean Sq NumDF DenDF F.value elim.num
## Age:IPVstatus:PovStat:Sex:CES1 0.0733 0.0733 1 74.83 1.3905 1
## Age:IPVstatus:PovStat:CES1 0.2073 0.2073 1 76.72 0.0890 2
## Age:PovStat:Sex:CES1 0.0034 0.0034 1 81.31 0.4199 3
## IPVstatus:PovStat:Sex:CES1 0.0422 0.0422 1 84.64 0.4761 4
## IPVstatus:PovStat:CES1 0.1542 0.1542 1 88.79 0.0660 5
## Age:PovStat:CES1 0.0126 0.0126 1 83.82 1.5658 6
## Age 2.0697 2.0697 1 96.86 28.7149 kept
## IPVstatus 0.1219 0.1219 1 87.50 18.0699 kept
## PovStat 1.3519 1.3519 1 83.89 28.8773 kept
## Sex 0.0677 0.0677 1 90.92 4.8506 kept
## CES1 0.0057 0.0057 1 89.77 3.4502 kept
## Age:IPVstatus 0.4919 0.4919 1 87.59 16.3846 kept
## Age:PovStat 0.1932 0.1932 1 89.03 8.2177 kept
## Age:Sex 0.0311 0.0311 1 83.35 5.3015 kept
## Age:CES1 0.0100 0.0100 1 89.31 0.4787 kept
## IPVstatus:PovStat 0.0381 0.0381 1 97.06 8.9974 kept
## IPVstatus:Sex 0.1985 0.1985 1 85.99 0.4014 kept
## IPVstatus:CES1 0.0140 0.0140 1 96.01 1.2870 kept
## PovStat:Sex 0.1329 0.1329 1 87.92 0.6062 kept
## PovStat:CES1 0.2606 0.2606 1 97.42 2.4621 kept
## Sex:CES1 0.3697 0.3697 1 97.01 1.7074 kept
## Age:IPVstatus:PovStat 0.4935 0.4935 1 97.79 10.0300 kept
## Age:IPVstatus:Sex 0.0000 0.0000 1 85.53 1.3564 kept
## Age:IPVstatus:CES1 0.0427 0.0427 1 93.11 0.4878 kept
## Age:PovStat:Sex 0.0174 0.0174 1 91.19 1.8130 kept
## Age:Sex:CES1 0.4308 0.4308 1 94.72 4.2872 kept
## IPVstatus:PovStat:Sex 0.0000 0.0000 1 91.00 4.8846 kept
## IPVstatus:Sex:CES1 0.1554 0.1554 1 93.03 10.9756 kept
## PovStat:Sex:CES1 0.1581 0.1581 1 95.31 4.5672 kept
## Age:IPVstatus:PovStat:Sex 0.2898 0.2898 1 87.40 6.0739 kept
## Age:IPVstatus:Sex:CES1 0.2693 0.2693 1 91.21 9.9724 kept
## Pr(>F)
## Age:IPVstatus:PovStat:Sex:CES1 0.2421
## Age:IPVstatus:PovStat:CES1 0.7663
## Age:PovStat:Sex:CES1 0.5188
## IPVstatus:PovStat:Sex:CES1 0.4921
## IPVstatus:PovStat:CES1 0.7979
## Age:PovStat:CES1 0.2143
## Age 0e+00
## IPVstatus 1e-04
## PovStat 0e+00
## Sex 0.0302
## CES1 0.0665
## Age:IPVstatus 0.0001
## Age:PovStat 0.0052
## Age:Sex 0.0238
## Age:CES1 0.4908
## IPVstatus:PovStat 0.0034
## IPVstatus:Sex 0.5280
## IPVstatus:CES1 0.2594
## PovStat:Sex 0.4383
## PovStat:CES1 0.1199
## Sex:CES1 0.1944
## Age:IPVstatus:PovStat 0.0021
## Age:IPVstatus:Sex 0.2474
## Age:IPVstatus:CES1 0.4866
## Age:PovStat:Sex 0.1815
## Age:Sex:CES1 0.0411
## IPVstatus:PovStat:Sex 0.0296
## IPVstatus:Sex:CES1 0.0013
## PovStat:Sex:CES1 0.0351
## Age:IPVstatus:PovStat:Sex 0.0157
## Age:IPVstatus:Sex:CES1 0.0022
##
## Least squares means:
## IPVstatus PovStat Sex CES1 Estimate
## IPVstatus 0 1.0 NA NA NA 3.41
## IPVstatus 1 2.0 NA NA NA 3.55
## PovStat Above NA 1.0 NA NA 3.33
## PovStat Below NA 2.0 NA NA 3.63
## Sex Women NA NA 2.0 NA 3.51
## Sex Men NA NA 1.0 NA 3.45
## CES1 0 NA NA NA 1.0 3.55
## CES1 1 NA NA NA 2.0 3.41
## IPVstatus:PovStat 0 Above 1.0 1.0 NA NA 3.30
## IPVstatus:PovStat 1 Above 2.0 1.0 NA NA 3.35
## IPVstatus:PovStat 0 Below 1.0 2.0 NA NA 3.52
## IPVstatus:PovStat 1 Below 2.0 2.0 NA NA 3.75
## IPVstatus:Sex 0 Women 1.0 NA 2.0 NA 3.45
## IPVstatus:Sex 1 Women 2.0 NA 2.0 NA 3.56
## IPVstatus:Sex 0 Men 1.0 NA 1.0 NA 3.37
## IPVstatus:Sex 1 Men 2.0 NA 1.0 NA 3.54
## IPVstatus:CES1 0 0 1.0 NA NA 1.0 3.45
## IPVstatus:CES1 1 0 2.0 NA NA 1.0 3.65
## IPVstatus:CES1 0 1 1.0 NA NA 2.0 3.37
## IPVstatus:CES1 1 1 2.0 NA NA 2.0 3.45
## PovStat:Sex Above Women NA 1.0 2.0 NA 3.37
## PovStat:Sex Below Women NA 2.0 2.0 NA 3.65
## PovStat:Sex Above Men NA 1.0 1.0 NA 3.28
## PovStat:Sex Below Men NA 2.0 1.0 NA 3.62
## PovStat:CES1 Above 0 NA 1.0 NA 1.0 3.33
## PovStat:CES1 Below 0 NA 2.0 NA 1.0 3.77
## PovStat:CES1 Above 1 NA 1.0 NA 2.0 3.32
## PovStat:CES1 Below 1 NA 2.0 NA 2.0 3.50
## Sex:CES1 Women 0 NA NA 2.0 1.0 3.56
## Sex:CES1 Men 0 NA NA 1.0 1.0 3.54
## Sex:CES1 Women 1 NA NA 2.0 2.0 3.46
## Sex:CES1 Men 1 NA NA 1.0 2.0 3.36
## IPVstatus:PovStat:Sex 0 Above Women 1.0 1.0 2.0 NA 3.39
## IPVstatus:PovStat:Sex 1 Above Women 2.0 1.0 2.0 NA 3.35
## IPVstatus:PovStat:Sex 0 Below Women 1.0 2.0 2.0 NA 3.52
## IPVstatus:PovStat:Sex 1 Below Women 2.0 2.0 2.0 NA 3.78
## IPVstatus:PovStat:Sex 0 Above Men 1.0 1.0 1.0 NA 3.21
## IPVstatus:PovStat:Sex 1 Above Men 2.0 1.0 1.0 NA 3.35
## IPVstatus:PovStat:Sex 0 Below Men 1.0 2.0 1.0 NA 3.52
## IPVstatus:PovStat:Sex 1 Below Men 2.0 2.0 1.0 NA 3.73
## IPVstatus:Sex:CES1 0 Women 0 1.0 NA 2.0 1.0 3.42
## IPVstatus:Sex:CES1 1 Women 0 2.0 NA 2.0 1.0 3.70
## IPVstatus:Sex:CES1 0 Men 0 1.0 NA 1.0 1.0 3.47
## IPVstatus:Sex:CES1 1 Men 0 2.0 NA 1.0 1.0 3.61
## IPVstatus:Sex:CES1 0 Women 1 1.0 NA 2.0 2.0 3.48
## IPVstatus:Sex:CES1 1 Women 1 2.0 NA 2.0 2.0 3.43
## IPVstatus:Sex:CES1 0 Men 1 1.0 NA 1.0 2.0 3.26
## IPVstatus:Sex:CES1 1 Men 1 2.0 NA 1.0 2.0 3.46
## PovStat:Sex:CES1 Above Women 0 NA 1.0 2.0 1.0 3.27
## PovStat:Sex:CES1 Below Women 0 NA 2.0 2.0 1.0 3.85
## PovStat:Sex:CES1 Above Men 0 NA 1.0 1.0 1.0 3.39
## PovStat:Sex:CES1 Below Men 0 NA 2.0 1.0 1.0 3.69
## PovStat:Sex:CES1 Above Women 1 NA 1.0 2.0 2.0 3.47
## PovStat:Sex:CES1 Below Women 1 NA 2.0 2.0 2.0 3.45
## PovStat:Sex:CES1 Above Men 1 NA 1.0 1.0 2.0 3.17
## PovStat:Sex:CES1 Below Men 1 NA 2.0 1.0 2.0 3.55
## Standard Error DF t-value Lower CI
## IPVstatus 0 0.0678 23.2 50.3 3.27
## IPVstatus 1 0.0812 42.7 43.7 3.39
## PovStat Above 0.0674 22.9 49.4 3.19
## PovStat Below 0.0820 41.5 44.3 3.47
## Sex Women 0.0700 25.5 50.1 3.37
## Sex Men 0.0855 46.3 40.4 3.28
## CES1 0 0.0760 33.5 46.7 3.40
## CES1 1 0.0730 30.0 46.7 3.26
## IPVstatus:PovStat 0 Above 0.0681 23.7 48.5 3.16
## IPVstatus:PovStat 1 Above 0.0907 58.7 37.0 3.17
## IPVstatus:PovStat 0 Below 0.0899 56.1 39.1 3.34
## IPVstatus:PovStat 1 Below 0.1089 75.0 34.5 3.54
## IPVstatus:Sex 0 Women 0.0715 27.2 48.3 3.31
## IPVstatus:Sex 1 Women 0.0904 55.8 39.4 3.38
## IPVstatus:Sex 0 Men 0.0911 52.6 37.0 3.18
## IPVstatus:Sex 1 Men 0.1149 84.3 30.8 3.31
## IPVstatus:CES1 0 0 0.0733 29.8 47.0 3.30
## IPVstatus:CES1 1 0 0.1106 78.8 33.0 3.43
## IPVstatus:CES1 0 1 0.0806 41.6 41.9 3.21
## IPVstatus:CES1 1 1 0.0885 52.6 39.0 3.27
## PovStat:Sex Above Women 0.0699 25.5 48.2 3.23
## PovStat:Sex Below Women 0.0907 55.3 40.2 3.47
## PovStat:Sex Above Men 0.0892 52.4 36.8 3.10
## PovStat:Sex Below Men 0.1131 80.0 32.0 3.40
## PovStat:CES1 Above 0 0.0915 55.4 36.4 3.15
## PovStat:CES1 Below 0 0.0931 55.6 40.5 3.58
## PovStat:CES1 Above 1 0.0749 31.2 44.3 3.17
## PovStat:CES1 Below 1 0.1063 74.1 32.9 3.29
## Sex:CES1 Women 0 0.0869 49.9 41.0 3.39
## Sex:CES1 Men 0 0.0994 67.1 35.6 3.34
## Sex:CES1 Women 1 0.0822 38.2 42.1 3.29
## Sex:CES1 Men 1 0.0987 66.0 34.1 3.17
## IPVstatus:PovStat:Sex 0 Above Women 0.0768 34.1 44.1 3.24
## IPVstatus:PovStat:Sex 1 Above Women 0.0849 49.4 39.4 3.18
## IPVstatus:PovStat:Sex 0 Below Women 0.1001 61.5 35.1 3.32
## IPVstatus:PovStat:Sex 1 Below Women 0.1380 94.5 27.4 3.51
## IPVstatus:PovStat:Sex 0 Above Men 0.0865 47.7 37.1 3.04
## IPVstatus:PovStat:Sex 1 Above Men 0.1426 98.3 23.5 3.07
## IPVstatus:PovStat:Sex 0 Below Men 0.1371 96.1 25.7 3.25
## IPVstatus:PovStat:Sex 1 Below Men 0.1538 95.6 24.2 3.42
## IPVstatus:Sex:CES1 0 Women 0 0.0804 42.1 42.6 3.26
## IPVstatus:Sex:CES1 1 Women 0 0.1284 90.4 28.8 3.44
## IPVstatus:Sex:CES1 0 Men 0 0.0960 59.7 36.1 3.28
## IPVstatus:Sex:CES1 1 Men 0 0.1661 99.5 21.8 3.29
## IPVstatus:Sex:CES1 0 Women 1 0.0916 51.9 38.0 3.30
## IPVstatus:Sex:CES1 1 Women 1 0.0986 63.5 34.8 3.24
## IPVstatus:Sex:CES1 0 Men 1 0.1177 83.9 27.7 3.03
## IPVstatus:Sex:CES1 1 Men 1 0.1298 91.2 26.7 3.21
## PovStat:Sex:CES1 Above Women 0 0.0827 43.9 39.6 3.11
## PovStat:Sex:CES1 Below Women 0 0.1233 89.5 31.2 3.60
## PovStat:Sex:CES1 Above Men 0 0.1401 96.2 24.2 3.11
## PovStat:Sex:CES1 Below Men 0 0.1153 80.0 32.0 3.46
## PovStat:Sex:CES1 Above Women 1 0.0830 41.6 41.8 3.30
## PovStat:Sex:CES1 Below Women 1 0.1074 72.7 32.1 3.23
## PovStat:Sex:CES1 Above Men 1 0.1056 62.4 30.1 2.96
## PovStat:Sex:CES1 Below Men 1 0.1682 99.3 21.1 3.22
## Upper CI p-value
## IPVstatus 0 3.55 <2e-16
## IPVstatus 1 3.72 <2e-16
## PovStat Above 3.47 <2e-16
## PovStat Below 3.80 <2e-16
## Sex Women 3.65 <2e-16
## Sex Men 3.62 <2e-16
## CES1 0 3.71 <2e-16
## CES1 1 3.56 <2e-16
## IPVstatus:PovStat 0 Above 3.44 <2e-16
## IPVstatus:PovStat 1 Above 3.53 <2e-16
## IPVstatus:PovStat 0 Below 3.70 <2e-16
## IPVstatus:PovStat 1 Below 3.97 <2e-16
## IPVstatus:Sex 0 Women 3.60 <2e-16
## IPVstatus:Sex 1 Women 3.75 <2e-16
## IPVstatus:Sex 0 Men 3.55 <2e-16
## IPVstatus:Sex 1 Men 3.77 <2e-16
## IPVstatus:CES1 0 0 3.60 <2e-16
## IPVstatus:CES1 1 0 3.88 <2e-16
## IPVstatus:CES1 0 1 3.54 <2e-16
## IPVstatus:CES1 1 1 3.63 <2e-16
## PovStat:Sex Above Women 3.51 <2e-16
## PovStat:Sex Below Women 3.83 <2e-16
## PovStat:Sex Above Men 3.46 <2e-16
## PovStat:Sex Below Men 3.85 <2e-16
## PovStat:CES1 Above 0 3.52 <2e-16
## PovStat:CES1 Below 0 3.96 <2e-16
## PovStat:CES1 Above 1 3.47 <2e-16
## PovStat:CES1 Below 1 3.71 <2e-16
## Sex:CES1 Women 0 3.73 <2e-16
## Sex:CES1 Men 0 3.74 <2e-16
## Sex:CES1 Women 1 3.62 <2e-16
## Sex:CES1 Men 1 3.56 <2e-16
## IPVstatus:PovStat:Sex 0 Above Women 3.55 <2e-16
## IPVstatus:PovStat:Sex 1 Above Women 3.52 <2e-16
## IPVstatus:PovStat:Sex 0 Below Women 3.72 <2e-16
## IPVstatus:PovStat:Sex 1 Below Women 4.05 <2e-16
## IPVstatus:PovStat:Sex 0 Above Men 3.39 <2e-16
## IPVstatus:PovStat:Sex 1 Above Men 3.64 <2e-16
## IPVstatus:PovStat:Sex 0 Below Men 3.79 <2e-16
## IPVstatus:PovStat:Sex 1 Below Men 4.03 <2e-16
## IPVstatus:Sex:CES1 0 Women 0 3.59 <2e-16
## IPVstatus:Sex:CES1 1 Women 0 3.95 <2e-16
## IPVstatus:Sex:CES1 0 Men 0 3.66 <2e-16
## IPVstatus:Sex:CES1 1 Men 0 3.94 <2e-16
## IPVstatus:Sex:CES1 0 Women 1 3.67 <2e-16
## IPVstatus:Sex:CES1 1 Women 1 3.63 <2e-16
## IPVstatus:Sex:CES1 0 Men 1 3.50 <2e-16
## IPVstatus:Sex:CES1 1 Men 1 3.72 <2e-16
## PovStat:Sex:CES1 Above Women 0 3.44 <2e-16
## PovStat:Sex:CES1 Below Women 0 4.09 <2e-16
## PovStat:Sex:CES1 Above Men 0 3.67 <2e-16
## PovStat:Sex:CES1 Below Men 0 3.92 <2e-16
## PovStat:Sex:CES1 Above Women 1 3.64 <2e-16
## PovStat:Sex:CES1 Below Women 1 3.66 <2e-16
## PovStat:Sex:CES1 Above Men 1 3.39 <2e-16
## PovStat:Sex:CES1 Below Men 1 3.89 <2e-16
##
## Differences of LSMEANS:
## Estimate Standard Error DF t-value
## IPVstatus 0-1 -0.1 0.0693 87.7 -2.05
## PovStat Above-Below -0.3 0.0703 94.4 -4.38
## Sex Women-Men 0.1 0.0826 96.3 0.68
## CES1 0-1 0.1 0.0680 93.0 2.06
## IPVstatus:PovStat 0 Above- 1 Above 0.0 0.0870 87.1 -0.57
## IPVstatus:PovStat 0 Above- 0 Below -0.2 0.0839 88.0 -2.57
## IPVstatus:PovStat 0 Above- 1 Below -0.5 0.1038 93.0 -4.34
## IPVstatus:PovStat 1 Above- 0 Below -0.2 0.0933 89.0 -1.78
## IPVstatus:PovStat 1 Above- 1 Below -0.4 0.1174 95.2 -3.41
## IPVstatus:PovStat 0 Below- 1 Below -0.2 0.1139 91.4 -2.06
## IPVstatus:Sex 0 Women- 1 Women -0.1 0.0835 90.5 -1.32
## IPVstatus:Sex 0 Women- 0 Men 0.1 0.0917 98.3 0.96
## IPVstatus:Sex 0 Women- 1 Men -0.1 0.1087 90.0 -0.79
## IPVstatus:Sex 1 Women- 0 Men 0.2 0.1069 96.0 1.86
## IPVstatus:Sex 1 Women- 1 Men 0.0 0.1279 92.5 0.19
## IPVstatus:Sex 0 Men- 1 Men -0.2 0.1173 90.8 -1.48
## IPVstatus:CES1 0 0- 1 0 -0.2 0.1100 92.9 -1.89
## IPVstatus:CES1 0 0- 0 1 0.1 0.0730 91.4 1.02
## IPVstatus:CES1 0 0- 1 1 0.0 0.0862 90.9 -0.02
## IPVstatus:CES1 1 0- 0 1 0.3 0.1068 90.2 2.65
## IPVstatus:CES1 1 0- 1 1 0.2 0.1172 95.4 1.76
## IPVstatus:CES1 0 1- 1 1 -0.1 0.0856 90.7 -0.90
## PovStat:Sex Above Women- Below Women -0.3 0.0813 91.3 -3.42
## PovStat:Sex Above Women- Above Men 0.1 0.0867 96.3 1.00
## PovStat:Sex Above Women- Below Men -0.3 0.1111 95.6 -2.26
## PovStat:Sex Below Women- Above Men 0.4 0.1057 95.4 3.45
## PovStat:Sex Below Women- Below Men 0.0 0.1230 92.3 0.21
## PovStat:Sex Above Men- Below Men -0.3 0.1107 92.4 -3.05
## PovStat:CES1 Above 0- Below 0 -0.4 0.1047 96.6 -4.18
## PovStat:CES1 Above 0- Above 1 0.0 0.0991 97.3 0.11
## PovStat:CES1 Above 0- Below 1 -0.2 0.1058 94.4 -1.59
## PovStat:CES1 Below 0- Above 1 0.4 0.0890 92.8 5.04
## PovStat:CES1 Below 0- Below 1 0.3 0.1141 94.6 2.36
## PovStat:CES1 Above 1- Below 1 -0.2 0.1119 95.9 -1.60
## Sex:CES1 Women 0- Men 0 0.0 0.1083 91.5 0.16
## Sex:CES1 Women 0- Women 1 0.1 0.0950 98.3 1.07
## Sex:CES1 Women 0- Men 1 0.2 0.1065 91.2 1.85
## Sex:CES1 Men 0- Women 1 0.1 0.1074 97.9 0.78
## Sex:CES1 Men 0- Men 1 0.2 0.0998 87.9 1.80
## Sex:CES1 Women 1- Men 1 0.1 0.1080 98.5 0.88
## Lower CI Upper CI p-value
## IPVstatus 0-1 -0.2800 -0.0046 0.043
## PovStat Above-Below -0.4477 -0.1686 <2e-16
## Sex Women-Men -0.1074 0.2204 0.495
## CES1 0-1 0.0054 0.2753 0.042
## IPVstatus:PovStat 0 Above- 1 Above -0.2227 0.1231 0.568
## IPVstatus:PovStat 0 Above- 0 Below -0.3825 -0.0489 0.012
## IPVstatus:PovStat 0 Above- 1 Below -0.6565 -0.2443 <2e-16
## IPVstatus:PovStat 1 Above- 0 Below -0.3512 0.0195 0.079
## IPVstatus:PovStat 1 Above- 1 Below -0.6336 -0.1676 9e-04
## IPVstatus:PovStat 0 Below- 1 Below -0.4610 -0.0085 0.042
## IPVstatus:Sex 0 Women- 1 Women -0.2764 0.0554 0.189
## IPVstatus:Sex 0 Women- 0 Men -0.0936 0.2703 0.338
## IPVstatus:Sex 0 Women- 1 Men -0.3018 0.1302 0.432
## IPVstatus:Sex 1 Women- 0 Men -0.0133 0.4109 0.066
## IPVstatus:Sex 1 Women- 1 Men -0.2293 0.2787 0.847
## IPVstatus:Sex 0 Men- 1 Men -0.4070 0.0588 0.141
## IPVstatus:CES1 0 0- 1 0 -0.4263 0.0107 0.062
## IPVstatus:CES1 0 0- 0 1 -0.0703 0.2199 0.308
## IPVstatus:CES1 0 0- 1 1 -0.1732 0.1692 0.982
## IPVstatus:CES1 1 0- 0 1 0.0704 0.4948 0.010
## IPVstatus:CES1 1 0- 1 1 -0.0268 0.4385 0.082
## IPVstatus:CES1 0 1- 1 1 -0.2468 0.0932 0.372
## PovStat:Sex Above Women- Below Women -0.4396 -0.1165 9e-04
## PovStat:Sex Above Women- Above Men -0.0854 0.2586 0.320
## PovStat:Sex Above Women- Below Men -0.4722 -0.0311 0.026
## PovStat:Sex Below Women- Above Men 0.1548 0.5745 8e-04
## PovStat:Sex Below Women- Below Men -0.2179 0.2707 0.830
## PovStat:Sex Above Men- Below Men -0.5581 -0.1183 0.003
## PovStat:CES1 Above 0- Below 0 -0.6454 -0.2298 1e-04
## PovStat:CES1 Above 0- Above 1 -0.1858 0.2076 0.913
## PovStat:CES1 Above 0- Below 1 -0.3779 0.0422 0.116
## PovStat:CES1 Below 0- Above 1 0.2717 0.6252 <2e-16
## PovStat:CES1 Below 0- Below 1 0.0432 0.4963 0.020
## PovStat:CES1 Above 1- Below 1 -0.4008 0.0434 0.114
## Sex:CES1 Women 0- Men 0 -0.1975 0.2328 0.871
## Sex:CES1 Women 0- Women 1 -0.0871 0.2900 0.288
## Sex:CES1 Women 0- Men 1 -0.0148 0.4085 0.068
## Sex:CES1 Men 0- Women 1 -0.1292 0.2969 0.437
## Sex:CES1 Men 0- Men 1 -0.0192 0.3775 0.076
## Sex:CES1 Women 1- Men 1 -0.1189 0.3097 0.379
##
## Final model:
## lme4::lmer(formula = logTrailsA ~ 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(logTrailsA ~ 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: logTrailsA ~ 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: 114.8
## Random effects:
## Groups Name Std.Dev. Corr
## HNDid (Intercept) 0.23993
## Age 0.00999 1.00
## subclass (Intercept) 0.23593
## Residual 0.19117
## Number of obs: 126, groups: HNDid, 63; subclass, 21
## Fixed Effects:
## (Intercept) Age
## 3.30551 0.01561
## IPVstatus1 PovStatBelow
## 0.54121 0.31640
## SexMen CES11
## 0.15118 0.45066
## Age:IPVstatus1 Age:PovStatBelow
## 0.05086 -0.01324
## Age:SexMen Age:CES11
## -0.00178 0.01273
## IPVstatus1:PovStatBelow IPVstatus1:SexMen
## 0.98871 -0.66280
## IPVstatus1:CES11 PovStatBelow:SexMen
## -1.03006 -0.08332
## PovStatBelow:CES11 SexMen:CES11
## -0.52057 -0.79333
## Age:IPVstatus1:PovStatBelow Age:IPVstatus1:SexMen
## 0.11065 -0.07175
## Age:IPVstatus1:CES11 Age:PovStatBelow:SexMen
## -0.09548 0.01042
## Age:SexMen:CES11 IPVstatus1:PovStatBelow:SexMen
## -0.02649 -0.64647
## IPVstatus1:SexMen:CES11 PovStatBelow:SexMen:CES11
## 1.56956 0.60334
## Age:IPVstatus1:PovStatBelow:SexMen Age:IPVstatus1:SexMen:CES11
## -0.07537 0.15571
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: logTrailsA ~ 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: 114.8
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 5.76e-02 0.23993
## Age 9.98e-05 0.00999 1.00
## subclass (Intercept) 5.57e-02 0.23593
## Residual 3.65e-02 0.19117
## Number of obs: 126, groups: HNDid, 63; subclass, 21
##
## Fixed effects:
## Estimate Std. Error df t value
## (Intercept) 3.30551 0.15427 20.50000 21.43
## Age 0.01561 0.01076 38.20000 1.45
## IPVstatus1 0.54121 0.27496 43.30000 1.97
## PovStatBelow 0.31640 0.22467 29.60000 1.41
## SexMen 0.15118 0.20136 16.50000 0.75
## CES11 0.45066 0.19599 25.90000 2.30
## Age:IPVstatus1 0.05086 0.02336 82.40000 2.18
## Age:PovStatBelow -0.01324 0.01615 32.60000 -0.82
## Age:SexMen -0.00178 0.01519 29.50000 -0.12
## Age:CES11 0.01273 0.01530 66.80000 0.83
## IPVstatus1:PovStatBelow 0.98871 0.42255 43.20000 2.34
## IPVstatus1:SexMen -0.66280 0.54993 32.00000 -1.21
## IPVstatus1:CES11 -1.03006 0.34485 38.20000 -2.99
## PovStatBelow:SexMen -0.08332 0.29852 22.00000 -0.28
## PovStatBelow:CES11 -0.52057 0.16620 18.20000 -3.13
## SexMen:CES11 -0.79333 0.34596 46.20000 -2.29
## Age:IPVstatus1:PovStatBelow 0.11065 0.03737 75.70000 2.96
## Age:IPVstatus1:SexMen -0.07175 0.03836 58.80000 -1.87
## Age:IPVstatus1:CES11 -0.09548 0.02960 75.80000 -3.23
## Age:PovStatBelow:SexMen 0.01042 0.02638 31.30000 0.39
## Age:SexMen:CES11 -0.02649 0.02998 80.40000 -0.88
## IPVstatus1:PovStatBelow:SexMen -0.64647 0.58127 34.70000 -1.11
## IPVstatus1:SexMen:CES11 1.56956 0.53900 22.70000 2.91
## PovStatBelow:SexMen:CES11 0.60334 0.35877 53.40000 1.68
## Age:IPVstatus1:PovStatBelow:SexMen -0.07537 0.05057 63.20000 -1.49
## Age:IPVstatus1:SexMen:CES11 0.15571 0.04565 63.90000 3.41
## Pr(>|t|)
## (Intercept) 1.6e-15
## Age 0.1549
## IPVstatus1 0.0555
## PovStatBelow 0.1695
## SexMen 0.4634
## CES11 0.0298
## Age:IPVstatus1 0.0323
## Age:PovStatBelow 0.4183
## Age:SexMen 0.9078
## Age:CES11 0.4081
## IPVstatus1:PovStatBelow 0.0240
## IPVstatus1:SexMen 0.2369
## IPVstatus1:CES11 0.0049
## PovStatBelow:SexMen 0.7828
## PovStatBelow:CES11 0.0057
## SexMen:CES11 0.0264
## Age:IPVstatus1:PovStatBelow 0.0041
## Age:IPVstatus1:SexMen 0.0664
## Age:IPVstatus1:CES11 0.0019
## Age:PovStatBelow:SexMen 0.6956
## Age:SexMen:CES11 0.3795
## IPVstatus1:PovStatBelow:SexMen 0.2737
## IPVstatus1:SexMen:CES11 0.0079
## PovStatBelow:SexMen:CES11 0.0985
## Age:IPVstatus1:PovStatBelow:SexMen 0.1411
## Age:IPVstatus1:SexMen:CES11 0.0011
##
## 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)
