Dan's Thesis Analyses - Tandem Stands
load(file="/Users/meganwilliams/Desktop/UMBC 2016/Dan's Thesis/DanThesisFinal.rdata")
library(lmerTest)
Tandem Stands - Model 1 (Main Effects)
mm1 = lmer(PhyPerfSecsTandem ~Age + Sex + EF + Race + BMI + PovStat + WRATtotal + (Age|HNDid), data=DanThesisFinal,REML=F)
summary(mm1)
## Linear mixed model fit by maximum likelihood ['merModLmerTest']
## Formula: PhyPerfSecsTandem ~ Age + Sex + EF + Race + BMI + PovStat + WRATtotal + (Age | HNDid)
## Data: DanThesisFinal
##
## AIC BIC logLik deviance
## 10763 10831 -5369 10739
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 2.2843 1.511
## Age 0.0152 0.123 1.00
## Residual 5.6685 2.381
## Number of obs: 2168, groups: HNDid, 1440
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 2.85e+01 5.37e-01 1.96e+03 53.11 < 2e-16
## Age -3.55e-02 7.69e-03 1.42e+03 -4.62 4.2e-06
## SexMen 8.90e-02 1.35e-01 2.01e+03 0.66 0.5088
## EF 4.11e-03 2.90e-02 2.08e+03 0.14 0.8872
## RaceAfrAm 3.46e-01 1.39e-01 2.02e+03 2.49 0.0128
## BMI -1.19e-02 8.61e-03 2.07e+03 -1.39 0.1658
## PovStatBelow -2.55e-01 1.40e-01 2.02e+03 -1.83 0.0679
## WRATtotal 2.61e-02 9.99e-03 1.99e+03 2.61 0.0091
##
## Correlation of Fixed Effects:
## (Intr) Age SexMen EF RcAfrA BMI PvSttB
## Age 0.034
## SexMen -0.208 -0.018
## EF 0.275 0.102 -0.042
## RaceAfrAm -0.269 0.030 -0.020 0.122
## BMI -0.489 -0.010 0.182 0.041 -0.033
## PovStatBelw -0.267 0.077 0.031 0.081 -0.053 0.063
## WRATtotal -0.831 0.009 0.022 -0.409 0.174 -0.016 0.180
Tandem Stands - Model 2 (Two-Way Interactions)
mm2 = lmer(PhyPerfSecsTandem ~(Age + Sex + EF + Race)^2 + BMI + PovStat + WRATtotal + (Age|HNDid), data=DanThesisFinal,REML=F)
summary(mm2)
## Linear mixed model fit by maximum likelihood ['merModLmerTest']
## Formula: PhyPerfSecsTandem ~ (Age + Sex + EF + Race)^2 + BMI + PovStat + WRATtotal + (Age | HNDid)
## Data: DanThesisFinal
##
## AIC BIC logLik deviance
## 10766 10868 -5365 10730
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 2.2575 1.503
## Age 0.0149 0.122 1.00
## Residual 5.6656 2.380
## Number of obs: 2168, groups: HNDid, 1440
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 2.83e+01 5.45e-01 1.94e+03 52.00 < 2e-16
## Age -6.79e-02 1.40e-02 1.46e+03 -4.84 1.4e-06
## SexMen 1.04e-01 2.18e-01 1.80e+03 0.48 0.6320
## EF 1.12e-02 4.75e-02 2.11e+03 0.24 0.8129
## RaceAfrAm 4.46e-01 1.88e-01 1.66e+03 2.37 0.0181
## BMI -1.07e-02 8.71e-03 2.08e+03 -1.22 0.2208
## PovStatBelow -2.52e-01 1.40e-01 2.02e+03 -1.80 0.0719
## WRATtotal 2.75e-02 1.00e-02 1.98e+03 2.74 0.0061
## Age:SexMen 1.20e-02 1.57e-02 1.45e+03 0.76 0.4447
## Age:EF 5.12e-03 2.96e-03 1.64e+03 1.73 0.0839
## Age:RaceAfrAm 4.24e-02 1.61e-02 1.46e+03 2.64 0.0083
## SexMen:EF -1.54e-02 5.29e-02 2.07e+03 -0.29 0.7708
## SexMen:RaceAfrAm 6.46e-02 2.78e-01 2.04e+03 0.23 0.8160
## EF:RaceAfrAm 2.88e-02 5.23e-02 2.07e+03 0.55 0.5817
##
## Correlation of Fixed Effects:
## (Intr) Age SexMen EF RcAfrA BMI PvSttB WRATtt Ag:SxM
## Age 0.103
## SexMen -0.223 -0.108
## EF 0.096 0.064 0.085
## RaceAfrAm -0.291 -0.171 0.459 0.163
## BMI -0.479 -0.024 0.039 0.096 -0.067
## PovStatBelw -0.266 0.018 0.025 0.050 -0.026 0.064
## WRATtotal -0.824 -0.024 0.048 -0.258 0.155 -0.017 0.179
## Age:SexMen -0.026 -0.434 0.199 -0.046 -0.018 -0.029 0.011 0.006
## Age:EF -0.033 -0.210 0.019 0.219 0.049 0.016 -0.042 0.029 -0.013
## Age:RcAfrAm -0.100 -0.711 0.016 -0.038 0.269 0.041 0.031 0.032 -0.004
## SexMen:EF 0.120 -0.031 -0.221 -0.459 -0.158 -0.083 -0.008 -0.054 0.119
## SxMn:RcAfrA 0.106 -0.005 -0.750 -0.074 -0.616 0.100 -0.004 -0.040 0.032
## EF:RacAfrAm 0.012 -0.092 0.019 -0.604 -0.046 -0.052 -0.012 0.067 0.010
## Age:EF Ag:RAA SxM:EF SM:RAA
## Age
## SexMen
## EF
## RaceAfrAm
## BMI
## PovStatBelw
## WRATtotal
## Age:SexMen
## Age:EF
## Age:RcAfrAm 0.212
## SexMen:EF -0.021 -0.004
## SxMn:RcAfrA 0.000 -0.012 0.219
## EF:RacAfrAm 0.015 0.111 -0.009 -0.049
Tandem Stands - Model 3 (Three-Way Interactions)
mm3 = lmer(PhyPerfSecsTandem~(Age + Sex + EF + Race)^3 + BMI + PovStat + WRATtotal + (Age|HNDid), data=DanThesisFinal,REML=F)
summary(mm3)
## Linear mixed model fit by maximum likelihood ['merModLmerTest']
## Formula: PhyPerfSecsTandem ~ (Age + Sex + EF + Race)^3 + BMI + PovStat + WRATtotal + (Age | HNDid)
## Data: DanThesisFinal
##
## AIC BIC logLik deviance
## 10772 10897 -5364 10728
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 2.2498 1.500
## Age 0.0148 0.122 1.00
## Residual 5.6665 2.380
## Number of obs: 2168, groups: HNDid, 1440
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 2.83e+01 5.46e-01 1.94e+03 51.91 < 2e-16
## Age -7.44e-02 1.68e-02 1.51e+03 -4.43 9.9e-06
## SexMen 1.45e-01 2.31e-01 1.40e+03 0.63 0.5299
## EF -4.13e-03 5.48e-02 2.04e+03 -0.08 0.9399
## RaceAfrAm 4.84e-01 1.93e-01 1.37e+03 2.51 0.0124
## BMI -1.04e-02 8.71e-03 2.08e+03 -1.20 0.2319
## PovStatBelow -2.55e-01 1.40e-01 2.01e+03 -1.83 0.0680
## WRATtotal 2.70e-02 1.00e-02 1.98e+03 2.70 0.0070
## Age:SexMen 3.22e-02 2.58e-02 1.54e+03 1.25 0.2111
## Age:EF 3.35e-03 5.16e-03 1.71e+03 0.65 0.5155
## Age:RaceAfrAm 5.45e-02 2.11e-02 1.44e+03 2.59 0.0098
## SexMen:EF 1.51e-02 8.08e-02 2.12e+03 0.19 0.8521
## SexMen:RaceAfrAm -2.17e-02 2.98e-01 1.36e+03 -0.07 0.9418
## EF:RaceAfrAm 5.49e-02 7.24e-02 2.10e+03 0.76 0.4478
## Age:SexMen:EF 1.35e-03 6.02e-03 1.67e+03 0.22 0.8224
## Age:SexMen:RaceAfrAm -3.45e-02 3.28e-02 1.49e+03 -1.05 0.2925
## Age:EF:RaceAfrAm 2.53e-03 5.98e-03 1.67e+03 0.42 0.6720
## SexMen:EF:RaceAfrAm -4.61e-02 1.06e-01 2.07e+03 -0.44 0.6632
##
## Correlation of Fixed Effects:
## (Intr) Age SexMen EF RcAfrA BMI PvSttB WRATtt Ag:SxM
## Age 0.097
## SexMen -0.229 -0.201
## EF 0.057 -0.015 0.172
## RaceAfrAm -0.293 -0.241 0.489 0.172
## BMI -0.479 -0.027 0.044 0.084 -0.059
## PovStatBelw -0.268 0.016 0.030 0.062 -0.024 0.065
## WRATtotal -0.822 -0.001 0.040 -0.218 0.144 -0.018 0.179
## Age:SexMen -0.034 -0.625 0.291 -0.006 0.154 0.000 0.004 -0.024
## Age:EF -0.020 -0.351 0.033 0.313 0.050 0.000 -0.015 0.012 0.157
## Age:RcAfrAm -0.091 -0.801 0.162 0.022 0.331 0.043 0.022 0.001 0.505
## SexMen:EF 0.121 -0.026 -0.282 -0.598 -0.147 -0.061 -0.036 -0.049 0.076
## SxMn:RcAfrA 0.116 0.161 -0.772 -0.132 -0.640 0.084 -0.007 -0.027 -0.237
## EF:RacAfrAm 0.043 -0.012 -0.133 -0.703 -0.070 -0.038 -0.030 0.041 0.039
## Age:SxMn:EF 0.011 0.169 -0.007 -0.155 -0.043 -0.003 -0.022 -0.002 -0.221
## Ag:SxMn:RAA 0.024 0.495 -0.229 -0.002 -0.209 -0.023 0.005 0.036 -0.790
## Ag:EF:RcAfA -0.006 0.232 -0.025 -0.210 0.005 0.017 0.001 0.007 -0.045
## SxMn:EF:RAA -0.057 0.056 0.212 0.428 0.056 0.008 0.035 0.017 -0.103
## Age:EF Ag:RAA SxM:EF SM:RAA EF:RAA A:SM:E A:SM:R A:EF:R
## Age
## SexMen
## EF
## RaceAfrAm
## BMI
## PovStatBelw
## WRATtotal
## Age:SexMen
## Age:EF
## Age:RcAfrAm 0.293
## SexMen:EF -0.109 0.010
## SxMn:RcAfrA -0.049 -0.218 0.224
## EF:RacAfrAm -0.171 0.039 0.418 0.056
## Age:SxMn:EF -0.482 -0.165 0.272 0.052 -0.024
## Ag:SxMn:RAA -0.130 -0.630 -0.056 0.324 -0.051 0.187
## Ag:EF:RcAfA -0.662 -0.153 -0.019 0.024 0.306 0.004 0.010
## SxMn:EF:RAA -0.016 -0.058 -0.706 -0.119 -0.621 -0.003 0.123 -0.008
Tandem Stands - Model 4 (Four-Way Interactions)
mm4 = lmer(PhyPerfSecsTandem~(Age + Sex + EF + Race)^4 + BMI + PovStat + WRATtotal + (Age|HNDid), data=DanThesisFinal,REML=F)
summary(mm4)
## Linear mixed model fit by maximum likelihood ['merModLmerTest']
## Formula: PhyPerfSecsTandem ~ (Age + Sex + EF + Race)^4 + BMI + PovStat + WRATtotal + (Age | HNDid)
## Data: DanThesisFinal
##
## AIC BIC logLik deviance
## 10773 10904 -5363 10727
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 2.2593 1.503
## Age 0.0148 0.122 1.00
## Residual 5.6557 2.378
## Number of obs: 2168, groups: HNDid, 1440
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 2.83e+01 5.46e-01 1.94e+03 51.91 < 2e-16
## Age -7.08e-02 1.71e-02 1.55e+03 -4.14 3.7e-05
## SexMen 1.43e-01 2.31e-01 1.40e+03 0.62 0.5355
## EF -1.50e-02 5.56e-02 1.96e+03 -0.27 0.7878
## RaceAfrAm 4.85e-01 1.93e-01 1.37e+03 2.51 0.0122
## BMI -1.05e-02 8.71e-03 2.08e+03 -1.21 0.2262
## PovStatBelow -2.51e-01 1.40e-01 2.01e+03 -1.80 0.0725
## WRATtotal 2.72e-02 1.00e-02 1.98e+03 2.71 0.0067
## Age:SexMen 2.50e-02 2.66e-02 1.56e+03 0.94 0.3477
## Age:EF 1.86e-04 5.90e-03 1.71e+03 0.03 0.9748
## Age:RaceAfrAm 5.15e-02 2.13e-02 1.47e+03 2.42 0.0155
## SexMen:EF 4.28e-02 8.46e-02 1.92e+03 0.51 0.6129
## SexMen:RaceAfrAm -3.34e-02 2.98e-01 1.37e+03 -0.11 0.9108
## EF:RaceAfrAm 7.54e-02 7.47e-02 1.93e+03 1.01 0.3127
## Age:SexMen:EF 8.98e-03 9.17e-03 1.76e+03 0.98 0.3272
## Age:SexMen:RaceAfrAm -2.76e-02 3.33e-02 1.51e+03 -0.83 0.4077
## Age:EF:RaceAfrAm 8.06e-03 7.80e-03 1.62e+03 1.03 0.3015
## SexMen:EF:RaceAfrAm -9.44e-02 1.14e-01 1.91e+03 -0.82 0.4095
## Age:SexMen:EF:RaceAfrAm -1.34e-02 1.22e-02 1.70e+03 -1.10 0.2702
##
## Correlation of Fixed Effects:
## (Intr) Age SexMen EF RcAfrA BMI PvSttB WRATtt Ag:SxM
## Age 0.095
## SexMen -0.229 -0.199
## EF 0.056 -0.048 0.171
## RaceAfrAm -0.293 -0.236 0.489 0.168
## BMI -0.479 -0.029 0.044 0.085 -0.059
## PovStatBelw -0.267 0.020 0.029 0.056 -0.024 0.064
## WRATtotal -0.822 0.002 0.039 -0.217 0.144 -0.018 0.179
## Age:SexMen -0.033 -0.642 0.285 0.037 0.148 0.003 -0.002 -0.027
## Age:EF -0.018 -0.394 0.033 0.355 0.041 0.007 -0.026 0.003 0.253
## Age:RcAfrAm -0.090 -0.804 0.162 0.044 0.329 0.045 0.019 -0.001 0.517
## SexMen:EF 0.116 0.033 -0.271 -0.614 -0.139 -0.062 -0.027 -0.042 -0.003
## SxMn:RcAfrA 0.116 0.151 -0.771 -0.123 -0.640 0.085 -0.008 -0.027 -0.221
## EF:RacAfrAm 0.041 0.036 -0.130 -0.714 -0.066 -0.040 -0.022 0.043 -0.024
## Age:SxMn:EF 0.007 0.253 -0.010 -0.232 -0.024 -0.012 0.005 0.009 -0.327
## Ag:SxMn:RAA 0.023 0.513 -0.228 -0.034 -0.205 -0.025 0.010 0.038 -0.798
## Ag:EF:RcAfA -0.005 0.297 -0.024 -0.272 0.008 0.004 0.017 0.014 -0.191
## SxMn:EF:RAA -0.052 -0.023 0.198 0.457 0.049 0.013 0.023 0.010 0.003
## A:SM:EF:RAA 0.000 -0.192 0.007 0.175 -0.006 0.014 -0.026 -0.014 0.246
## Age:EF Ag:RAA SxM:EF SM:RAA EF:RAA Ag:SM:EF A:SM:R A:EF:R
## Age
## SexMen
## EF
## RaceAfrAm
## BMI
## PovStatBelw
## WRATtotal
## Age:SexMen
## Age:EF
## Age:RcAfrAm 0.316
## SexMen:EF -0.235 -0.028
## SxMn:RcAfrA -0.025 -0.212 0.203
## EF:RacAfrAm -0.266 0.006 0.461 0.045
## Age:SxMn:EF -0.643 -0.204 0.395 0.007 0.172
## Ag:SxMn:RAA -0.202 -0.638 0.003 0.312 -0.003 0.261
## Ag:EF:RcAfA -0.756 -0.198 0.177 -0.005 0.387 0.486 0.127
## SxMn:EF:RAA 0.173 -0.004 -0.736 -0.096 -0.651 -0.292 0.041 -0.252
## A:SM:EF:RAA 0.485 0.128 -0.297 0.036 -0.249 -0.754 -0.186 -0.642
## SM:EF:
## Age
## SexMen
## EF
## RaceAfrAm
## BMI
## PovStatBelw
## WRATtotal
## Age:SexMen
## Age:EF
## Age:RcAfrAm
## SexMen:EF
## SxMn:RcAfrA
## EF:RacAfrAm
## Age:SxMn:EF
## Ag:SxMn:RAA
## Ag:EF:RcAfA
## SxMn:EF:RAA
## A:SM:EF:RAA 0.384
Tandem Stands - Model 5 (Retaining significant higher-order terms)
mm5 = lmer(PhyPerfSecsTandem ~Age + Sex + EF + Race + BMI + PovStat + WRATtotal + (Age|HNDid) + Age*Race, data=DanThesisFinal,REML=F)
summary(mm5)
## Linear mixed model fit by maximum likelihood ['merModLmerTest']
## Formula: PhyPerfSecsTandem ~ Age + Sex + EF + Race + BMI + PovStat + WRATtotal + (Age | HNDid) + Age * Race
## Data: DanThesisFinal
##
## AIC BIC logLik deviance
## 10760 10834 -5367 10734
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 2.2753 1.508
## Age 0.0151 0.123 1.00
## Residual 5.6615 2.379
## Number of obs: 2168, groups: HNDid, 1440
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 2.84e+01 5.39e-01 1.96e+03 52.68 < 2e-16
## Age -5.72e-02 1.22e-02 1.47e+03 -4.68 3.1e-06
## SexMen 9.20e-02 1.35e-01 2.01e+03 0.68 0.4941
## EF 2.18e-03 2.89e-02 2.08e+03 0.08 0.9398
## RaceAfrAm 4.62e-01 1.48e-01 1.37e+03 3.13 0.0018
## BMI -1.10e-02 8.61e-03 2.07e+03 -1.28 0.1999
## PovStatBelow -2.42e-01 1.40e-01 2.02e+03 -1.73 0.0836
## WRATtotal 2.65e-02 9.98e-03 1.99e+03 2.65 0.0080
## Age:RaceAfrAm 3.57e-02 1.56e-02 1.43e+03 2.28 0.0225
##
## Correlation of Fixed Effects:
## (Intr) Age SexMen EF RcAfrA BMI PvSttB WRATtt
## Age 0.098
## SexMen -0.208 -0.019
## EF 0.277 0.088 -0.042
## RaceAfrAm -0.285 -0.250 -0.016 0.104
## BMI -0.490 -0.042 0.182 0.039 -0.015
## PovStatBelw -0.269 0.015 0.031 0.079 -0.035 0.065
## WRATtotal -0.829 -0.009 0.022 -0.409 0.170 -0.015 0.181
## Age:RcAfrAm -0.098 -0.778 0.010 -0.030 0.344 0.046 0.043 0.019