Trails A Regression Model 1
load("~/Desktop/Research/IPV and Cognition Paper/R Commands and Scripts/LifeIPVandCogData.rda")
library(lme4)
library(lmerTest)
library(compute.es)
(mm1 = lmer(logTrailsA~(Age + LifeIPV + Sex + PovStat)^4 + (Age|HNDid) + (1|subclass) , data=LifeIPVandCogData))
(st = step(mm1))
Trails A Final Regression Model 1
mm1 = lmer(logTrailsA~Age + Sex + PovStat + Sex:PovStat + (Age|HNDid) + (1|subclass) , data=LifeIPVandCogData)
summary(mm1)
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: logTrailsA ~ Age + Sex + PovStat + Sex:PovStat + (Age | HNDid) + (1 | subclass)
## Data: LifeIPVandCogData
##
## REML criterion at convergence: 619.1
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 0.058549 0.2420
## Age 0.000195 0.0140 0.71
## subclass (Intercept) 0.001859 0.0431
## Residual 0.110865 0.3330
## Number of obs: 594, groups: HNDid, 297; subclass, 99
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 3.43e+00 3.31e-02 2.40e+02 103.59 < 2e-16
## Age 1.30e-02 2.33e-03 1.82e+02 5.56 9.5e-08
## SexMen 4.64e-02 5.01e-02 2.60e+02 0.93 0.3553
## PovStatBelow 3.76e-02 4.83e-02 2.72e+02 0.78 0.4375
## SexMen:PovStatBelow 2.23e-01 8.37e-02 2.71e+02 2.66 0.0082
##
## Correlation of Fixed Effects:
## (Intr) Age SexMen PvSttB
## Age 0.335
## SexMen -0.567 0.037
## PovStatBelw -0.603 -0.013 0.396
## SxMn:PvSttB 0.341 -0.018 -0.598 -0.578
Trails A Regression Model 2 (with CES)
(mm2 = lmer(logTrailsA~(Age + LifeIPV + PovStat + Sex + CESD)^5 + (Age|HNDid) + (1|subclass) , data=LifeIPVandCogData))
(st=step(mm2))
Trails A Final Regression Model 2 (with CES)
mm2=lmer(logTrailsA~ Age + LifeIPV + PovStat + Sex + CESD + Age:PovStat + Age:Sex + Age:CESD + LifeIPV:PovStat + LifeIPV:Sex + LifeIPV:CESD + PovStat:Sex + PovStat:CESD + Sex:CESD + Age:PovStat:Sex + Age:Sex:CESD + LifeIPV:PovStat:Sex + LifeIPV:PovStat:CESD + LifeIPV:Sex:CESD + PovStat:Sex:CESD + LifeIPV:PovStat:Sex:CESD +(Age|HNDid) + (1|subclass),data=LifeIPVandCogData)
summary(mm2)
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: logTrailsA ~ Age + LifeIPV + PovStat + Sex + CESD + Age:PovStat + Age:Sex + Age:CESD + LifeIPV:PovStat + LifeIPV:Sex + LifeIPV:CESD + PovStat:Sex + PovStat:CESD + Sex:CESD + Age:PovStat:Sex + Age:Sex:CESD + LifeIPV:PovStat:Sex + LifeIPV:PovStat:CESD + LifeIPV:Sex:CESD + PovStat:Sex:CESD + LifeIPV:PovStat:Sex:CESD + (Age | HNDid) + (1 | subclass)
## Data: LifeIPVandCogData
##
## REML criterion at convergence: 665.9
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 0.061764 0.2485
## Age 0.000214 0.0146 0.63
## subclass (Intercept) 0.002287 0.0478
## Residual 0.104649 0.3235
## Number of obs: 594, groups: HNDid, 297; subclass, 99
##
## Fixed effects:
## Estimate Std. Error df t value
## (Intercept) 3.39567 0.05363 314.00000 63.31
## Age 0.01003 0.00439 210.00000 2.29
## LifeIPV1 -0.01544 0.09040 382.00000 -0.17
## PovStatBelow 0.04721 0.08883 340.00000 0.53
## SexMen 0.04472 0.07631 274.00000 0.59
## CESD1 0.03327 0.07688 568.00000 0.43
## Age:PovStatBelow 0.00759 0.00595 192.00000 1.28
## Age:SexMen 0.00151 0.00714 246.00000 0.21
## Age:CESD1 -0.00248 0.00530 327.00000 -0.47
## LifeIPV1:PovStatBelow 0.15867 0.14034 399.00000 1.13
## LifeIPV1:SexMen 0.20059 0.26585 176.00000 0.75
## LifeIPV1:CESD1 0.01760 0.11763 520.00000 0.15
## PovStatBelow:SexMen 0.24469 0.13338 312.00000 1.83
## PovStatBelow:CESD1 -0.02720 0.12105 467.00000 -0.22
## SexMen:CESD1 0.26842 0.12648 534.00000 2.12
## Age:PovStatBelow:SexMen -0.02090 0.01061 212.00000 -1.97
## Age:SexMen:CESD1 0.02830 0.00955 382.00000 2.96
## LifeIPV1:PovStatBelow:SexMen -0.45977 0.34912 263.00000 -1.32
## LifeIPV1:PovStatBelow:CESD1 -0.11262 0.17951 538.00000 -0.63
## LifeIPV1:SexMen:CESD1 -0.42092 0.31468 292.00000 -1.34
## PovStatBelow:SexMen:CESD1 -0.36350 0.19024 531.00000 -1.91
## LifeIPV1:PovStatBelow:SexMen:CESD1 0.89285 0.41770 404.00000 2.14
## Pr(>|t|)
## (Intercept) <2e-16
## Age 0.0233
## LifeIPV1 0.8645
## PovStatBelow 0.5955
## SexMen 0.5584
## CESD1 0.6654
## Age:PovStatBelow 0.2038
## Age:SexMen 0.8330
## Age:CESD1 0.6402
## LifeIPV1:PovStatBelow 0.2589
## LifeIPV1:SexMen 0.4515
## LifeIPV1:CESD1 0.8811
## PovStatBelow:SexMen 0.0675
## PovStatBelow:CESD1 0.8223
## SexMen:CESD1 0.0343
## Age:PovStatBelow:SexMen 0.0503
## Age:SexMen:CESD1 0.0032
## LifeIPV1:PovStatBelow:SexMen 0.1890
## LifeIPV1:PovStatBelow:CESD1 0.5307
## LifeIPV1:SexMen:CESD1 0.1821
## PovStatBelow:SexMen:CESD1 0.0566
## LifeIPV1:PovStatBelow:SexMen:CESD1 0.0332
Effect Size
es = tes(t=2.14, n.1=99, n.2=198, level=95, dig=2, id=HNDid, data=LifeIPVandCogData)
summary(es)
## id N.total n.1 n.2
## Min. :8.03e+09 Min. :297 Min. :99 Min. :198
## 1st Qu.:8.17e+09 1st Qu.:297 1st Qu.:99 1st Qu.:198
## Median :8.19e+09 Median :297 Median :99 Median :198
## Mean :8.19e+09 Mean :297 Mean :99 Mean :198
## 3rd Qu.:8.21e+09 3rd Qu.:297 3rd Qu.:99 3rd Qu.:198
## Max. :8.22e+09 Max. :297 Max. :99 Max. :198
## d var.d l.d u.d
## Min. :0.26 Min. :0.02 Min. :0.02 Min. :0.51
## 1st Qu.:0.26 1st Qu.:0.02 1st Qu.:0.02 1st Qu.:0.51
## Median :0.26 Median :0.02 Median :0.02 Median :0.51
## Mean :0.26 Mean :0.02 Mean :0.02 Mean :0.51
## 3rd Qu.:0.26 3rd Qu.:0.02 3rd Qu.:0.02 3rd Qu.:0.51
## Max. :0.26 Max. :0.02 Max. :0.02 Max. :0.51
## U3.d cl.d cliffs.d pval.d
## Min. :60.4 Min. :57.4 Min. :0.15 Min. :0.03
## 1st Qu.:60.4 1st Qu.:57.4 1st Qu.:0.15 1st Qu.:0.03
## Median :60.4 Median :57.4 Median :0.15 Median :0.03
## Mean :60.4 Mean :57.4 Mean :0.15 Mean :0.03
## 3rd Qu.:60.4 3rd Qu.:57.4 3rd Qu.:0.15 3rd Qu.:0.03
## Max. :60.4 Max. :57.4 Max. :0.15 Max. :0.03
## g var.g l.g u.g
## Min. :0.26 Min. :0.02 Min. :0.02 Min. :0.51
## 1st Qu.:0.26 1st Qu.:0.02 1st Qu.:0.02 1st Qu.:0.51
## Median :0.26 Median :0.02 Median :0.02 Median :0.51
## Mean :0.26 Mean :0.02 Mean :0.02 Mean :0.51
## 3rd Qu.:0.26 3rd Qu.:0.02 3rd Qu.:0.02 3rd Qu.:0.51
## Max. :0.26 Max. :0.02 Max. :0.02 Max. :0.51
## U3.g cl.g pval.g r var.r
## Min. :60.4 Min. :57.4 Min. :0.03 Min. :0.12 Min. :0
## 1st Qu.:60.4 1st Qu.:57.4 1st Qu.:0.03 1st Qu.:0.12 1st Qu.:0
## Median :60.4 Median :57.4 Median :0.03 Median :0.12 Median :0
## Mean :60.4 Mean :57.4 Mean :0.03 Mean :0.12 Mean :0
## 3rd Qu.:60.4 3rd Qu.:57.4 3rd Qu.:0.03 3rd Qu.:0.12 3rd Qu.:0
## Max. :60.4 Max. :57.4 Max. :0.03 Max. :0.12 Max. :0
## l.r u.r pval.r fisher.z var.z
## Min. :0.01 Min. :0.23 Min. :0.03 Min. :0.12 Min. :0
## 1st Qu.:0.01 1st Qu.:0.23 1st Qu.:0.03 1st Qu.:0.12 1st Qu.:0
## Median :0.01 Median :0.23 Median :0.03 Median :0.12 Median :0
## Mean :0.01 Mean :0.23 Mean :0.03 Mean :0.12 Mean :0
## 3rd Qu.:0.01 3rd Qu.:0.23 3rd Qu.:0.03 3rd Qu.:0.12 3rd Qu.:0
## Max. :0.01 Max. :0.23 Max. :0.03 Max. :0.12 Max. :0
## l.z u.z OR l.or
## Min. :0.01 Min. :0.24 Min. :1.61 Min. :1.04
## 1st Qu.:0.01 1st Qu.:0.24 1st Qu.:1.61 1st Qu.:1.04
## Median :0.01 Median :0.24 Median :1.61 Median :1.04
## Mean :0.01 Mean :0.24 Mean :1.61 Mean :1.04
## 3rd Qu.:0.01 3rd Qu.:0.24 3rd Qu.:1.61 3rd Qu.:1.04
## Max. :0.01 Max. :0.24 Max. :1.61 Max. :1.04
## u.or pval.or lOR l.lor
## Min. :2.51 Min. :0.03 Min. :0.48 Min. :0.04
## 1st Qu.:2.51 1st Qu.:0.03 1st Qu.:0.48 1st Qu.:0.04
## Median :2.51 Median :0.03 Median :0.48 Median :0.04
## Mean :2.51 Mean :0.03 Mean :0.48 Mean :0.04
## 3rd Qu.:2.51 3rd Qu.:0.03 3rd Qu.:0.48 3rd Qu.:0.04
## Max. :2.51 Max. :0.03 Max. :0.48 Max. :0.04
## u.lor pval.lor NNT
## Min. :0.92 Min. :0.03 Min. :12.3
## 1st Qu.:0.92 1st Qu.:0.03 1st Qu.:12.3
## Median :0.92 Median :0.03 Median :12.3
## Mean :0.92 Mean :0.03 Mean :12.3
## 3rd Qu.:0.92 3rd Qu.:0.03 3rd Qu.:12.3
## Max. :0.92 Max. :0.03 Max. :12.3
Trail Making Test A
LifeIPV x Sex x Poverty x CESD Interaction Plot
library(zUtil)
library(zStat)
#Stratify by Gender
LifeIPVWomen <- LifeIPVandCogData[ which(LifeIPVandCogData$Sex=='Women'),]
LifeIPVMen <- LifeIPVandCogData[ which(LifeIPVandCogData$Sex=='Men'),]
pAge = seq(-20,20)
##Women
pCESD=(c(0,1))
hatIPVcog1 = zMixHat(LifeIPVWomen, mm2, vary = "CESD=pCESD,LifeIPV=zQ(0,1),PovStat=zQ(Above,Below)",fixedCov=c("Sex","Age"))
par(mar=c(4,4,.5,2),las = 1, lwd = 2)
HNDcolors = HNDpltColors()
with(hatIPVcog1[hatIPVcog1$PovStat == "Above" & hatIPVcog1$LifeIPV == "0", ], plot(pCESD,hat, lty = 1, col = "black", type = "l",ylim = c(3,4), ylab = "log(Trails A)", xlab = "Depression",xaxt = "n"))
with(hatIPVcog1[hatIPVcog1$PovStat == "Above" & hatIPVcog1$LifeIPV == "1", ], lines(pCESD,hat, lty = 2, col = "black"))
with(hatIPVcog1[hatIPVcog1$PovStat == "Below" & hatIPVcog1$LifeIPV == "0", ], lines(pCESD,hat, lty = 1, col = "red", typ = "l", ylim = c(3,4), ylab = "Trails A", xlab = "Depression",xaxt = "n"))
with(hatIPVcog1[hatIPVcog1$PovStat == "Below" & hatIPVcog1$LifeIPV == "1", ], lines(pCESD,hat, lty = 2, col = "red"))
axis(1,at=c(0,1),labels = c("0","1"))
legend(0,4, zQ(NoIPV,IPV), lty = 1:2, col = "black",cex=.90,bty="n")
text(0,3.1,"Above Poverty in Black", adj = c(0,0), col = "black",cex=.95)
text(0,3,"Below Poverty in Red", adj = c(0,0), col = "red",cex=.95)
text(.8,3.9,"WOMEN", adj = c(0,0), col = "red",cex=.95)
##Men
pCESD=(c(0,1))
hatIPVcog1 = zMixHat(LifeIPVMen, mm2, vary = "CESD=pCESD,LifeIPV=zQ(0,1),PovStat=zQ(Above,Below)",fixedCov=c("Sex","Age"))
par(mar=c(4,4,.5,2),las = 1, lwd = 2)
HNDcolors = HNDpltColors()
with(hatIPVcog1[hatIPVcog1$PovStat == "Above" & hatIPVcog1$LifeIPV == "0", ], plot(pCESD,hat, lty = 1, col = "black", type = "l",ylim = c(3,4), ylab = "log(Trails A)", xlab = "Depression",xaxt = "n"))
with(hatIPVcog1[hatIPVcog1$PovStat == "Above" & hatIPVcog1$LifeIPV == "1", ], lines(pCESD,hat, lty = 2, col = "black"))
with(hatIPVcog1[hatIPVcog1$PovStat == "Below" & hatIPVcog1$LifeIPV == "0", ], lines(pCESD,hat, lty = 1, col = "red", typ = "l", ylim = c(3,4), ylab = "Trails A", xlab = "Depression",xaxt = "n"))
with(hatIPVcog1[hatIPVcog1$PovStat == "Below" & hatIPVcog1$LifeIPV == "1", ], lines(pCESD,hat, lty = 2, col = "red"))
axis(1,at=c(0,1),labels = c("0","1"))
legend(0,4, zQ(NoIPV,IPV), lty = 1:2, col = "black",cex=.90,bty="n")
text(0,3.1,"Above Poverty in Black", adj = c(0,0), col = "black",cex=.95)
text(0,3,"Below Poverty in Red", adj = c(0,0), col = "red",cex=.95)
text(.8,3.9,"MEN", adj = c(0,0), col = "blue",cex=.95)
Trails B Regression Model 1
(mm1 = lmer(logTrailsB~(Age + LifeIPV + Sex + PovStat)^4 + (Age|HNDid) + (1|subclass) , data=LifeIPVandCogData))
(st = step(mm1))
Trails B Final Regression Model 1
mm1 = lmer(logTrailsB~Age + LifeIPV + Sex + PovStat + Age:PovStat + LifeIPV:Sex + (Age|HNDid) + (1|subclass) , data=LifeIPVandCogData)
summary(mm1)
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: logTrailsB ~ Age + LifeIPV + Sex + PovStat + Age:PovStat + LifeIPV:Sex + (Age | HNDid) + (1 | subclass)
## Data: LifeIPVandCogData
##
## REML criterion at convergence: 1117
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 0.381645 0.6178
## Age 0.000168 0.0130 0.33
## subclass (Intercept) 0.006381 0.0799
## Residual 0.145953 0.3820
## Number of obs: 594, groups: HNDid, 297; subclass, 99
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 4.56e+00 7.45e-02 2.56e+02 61.21 <2e-16
## Age 9.60e-03 5.11e-03 2.25e+02 1.88 0.0617
## LifeIPV1 1.02e-03 9.81e-02 2.24e+02 0.01 0.9917
## SexMen 4.76e-02 9.65e-02 2.90e+02 0.49 0.6218
## PovStatBelow 2.22e-01 8.60e-02 2.19e+02 2.58 0.0105
## Age:PovStatBelow 2.09e-02 7.97e-03 1.86e+02 2.62 0.0095
## LifeIPV1:SexMen 4.41e-01 2.16e-01 2.80e+02 2.04 0.0421
##
## Correlation of Fixed Effects:
## (Intr) Age LfIPV1 SexMen PvSttB Ag:PSB
## Age 0.271
## LifeIPV1 -0.553 -0.047
## SexMen -0.605 0.003 0.437
## PovStatBelw -0.439 -0.222 -0.071 0.053
## Ag:PvSttBlw -0.144 -0.635 -0.029 -0.022 0.328
## LfIPV1:SxMn 0.301 0.053 -0.453 -0.447 -0.065 0.000
Trails B Regression Model 2 (with CES)
mm2 = lmer(logTrailsB~(Age + LifeIPV + PovStat + Sex + CESD)^5 + (Age|HNDid) + (1|subclass) , data=LifeIPVandCogData)
(st=step(mm2))
Trails B Final Regression Model 2 (with CES)
mm2=lmer(logTrailsB~Age + LifeIPV + Sex + PovStat + Age:PovStat + LifeIPV:Sex + (Age|HNDid) + (1|subclass) , data=LifeIPVandCogData)
summary(mm2)
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: logTrailsB ~ Age + LifeIPV + Sex + PovStat + Age:PovStat + LifeIPV:Sex + (Age | HNDid) + (1 | subclass)
## Data: LifeIPVandCogData
##
## REML criterion at convergence: 1117
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 0.381645 0.6178
## Age 0.000168 0.0130 0.33
## subclass (Intercept) 0.006381 0.0799
## Residual 0.145953 0.3820
## Number of obs: 594, groups: HNDid, 297; subclass, 99
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 4.56e+00 7.45e-02 2.56e+02 61.21 <2e-16
## Age 9.60e-03 5.11e-03 2.25e+02 1.88 0.0617
## LifeIPV1 1.02e-03 9.81e-02 2.24e+02 0.01 0.9917
## SexMen 4.76e-02 9.65e-02 2.90e+02 0.49 0.6218
## PovStatBelow 2.22e-01 8.60e-02 2.19e+02 2.58 0.0105
## Age:PovStatBelow 2.09e-02 7.97e-03 1.86e+02 2.62 0.0095
## LifeIPV1:SexMen 4.41e-01 2.16e-01 2.80e+02 2.04 0.0421
##
## Correlation of Fixed Effects:
## (Intr) Age LfIPV1 SexMen PvSttB Ag:PSB
## Age 0.271
## LifeIPV1 -0.553 -0.047
## SexMen -0.605 0.003 0.437
## PovStatBelw -0.439 -0.222 -0.071 0.053
## Ag:PvSttBlw -0.144 -0.635 -0.029 -0.022 0.328
## LfIPV1:SxMn 0.301 0.053 -0.453 -0.447 -0.065 0.000
Effect Size
es = tes(t=2.04, n.1=99, n.2=198, level=95, dig=2, id=HNDid, data=LifeIPVandCogData)
summary(es)
## id N.total n.1 n.2
## Min. :8.03e+09 Min. :297 Min. :99 Min. :198
## 1st Qu.:8.17e+09 1st Qu.:297 1st Qu.:99 1st Qu.:198
## Median :8.19e+09 Median :297 Median :99 Median :198
## Mean :8.19e+09 Mean :297 Mean :99 Mean :198
## 3rd Qu.:8.21e+09 3rd Qu.:297 3rd Qu.:99 3rd Qu.:198
## Max. :8.22e+09 Max. :297 Max. :99 Max. :198
## d var.d l.d u.d
## Min. :0.25 Min. :0.02 Min. :0.01 Min. :0.49
## 1st Qu.:0.25 1st Qu.:0.02 1st Qu.:0.01 1st Qu.:0.49
## Median :0.25 Median :0.02 Median :0.01 Median :0.49
## Mean :0.25 Mean :0.02 Mean :0.01 Mean :0.49
## 3rd Qu.:0.25 3rd Qu.:0.02 3rd Qu.:0.01 3rd Qu.:0.49
## Max. :0.25 Max. :0.02 Max. :0.01 Max. :0.49
## U3.d cl.d cliffs.d pval.d g
## Min. :59.9 Min. :57 Min. :0.14 Min. :0.04 Min. :0.25
## 1st Qu.:59.9 1st Qu.:57 1st Qu.:0.14 1st Qu.:0.04 1st Qu.:0.25
## Median :59.9 Median :57 Median :0.14 Median :0.04 Median :0.25
## Mean :59.9 Mean :57 Mean :0.14 Mean :0.04 Mean :0.25
## 3rd Qu.:59.9 3rd Qu.:57 3rd Qu.:0.14 3rd Qu.:0.04 3rd Qu.:0.25
## Max. :59.9 Max. :57 Max. :0.14 Max. :0.04 Max. :0.25
## var.g l.g u.g U3.g cl.g
## Min. :0.02 Min. :0.01 Min. :0.49 Min. :59.9 Min. :57
## 1st Qu.:0.02 1st Qu.:0.01 1st Qu.:0.49 1st Qu.:59.9 1st Qu.:57
## Median :0.02 Median :0.01 Median :0.49 Median :59.9 Median :57
## Mean :0.02 Mean :0.01 Mean :0.49 Mean :59.9 Mean :57
## 3rd Qu.:0.02 3rd Qu.:0.01 3rd Qu.:0.49 3rd Qu.:59.9 3rd Qu.:57
## Max. :0.02 Max. :0.01 Max. :0.49 Max. :59.9 Max. :57
## pval.g r var.r l.r u.r
## Min. :0.04 Min. :0.12 Min. :0 Min. :0 Min. :0.23
## 1st Qu.:0.04 1st Qu.:0.12 1st Qu.:0 1st Qu.:0 1st Qu.:0.23
## Median :0.04 Median :0.12 Median :0 Median :0 Median :0.23
## Mean :0.04 Mean :0.12 Mean :0 Mean :0 Mean :0.23
## 3rd Qu.:0.04 3rd Qu.:0.12 3rd Qu.:0 3rd Qu.:0 3rd Qu.:0.23
## Max. :0.04 Max. :0.12 Max. :0 Max. :0 Max. :0.23
## pval.r fisher.z var.z l.z u.z
## Min. :0.04 Min. :0.12 Min. :0 Min. :0 Min. :0.23
## 1st Qu.:0.04 1st Qu.:0.12 1st Qu.:0 1st Qu.:0 1st Qu.:0.23
## Median :0.04 Median :0.12 Median :0 Median :0 Median :0.23
## Mean :0.04 Mean :0.12 Mean :0 Mean :0 Mean :0.23
## 3rd Qu.:0.04 3rd Qu.:0.12 3rd Qu.:0 3rd Qu.:0 3rd Qu.:0.23
## Max. :0.04 Max. :0.12 Max. :0 Max. :0 Max. :0.23
## OR l.or u.or pval.or
## Min. :1.58 Min. :1.01 Min. :2.45 Min. :0.04
## 1st Qu.:1.58 1st Qu.:1.01 1st Qu.:2.45 1st Qu.:0.04
## Median :1.58 Median :1.01 Median :2.45 Median :0.04
## Mean :1.58 Mean :1.01 Mean :2.45 Mean :0.04
## 3rd Qu.:1.58 3rd Qu.:1.01 3rd Qu.:2.45 3rd Qu.:0.04
## Max. :1.58 Max. :1.01 Max. :2.45 Max. :0.04
## lOR l.lor u.lor pval.lor NNT
## Min. :0.46 Min. :0.01 Min. :0.9 Min. :0.04 Min. :12.9
## 1st Qu.:0.46 1st Qu.:0.01 1st Qu.:0.9 1st Qu.:0.04 1st Qu.:12.9
## Median :0.46 Median :0.01 Median :0.9 Median :0.04 Median :12.9
## Mean :0.46 Mean :0.01 Mean :0.9 Mean :0.04 Mean :12.9
## 3rd Qu.:0.46 3rd Qu.:0.01 3rd Qu.:0.9 3rd Qu.:0.04 3rd Qu.:12.9
## Max. :0.46 Max. :0.01 Max. :0.9 Max. :0.04 Max. :12.9
Trail Making Test B
LifeIPV x Sex Interaction Plot
FluencyWord Regression Model 1
(mm1 = lmer(FluencyWord~(Age + LifeIPV + Sex + PovStat)^4 + (Age|HNDid) + (1|subclass) , data=LifeIPVandCogData))
(st = step(mm1))
FluencyWord Final Regression Model 1
mm1 = lmer(FluencyWord~Age + LifeIPV + Sex + PovStat + Age:LifeIPV + LifeIPV:Sex + LifeIPV:PovStat + Sex:PovStat + LifeIPV:Sex:PovStat + (Age|HNDid) + (1|subclass) , data=LifeIPVandCogData)
summary(mm1)
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: FluencyWord ~ Age + LifeIPV + Sex + PovStat + Age:LifeIPV + LifeIPV:Sex + LifeIPV:PovStat + Sex:PovStat + LifeIPV:Sex:PovStat + (Age | HNDid) + (1 | subclass)
## Data: LifeIPVandCogData
##
## REML criterion at convergence: 3528
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 21.39050 4.6250
## Age 0.00408 0.0639 -0.10
## subclass (Intercept) 1.86022 1.3639
## Residual 9.20251 3.0336
## Number of obs: 594, groups: HNDid, 297; subclass, 99
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 19.3854 0.6594 275.0000 29.40 <2e-16
## Age -0.1047 0.0373 190.8000 -2.81 0.0055
## LifeIPV1 -0.2725 1.0046 233.2000 -0.27 0.7865
## SexMen 0.2943 0.9314 282.0000 0.32 0.7523
## PovStatBelow -1.6112 1.0432 270.6000 -1.54 0.1236
## Age:LifeIPV1 0.1695 0.0643 162.1000 2.64 0.0092
## LifeIPV1:SexMen 5.3653 2.5053 254.8000 2.14 0.0332
## LifeIPV1:PovStatBelow 1.0322 1.5517 273.0000 0.67 0.5065
## SexMen:PovStatBelow 1.6506 1.6058 277.4000 1.03 0.3049
## LifeIPV1:SexMen:PovStatBelow -8.4534 3.4918 269.8000 -2.42 0.0161
##
## Correlation of Fixed Effects:
## (Intr) Age LfIPV1 SexMen PvSttB A:LIPV LfIPV1:SM LIPV1:P
## Age 0.169
## LifeIPV1 -0.625 -0.103
## SexMen -0.670 0.008 0.438
## PovStatBelw -0.590 0.020 0.385 0.420
## Age:LifIPV1 -0.095 -0.556 0.119 0.000 -0.014
## LfIPV1:SxMn 0.243 -0.005 -0.392 -0.368 -0.151 0.098
## LfIPV1:PvSB 0.395 -0.019 -0.650 -0.281 -0.670 -0.010 0.255
## SxMn:PvSttB 0.381 -0.034 -0.248 -0.577 -0.650 0.019 0.218 0.433
## LIPV1:SM:PS -0.174 0.013 0.283 0.268 0.294 -0.041 -0.716 -0.439
## SM:PSB
## Age
## LifeIPV1
## SexMen
## PovStatBelw
## Age:LifIPV1
## LfIPV1:SxMn
## LfIPV1:PvSB
## SxMn:PvSttB
## LIPV1:SM:PS -0.463
FluencyWord Regression Model 2 (with CES)
(mm2 = lmer(FluencyWord~(Age + LifeIPV + PovStat + Sex + CESD)^5 + (Age|HNDid) + (1|subclass) , data=LifeIPVandCogData))
(st=step(mm2))
FluencyWord Final Regression Model 2 (with CES)
mm2=lmer(FluencyWord~Age + LifeIPV + Sex + PovStat + Age:LifeIPV + LifeIPV:PovStat + LifeIPV:Sex + PovStat:Sex + LifeIPV:PovStat:Sex + (Age|HNDid) + (1|subclass) , data=LifeIPVandCogData)
summary(mm2)
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: FluencyWord ~ Age + LifeIPV + Sex + PovStat + Age:LifeIPV + LifeIPV:PovStat + LifeIPV:Sex + PovStat:Sex + LifeIPV:PovStat:Sex + (Age | HNDid) + (1 | subclass)
## Data: LifeIPVandCogData
##
## REML criterion at convergence: 3528
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 21.39050 4.6250
## Age 0.00408 0.0639 -0.10
## subclass (Intercept) 1.86022 1.3639
## Residual 9.20251 3.0336
## Number of obs: 594, groups: HNDid, 297; subclass, 99
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 19.3854 0.6594 275.0000 29.40 <2e-16
## Age -0.1047 0.0373 190.8000 -2.81 0.0055
## LifeIPV1 -0.2725 1.0046 233.2000 -0.27 0.7865
## SexMen 0.2943 0.9314 282.0000 0.32 0.7523
## PovStatBelow -1.6112 1.0432 270.6000 -1.54 0.1236
## Age:LifeIPV1 0.1695 0.0643 162.1000 2.64 0.0092
## LifeIPV1:PovStatBelow 1.0322 1.5517 273.0000 0.67 0.5065
## LifeIPV1:SexMen 5.3653 2.5053 254.8000 2.14 0.0332
## SexMen:PovStatBelow 1.6506 1.6058 277.4000 1.03 0.3049
## LifeIPV1:SexMen:PovStatBelow -8.4534 3.4918 269.8000 -2.42 0.0161
##
## Correlation of Fixed Effects:
## (Intr) Age LfIPV1 SexMen PvSttB A:LIPV LIPV1:P LfIPV1:SM
## Age 0.169
## LifeIPV1 -0.625 -0.103
## SexMen -0.670 0.008 0.438
## PovStatBelw -0.590 0.020 0.385 0.420
## Age:LifIPV1 -0.095 -0.556 0.119 0.000 -0.014
## LfIPV1:PvSB 0.395 -0.019 -0.650 -0.281 -0.670 -0.010
## LfIPV1:SxMn 0.243 -0.005 -0.392 -0.368 -0.151 0.098 0.255
## SxMn:PvSttB 0.381 -0.034 -0.248 -0.577 -0.650 0.019 0.433 0.218
## LIPV1:SM:PS -0.174 0.013 0.283 0.268 0.294 -0.041 -0.439 -0.716
## SM:PSB
## Age
## LifeIPV1
## SexMen
## PovStatBelw
## Age:LifIPV1
## LfIPV1:PvSB
## LfIPV1:SxMn
## SxMn:PvSttB
## LIPV1:SM:PS -0.463
Effect Size
Age:Life
Est. = 0.16949883 SE = 0.06428640 df = 162.13299000 t = 2.63662 p = 0.0091877
es = tes(t=2.64, n.1=99, n.2=198, level=95, dig=2, id=HNDid, data=LifeIPVandCogData)
summary(es)
## id N.total n.1 n.2
## Min. :8.03e+09 Min. :297 Min. :99 Min. :198
## 1st Qu.:8.17e+09 1st Qu.:297 1st Qu.:99 1st Qu.:198
## Median :8.19e+09 Median :297 Median :99 Median :198
## Mean :8.19e+09 Mean :297 Mean :99 Mean :198
## 3rd Qu.:8.21e+09 3rd Qu.:297 3rd Qu.:99 3rd Qu.:198
## Max. :8.22e+09 Max. :297 Max. :99 Max. :198
## d var.d l.d u.d
## Min. :0.32 Min. :0.02 Min. :0.08 Min. :0.57
## 1st Qu.:0.32 1st Qu.:0.02 1st Qu.:0.08 1st Qu.:0.57
## Median :0.32 Median :0.02 Median :0.08 Median :0.57
## Mean :0.32 Mean :0.02 Mean :0.08 Mean :0.57
## 3rd Qu.:0.32 3rd Qu.:0.02 3rd Qu.:0.08 3rd Qu.:0.57
## Max. :0.32 Max. :0.02 Max. :0.08 Max. :0.57
## U3.d cl.d cliffs.d pval.d
## Min. :62.7 Min. :59.1 Min. :0.18 Min. :0.01
## 1st Qu.:62.7 1st Qu.:59.1 1st Qu.:0.18 1st Qu.:0.01
## Median :62.7 Median :59.1 Median :0.18 Median :0.01
## Mean :62.7 Mean :59.1 Mean :0.18 Mean :0.01
## 3rd Qu.:62.7 3rd Qu.:59.1 3rd Qu.:0.18 3rd Qu.:0.01
## Max. :62.7 Max. :59.1 Max. :0.18 Max. :0.01
## g var.g l.g u.g
## Min. :0.32 Min. :0.02 Min. :0.08 Min. :0.57
## 1st Qu.:0.32 1st Qu.:0.02 1st Qu.:0.08 1st Qu.:0.57
## Median :0.32 Median :0.02 Median :0.08 Median :0.57
## Mean :0.32 Mean :0.02 Mean :0.08 Mean :0.57
## 3rd Qu.:0.32 3rd Qu.:0.02 3rd Qu.:0.08 3rd Qu.:0.57
## Max. :0.32 Max. :0.02 Max. :0.08 Max. :0.57
## U3.g cl.g pval.g r var.r
## Min. :62.7 Min. :59.1 Min. :0.01 Min. :0.15 Min. :0
## 1st Qu.:62.7 1st Qu.:59.1 1st Qu.:0.01 1st Qu.:0.15 1st Qu.:0
## Median :62.7 Median :59.1 Median :0.01 Median :0.15 Median :0
## Mean :62.7 Mean :59.1 Mean :0.01 Mean :0.15 Mean :0
## 3rd Qu.:62.7 3rd Qu.:59.1 3rd Qu.:0.01 3rd Qu.:0.15 3rd Qu.:0
## Max. :62.7 Max. :59.1 Max. :0.01 Max. :0.15 Max. :0
## l.r u.r pval.r fisher.z var.z
## Min. :0.04 Min. :0.26 Min. :0.01 Min. :0.15 Min. :0
## 1st Qu.:0.04 1st Qu.:0.26 1st Qu.:0.01 1st Qu.:0.15 1st Qu.:0
## Median :0.04 Median :0.26 Median :0.01 Median :0.15 Median :0
## Mean :0.04 Mean :0.26 Mean :0.01 Mean :0.15 Mean :0
## 3rd Qu.:0.04 3rd Qu.:0.26 3rd Qu.:0.01 3rd Qu.:0.15 3rd Qu.:0
## Max. :0.04 Max. :0.26 Max. :0.01 Max. :0.15 Max. :0
## l.z u.z OR l.or u.or
## Min. :0.04 Min. :0.27 Min. :1.8 Min. :1.16 Min. :2.8
## 1st Qu.:0.04 1st Qu.:0.27 1st Qu.:1.8 1st Qu.:1.16 1st Qu.:2.8
## Median :0.04 Median :0.27 Median :1.8 Median :1.16 Median :2.8
## Mean :0.04 Mean :0.27 Mean :1.8 Mean :1.16 Mean :2.8
## 3rd Qu.:0.04 3rd Qu.:0.27 3rd Qu.:1.8 3rd Qu.:1.16 3rd Qu.:2.8
## Max. :0.04 Max. :0.27 Max. :1.8 Max. :1.16 Max. :2.8
## pval.or lOR l.lor u.lor
## Min. :0.01 Min. :0.59 Min. :0.15 Min. :1.03
## 1st Qu.:0.01 1st Qu.:0.59 1st Qu.:0.15 1st Qu.:1.03
## Median :0.01 Median :0.59 Median :0.15 Median :1.03
## Mean :0.01 Mean :0.59 Mean :0.15 Mean :1.03
## 3rd Qu.:0.01 3rd Qu.:0.59 3rd Qu.:0.15 3rd Qu.:1.03
## Max. :0.01 Max. :0.59 Max. :0.15 Max. :1.03
## pval.lor NNT
## Min. :0.01 Min. :9.74
## 1st Qu.:0.01 1st Qu.:9.74
## Median :0.01 Median :9.74
## Mean :0.01 Mean :9.74
## 3rd Qu.:0.01 3rd Qu.:9.74
## Max. :0.01 Max. :9.74
LifeIPV:Sex:PovStat
Est. = -8.45335593 SE = 3.49175509 df = 269.80132000 t = -2.42095 p = 0.0161404
es = tes(t=-2.42, n.1=99, n.2=198, level=95, dig=2, id=HNDid, data=LifeIPVandCogData)
summary(es)
## id N.total n.1 n.2
## Min. :8.03e+09 Min. :297 Min. :99 Min. :198
## 1st Qu.:8.17e+09 1st Qu.:297 1st Qu.:99 1st Qu.:198
## Median :8.19e+09 Median :297 Median :99 Median :198
## Mean :8.19e+09 Mean :297 Mean :99 Mean :198
## 3rd Qu.:8.21e+09 3rd Qu.:297 3rd Qu.:99 3rd Qu.:198
## Max. :8.22e+09 Max. :297 Max. :99 Max. :198
## d var.d l.d u.d
## Min. :-0.3 Min. :0.02 Min. :-0.54 Min. :-0.05
## 1st Qu.:-0.3 1st Qu.:0.02 1st Qu.:-0.54 1st Qu.:-0.05
## Median :-0.3 Median :0.02 Median :-0.54 Median :-0.05
## Mean :-0.3 Mean :0.02 Mean :-0.54 Mean :-0.05
## 3rd Qu.:-0.3 3rd Qu.:0.02 3rd Qu.:-0.54 3rd Qu.:-0.05
## Max. :-0.3 Max. :0.02 Max. :-0.54 Max. :-0.05
## U3.d cl.d cliffs.d pval.d
## Min. :38.3 Min. :41.7 Min. :-0.17 Min. :0.02
## 1st Qu.:38.3 1st Qu.:41.7 1st Qu.:-0.17 1st Qu.:0.02
## Median :38.3 Median :41.7 Median :-0.17 Median :0.02
## Mean :38.3 Mean :41.7 Mean :-0.17 Mean :0.02
## 3rd Qu.:38.3 3rd Qu.:41.7 3rd Qu.:-0.17 3rd Qu.:0.02
## Max. :38.3 Max. :41.7 Max. :-0.17 Max. :0.02
## g var.g l.g u.g
## Min. :-0.3 Min. :0.02 Min. :-0.54 Min. :-0.05
## 1st Qu.:-0.3 1st Qu.:0.02 1st Qu.:-0.54 1st Qu.:-0.05
## Median :-0.3 Median :0.02 Median :-0.54 Median :-0.05
## Mean :-0.3 Mean :0.02 Mean :-0.54 Mean :-0.05
## 3rd Qu.:-0.3 3rd Qu.:0.02 3rd Qu.:-0.54 3rd Qu.:-0.05
## Max. :-0.3 Max. :0.02 Max. :-0.54 Max. :-0.05
## U3.g cl.g pval.g r var.r
## Min. :38.3 Min. :41.7 Min. :0.02 Min. :0.14 Min. :0
## 1st Qu.:38.3 1st Qu.:41.7 1st Qu.:0.02 1st Qu.:0.14 1st Qu.:0
## Median :38.3 Median :41.7 Median :0.02 Median :0.14 Median :0
## Mean :38.3 Mean :41.7 Mean :0.02 Mean :0.14 Mean :0
## 3rd Qu.:38.3 3rd Qu.:41.7 3rd Qu.:0.02 3rd Qu.:0.14 3rd Qu.:0
## Max. :38.3 Max. :41.7 Max. :0.02 Max. :0.14 Max. :0
## l.r u.r pval.r fisher.z var.z
## Min. :0.03 Min. :0.25 Min. :0.02 Min. :0.14 Min. :0
## 1st Qu.:0.03 1st Qu.:0.25 1st Qu.:0.02 1st Qu.:0.14 1st Qu.:0
## Median :0.03 Median :0.25 Median :0.02 Median :0.14 Median :0
## Mean :0.03 Mean :0.25 Mean :0.02 Mean :0.14 Mean :0
## 3rd Qu.:0.03 3rd Qu.:0.25 3rd Qu.:0.02 3rd Qu.:0.14 3rd Qu.:0
## Max. :0.03 Max. :0.25 Max. :0.02 Max. :0.14 Max. :0
## l.z u.z OR l.or
## Min. :0.03 Min. :0.26 Min. :0.58 Min. :0.37
## 1st Qu.:0.03 1st Qu.:0.26 1st Qu.:0.58 1st Qu.:0.37
## Median :0.03 Median :0.26 Median :0.58 Median :0.37
## Mean :0.03 Mean :0.26 Mean :0.58 Mean :0.37
## 3rd Qu.:0.03 3rd Qu.:0.26 3rd Qu.:0.58 3rd Qu.:0.37
## Max. :0.03 Max. :0.26 Max. :0.58 Max. :0.37
## u.or pval.or lOR l.lor
## Min. :0.91 Min. :0.02 Min. :-0.54 Min. :-0.98
## 1st Qu.:0.91 1st Qu.:0.02 1st Qu.:-0.54 1st Qu.:-0.98
## Median :0.91 Median :0.02 Median :-0.54 Median :-0.98
## Mean :0.91 Mean :0.02 Mean :-0.54 Mean :-0.98
## 3rd Qu.:0.91 3rd Qu.:0.02 3rd Qu.:-0.54 3rd Qu.:-0.98
## Max. :0.91 Max. :0.02 Max. :-0.54 Max. :-0.98
## u.lor pval.lor NNT
## Min. :-0.1 Min. :0.02 Min. :-13.8
## 1st Qu.:-0.1 1st Qu.:0.02 1st Qu.:-13.8
## Median :-0.1 Median :0.02 Median :-13.8
## Mean :-0.1 Mean :0.02 Mean :-13.8
## 3rd Qu.:-0.1 3rd Qu.:0.02 3rd Qu.:-13.8
## Max. :-0.1 Max. :0.02 Max. :-13.8
Animal Naming Test
Age x LifeIPV Interaction Plot
Animal Naming Test
LifeIPV x PovStat x Sex Interaction Plot
