## 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
Trails B Regression Model 1
(mm1 = lmer(TrailsBtestSec ~ (Sex + CES1)^2 + (Age | HNDid), data = AllPsDATA))
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: TrailsBtestSec ~ (Sex + CES1)^2 + (Age | HNDid)
## Data: AllPsDATA
## REML criterion at convergence: 58683
## Random effects:
## Groups Name Std.Dev. Corr
## HNDid (Intercept) 130.508
## Age 0.876 1.00
## Residual 77.877
## Number of obs: 4822, groups: HNDid, 1217
## Fixed Effects:
## (Intercept) SexMen CES11 SexMen:CES11
## 121.76 17.70 33.38 -1.96
(st = step(mm1))
## Random term (Age | HNDid) was eliminated because of having correlation +-1 or NaN
##
## 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) 3264 1 kept < 1e-07
##
## Fixed effects:
## Sum Sq Mean Sq NumDF DenDF F.value elim.num Pr(>F)
## Sex:CES1 261.2 261.2 1 1195 0.0431 1 0.8355
## Sex 21229.1 21229.1 1 1196 4.8832 kept 0.0273
## CES1 101238.6 101238.6 1 1196 16.7263 kept 0
##
## Least squares means:
## Sex CES1 Estimate Standard Error DF t-value Lower CI
## Sex Women 2 NA 140.92 5.16 1196 27.3 131
## Sex Men 1 NA 158.53 6.19 1196 25.6 146
## CES1 0 NA 1 133.28 5.05 1195 26.4 123
## CES1 1 NA 2 166.17 6.32 1197 26.3 154
## Upper CI p-value
## Sex Women 151 <2e-16
## Sex Men 171 <2e-16
## CES1 0 143 <2e-16
## CES1 1 179 <2e-16
##
## Differences of LSMEANS:
## Estimate Standard Error DF t-value Lower CI Upper CI
## Sex Women-Men -17.6 7.97 1195.5 -2.21 -33.3 -1.98
## CES1 0-1 -32.9 8.04 1196.1 -4.09 -48.7 -17.11
## p-value
## Sex Women-Men 0.03
## CES1 0-1 <2e-16
##
## Final model:
## lme4::lmer(formula = TrailsBtestSec ~ Sex + CES1 + (1 | HNDid),
## data = AllPsDATA, REML = reml, contrasts = l)
summary(mm1)
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: TrailsBtestSec ~ (Sex + CES1)^2 + (Age | HNDid)
## Data: AllPsDATA
##
## REML criterion at convergence: 58683
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 1.70e+04 130.508
## Age 7.68e-01 0.876 1.00
## Residual 6.06e+03 77.877
## Number of obs: 4822, groups: HNDid, 1217
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 121.76 6.75 1177.00 18.04 <2e-16
## SexMen 17.70 10.05 1176.00 1.76 0.0784
## CES11 33.38 10.25 1151.00 3.26 0.0012
## SexMen:CES11 -1.96 16.23 1154.00 -0.12 0.9040
##
## Correlation of Fixed Effects:
## (Intr) SexMen CES11
## SexMen -0.672
## CES11 -0.659 0.443
## SexMn:CES11 0.416 -0.619 -0.631
Trails B Regression Model 2
(mm2 = lmer(TrailsBtestSec ~ (Age + Sex + CES1 + Race)^4 + (Age | HNDid), data = AllPsDATA))
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: TrailsBtestSec ~ (Age + Sex + CES1 + Race)^4 + (Age | HNDid)
## Data: AllPsDATA
## REML criterion at convergence: 58548
## Random effects:
## Groups Name Std.Dev. Corr
## HNDid (Intercept) 124.759
## Age 0.753 1.00
## Residual 78.072
## Number of obs: 4822, groups: HNDid, 1217
## Fixed Effects:
## (Intercept) Age
## 87.1194 0.8000
## SexMen CES11
## 32.1114 33.5799
## RaceAfrAm Age:SexMen
## 64.6406 1.2380
## Age:CES11 Age:RaceAfrAm
## 0.0861 0.8495
## SexMen:CES11 SexMen:RaceAfrAm
## -10.9088 -23.1187
## CES11:RaceAfrAm Age:SexMen:CES11
## 12.9033 0.2467
## Age:SexMen:RaceAfrAm Age:CES11:RaceAfrAm
## -0.9019 1.0300
## SexMen:CES11:RaceAfrAm Age:SexMen:CES11:RaceAfrAm
## 6.7023 -2.0882
(st = step(mm2))
## Random term (Age | HNDid) was eliminated because of having correlation +-1 or NaN
##
## 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) 3002 1 kept < 1e-07
##
## Fixed effects:
## Sum Sq Mean Sq NumDF DenDF F.value elim.num Pr(>F)
## Age:Sex:CES1:Race 4313.3 4313.3 1 3230 0.7083 1 0.4001
## Age:CES1:Race 245.1 245.1 1 3168 0.0407 2 0.8402
## Sex:CES1:Race 285.5 285.5 1 1179 0.0464 3 0.8294
## Age:Sex:CES1 2771.6 2771.6 1 3249 0.5355 4 0.4643
## Age:CES1 386.8 386.8 1 3190 0.1444 5 0.7040
## Sex:CES1 874.3 874.3 1 1179 0.1754 6 0.6755
## CES1:Race 5901.8 5901.8 1 1181 0.9686 7 0.3252
## Age:Sex:Race 12279.3 12279.3 1 3205 1.8927 8 0.1690
## Age:Sex 3156.5 3156.5 1 3217 0.4474 9 0.5036
## Age:Race 5570.8 5570.8 1 3160 0.8812 10 0.3479
## Sex:Race 11938.6 11938.6 1 1181 1.9999 11 0.1576
## Age 142716.5 142716.5 1 3186 28.9689 kept <1e-07
## Sex 23230.5 23230.5 1 1182 5.2925 kept 0.0216
## CES1 124342.4 124342.4 1 1185 23.2162 kept 0
## Race 399006.1 399006.1 1 1183 65.4065 kept <1e-07
##
## Least squares means:
## Sex CES1 Race Estimate Standard Error DF t-value Lower CI
## Sex Women 2 NA NA 138.93 4.96 1182 28.0 129
## Sex Men 1 NA NA 156.54 5.95 1182 26.3 145
## CES1 0 NA 1 NA 129.10 4.87 1182 26.5 120
## CES1 1 NA 2 NA 166.38 6.08 1185 27.4 154
## Race White NA NA 2 117.18 5.63 1180 20.8 106
## Race AfrAm NA NA 1 178.29 5.25 1186 33.9 168
## Upper CI p-value
## Sex Women 149 <2e-16
## Sex Men 168 <2e-16
## CES1 0 139 <2e-16
## CES1 1 178 <2e-16
## Race White 128 <2e-16
## Race AfrAm 189 <2e-16
##
## Differences of LSMEANS:
## Estimate Standard Error DF t-value Lower CI Upper CI
## Sex Women-Men -17.6 7.65 1181.7 -2.30 -32.6 -2.59
## CES1 0-1 -37.3 7.74 1185.4 -4.82 -52.5 -22.10
## Race White-AfrAm -61.1 7.56 1183.3 -8.09 -75.9 -46.29
## p-value
## Sex Women-Men 0.02
## CES1 0-1 <2e-16
## Race White-AfrAm <2e-16
##
## Final model:
## lme4::lmer(formula = TrailsBtestSec ~ Age + Sex + CES1 + Race +
## (1 | HNDid), data = AllPsDATA, REML = reml, contrasts = l)
Re-run the suggested final Model 2
(mm2 = lmer(TrailsBtestSec ~ Age + Sex + CES1 + Race + (Age | HNDid), data = AllPsDATA))
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: TrailsBtestSec ~ Age + Sex + CES1 + Race + (Age | HNDid)
## Data: AllPsDATA
## REML criterion at convergence: 58599
## Random effects:
## Groups Name Std.Dev. Corr
## HNDid (Intercept) 124.8
## Age 0.8 1.00
## Residual 78.0
## Number of obs: 4822, groups: HNDid, 1217
## Fixed Effects:
## (Intercept) Age SexMen CES11 RaceAfrAm
## 90.7 1.7 16.8 36.9 60.1
summary(mm2)
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: TrailsBtestSec ~ Age + Sex + CES1 + Race + (Age | HNDid)
## Data: AllPsDATA
##
## REML criterion at convergence: 58599
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 15580.49 124.8
## Age 0.64 0.8 1.00
## Residual 6087.38 78.0
## Number of obs: 4822, groups: HNDid, 1217
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 90.703 7.222 1173.000 12.56 < 2e-16
## Age 1.696 0.316 3113.000 5.36 8.8e-08
## SexMen 16.813 7.584 1154.000 2.22 0.027
## CES11 36.920 7.653 1146.000 4.82 1.6e-06
## RaceAfrAm 60.080 7.491 1158.000 8.02 2.7e-15
##
## Correlation of Fixed Effects:
## (Intr) Age SexMen CES11
## Age 0.024
## SexMen -0.473 0.000
## CES11 -0.478 0.051 0.085
## RaceAfrAm -0.575 0.045 -0.003 0.039
Trails B Regression Model 3
(mm3 = lmer(TrailsBtestSec ~ (Age + Sex + CES1 + Race + IPVincid)^5 + (Age |
HNDid), data = AllPsDATA))
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: TrailsBtestSec ~ (Age + Sex + CES1 + Race + IPVincid)^5 + (Age | HNDid)
## Data: AllPsDATA
## REML criterion at convergence: 58384
## Random effects:
## Groups Name Std.Dev. Corr
## HNDid (Intercept) 123.85
## Age 0.79 1.00
## Residual 78.15
## Number of obs: 4822, groups: HNDid, 1217
## Fixed Effects:
## (Intercept)
## 87.1807
## Age
## 0.8039
## SexMen
## 29.2018
## CES11
## 34.8687
## RaceAfrAm
## 65.7618
## IPVincid1
## -20.1807
## Age:SexMen
## 1.1573
## Age:CES11
## 0.0373
## Age:RaceAfrAm
## 0.7900
## Age:IPVincid1
## -0.8039
## SexMen:CES11
## -8.7396
## SexMen:RaceAfrAm
## -23.3510
## SexMen:IPVincid1
## 335.1964
## CES11:RaceAfrAm
## 12.6523
## CES11:IPVincid1
## -43.0578
## RaceAfrAm:IPVincid1
## -53.0790
## Age:SexMen:CES11
## 0.3554
## Age:SexMen:RaceAfrAm
## -0.4943
## Age:SexMen:IPVincid1
## 16.9215
## Age:CES11:RaceAfrAm
## 1.1417
## Age:CES11:IPVincid1
## -0.6422
## Age:RaceAfrAm:IPVincid1
## 0.2228
## SexMen:CES11:RaceAfrAm
## 2.1571
## SexMen:CES11:IPVincid1
## -344.9696
## SexMen:RaceAfrAm:IPVincid1
## 92.2699
## CES11:RaceAfrAm:IPVincid1
## -2.9857
## Age:SexMen:CES11:RaceAfrAm
## -2.3550
## Age:SexMen:CES11:IPVincid1
## -17.8292
## Age:SexMen:RaceAfrAm:IPVincid1
## -18.5973
## Age:CES11:RaceAfrAm:IPVincid1
## 0.2294
## SexMen:CES11:RaceAfrAm:IPVincid1
## 91.5900
## Age:SexMen:CES11:RaceAfrAm:IPVincid1
## 17.6067
(st = step(mm3))
## Random term (Age | HNDid) was eliminated because of having correlation +-1 or NaN
##
## 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) 2976 1 kept < 1e-07
##
## Fixed effects:
## Sum Sq Mean Sq NumDF DenDF F.value
## Age:Sex:CES1:Race:IPVincid 952.46 952.46 1 3992 0.1561
## Age:CES1:Race:IPVincid 460.22 460.22 1 4501 0.0383
## Sex:CES1:Race:IPVincid 560.24 560.24 1 2182 0.1290
## CES1:Race:IPVincid 180.33 180.33 1 1567 0.0652
## Age:Sex:CES1:IPVincid 4283.78 4283.78 1 4768 0.0534
## Age:CES1:IPVincid 4382.60 4382.60 1 4795 0.0354
## Age:Sex:CES1:Race 4852.89 4852.89 1 3177 0.7864
## Sex:CES1:Race 35.45 35.45 1 1175 0.0037
## Age:CES1:Race 421.13 421.13 1 3124 0.0666
## Age:Sex:CES1 2607.83 2607.83 1 3213 0.5664
## Age:CES1 402.33 402.33 1 3159 0.1295
## CES1:Race 3869.60 3869.60 1 1175 0.7117
## Age:Sex:Race:IPVincid 2428.32 2428.32 1 4464 1.8404
## Age:Sex:IPVincid 340.79 340.79 1 4558 0.3902
## Age:Race:IPVincid 9383.89 9383.89 1 4573 1.4221
## Age:IPVincid 96.20 96.20 1 4601 0.0718
## Age:Sex:Race 13157.12 13157.12 1 3193 1.6658
## Age:Sex 3225.71 3225.71 1 3211 0.8555
## Age:Race 5669.99 5669.99 1 3150 0.9896
## Sex:Race:IPVincid 14041.98 14041.98 1 1217 2.3809
## Race:IPVincid 1936.74 1936.74 1 1199 0.0204
## Sex:Race 12181.34 12181.34 1 1178 1.9384
## Sex:CES1:IPVincid 17739.15 17739.15 1 1228 2.9828
## Sex:CES1 931.26 931.26 1 1180 0.4288
## CES1:IPVincid 18168.40 18168.40 1 1212 3.2774
## Age 145630.00 145630.00 1 3187 29.0835
## Sex 23553.96 23553.96 1 1201 13.3453
## CES1 126166.37 126166.37 1 1184 22.2947
## Race 404656.30 404656.30 1 1182 65.0457
## IPVincid 984.54 984.54 1 1206 0.7688
## Sex:IPVincid 67532.93 67532.93 1 1201 9.9915
## elim.num Pr(>F)
## Age:Sex:CES1:Race:IPVincid 1 0.6928
## Age:CES1:Race:IPVincid 2 0.8449
## Sex:CES1:Race:IPVincid 3 0.7195
## CES1:Race:IPVincid 4 0.7985
## Age:Sex:CES1:IPVincid 5 0.8173
## Age:CES1:IPVincid 6 0.8508
## Age:Sex:CES1:Race 7 0.3752
## Sex:CES1:Race 8 0.9512
## Age:CES1:Race 9 0.7963
## Age:Sex:CES1 10 0.4517
## Age:CES1 11 0.7189
## CES1:Race 12 0.3990
## Age:Sex:Race:IPVincid 13 0.1750
## Age:Sex:IPVincid 14 0.5322
## Age:Race:IPVincid 15 0.2331
## Age:IPVincid 16 0.7887
## Age:Sex:Race 17 0.1969
## Age:Sex 18 0.3551
## Age:Race 19 0.3199
## Sex:Race:IPVincid 20 0.1231
## Race:IPVincid 21 0.8863
## Sex:Race 22 0.1641
## Sex:CES1:IPVincid 23 0.0844
## Sex:CES1 24 0.5127
## CES1:IPVincid 25 0.0705
## Age kept <1e-07
## Sex kept 0.0003
## CES1 kept 0
## Race kept <1e-07
## IPVincid kept 0.3808
## Sex:IPVincid kept 0.0016
##
## Least squares means:
## Sex CES1 Race IPVincid Estimate Standard Error DF
## Sex Women 2 NA NA NA 106.7 19.06 1181
## Sex Men 1 NA NA NA 213.4 22.21 1220
## CES1 0 NA 1 NA NA 141.8 15.16 1204
## CES1 1 NA 2 NA NA 178.3 15.15 1204
## Race White NA NA 2 NA 129.7 15.28 1203
## Race AfrAm NA NA 1 NA 190.4 14.99 1206
## IPVincid 0 NA NA NA 1 147.2 3.95 1181
## IPVincid 1 NA NA NA 2 172.9 29.07 1206
## Sex:IPVincid Women 0 2 NA NA 1 140.0 4.99 1181
## Sex:IPVincid Men 0 1 NA NA 1 154.4 5.99 1180
## Sex:IPVincid Women 1 2 NA NA 2 73.4 37.80 1181
## Sex:IPVincid Men 1 1 NA NA 2 272.5 44.04 1221
## t-value Lower CI Upper CI p-value
## Sex Women 5.60 69.293 144 <2e-16
## Sex Men 9.61 169.872 257 <2e-16
## CES1 0 9.35 112.068 172 <2e-16
## CES1 1 11.77 148.593 208 <2e-16
## Race White 8.49 99.710 160 <2e-16
## Race AfrAm 12.71 161.049 220 <2e-16
## IPVincid 0 37.30 139.454 155 <2e-16
## IPVincid 1 5.95 115.910 230 <2e-16
## Sex:IPVincid Women 0 28.05 130.223 150 <2e-16
## Sex:IPVincid Men 0 25.78 142.627 166 <2e-16
## Sex:IPVincid Women 1 1.94 -0.801 148 0.0525
## Sex:IPVincid Men 1 6.19 186.107 359 <2e-16
##
## Differences of LSMEANS:
## Estimate Standard Error DF t-value
## Sex Women-Men -106.7 29.22 1200.9 -3.65
## CES1 0-1 -36.5 7.73 1184.5 -4.72
## Race White-AfrAm -60.8 7.53 1182.3 -8.07
## IPVincid 0-1 -25.7 29.36 1205.9 -0.88
## Sex:IPVincid Women 0- Men 0 -14.4 7.70 1180.4 -1.86
## Sex:IPVincid Women 0- Women 1 66.7 38.14 1181.2 1.75
## Sex:IPVincid Women 0- Men 1 -132.5 44.33 1220.4 -2.99
## Sex:IPVincid Men 0- Women 1 81.0 38.29 1181.1 2.12
## Sex:IPVincid Men 0- Men 1 -118.1 44.47 1220.2 -2.66
## Sex:IPVincid Women 1- Men 1 -199.1 57.94 1201.3 -3.44
## Lower CI Upper CI p-value
## Sex Women-Men -164.08 -49.419 3e-04
## CES1 0-1 -51.65 -21.329 <2e-16
## Race White-AfrAm -75.55 -45.985 <2e-16
## IPVincid 0-1 -83.33 31.853 0.381
## Sex:IPVincid Women 0- Men 0 -29.46 0.748 0.062
## Sex:IPVincid Women 0- Women 1 -8.17 141.476 0.081
## Sex:IPVincid Women 0- Men 1 -219.46 -45.514 0.003
## Sex:IPVincid Men 0- Women 1 5.89 156.129 0.035
## Sex:IPVincid Men 0- Men 1 -205.39 -30.880 0.008
## Sex:IPVincid Women 1- Men 1 -312.82 -85.464 6e-04
##
## Final model:
## lme4::lmer(formula = TrailsBtestSec ~ Age + Sex + CES1 + Race +
## IPVincid + (1 | HNDid) + Sex:IPVincid, data = AllPsDATA,
## REML = reml, contrasts = l)
Re-run suggested final model 3
(mm3 = lmer(TrailsBtestSec ~ Age + Sex + CES1 + Race + IPVincid + (Age | HNDid) +
Sex:IPVincid, data = AllPsDATA))
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: TrailsBtestSec ~ Age + Sex + CES1 + Race + IPVincid + (Age | HNDid) + Sex:IPVincid
## Data: AllPsDATA
## REML criterion at convergence: 58570
## Random effects:
## Groups Name Std.Dev. Corr
## HNDid (Intercept) 124.35
## Age 0.83 1.00
## Residual 78.02
## Number of obs: 4822, groups: HNDid, 1217
## Fixed Effects:
## (Intercept) Age SexMen CES11
## 92.5 1.7 13.3 36.1
## RaceAfrAm IPVincid1 SexMen:IPVincid1
## 59.7 -65.3 184.3
summary(mm3)
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: TrailsBtestSec ~ Age + Sex + CES1 + Race + IPVincid + (Age | HNDid) + Sex:IPVincid
## Data: AllPsDATA
##
## REML criterion at convergence: 58570
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 1.55e+04 124.35
## Age 6.89e-01 0.83 1.00
## Residual 6.09e+03 78.02
## Number of obs: 4822, groups: HNDid, 1217
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 92.478 7.220 1172.000 12.81 < 2e-16
## Age 1.702 0.316 3105.000 5.38 7.9e-08
## SexMen 13.295 7.627 1154.000 1.74 0.082
## CES11 36.084 7.638 1142.000 4.72 2.6e-06
## RaceAfrAm 59.674 7.464 1154.000 8.00 3.1e-15
## IPVincid1 -65.343 36.254 752.000 -1.80 0.072
## SexMen:IPVincid1 184.326 55.841 810.000 3.30 0.001
##
## Correlation of Fixed Effects:
## (Intr) Age SexMen CES11 RcAfrA IPVnc1
## Age 0.024
## SexMen -0.477 0.000
## CES11 -0.476 0.047 0.088
## RaceAfrAm -0.573 0.044 -0.001 0.041
## IPVincid1 -0.072 0.054 0.084 -0.028 -0.007
## SxMn:IPVnc1 0.078 -0.007 -0.138 -0.025 -0.012 -0.646
Trails B Regression Model 4
(mm4 = lmer(TrailsBtestSec ~ (Age + Sex + CES1 + Race + IPVlife)^5 + (Age |
HNDid), data = AllPsDATA))
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: TrailsBtestSec ~ (Age + Sex + CES1 + Race + IPVlife)^5 + (Age | HNDid)
## Data: AllPsDATA
## REML criterion at convergence: 58418
## Random effects:
## Groups Name Std.Dev. Corr
## HNDid (Intercept) 124.478
## Age 0.807 1.00
## Residual 78.139
## Number of obs: 4822, groups: HNDid, 1217
## Fixed Effects:
## (Intercept) Age
## 85.4669 0.8547
## SexMen CES11
## 30.6436 33.1425
## RaceAfrAm IPVlife1
## 65.9285 20.8633
## Age:SexMen Age:CES11
## 1.0829 0.0643
## Age:RaceAfrAm Age:IPVlife1
## 0.8939 -0.7278
## SexMen:CES11 SexMen:RaceAfrAm
## -5.7430 -23.0316
## SexMen:IPVlife1 CES11:RaceAfrAm
## 106.0942 20.4702
## CES11:IPVlife1 RaceAfrAm:IPVlife1
## -11.1286 -20.1423
## Age:SexMen:CES11 Age:SexMen:RaceAfrAm
## 0.3856 -0.5843
## Age:SexMen:IPVlife1 Age:CES11:RaceAfrAm
## 6.0816 1.0380
## Age:CES11:IPVlife1 Age:RaceAfrAm:IPVlife1
## 0.6525 -0.8982
## SexMen:CES11:RaceAfrAm SexMen:CES11:IPVlife1
## -6.6556 -190.1380
## SexMen:RaceAfrAm:IPVlife1 CES11:RaceAfrAm:IPVlife1
## -123.1354 -29.9952
## Age:SexMen:CES11:RaceAfrAm Age:SexMen:CES11:IPVlife1
## -2.2713 -8.4119
## Age:SexMen:RaceAfrAm:IPVlife1 Age:CES11:RaceAfrAm:IPVlife1
## -22.2136 0.5845
## SexMen:CES11:RaceAfrAm:IPVlife1 Age:SexMen:CES11:RaceAfrAm:IPVlife1
## 270.5640 22.8519
(st = step(mm4))
## Random term (Age | HNDid) was eliminated because of having correlation +-1 or NaN
##
## 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) 2993 1 kept < 1e-07
##
## Fixed effects:
## Sum Sq Mean Sq NumDF DenDF F.value elim.num
## Age:Sex:CES1:Race:IPVlife 9.329e+03 9.329e+03 1 4059 1.5295 1
## Age:Sex:CES1:IPVlife 1.022e+03 1.022e+03 1 4080 0.1554 2
## Age:CES1:Race:IPVlife 9.293e+02 9.293e+02 1 3417 0.2541 3
## Age:CES1:IPVlife 9.579e+02 9.579e+02 1 3405 0.2753 4
## Sex:CES1:Race:IPVlife 4.277e+03 4.277e+03 1 1246 0.4278 5
## CES1:Race:IPVlife 5.088e+02 5.088e+02 1 1184 0.0844 6
## Age:Sex:CES1:Race 4.737e+03 4.737e+03 1 3196 0.6177 7
## Sex:CES1:Race 4.885e+00 4.885e+00 1 1173 0.0006 8
## Age:CES1:Race 4.094e+02 4.094e+02 1 3151 0.0700 9
## Age:Sex:CES1 3.215e+03 3.215e+03 1 3218 0.5042 10
## Age:CES1 3.908e+02 3.908e+02 1 3175 0.2412 11
## CES1:Race 4.530e+03 4.530e+03 1 1174 1.0324 12
## Age:Sex:Race:IPVlife 6.829e+03 6.829e+03 1 3918 1.3289 13
## Age:Sex:IPVlife 1.383e+03 1.383e+03 1 4132 0.0851 14
## Age:Race:IPVlife 1.834e+03 1.834e+03 1 3364 0.1468 15
## Age:IPVlife 2.603e+03 2.603e+03 1 3340 0.1527 16
## Age:Sex:Race 1.189e+04 1.189e+04 1 3197 1.6647 17
## Age:Sex 3.208e+03 3.208e+03 1 3215 0.5751 18
## Age:Race 5.644e+03 5.644e+03 1 3150 0.9967 19
## Sex:Race:IPVlife 8.718e+03 8.718e+03 1 1202 1.8003 20
## Race:IPVlife 6.723e+03 6.723e+03 1 1181 0.4630 21
## Sex:CES1:IPVlife 6.664e+03 6.664e+03 1 1206 1.5545 22
## Sex:CES1 7.876e+02 7.876e+02 1 1179 0.7187 23
## CES1:IPVlife 1.276e+04 1.276e+04 1 1183 1.7899 24
## Sex:Race 1.240e+04 1.240e+04 1 1179 2.1254 25
## Sex:IPVlife 1.806e+04 1.806e+04 1 1195 2.6562 26
## IPVlife 1.514e+03 1.514e+03 1 1186 0.2467 27
## Age 1.436e+05 1.436e+05 1 3186 28.9689 kept
## Sex 2.333e+04 2.333e+04 1 1182 5.2925 kept
## CES1 1.249e+05 1.249e+05 1 1185 23.2162 kept
## Race 4.008e+05 4.008e+05 1 1183 65.4065 kept
## Pr(>F)
## Age:Sex:CES1:Race:IPVlife 0.2163
## Age:Sex:CES1:IPVlife 0.6934
## Age:CES1:Race:IPVlife 0.6142
## Age:CES1:IPVlife 0.5998
## Sex:CES1:Race:IPVlife 0.5132
## CES1:Race:IPVlife 0.7714
## Age:Sex:CES1:Race 0.4320
## Sex:CES1:Race 0.9812
## Age:CES1:Race 0.7913
## Age:Sex:CES1 0.4777
## Age:CES1 0.6234
## CES1:Race 0.3098
## Age:Sex:Race:IPVlife 0.2491
## Age:Sex:IPVlife 0.7706
## Age:Race:IPVlife 0.7016
## Age:IPVlife 0.6960
## Age:Sex:Race 0.1971
## Age:Sex 0.4483
## Age:Race 0.3182
## Sex:Race:IPVlife 0.1799
## Race:IPVlife 0.4964
## Sex:CES1:IPVlife 0.2127
## Sex:CES1 0.3967
## CES1:IPVlife 0.1812
## Sex:Race 0.1451
## Sex:IPVlife 0.1034
## IPVlife 0.6195
## Age <1e-07
## Sex 0.0216
## CES1 0
## Race <1e-07
##
## Least squares means:
## Sex CES1 Race Estimate Standard Error DF t-value Lower CI
## Sex Women 2 NA NA 138.93 4.96 1182 28.0 129
## Sex Men 1 NA NA 156.54 5.95 1182 26.3 145
## CES1 0 NA 1 NA 129.10 4.87 1182 26.5 120
## CES1 1 NA 2 NA 166.38 6.08 1185 27.4 154
## Race White NA NA 2 117.18 5.63 1180 20.8 106
## Race AfrAm NA NA 1 178.29 5.25 1186 33.9 168
## Upper CI p-value
## Sex Women 149 <2e-16
## Sex Men 168 <2e-16
## CES1 0 139 <2e-16
## CES1 1 178 <2e-16
## Race White 128 <2e-16
## Race AfrAm 189 <2e-16
##
## Differences of LSMEANS:
## Estimate Standard Error DF t-value Lower CI Upper CI
## Sex Women-Men -17.6 7.65 1181.7 -2.30 -32.6 -2.59
## CES1 0-1 -37.3 7.74 1185.4 -4.82 -52.5 -22.10
## Race White-AfrAm -61.1 7.56 1183.3 -8.09 -75.9 -46.29
## p-value
## Sex Women-Men 0.02
## CES1 0-1 <2e-16
## Race White-AfrAm <2e-16
##
## Final model:
## lme4::lmer(formula = TrailsBtestSec ~ Age + Sex + CES1 + Race +
## (1 | HNDid), data = AllPsDATA, REML = reml, contrasts = l)
Re-run suggested model 4
(mm4 = lmer(TrailsBtestSec ~ Age + Sex + CES1 + Race + (Age | HNDid), data = AllPsDATA))
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: TrailsBtestSec ~ Age + Sex + CES1 + Race + (Age | HNDid)
## Data: AllPsDATA
## REML criterion at convergence: 58599
## Random effects:
## Groups Name Std.Dev. Corr
## HNDid (Intercept) 124.8
## Age 0.8 1.00
## Residual 78.0
## Number of obs: 4822, groups: HNDid, 1217
## Fixed Effects:
## (Intercept) Age SexMen CES11 RaceAfrAm
## 90.7 1.7 16.8 36.9 60.1
summary(mm4)
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: TrailsBtestSec ~ Age + Sex + CES1 + Race + (Age | HNDid)
## Data: AllPsDATA
##
## REML criterion at convergence: 58599
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 15580.49 124.8
## Age 0.64 0.8 1.00
## Residual 6087.38 78.0
## Number of obs: 4822, groups: HNDid, 1217
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 90.703 7.222 1173.000 12.56 < 2e-16
## Age 1.696 0.316 3113.000 5.36 8.8e-08
## SexMen 16.813 7.584 1154.000 2.22 0.027
## CES11 36.920 7.653 1146.000 4.82 1.6e-06
## RaceAfrAm 60.080 7.491 1158.000 8.02 2.7e-15
##
## Correlation of Fixed Effects:
## (Intr) Age SexMen CES11
## Age 0.024
## SexMen -0.473 0.000
## CES11 -0.478 0.051 0.085
## RaceAfrAm -0.575 0.045 -0.003 0.039