Clock Total 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: ClockTotal ~ (Age + IPVstatus + Sex + PovStat)^4 + (Age | HNDid) + (1 | subclass)
## Data: IPVandCognitionDataSet2
## REML criterion at convergence: 394.7
## Random effects:
## Groups Name Std.Dev. Corr
## HNDid (Intercept) 0.6582
## Age 0.0401 1.00
## subclass (Intercept) 0.5111
## Residual 0.9210
## Number of obs: 126, groups: HNDid, 63; subclass, 21
## Fixed Effects:
## (Intercept) Age
## 8.2669 -0.0190
## IPVstatus1 SexMen
## 0.6067 1.2351
## PovStatBelow Age:IPVstatus1
## 0.7636 -0.0165
## Age:SexMen Age:PovStatBelow
## 0.0429 0.0692
## IPVstatus1:SexMen IPVstatus1:PovStatBelow
## -1.7680 -2.3647
## SexMen:PovStatBelow Age:IPVstatus1:SexMen
## -1.7334 -0.0523
## Age:IPVstatus1:PovStatBelow Age:SexMen:PovStatBelow
## -0.2171 -0.2327
## IPVstatus1:SexMen:PovStatBelow Age:IPVstatus1:SexMen:PovStatBelow
## 4.2674 0.4593
## 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.0000
## (1 | subclass) 7.33 1 kept 0.0068
##
## Fixed effects:
## Sum Sq Mean Sq NumDF DenDF F.value elim.num
## Age 0.0170 0.0170 1 94.91 2.9501 kept
## IPVstatus 1.6095 1.6095 1 101.98 0.5630 kept
## Sex 2.0890 2.0890 1 107.34 1.8703 kept
## PovStat 0.3506 0.3506 1 108.14 0.7424 kept
## Age:IPVstatus 0.0115 0.0115 1 106.43 0.0622 kept
## Age:Sex 1.8739 1.8739 1 105.82 0.0000 kept
## Age:PovStat 0.0038 0.0038 1 109.98 2.1384 kept
## IPVstatus:Sex 2.9422 2.9422 1 106.56 0.1263 kept
## IPVstatus:PovStat 0.2890 0.2890 1 104.40 0.0043 kept
## Sex:PovStat 0.3687 0.3687 1 108.02 0.1028 kept
## Age:IPVstatus:Sex 2.5876 2.5876 1 109.50 5.0828 kept
## Age:IPVstatus:PovStat 0.5338 0.5338 1 107.75 0.5843 kept
## Age:Sex:PovStat 11.3140 11.3140 1 109.56 0.2061 kept
## IPVstatus:Sex:PovStat 0.6715 0.6715 1 109.98 8.9214 kept
## Age:IPVstatus:Sex:PovStat 10.8651 10.8651 1 108.15 11.1195 kept
## Pr(>F)
## Age 0.0891
## IPVstatus 0.4548
## Sex 0.1743
## PovStat 0.3908
## Age:IPVstatus 0.8036
## Age:Sex 0.9986
## Age:PovStat 0.1465
## IPVstatus:Sex 0.7230
## IPVstatus:PovStat 0.9478
## Sex:PovStat 0.7491
## Age:IPVstatus:Sex 0.0262
## Age:IPVstatus:PovStat 0.4463
## Age:Sex:PovStat 0.6508
## IPVstatus:Sex:PovStat 0.0035
## Age:IPVstatus:Sex:PovStat 0.0012
##
## Least squares means:
## IPVstatus Sex PovStat Estimate
## IPVstatus 0 1.0 NA NA 9.078
## IPVstatus 1 2.0 NA NA 8.876
## Sex Women NA 2.0 NA 8.700
## Sex Men NA 1.0 NA 9.252
## PovStat Above NA NA 1.0 8.913
## PovStat Below NA NA 2.0 9.040
## IPVstatus:Sex 0 Women 1.0 2.0 NA 8.521
## IPVstatus:Sex 1 Women 2.0 2.0 NA 8.879
## IPVstatus:Sex 0 Men 1.0 1.0 NA 9.633
## IPVstatus:Sex 1 Men 2.0 1.0 NA 8.872
## IPVstatus:PovStat 0 Above 1.0 NA 1.0 8.905
## IPVstatus:PovStat 1 Above 2.0 NA 1.0 8.920
## IPVstatus:PovStat 0 Below 1.0 NA 2.0 9.250
## IPVstatus:PovStat 1 Below 2.0 NA 2.0 8.831
## Sex:PovStat Women Above NA 2.0 1.0 8.768
## Sex:PovStat Men Above NA 1.0 1.0 9.058
## Sex:PovStat Women Below NA 2.0 2.0 8.633
## Sex:PovStat Men Below NA 1.0 2.0 9.447
## IPVstatus:Sex:PovStat 0 Women Above 1.0 2.0 1.0 8.400
## IPVstatus:Sex:PovStat 1 Women Above 2.0 2.0 1.0 9.135
## IPVstatus:Sex:PovStat 0 Men Above 1.0 1.0 1.0 9.411
## IPVstatus:Sex:PovStat 1 Men Above 2.0 1.0 1.0 8.705
## IPVstatus:Sex:PovStat 0 Women Below 1.0 2.0 2.0 8.643
## IPVstatus:Sex:PovStat 1 Women Below 2.0 2.0 2.0 8.624
## IPVstatus:Sex:PovStat 0 Men Below 1.0 1.0 2.0 9.856
## IPVstatus:Sex:PovStat 1 Men Below 2.0 1.0 2.0 9.038
## Standard Error DF t-value Lower CI
## IPVstatus 0 0.186 38.8 48.7 8.70
## IPVstatus 1 0.215 60.2 41.3 8.45
## Sex Women 0.193 38.5 45.1 8.31
## Sex Men 0.230 48.1 40.3 8.79
## PovStat Above 0.170 29.0 52.5 8.57
## PovStat Below 0.238 61.8 38.0 8.56
## IPVstatus:Sex 0 Women 0.206 47.6 41.3 8.11
## IPVstatus:Sex 1 Women 0.282 80.6 31.5 8.32
## IPVstatus:Sex 0 Men 0.285 72.7 33.8 9.07
## IPVstatus:Sex 1 Men 0.309 83.4 28.7 8.26
## IPVstatus:PovStat 0 Above 0.181 34.8 49.2 8.54
## IPVstatus:PovStat 1 Above 0.240 71.5 37.1 8.44
## IPVstatus:PovStat 0 Below 0.292 86.6 31.7 8.67
## IPVstatus:PovStat 1 Below 0.339 93.2 26.1 8.16
## Sex:PovStat Women Above 0.198 41.1 44.4 8.37
## Sex:PovStat Men Above 0.243 57.6 37.3 8.57
## Sex:PovStat Women Below 0.292 84.0 29.6 8.05
## Sex:PovStat Men Below 0.346 89.2 27.3 8.76
## IPVstatus:Sex:PovStat 0 Women Above 0.233 56.2 36.0 7.93
## IPVstatus:Sex:PovStat 1 Women Above 0.275 81.8 33.2 8.59
## IPVstatus:Sex:PovStat 0 Men Above 0.237 66.0 39.6 8.94
## IPVstatus:Sex:PovStat 1 Men Above 0.379 97.9 22.9 7.95
## IPVstatus:Sex:PovStat 0 Women Below 0.323 74.4 26.7 8.00
## IPVstatus:Sex:PovStat 1 Women Below 0.478 102.3 18.1 7.68
## IPVstatus:Sex:PovStat 0 Men Below 0.484 104.3 20.4 8.90
## IPVstatus:Sex:PovStat 1 Men Below 0.474 97.8 19.1 8.10
## Upper CI p-value
## IPVstatus 0 9.45 <2e-16
## IPVstatus 1 9.30 <2e-16
## Sex Women 9.09 <2e-16
## Sex Men 9.71 <2e-16
## PovStat Above 9.26 <2e-16
## PovStat Below 9.52 <2e-16
## IPVstatus:Sex 0 Women 8.94 <2e-16
## IPVstatus:Sex 1 Women 9.44 <2e-16
## IPVstatus:Sex 0 Men 10.20 <2e-16
## IPVstatus:Sex 1 Men 9.49 <2e-16
## IPVstatus:PovStat 0 Above 9.27 <2e-16
## IPVstatus:PovStat 1 Above 9.40 <2e-16
## IPVstatus:PovStat 0 Below 9.83 <2e-16
## IPVstatus:PovStat 1 Below 9.50 <2e-16
## Sex:PovStat Women Above 9.17 <2e-16
## Sex:PovStat Men Above 9.54 <2e-16
## Sex:PovStat Women Below 9.21 <2e-16
## Sex:PovStat Men Below 10.13 <2e-16
## IPVstatus:Sex:PovStat 0 Women Above 8.87 <2e-16
## IPVstatus:Sex:PovStat 1 Women Above 9.68 <2e-16
## IPVstatus:Sex:PovStat 0 Men Above 9.88 <2e-16
## IPVstatus:Sex:PovStat 1 Men Above 9.46 <2e-16
## IPVstatus:Sex:PovStat 0 Women Below 9.29 <2e-16
## IPVstatus:Sex:PovStat 1 Women Below 9.57 <2e-16
## IPVstatus:Sex:PovStat 0 Men Below 10.82 <2e-16
## IPVstatus:Sex:PovStat 1 Men Below 9.98 <2e-16
##
## Differences of LSMEANS:
## Estimate Standard Error DF
## IPVstatus 0-1 0.2 0.2335 99.8
## Sex Women-Men -0.6 0.2697 97.6
## PovStat Above-Below -0.1 0.2517 110.0
## IPVstatus:Sex 0 Women- 1 Women -0.4 0.3086 109.0
## IPVstatus:Sex 0 Women- 0 Men -1.1 0.3297 104.0
## IPVstatus:Sex 0 Women- 1 Men -0.4 0.3302 109.0
## IPVstatus:Sex 1 Women- 0 Men -0.8 0.3814 106.3
## IPVstatus:Sex 1 Women- 1 Men 0.0 0.4065 104.9
## IPVstatus:Sex 0 Men- 1 Men 0.8 0.3771 108.7
## IPVstatus:PovStat 0 Above- 1 Above 0.0 0.2569 107.0
## IPVstatus:PovStat 0 Above- 0 Below -0.3 0.3113 109.9
## IPVstatus:PovStat 0 Above- 1 Below 0.1 0.3486 109.4
## IPVstatus:PovStat 1 Above- 0 Below -0.3 0.3378 105.1
## IPVstatus:PovStat 1 Above- 1 Below 0.1 0.4010 109.1
## IPVstatus:PovStat 0 Below- 1 Below 0.4 0.4172 109.3
## Sex:PovStat Women Above- Men Above -0.3 0.2845 102.8
## Sex:PovStat Women Above- Women Below 0.1 0.3149 109.5
## Sex:PovStat Women Above- Men Below -0.7 0.3791 105.4
## Sex:PovStat Men Above- Women Below 0.4 0.3584 106.5
## Sex:PovStat Men Above- Men Below -0.4 0.3827 109.3
## Sex:PovStat Women Below- Men Below -0.8 0.4284 109.6
## t-value Lower CI Upper CI p-value
## IPVstatus 0-1 0.87 -0.2612 0.6652 0.389
## Sex Women-Men -2.05 -1.0873 -0.0168 0.043
## PovStat Above-Below -0.51 -0.6262 0.3713 0.614
## IPVstatus:Sex 0 Women- 1 Women -1.16 -0.9696 0.2539 0.249
## IPVstatus:Sex 0 Women- 0 Men -3.37 -1.7657 -0.4580 0.001
## IPVstatus:Sex 0 Women- 1 Men -1.06 -1.0045 0.3043 0.291
## IPVstatus:Sex 1 Women- 0 Men -1.98 -1.5102 0.0022 0.051
## IPVstatus:Sex 1 Women- 1 Men 0.02 -0.7983 0.8139 0.985
## IPVstatus:Sex 0 Men- 1 Men 2.02 0.0143 1.5093 0.046
## IPVstatus:PovStat 0 Above- 1 Above -0.06 -0.5242 0.4943 0.954
## IPVstatus:PovStat 0 Above- 0 Below -1.11 -0.9612 0.2725 0.271
## IPVstatus:PovStat 0 Above- 1 Below 0.21 -0.6165 0.7655 0.831
## IPVstatus:PovStat 1 Above- 0 Below -0.98 -0.9992 0.3404 0.332
## IPVstatus:PovStat 1 Above- 1 Below 0.22 -0.7053 0.8842 0.824
## IPVstatus:PovStat 0 Below- 1 Below 1.00 -0.4079 1.2456 0.318
## Sex:PovStat Women Above- Men Above -1.02 -0.8546 0.2739 0.310
## Sex:PovStat Women Above- Women Below 0.43 -0.4899 0.7584 0.671
## Sex:PovStat Women Above- Men Below -1.79 -1.4312 0.0723 0.076
## Sex:PovStat Men Above- Women Below 1.18 -0.2859 1.1351 0.239
## Sex:PovStat Men Above- Men Below -1.02 -1.1475 0.3693 0.311
## Sex:PovStat Women Below- Men Below -1.90 -1.6627 0.0352 0.060
##
## Final model:
## lme4::lmer(formula = ClockTotal ~ Age + IPVstatus + Sex + PovStat +
## (1 | subclass) + Age:IPVstatus + Age:Sex + Age:PovStat +
## IPVstatus:Sex + IPVstatus:PovStat + Sex:PovStat + Age:IPVstatus:Sex +
## Age:IPVstatus:PovStat + Age:Sex:PovStat + IPVstatus:Sex:PovStat +
## Age:IPVstatus:Sex:PovStat, data = IPVandCognitionDataSet2,
## REML = reml, contrasts = l)
Re-run suggested final Model 1
(mm1 = lmer(ClockTotal ~ Age + IPVstatus + Sex + PovStat + (Age | HNDid) + (1 |
subclass) + Age:IPVstatus + Age:Sex + Age:PovStat + IPVstatus:Sex + IPVstatus:PovStat +
Sex:PovStat + Age:IPVstatus:Sex + Age:IPVstatus:PovStat + Age:Sex:PovStat +
IPVstatus:Sex:PovStat + Age:IPVstatus:Sex:PovStat, data = IPVandCognitionDataSet2))
## Warning: number of observations <= rank(Z); variance-covariance matrix
## will be unidentifiable
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: ClockTotal ~ Age + IPVstatus + Sex + PovStat + (Age | HNDid) + (1 | subclass) + Age:IPVstatus + Age:Sex + Age:PovStat + IPVstatus:Sex + IPVstatus:PovStat + Sex:PovStat + Age:IPVstatus:Sex + Age:IPVstatus:PovStat + Age:Sex:PovStat + IPVstatus:Sex:PovStat + Age:IPVstatus:Sex:PovStat
## Data: IPVandCognitionDataSet2
## REML criterion at convergence: 394.7
## Random effects:
## Groups Name Std.Dev. Corr
## HNDid (Intercept) 0.6582
## Age 0.0401 1.00
## subclass (Intercept) 0.5111
## Residual 0.9210
## Number of obs: 126, groups: HNDid, 63; subclass, 21
## Fixed Effects:
## (Intercept) Age
## 8.2669 -0.0190
## IPVstatus1 SexMen
## 0.6067 1.2351
## PovStatBelow Age:IPVstatus1
## 0.7636 -0.0165
## Age:SexMen Age:PovStatBelow
## 0.0429 0.0692
## IPVstatus1:SexMen IPVstatus1:PovStatBelow
## -1.7680 -2.3647
## SexMen:PovStatBelow Age:IPVstatus1:SexMen
## -1.7334 -0.0523
## Age:IPVstatus1:PovStatBelow Age:SexMen:PovStatBelow
## -0.2171 -0.2327
## IPVstatus1:SexMen:PovStatBelow Age:IPVstatus1:SexMen:PovStatBelow
## 4.2674 0.4593
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: ClockTotal ~ Age + IPVstatus + Sex + PovStat + (Age | HNDid) + (1 | subclass) + Age:IPVstatus + Age:Sex + Age:PovStat + IPVstatus:Sex + IPVstatus:PovStat + Sex:PovStat + Age:IPVstatus:Sex + Age:IPVstatus:PovStat + Age:Sex:PovStat + IPVstatus:Sex:PovStat + Age:IPVstatus:Sex:PovStat
## Data: IPVandCognitionDataSet2
##
## REML criterion at convergence: 394.7
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 0.43322 0.6582
## Age 0.00161 0.0401 1.00
## subclass (Intercept) 0.26121 0.5111
## Residual 0.84823 0.9210
## Number of obs: 126, groups: HNDid, 63; subclass, 21
##
## Fixed effects:
## Estimate Std. Error df t value
## (Intercept) 8.2669 0.4217 21.9000 19.60
## Age -0.0190 0.0333 33.2000 -0.57
## IPVstatus1 0.6067 0.6474 21.2000 0.94
## SexMen 1.2351 0.6299 20.2000 1.96
## PovStatBelow 0.7636 0.7021 11.2000 1.09
## Age:IPVstatus1 -0.0165 0.0597 40.6000 -0.28
## Age:SexMen 0.0429 0.0530 29.6000 0.81
## Age:PovStatBelow 0.0692 0.0568 12.1000 1.22
## IPVstatus1:SexMen -1.7680 1.0747 25.8000 -1.65
## IPVstatus1:PovStatBelow -2.3647 1.4873 35.1000 -1.59
## SexMen:PovStatBelow -1.7334 0.9367 11.9000 -1.85
## Age:IPVstatus1:SexMen -0.0523 0.0978 49.0000 -0.54
## Age:IPVstatus1:PovStatBelow -0.2171 0.1454 67.0000 -1.49
## Age:SexMen:PovStatBelow -0.2327 0.1001 15.8000 -2.32
## IPVstatus1:SexMen:PovStatBelow 4.2674 1.8483 27.0000 2.31
## Age:IPVstatus1:SexMen:PovStatBelow 0.4593 0.1868 46.5000 2.46
## Pr(>|t|)
## (Intercept) 2.2e-15
## Age 0.573
## IPVstatus1 0.359
## SexMen 0.064
## PovStatBelow 0.300
## Age:IPVstatus1 0.784
## Age:SexMen 0.425
## Age:PovStatBelow 0.246
## IPVstatus1:SexMen 0.112
## IPVstatus1:PovStatBelow 0.121
## SexMen:PovStatBelow 0.089
## Age:IPVstatus1:SexMen 0.595
## Age:IPVstatus1:PovStatBelow 0.140
## Age:SexMen:PovStatBelow 0.034
## IPVstatus1:SexMen:PovStatBelow 0.029
## Age:IPVstatus1:SexMen:PovStatBelow 0.018
##
## Correlation of Fixed Effects:
## (Intr) Age IPVst1 SexMen PvSttB Ag:IPV1 Ag:SxM Ag:PSB
## Age 0.828
## IPVstatus1 -0.593 -0.507
## SexMen -0.624 -0.523 0.387
## PovStatBelw -0.590 -0.506 0.402 0.384
## Ag:IPVstts1 -0.443 -0.503 0.853 0.280 0.306
## Age:SexMen -0.497 -0.572 0.304 0.880 0.293 0.301
## Ag:PvSttBlw -0.482 -0.575 0.331 0.302 0.823 0.333 0.325
## IPVstts1:SM 0.351 0.309 -0.602 -0.529 -0.233 -0.513 -0.473 -0.202
## IPVstt1:PSB 0.329 0.264 -0.501 -0.219 -0.523 -0.414 -0.154 -0.419
## SxMn:PvSttB 0.440 0.369 -0.268 -0.671 -0.735 -0.192 -0.582 -0.590
## Ag:IPVs1:SM 0.271 0.320 -0.520 -0.448 -0.178 -0.620 -0.518 -0.209
## Ag:IPV1:PSB 0.218 0.240 -0.387 -0.133 -0.355 -0.445 -0.129 -0.418
## Ag:SxMn:PSB 0.264 0.300 -0.142 -0.485 -0.424 -0.144 -0.537 -0.516
## IPV1:SM:PSB -0.259 -0.209 0.383 0.351 0.417 0.313 0.290 0.327
## A:IPV1:SM:P -0.166 -0.189 0.277 0.263 0.252 0.329 0.281 0.332
## IPVs1:SM IPV1:P SM:PSB Ag:IPV1:SM A:IPV1:P A:SM:P IPV1:SM:
## Age
## IPVstatus1
## SexMen
## PovStatBelw
## Ag:IPVstts1
## Age:SexMen
## Ag:PvSttBlw
## IPVstts1:SM
## IPVstt1:PSB 0.297
## SxMn:PvSttB 0.366 0.371
## Ag:IPVs1:SM 0.866 0.246 0.304
## Ag:IPV1:PSB 0.231 0.879 0.240 0.269
## Ag:SxMn:PSB 0.231 0.201 0.654 0.262 0.197
## IPV1:SM:PSB -0.610 -0.795 -0.527 -0.515 -0.695 -0.315
## A:IPV1:SM:P -0.454 -0.662 -0.354 -0.521 -0.770 -0.539 0.776
plot(st)
plot(mm1)
Clock Total Regression Model 2 (with CES)
load("~/Desktop/Megan/Research/IPV and Cognition Paper/IPV R Output/IPVandCognitionDataSet2.rda")
library(lme4)
library(lmerTest)
(mm2 = lmer(ClockTotal ~ (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: ClockTotal ~ (Age + IPVstatus + Sex + PovStat + CES1)^5 + (Age | HNDid) + (1 | subclass)
## Data: IPVandCognitionDataSet2
## REML criterion at convergence: 387.4
## Random effects:
## Groups Name Std.Dev. Corr
## HNDid (Intercept) 0.4395
## Age 0.0252 1.00
## subclass (Intercept) 0.6098
## Residual 0.9508
## Number of obs: 126, groups: HNDid, 63; subclass, 21
## Fixed Effects:
## (Intercept)
## 8.2464
## Age
## -0.0294
## IPVstatus1
## -0.8196
## SexMen
## 1.5584
## PovStatBelow
## 0.5721
## CES11
## 0.2792
## Age:IPVstatus1
## -0.1344
## Age:SexMen
## 0.0765
## Age:PovStatBelow
## 0.0723
## Age:CES11
## 0.0421
## IPVstatus1:SexMen
## -0.0143
## IPVstatus1:PovStatBelow
## -3.1305
## IPVstatus1:CES11
## 1.7560
## SexMen:PovStatBelow
## -1.6411
## SexMen:CES11
## -1.7681
## PovStatBelow:CES11
## 0.0553
## Age:IPVstatus1:SexMen
## 0.0873
## Age:IPVstatus1:PovStatBelow
## -0.3085
## Age:IPVstatus1:CES11
## 0.1348
## Age:SexMen:PovStatBelow
## -0.3267
## Age:SexMen:CES11
## -0.1572
## Age:PovStatBelow:CES11
## -0.0308
## IPVstatus1:SexMen:PovStatBelow
## 4.9071
## IPVstatus1:SexMen:CES11
## -1.0320
## IPVstatus1:PovStatBelow:CES11
## 1.0201
## SexMen:PovStatBelow:CES11
## 1.5830
## Age:IPVstatus1:SexMen:PovStatBelow
## 0.6091
## Age:IPVstatus1:SexMen:CES11
## -0.0886
## Age:IPVstatus1:PovStatBelow:CES11
## 0.1546
## Age:SexMen:PovStatBelow:CES11
## 0.5336
## IPVstatus1:SexMen:PovStatBelow:CES11
## -2.6328
## Age:IPVstatus1:SexMen:PovStatBelow:CES11
## -0.4914
(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.0000
## (1 | subclass) 10.24 1 kept 0.0014
##
## Fixed effects:
## Sum Sq Mean Sq NumDF DenDF F.value
## Age:IPVstatus:Sex:PovStat:CES1 0.7407 0.7407 1 83.38 0.7770
## Age:IPVstatus:PovStat:CES1 0.0625 0.0625 1 84.92 0.0119
## IPVstatus:Sex:PovStat:CES1 0.3459 0.3459 1 91.81 0.3769
## IPVstatus:PovStat:CES1 0.5225 0.5225 1 96.98 0.0034
## Age:IPVstatus:Sex:CES1 1.1427 1.1427 1 97.61 1.9237
## Age:IPVstatus:CES1 0.5360 0.5360 1 98.68 0.3752
## IPVstatus:Sex:CES1 0.4127 0.4127 1 99.54 0.4045
## IPVstatus:CES1 0.0296 0.0296 1 99.45 0.1986
## Age:Sex:PovStat:CES1 2.9045 2.9045 1 97.37 3.8139
## Age:PovStat:CES1 1.1581 1.1581 1 97.66 0.6533
## Sex:PovStat:CES1 0.0555 0.0555 1 103.16 0.3543
## PovStat:CES1 0.1013 0.1013 1 103.16 0.0129
## Age:Sex:CES1 2.9691 2.9691 1 105.76 1.8343
## Age:CES1 3.4250 3.4250 1 105.47 0.5689
## Sex:CES1 4.3330 4.3330 1 106.51 3.1836
## CES1 0.3324 0.3324 1 107.93 0.3860
## Age 0.0000 0.0000 1 94.91 2.9501
## IPVstatus 1.5929 1.5929 1 101.98 0.5630
## Sex 1.8193 1.8193 1 107.34 1.8703
## PovStat 0.2899 0.2899 1 108.14 0.7424
## Age:IPVstatus 0.0390 0.0390 1 106.43 0.0622
## Age:Sex 2.2565 2.2565 1 105.82 0.0000
## Age:PovStat 0.0002 0.0002 1 109.98 2.1384
## IPVstatus:Sex 3.1134 3.1134 1 106.56 0.1263
## IPVstatus:PovStat 0.3558 0.3558 1 104.40 0.0043
## Sex:PovStat 0.2897 0.2897 1 108.02 0.1028
## Age:IPVstatus:Sex 3.8198 3.8198 1 109.50 5.0828
## Age:IPVstatus:PovStat 0.8921 0.8921 1 107.75 0.5843
## Age:Sex:PovStat 6.8599 6.8599 1 109.56 0.2061
## IPVstatus:Sex:PovStat 0.0059 0.0059 1 109.98 8.9214
## Age:IPVstatus:Sex:PovStat 7.9306 7.9306 1 108.15 11.1195
## elim.num Pr(>F)
## Age:IPVstatus:Sex:PovStat:CES1 1 0.3806
## Age:IPVstatus:PovStat:CES1 2 0.9135
## IPVstatus:Sex:PovStat:CES1 3 0.5408
## IPVstatus:PovStat:CES1 4 0.9539
## Age:IPVstatus:Sex:CES1 5 0.1686
## Age:IPVstatus:CES1 6 0.5416
## IPVstatus:Sex:CES1 7 0.5262
## IPVstatus:CES1 8 0.6568
## Age:Sex:PovStat:CES1 9 0.0537
## Age:PovStat:CES1 10 0.4209
## Sex:PovStat:CES1 11 0.5530
## PovStat:CES1 12 0.9097
## Age:Sex:CES1 13 0.1785
## Age:CES1 14 0.4524
## Sex:CES1 15 0.0772
## CES1 16 0.5357
## Age kept 0.0891
## IPVstatus kept 0.4548
## Sex kept 0.1743
## PovStat kept 0.3908
## Age:IPVstatus kept 0.8036
## Age:Sex kept 0.9986
## Age:PovStat kept 0.1465
## IPVstatus:Sex kept 0.7230
## IPVstatus:PovStat kept 0.9478
## Sex:PovStat kept 0.7491
## Age:IPVstatus:Sex kept 0.0262
## Age:IPVstatus:PovStat kept 0.4463
## Age:Sex:PovStat kept 0.6508
## IPVstatus:Sex:PovStat kept 0.0035
## Age:IPVstatus:Sex:PovStat kept 0.0012
##
## Least squares means:
## IPVstatus Sex PovStat Estimate
## IPVstatus 0 1.0 NA NA 9.078
## IPVstatus 1 2.0 NA NA 8.876
## Sex Women NA 2.0 NA 8.700
## Sex Men NA 1.0 NA 9.252
## PovStat Above NA NA 1.0 8.913
## PovStat Below NA NA 2.0 9.040
## IPVstatus:Sex 0 Women 1.0 2.0 NA 8.521
## IPVstatus:Sex 1 Women 2.0 2.0 NA 8.879
## IPVstatus:Sex 0 Men 1.0 1.0 NA 9.633
## IPVstatus:Sex 1 Men 2.0 1.0 NA 8.872
## IPVstatus:PovStat 0 Above 1.0 NA 1.0 8.905
## IPVstatus:PovStat 1 Above 2.0 NA 1.0 8.920
## IPVstatus:PovStat 0 Below 1.0 NA 2.0 9.250
## IPVstatus:PovStat 1 Below 2.0 NA 2.0 8.831
## Sex:PovStat Women Above NA 2.0 1.0 8.768
## Sex:PovStat Men Above NA 1.0 1.0 9.058
## Sex:PovStat Women Below NA 2.0 2.0 8.633
## Sex:PovStat Men Below NA 1.0 2.0 9.447
## IPVstatus:Sex:PovStat 0 Women Above 1.0 2.0 1.0 8.400
## IPVstatus:Sex:PovStat 1 Women Above 2.0 2.0 1.0 9.135
## IPVstatus:Sex:PovStat 0 Men Above 1.0 1.0 1.0 9.411
## IPVstatus:Sex:PovStat 1 Men Above 2.0 1.0 1.0 8.705
## IPVstatus:Sex:PovStat 0 Women Below 1.0 2.0 2.0 8.643
## IPVstatus:Sex:PovStat 1 Women Below 2.0 2.0 2.0 8.624
## IPVstatus:Sex:PovStat 0 Men Below 1.0 1.0 2.0 9.856
## IPVstatus:Sex:PovStat 1 Men Below 2.0 1.0 2.0 9.038
## Standard Error DF t-value Lower CI
## IPVstatus 0 0.186 38.8 48.7 8.70
## IPVstatus 1 0.215 60.2 41.3 8.45
## Sex Women 0.193 38.5 45.1 8.31
## Sex Men 0.230 48.1 40.3 8.79
## PovStat Above 0.170 29.0 52.5 8.57
## PovStat Below 0.238 61.8 38.0 8.56
## IPVstatus:Sex 0 Women 0.206 47.6 41.3 8.11
## IPVstatus:Sex 1 Women 0.282 80.6 31.5 8.32
## IPVstatus:Sex 0 Men 0.285 72.7 33.8 9.07
## IPVstatus:Sex 1 Men 0.309 83.4 28.7 8.26
## IPVstatus:PovStat 0 Above 0.181 34.8 49.2 8.54
## IPVstatus:PovStat 1 Above 0.240 71.5 37.1 8.44
## IPVstatus:PovStat 0 Below 0.292 86.6 31.7 8.67
## IPVstatus:PovStat 1 Below 0.339 93.2 26.1 8.16
## Sex:PovStat Women Above 0.198 41.1 44.4 8.37
## Sex:PovStat Men Above 0.243 57.6 37.3 8.57
## Sex:PovStat Women Below 0.292 84.0 29.6 8.05
## Sex:PovStat Men Below 0.346 89.2 27.3 8.76
## IPVstatus:Sex:PovStat 0 Women Above 0.233 56.2 36.0 7.93
## IPVstatus:Sex:PovStat 1 Women Above 0.275 81.8 33.2 8.59
## IPVstatus:Sex:PovStat 0 Men Above 0.237 66.0 39.6 8.94
## IPVstatus:Sex:PovStat 1 Men Above 0.379 97.9 22.9 7.95
## IPVstatus:Sex:PovStat 0 Women Below 0.323 74.4 26.7 8.00
## IPVstatus:Sex:PovStat 1 Women Below 0.478 102.3 18.1 7.68
## IPVstatus:Sex:PovStat 0 Men Below 0.484 104.3 20.4 8.90
## IPVstatus:Sex:PovStat 1 Men Below 0.474 97.8 19.1 8.10
## Upper CI p-value
## IPVstatus 0 9.45 <2e-16
## IPVstatus 1 9.30 <2e-16
## Sex Women 9.09 <2e-16
## Sex Men 9.71 <2e-16
## PovStat Above 9.26 <2e-16
## PovStat Below 9.52 <2e-16
## IPVstatus:Sex 0 Women 8.94 <2e-16
## IPVstatus:Sex 1 Women 9.44 <2e-16
## IPVstatus:Sex 0 Men 10.20 <2e-16
## IPVstatus:Sex 1 Men 9.49 <2e-16
## IPVstatus:PovStat 0 Above 9.27 <2e-16
## IPVstatus:PovStat 1 Above 9.40 <2e-16
## IPVstatus:PovStat 0 Below 9.83 <2e-16
## IPVstatus:PovStat 1 Below 9.50 <2e-16
## Sex:PovStat Women Above 9.17 <2e-16
## Sex:PovStat Men Above 9.54 <2e-16
## Sex:PovStat Women Below 9.21 <2e-16
## Sex:PovStat Men Below 10.13 <2e-16
## IPVstatus:Sex:PovStat 0 Women Above 8.87 <2e-16
## IPVstatus:Sex:PovStat 1 Women Above 9.68 <2e-16
## IPVstatus:Sex:PovStat 0 Men Above 9.88 <2e-16
## IPVstatus:Sex:PovStat 1 Men Above 9.46 <2e-16
## IPVstatus:Sex:PovStat 0 Women Below 9.29 <2e-16
## IPVstatus:Sex:PovStat 1 Women Below 9.57 <2e-16
## IPVstatus:Sex:PovStat 0 Men Below 10.82 <2e-16
## IPVstatus:Sex:PovStat 1 Men Below 9.98 <2e-16
##
## Differences of LSMEANS:
## Estimate Standard Error DF
## IPVstatus 0-1 0.2 0.2335 99.8
## Sex Women-Men -0.6 0.2697 97.6
## PovStat Above-Below -0.1 0.2517 110.0
## IPVstatus:Sex 0 Women- 1 Women -0.4 0.3086 109.0
## IPVstatus:Sex 0 Women- 0 Men -1.1 0.3297 104.0
## IPVstatus:Sex 0 Women- 1 Men -0.4 0.3302 109.0
## IPVstatus:Sex 1 Women- 0 Men -0.8 0.3814 106.3
## IPVstatus:Sex 1 Women- 1 Men 0.0 0.4065 104.9
## IPVstatus:Sex 0 Men- 1 Men 0.8 0.3771 108.7
## IPVstatus:PovStat 0 Above- 1 Above 0.0 0.2569 107.0
## IPVstatus:PovStat 0 Above- 0 Below -0.3 0.3113 109.9
## IPVstatus:PovStat 0 Above- 1 Below 0.1 0.3486 109.4
## IPVstatus:PovStat 1 Above- 0 Below -0.3 0.3378 105.1
## IPVstatus:PovStat 1 Above- 1 Below 0.1 0.4010 109.1
## IPVstatus:PovStat 0 Below- 1 Below 0.4 0.4172 109.3
## Sex:PovStat Women Above- Men Above -0.3 0.2845 102.8
## Sex:PovStat Women Above- Women Below 0.1 0.3149 109.5
## Sex:PovStat Women Above- Men Below -0.7 0.3791 105.4
## Sex:PovStat Men Above- Women Below 0.4 0.3584 106.5
## Sex:PovStat Men Above- Men Below -0.4 0.3827 109.3
## Sex:PovStat Women Below- Men Below -0.8 0.4284 109.6
## t-value Lower CI Upper CI p-value
## IPVstatus 0-1 0.87 -0.2612 0.6652 0.389
## Sex Women-Men -2.05 -1.0873 -0.0168 0.043
## PovStat Above-Below -0.51 -0.6262 0.3713 0.614
## IPVstatus:Sex 0 Women- 1 Women -1.16 -0.9696 0.2539 0.249
## IPVstatus:Sex 0 Women- 0 Men -3.37 -1.7657 -0.4580 0.001
## IPVstatus:Sex 0 Women- 1 Men -1.06 -1.0045 0.3043 0.291
## IPVstatus:Sex 1 Women- 0 Men -1.98 -1.5102 0.0022 0.051
## IPVstatus:Sex 1 Women- 1 Men 0.02 -0.7983 0.8139 0.985
## IPVstatus:Sex 0 Men- 1 Men 2.02 0.0143 1.5093 0.046
## IPVstatus:PovStat 0 Above- 1 Above -0.06 -0.5242 0.4943 0.954
## IPVstatus:PovStat 0 Above- 0 Below -1.11 -0.9612 0.2725 0.271
## IPVstatus:PovStat 0 Above- 1 Below 0.21 -0.6165 0.7655 0.831
## IPVstatus:PovStat 1 Above- 0 Below -0.98 -0.9992 0.3404 0.332
## IPVstatus:PovStat 1 Above- 1 Below 0.22 -0.7053 0.8842 0.824
## IPVstatus:PovStat 0 Below- 1 Below 1.00 -0.4079 1.2456 0.318
## Sex:PovStat Women Above- Men Above -1.02 -0.8546 0.2739 0.310
## Sex:PovStat Women Above- Women Below 0.43 -0.4899 0.7584 0.671
## Sex:PovStat Women Above- Men Below -1.79 -1.4312 0.0723 0.076
## Sex:PovStat Men Above- Women Below 1.18 -0.2859 1.1351 0.239
## Sex:PovStat Men Above- Men Below -1.02 -1.1475 0.3693 0.311
## Sex:PovStat Women Below- Men Below -1.90 -1.6627 0.0352 0.060
##
## Final model:
## lme4::lmer(formula = ClockTotal ~ Age + IPVstatus + Sex + PovStat +
## (1 | subclass) + Age:IPVstatus + Age:Sex + Age:PovStat +
## IPVstatus:Sex + IPVstatus:PovStat + Sex:PovStat + Age:IPVstatus:Sex +
## Age:IPVstatus:PovStat + Age:Sex:PovStat + IPVstatus:Sex:PovStat +
## Age:IPVstatus:Sex:PovStat, data = IPVandCognitionDataSet2,
## REML = reml, contrasts = l)
Re-run suggested final Model 2
(mm2 = lmer(ClockTotal ~ Age + IPVstatus + Sex + PovStat + (Age | HNDid) + (1 |
subclass) + Age:IPVstatus + Age:Sex + Age:PovStat + IPVstatus:Sex + IPVstatus:PovStat +
Sex:PovStat + Age:IPVstatus:Sex + Age:IPVstatus:PovStat + Age:Sex:PovStat +
IPVstatus:Sex:PovStat + Age:IPVstatus:Sex:PovStat, data = IPVandCognitionDataSet2))
## Warning: number of observations <= rank(Z); variance-covariance matrix
## will be unidentifiable
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: ClockTotal ~ Age + IPVstatus + Sex + PovStat + (Age | HNDid) + (1 | subclass) + Age:IPVstatus + Age:Sex + Age:PovStat + IPVstatus:Sex + IPVstatus:PovStat + Sex:PovStat + Age:IPVstatus:Sex + Age:IPVstatus:PovStat + Age:Sex:PovStat + IPVstatus:Sex:PovStat + Age:IPVstatus:Sex:PovStat
## Data: IPVandCognitionDataSet2
## REML criterion at convergence: 394.7
## Random effects:
## Groups Name Std.Dev. Corr
## HNDid (Intercept) 0.6582
## Age 0.0401 1.00
## subclass (Intercept) 0.5111
## Residual 0.9210
## Number of obs: 126, groups: HNDid, 63; subclass, 21
## Fixed Effects:
## (Intercept) Age
## 8.2669 -0.0190
## IPVstatus1 SexMen
## 0.6067 1.2351
## PovStatBelow Age:IPVstatus1
## 0.7636 -0.0165
## Age:SexMen Age:PovStatBelow
## 0.0429 0.0692
## IPVstatus1:SexMen IPVstatus1:PovStatBelow
## -1.7680 -2.3647
## SexMen:PovStatBelow Age:IPVstatus1:SexMen
## -1.7334 -0.0523
## Age:IPVstatus1:PovStatBelow Age:SexMen:PovStatBelow
## -0.2171 -0.2327
## IPVstatus1:SexMen:PovStatBelow Age:IPVstatus1:SexMen:PovStatBelow
## 4.2674 0.4593
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: ClockTotal ~ Age + IPVstatus + Sex + PovStat + (Age | HNDid) + (1 | subclass) + Age:IPVstatus + Age:Sex + Age:PovStat + IPVstatus:Sex + IPVstatus:PovStat + Sex:PovStat + Age:IPVstatus:Sex + Age:IPVstatus:PovStat + Age:Sex:PovStat + IPVstatus:Sex:PovStat + Age:IPVstatus:Sex:PovStat
## Data: IPVandCognitionDataSet2
##
## REML criterion at convergence: 394.7
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 0.43322 0.6582
## Age 0.00161 0.0401 1.00
## subclass (Intercept) 0.26121 0.5111
## Residual 0.84823 0.9210
## Number of obs: 126, groups: HNDid, 63; subclass, 21
##
## Fixed effects:
## Estimate Std. Error df t value
## (Intercept) 8.2669 0.4217 21.9000 19.60
## Age -0.0190 0.0333 33.2000 -0.57
## IPVstatus1 0.6067 0.6474 21.2000 0.94
## SexMen 1.2351 0.6299 20.2000 1.96
## PovStatBelow 0.7636 0.7021 11.2000 1.09
## Age:IPVstatus1 -0.0165 0.0597 40.6000 -0.28
## Age:SexMen 0.0429 0.0530 29.6000 0.81
## Age:PovStatBelow 0.0692 0.0568 12.1000 1.22
## IPVstatus1:SexMen -1.7680 1.0747 25.8000 -1.65
## IPVstatus1:PovStatBelow -2.3647 1.4873 35.1000 -1.59
## SexMen:PovStatBelow -1.7334 0.9367 11.9000 -1.85
## Age:IPVstatus1:SexMen -0.0523 0.0978 49.0000 -0.54
## Age:IPVstatus1:PovStatBelow -0.2171 0.1454 67.0000 -1.49
## Age:SexMen:PovStatBelow -0.2327 0.1001 15.8000 -2.32
## IPVstatus1:SexMen:PovStatBelow 4.2674 1.8483 27.0000 2.31
## Age:IPVstatus1:SexMen:PovStatBelow 0.4593 0.1868 46.5000 2.46
## Pr(>|t|)
## (Intercept) 2.2e-15
## Age 0.573
## IPVstatus1 0.359
## SexMen 0.064
## PovStatBelow 0.300
## Age:IPVstatus1 0.784
## Age:SexMen 0.425
## Age:PovStatBelow 0.246
## IPVstatus1:SexMen 0.112
## IPVstatus1:PovStatBelow 0.121
## SexMen:PovStatBelow 0.089
## Age:IPVstatus1:SexMen 0.595
## Age:IPVstatus1:PovStatBelow 0.140
## Age:SexMen:PovStatBelow 0.034
## IPVstatus1:SexMen:PovStatBelow 0.029
## Age:IPVstatus1:SexMen:PovStatBelow 0.018
##
## Correlation of Fixed Effects:
## (Intr) Age IPVst1 SexMen PvSttB Ag:IPV1 Ag:SxM Ag:PSB
## Age 0.828
## IPVstatus1 -0.593 -0.507
## SexMen -0.624 -0.523 0.387
## PovStatBelw -0.590 -0.506 0.402 0.384
## Ag:IPVstts1 -0.443 -0.503 0.853 0.280 0.306
## Age:SexMen -0.497 -0.572 0.304 0.880 0.293 0.301
## Ag:PvSttBlw -0.482 -0.575 0.331 0.302 0.823 0.333 0.325
## IPVstts1:SM 0.351 0.309 -0.602 -0.529 -0.233 -0.513 -0.473 -0.202
## IPVstt1:PSB 0.329 0.264 -0.501 -0.219 -0.523 -0.414 -0.154 -0.419
## SxMn:PvSttB 0.440 0.369 -0.268 -0.671 -0.735 -0.192 -0.582 -0.590
## Ag:IPVs1:SM 0.271 0.320 -0.520 -0.448 -0.178 -0.620 -0.518 -0.209
## Ag:IPV1:PSB 0.218 0.240 -0.387 -0.133 -0.355 -0.445 -0.129 -0.418
## Ag:SxMn:PSB 0.264 0.300 -0.142 -0.485 -0.424 -0.144 -0.537 -0.516
## IPV1:SM:PSB -0.259 -0.209 0.383 0.351 0.417 0.313 0.290 0.327
## A:IPV1:SM:P -0.166 -0.189 0.277 0.263 0.252 0.329 0.281 0.332
## IPVs1:SM IPV1:P SM:PSB Ag:IPV1:SM A:IPV1:P A:SM:P IPV1:SM:
## Age
## IPVstatus1
## SexMen
## PovStatBelw
## Ag:IPVstts1
## Age:SexMen
## Ag:PvSttBlw
## IPVstts1:SM
## IPVstt1:PSB 0.297
## SxMn:PvSttB 0.366 0.371
## Ag:IPVs1:SM 0.866 0.246 0.304
## Ag:IPV1:PSB 0.231 0.879 0.240 0.269
## Ag:SxMn:PSB 0.231 0.201 0.654 0.262 0.197
## IPV1:SM:PSB -0.610 -0.795 -0.527 -0.515 -0.695 -0.315
## A:IPV1:SM:P -0.454 -0.662 -0.354 -0.521 -0.770 -0.539 0.776
plot(st)
plot(mm2)