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 + TrailsAtestSec)^5 + (Age | HNDid) + (1 | subclass)
## Data: IPVandCognitionDataSet2
## REML criterion at convergence: 1350
## Random effects:
## Groups Name Std.Dev. Corr
## HNDid (Intercept) 1.91e+02
## Age 8.77e+00 1.00
## subclass (Intercept) 2.75e-04
## Residual 6.34e+01
## Number of obs: 126, groups: HNDid, 63; subclass, 21
## Fixed Effects:
## (Intercept)
## 118.539
## Age
## 12.029
## IPVstatus1
## -94.829
## SexMen
## -94.751
## PovStatBelow
## -10.837
## TrailsAtestSec
## 0.371
## Age:IPVstatus1
## -13.312
## Age:SexMen
## -15.241
## Age:PovStatBelow
## 0.967
## Age:TrailsAtestSec
## -0.418
## IPVstatus1:SexMen
## 201.471
## IPVstatus1:PovStatBelow
## -16.152
## IPVstatus1:TrailsAtestSec
## 1.127
## SexMen:PovStatBelow
## -86.081
## SexMen:TrailsAtestSec
## 1.332
## PovStatBelow:TrailsAtestSec
## 3.965
## Age:IPVstatus1:SexMen
## 41.174
## Age:IPVstatus1:PovStatBelow
## -5.980
## Age:IPVstatus1:TrailsAtestSec
## 0.459
## Age:SexMen:PovStatBelow
## -14.661
## Age:SexMen:TrailsAtestSec
## 0.536
## Age:PovStatBelow:TrailsAtestSec
## 0.371
## IPVstatus1:SexMen:PovStatBelow
## 35.670
## IPVstatus1:SexMen:TrailsAtestSec
## 0.258
## IPVstatus1:PovStatBelow:TrailsAtestSec
## -2.887
## SexMen:PovStatBelow:TrailsAtestSec
## 2.090
## Age:IPVstatus1:SexMen:PovStatBelow
## -18.324
## Age:IPVstatus1:SexMen:TrailsAtestSec
## -1.610
## Age:IPVstatus1:PovStatBelow:TrailsAtestSec
## -0.116
## Age:SexMen:PovStatBelow:TrailsAtestSec
## 0.268
## IPVstatus1:SexMen:PovStatBelow:TrailsAtestSec
## -1.712
## Age:IPVstatus1:SexMen:PovStatBelow:TrailsAtestSec
## 1.131
## 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 term (1 | subclass) was eliminated because of standard deviation being equal to 0
##
## Random effects:
## Chi.sq Chi.DF elim.num p.value
## (1 | HNDid) 35.71 1 kept < 1e-07
##
## Fixed effects:
## Sum Sq Mean Sq NumDF DenDF
## Age:IPVstatus:Sex:PovStat:TrailsAtestSec 3289.8 3289.8 1 71.26
## Age:IPVstatus:Sex:PovStat 1190.1 1190.1 1 90.86
## IPVstatus:Sex:PovStat:TrailsAtestSec 2171.8 2171.8 1 87.48
## Age:Sex:PovStat:TrailsAtestSec 1695.1 1695.1 1 69.09
## Age:Sex:PovStat 12535.4 12535.4 1 79.41
## Sex:PovStat:TrailsAtestSec 1973.7 1973.7 1 89.52
## IPVstatus:Sex:PovStat 2826.4 2826.4 1 57.01
## Sex:PovStat 244.5 244.5 1 62.25
## Age:IPVstatus:Sex:TrailsAtestSec 2455.1 2455.1 1 86.94
## IPVstatus:Sex:TrailsAtestSec 754.7 754.7 1 100.98
## Age:IPVstatus:Sex 213.5 213.5 1 102.24
## Age:IPVstatus:PovStat:TrailsAtestSec 8880.3 8880.3 1 71.14
## Age:IPVstatus:PovStat 137.2 137.2 1 102.15
## Age:PovStat:TrailsAtestSec 122.9 122.9 1 77.55
## Age:IPVstatus:TrailsAtestSec 5504.1 5504.1 1 77.19
## Age:IPVstatus 318.5 318.5 1 102.83
## IPVstatus:PovStat:TrailsAtestSec 8829.0 8829.0 1 88.72
## PovStat:TrailsAtestSec 9941.0 9941.0 1 83.46
## IPVstatus:PovStat 4923.2 4923.2 1 54.64
## IPVstatus:TrailsAtestSec 1897.6 1897.6 1 102.03
## Age 42572.1 42572.1 1 84.88
## IPVstatus 3895.0 3895.0 1 53.79
## Sex 14129.0 14129.0 1 110.53
## PovStat 11066.9 11066.9 1 92.31
## TrailsAtestSec 46335.9 46335.9 1 91.39
## Age:Sex 6893.8 6893.8 1 89.24
## Age:PovStat 43771.7 43771.7 1 112.29
## Age:TrailsAtestSec 207.5 207.5 1 90.39
## IPVstatus:Sex 29813.7 29813.7 1 53.28
## Sex:TrailsAtestSec 5359.5 5359.5 1 104.39
## Age:Sex:TrailsAtestSec 8093.8 8093.8 1 98.73
## F.value elim.num Pr(>F)
## Age:IPVstatus:Sex:PovStat:TrailsAtestSec 0.7208 1 0.3987
## Age:IPVstatus:Sex:PovStat 0.4122 2 0.5225
## IPVstatus:Sex:PovStat:TrailsAtestSec 0.2970 3 0.5871
## Age:Sex:PovStat:TrailsAtestSec 0.4594 4 0.5002
## Age:Sex:PovStat 0.0145 5 0.9043
## Sex:PovStat:TrailsAtestSec 0.1307 6 0.7186
## IPVstatus:Sex:PovStat 0.5283 7 0.4703
## Sex:PovStat 0.0810 8 0.7769
## Age:IPVstatus:Sex:TrailsAtestSec 1.6923 9 0.1967
## IPVstatus:Sex:TrailsAtestSec 0.0009 10 0.9758
## Age:IPVstatus:Sex 0.6850 11 0.4098
## Age:IPVstatus:PovStat:TrailsAtestSec 1.7470 12 0.1905
## Age:IPVstatus:PovStat 0.0092 13 0.9239
## Age:PovStat:TrailsAtestSec 0.0694 14 0.7930
## Age:IPVstatus:TrailsAtestSec 0.7633 15 0.3850
## Age:IPVstatus 0.3901 16 0.5336
## IPVstatus:PovStat:TrailsAtestSec 1.6780 17 0.1986
## PovStat:TrailsAtestSec 0.8934 18 0.3473
## IPVstatus:PovStat 0.9682 19 0.3295
## IPVstatus:TrailsAtestSec 1.9835 20 0.1621
## Age 0.2666 kept 0.6069
## IPVstatus 1.0214 kept 0.3167
## Sex 0.7159 kept 0.3993
## PovStat 2.4795 kept 0.1188
## TrailsAtestSec 4.6028 kept 0.0346
## Age:Sex 3.3443 kept 0.0708
## Age:PovStat 5.8519 kept 0.0172
## Age:TrailsAtestSec 0.1210 kept 0.7287
## IPVstatus:Sex 6.2175 kept 0.0158
## Sex:TrailsAtestSec 4.0258 kept 0.0474
## Age:Sex:TrailsAtestSec 5.3194 kept 0.0232
##
## Least squares means:
## IPVstatus Sex PovStat Estimate Standard Error DF
## IPVstatus 0 1.0 NA NA 116.5 22.9 58.0
## IPVstatus 1 2.0 NA NA 153.7 31.2 56.4
## Sex Women NA 2.0 NA 108.8 26.5 59.4
## Sex Men NA 1.0 NA 161.4 29.0 59.5
## PovStat Above NA NA 1.0 134.3 22.3 56.0
## PovStat Below NA NA 2.0 135.9 32.9 61.3
## IPVstatus:Sex 0 Women 1.0 2.0 NA 136.1 30.7 59.3
## IPVstatus:Sex 1 Women 2.0 2.0 NA 81.5 40.5 54.5
## IPVstatus:Sex 0 Men 1.0 1.0 NA 96.9 32.3 57.8
## IPVstatus:Sex 1 Men 2.0 1.0 NA 225.8 46.5 56.6
## t-value Lower CI Upper CI p-value
## IPVstatus 0 5.09 70.687 162 <2e-16
## IPVstatus 1 4.92 91.135 216 <2e-16
## Sex Women 4.10 55.763 162 0.0001
## Sex Men 5.56 103.284 219 <2e-16
## PovStat Above 6.02 89.618 179 <2e-16
## PovStat Below 4.13 70.043 202 0.0001
## IPVstatus:Sex 0 Women 4.43 74.678 198 <2e-16
## IPVstatus:Sex 1 Women 2.01 0.252 163 0.0493
## IPVstatus:Sex 0 Men 3.00 32.288 162 0.0040
## IPVstatus:Sex 1 Men 4.85 132.585 319 <2e-16
##
## Differences of LSMEANS:
## Estimate Standard Error DF t-value
## IPVstatus 0-1 -37.1 36.73 53.8 -1.01
## Sex Women-Men -52.6 37.99 59.1 -1.38
## PovStat Above-Below -1.6 38.93 59.4 -0.04
## IPVstatus:Sex 0 Women- 1 Women 54.6 48.57 52.6 1.12
## IPVstatus:Sex 0 Women- 0 Men 39.2 43.29 59.1 0.91
## IPVstatus:Sex 0 Women- 1 Men -89.7 55.15 56.7 -1.63
## IPVstatus:Sex 1 Women- 0 Men -15.5 50.43 56.2 -0.31
## IPVstatus:Sex 1 Women- 1 Men -144.3 61.00 54.8 -2.37
## IPVstatus:Sex 0 Men- 1 Men -128.9 55.20 54.3 -2.33
## Lower CI Upper CI p-value
## IPVstatus 0-1 -110.8 36.5 0.32
## Sex Women-Men -128.6 23.4 0.17
## PovStat Above-Below -79.5 76.2 0.97
## IPVstatus:Sex 0 Women- 1 Women -42.8 152.1 0.27
## IPVstatus:Sex 0 Women- 0 Men -47.4 125.8 0.37
## IPVstatus:Sex 0 Women- 1 Men -200.1 20.8 0.11
## IPVstatus:Sex 1 Women- 0 Men -116.5 85.6 0.76
## IPVstatus:Sex 1 Women- 1 Men -266.6 -22.1 0.02
## IPVstatus:Sex 0 Men- 1 Men -239.5 -18.2 0.02
##
## Final model:
## lme4::lmer(formula = TrailsBtestSec ~ Age + IPVstatus + Sex +
## PovStat + TrailsAtestSec + (1 | HNDid) + Age:Sex + Age:PovStat +
## Age:TrailsAtestSec + IPVstatus:Sex + Sex:TrailsAtestSec +
## Age:Sex:TrailsAtestSec, data = IPVandCognitionDataSet2, REML = reml,
## contrasts = l)
Re-run the suggested final Model 1
(mm1 = lmer(TrailsBtestSec ~ Age + IPVstatus + Sex + PovStat + TrailsAtestSec +
(Age | HNDid) + (1 | subclass) + Age:Sex + Age:PovStat + Age:TrailsAtestSec +
IPVstatus:Sex + Sex:TrailsAtestSec + Age:Sex:TrailsAtestSec, 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 + TrailsAtestSec + (Age | HNDid) + (1 | subclass) + Age:Sex + Age:PovStat + Age:TrailsAtestSec + IPVstatus:Sex + Sex:TrailsAtestSec + Age:Sex:TrailsAtestSec
## Data: IPVandCognitionDataSet2
## REML criterion at convergence: 1467
## Random effects:
## Groups Name Std.Dev. Corr
## HNDid (Intercept) 1.60e+02
## Age 4.46e+00 1.00
## subclass (Intercept) 3.83e-04
## Residual 6.19e+01
## Number of obs: 126, groups: HNDid, 63; subclass, 21
## Fixed Effects:
## (Intercept) Age
## 106.168 5.367
## IPVstatus1 SexMen
## -50.549 -203.514
## PovStatBelow TrailsAtestSec
## 91.116 0.546
## Age:SexMen Age:PovStatBelow
## -16.392 12.121
## Age:TrailsAtestSec IPVstatus1:SexMen
## -0.216 218.183
## SexMen:TrailsAtestSec Age:SexMen:TrailsAtestSec
## 4.651 0.502
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 + TrailsAtestSec + (Age | HNDid) + (1 | subclass) + Age:Sex + Age:PovStat + Age:TrailsAtestSec + IPVstatus:Sex + Sex:TrailsAtestSec + Age:Sex:TrailsAtestSec
## Data: IPVandCognitionDataSet2
##
## REML criterion at convergence: 1467
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 2.55e+04 1.60e+02
## Age 1.99e+01 4.46e+00 1.00
## subclass (Intercept) 1.47e-07 3.83e-04
## Residual 3.83e+03 6.19e+01
## Number of obs: 126, groups: HNDid, 63; subclass, 21
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 106.168 71.235 100.000 1.49 0.1393
## Age 5.367 5.500 82.100 0.98 0.3320
## IPVstatus1 -50.549 45.194 51.700 -1.12 0.2685
## SexMen -203.514 105.570 92.700 -1.93 0.0569
## PovStatBelow 91.116 57.941 47.000 1.57 0.1225
## TrailsAtestSec 0.546 1.573 81.300 0.35 0.7296
## Age:SexMen -16.392 8.573 84.400 -1.91 0.0593
## Age:PovStatBelow 12.121 4.418 100.800 2.74 0.0072
## Age:TrailsAtestSec -0.216 0.157 82.800 -1.37 0.1735
## IPVstatus1:SexMen 218.183 69.254 52.400 3.15 0.0027
## SexMen:TrailsAtestSec 4.651 2.372 82.200 1.96 0.0533
## Age:SexMen:TrailsAtestSec 0.502 0.233 88.900 2.16 0.0337
##
## Correlation of Fixed Effects:
## (Intr) Age IPVst1 SexMen PvSttB TrlsAS Ag:SxM Ag:PSB Ag:TAS
## Age 0.782
## IPVstatus1 -0.303 -0.066
## SexMen -0.667 -0.524 0.203
## PovStatBelw -0.132 -0.061 0.033 0.041
## TralsAtstSc -0.778 -0.691 0.041 0.531 -0.157
## Age:SexMen -0.497 -0.641 0.042 0.698 0.082 0.424
## Ag:PvSttBlw -0.030 -0.023 0.009 -0.011 0.776 -0.212 0.105
## Ag:TrlsAtsS -0.544 -0.842 -0.010 0.371 -0.139 0.724 0.515 -0.240
## IPVstts1:SM 0.210 0.048 -0.655 -0.196 -0.057 -0.030 -0.048 0.005 -0.001
## SxMn:TrlsAS 0.541 0.470 -0.033 -0.838 -0.124 -0.620 -0.601 -0.067 -0.438
## Ag:SxMn:TAS 0.377 0.577 0.004 -0.476 -0.147 -0.425 -0.884 -0.147 -0.603
## IPV1:S SM:TAS
## Age
## IPVstatus1
## SexMen
## PovStatBelw
## TralsAtstSc
## Age:SexMen
## Ag:PvSttBlw
## Ag:TrlsAtsS
## IPVstts1:SM
## SxMn:TrlsAS -0.036
## Ag:SxMn:TAS 0.046 0.542
plot(st)
plot(mm1)
