Dan's Thesis Analyses - Single Leg Stand (Mean)
load(file="/Users/meganwilliams/Desktop/UMBC 2016/Dan's Thesis/DanThesisFinal.rdata")
library(lmerTest)
Single Leg Stand (Mean) - Model 1 (Main Effects)
mm1 = lmer(PhyPerfSingleMN~Age + Sex + EF + Race + BMI + PovStat + WRATtotal + (Age|HNDid), data=DanThesisFinal,REML=F)
summary(mm1)
## Linear mixed model fit by maximum likelihood ['merModLmerTest']
## Formula: PhyPerfSingleMN ~ Age + Sex + EF + Race + BMI + PovStat + WRATtotal + (Age | HNDid)
## Data: DanThesisFinal
##
## AIC BIC logLik deviance
## 8209 8273 -4092 8185
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 6.0891 2.468
## Age 0.0297 0.172 1.00
## Residual 6.0614 2.462
## Number of obs: 1569, groups: HNDid, 1112
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 29.3511 0.7345 1325.0000 39.96 < 2e-16
## Age -0.1003 0.0109 1119.0000 -9.21 < 2e-16
## SexMen 0.1947 0.1750 1345.0000 1.11 0.266
## EF 0.0774 0.0380 1391.0000 2.04 0.042
## RaceAfrAm 0.4032 0.1852 1375.0000 2.18 0.030
## BMI -0.0856 0.0122 1330.0000 -7.03 3.4e-12
## PovStatBelow -0.3005 0.1869 1374.0000 -1.61 0.108
## WRATtotal 0.0288 0.0135 1391.0000 2.13 0.033
##
## Correlation of Fixed Effects:
## (Intr) Age SexMen EF RcAfrA BMI PvSttB
## Age 0.069
## SexMen -0.179 -0.056
## EF 0.260 0.092 -0.045
## RaceAfrAm -0.283 0.019 -0.049 0.135
## BMI -0.493 0.028 0.147 0.046 -0.030
## PovStatBelw -0.267 0.085 0.032 0.108 -0.065 0.085
## WRATtotal -0.831 -0.001 0.005 -0.406 0.200 -0.008 0.178
Single Leg Stand (Mean) - Model 2 (Two-Way Interactions)
mm2 = lmer(PhyPerfSingleMN~(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: PhyPerfSingleMN ~ (Age + Sex + EF + Race)^2 + BMI + PovStat + WRATtotal + (Age | HNDid)
## Data: DanThesisFinal
##
## AIC BIC logLik deviance
## 8205 8302 -4085 8169
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 6.0273 2.455
## Age 0.0294 0.171 1.00
## Residual 6.0056 2.451
## Number of obs: 1569, groups: HNDid, 1112
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 2.88e+01 7.55e-01 1.33e+03 38.07 < 2e-16
## Age -1.45e-01 2.02e-02 1.14e+03 -7.19 1.2e-12
## SexMen 5.27e-01 3.07e-01 1.24e+03 1.71 0.08675
## EF 1.98e-01 6.77e-02 1.50e+03 2.93 0.00344
## RaceAfrAm 9.67e-01 2.87e-01 1.13e+03 3.37 0.00077
## BMI -8.43e-02 1.23e-02 1.33e+03 -6.85 1.1e-11
## PovStatBelow -3.00e-01 1.87e-01 1.37e+03 -1.61 0.10841
## WRATtotal 3.28e-02 1.35e-02 1.38e+03 2.43 0.01524
## Age:SexMen 4.35e-03 2.18e-02 1.12e+03 0.20 0.84191
## Age:EF 1.12e-02 4.31e-03 1.17e+03 2.61 0.00922
## Age:RaceAfrAm 6.28e-02 2.26e-02 1.12e+03 2.77 0.00566
## SexMen:EF -1.27e-01 6.88e-02 1.36e+03 -1.85 0.06420
## SexMen:RaceAfrAm -3.80e-01 3.65e-01 1.39e+03 -1.04 0.29888
## EF:RaceAfrAm 8.44e-03 6.83e-02 1.36e+03 0.12 0.90158
##
## Correlation of Fixed Effects:
## (Intr) Age SexMen EF RcAfrA BMI PvSttB WRATtt Ag:SxM
## Age 0.178
## SexMen -0.246 -0.242
## EF 0.031 -0.017 0.109
## RaceAfrAm -0.335 -0.339 0.415 0.206
## BMI -0.473 0.036 0.017 0.102 -0.062
## PovStatBelw -0.267 0.017 0.043 0.055 -0.012 0.084
## WRATtotal -0.820 -0.043 0.063 -0.212 0.181 -0.011 0.176
## Age:SexMen -0.075 -0.470 0.388 -0.034 -0.008 -0.040 0.024 0.019
## Age:EF -0.063 -0.233 0.033 0.439 0.115 0.002 -0.052 0.038 0.008
## Age:RcAfrAm -0.145 -0.696 0.049 0.029 0.496 -0.005 0.032 0.042 0.000
## SexMen:EF 0.145 -0.017 -0.257 -0.467 -0.174 -0.105 -0.006 -0.059 0.093
## SxMn:RcAfrA 0.139 0.036 -0.707 -0.085 -0.594 0.091 -0.020 -0.067 0.024
## EF:RacAfrAm 0.031 -0.070 0.034 -0.542 -0.063 -0.051 -0.023 0.052 0.022
## Age:EF Ag:RAA SxM:EF SM:RAA
## Age
## SexMen
## EF
## RaceAfrAm
## BMI
## PovStatBelw
## WRATtotal
## Age:SexMen
## Age:EF
## Age:RcAfrAm 0.216
## SexMen:EF -0.077 -0.003
## SxMn:RcAfrA -0.001 -0.059 0.253
## EF:RacAfrAm -0.024 0.080 -0.047 -0.069
Single Leg Stand (Mean) - Model 3 (Three-Way Interactions)
mm3 = lmer(PhyPerfSingleMN~(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: PhyPerfSingleMN ~ (Age + Sex + EF + Race)^3 + BMI + PovStat + WRATtotal + (Age | HNDid)
## Data: DanThesisFinal
##
## AIC BIC logLik deviance
## 8212 8330 -4084 8168
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 6.0019 2.450
## Age 0.0296 0.172 1.00
## Residual 6.0046 2.450
## Number of obs: 1569, groups: HNDid, 1112
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 2.88e+01 7.57e-01 1.33e+03 38.05 < 2e-16
## Age -1.44e-01 2.51e-02 1.15e+03 -5.75 1.1e-08
## SexMen 4.12e-01 3.56e-01 8.72e+02 1.16 0.2479
## EF 1.50e-01 8.66e-02 1.42e+03 1.73 0.0839
## RaceAfrAm 9.23e-01 3.16e-01 8.75e+02 2.92 0.0036
## BMI -8.39e-02 1.23e-02 1.33e+03 -6.82 1.3e-11
## PovStatBelow -3.08e-01 1.87e-01 1.36e+03 -1.65 0.0996
## WRATtotal 3.27e-02 1.35e-02 1.38e+03 2.42 0.0156
## Age:SexMen 3.33e-03 3.61e-02 1.17e+03 0.09 0.9264
## Age:EF 9.76e-03 7.84e-03 1.17e+03 1.24 0.2134
## Age:RaceAfrAm 6.19e-02 3.14e-02 1.09e+03 1.97 0.0487
## SexMen:EF -8.36e-03 1.15e-01 1.52e+03 -0.07 0.9423
## SexMen:RaceAfrAm -2.73e-01 4.57e-01 8.49e+02 -0.60 0.5507
## EF:RaceAfrAm 7.67e-02 1.11e-01 1.50e+03 0.69 0.4896
## Age:SexMen:EF 4.88e-03 8.67e-03 1.16e+03 0.56 0.5739
## Age:SexMen:RaceAfrAm -8.05e-04 4.59e-02 1.13e+03 -0.02 0.9860
## Age:EF:RaceAfrAm -9.75e-04 8.82e-03 1.17e+03 -0.11 0.9120
## SexMen:EF:RaceAfrAm -1.65e-01 1.37e-01 1.35e+03 -1.20 0.2293
##
## Correlation of Fixed Effects:
## (Intr) Age SexMen EF RcAfrA BMI PvSttB WRATtt Ag:SxM
## Age 0.174
## SexMen -0.242 -0.408
## EF -0.021 -0.190 0.233
## RaceAfrAm -0.326 -0.483 0.530 0.261
## BMI -0.469 0.045 -0.001 0.065 -0.069
## PovStatBelw -0.267 0.025 0.029 0.062 -0.020 0.083
## WRATtotal -0.814 0.001 0.027 -0.166 0.136 -0.009 0.178
## Age:SexMen -0.075 -0.667 0.545 0.096 0.321 -0.048 -0.009 -0.044
## Age:EF -0.077 -0.414 0.131 0.609 0.195 -0.014 -0.012 0.012 0.216
## Age:RcAfrAm -0.137 -0.805 0.331 0.166 0.610 -0.023 0.005 -0.015 0.541
## SexMen:EF 0.127 0.062 -0.328 -0.599 -0.193 -0.041 -0.033 -0.039 -0.030
## SxMn:RcAfrA 0.142 0.331 -0.784 -0.187 -0.682 0.091 -0.002 -0.015 -0.439
## EF:RacAfrAm 0.077 0.095 -0.152 -0.677 -0.135 -0.024 -0.030 0.025 -0.007
## Age:SxMn:EF 0.017 0.205 -0.117 -0.307 -0.115 0.028 -0.027 0.012 -0.240
## Ag:SxMn:RAA 0.033 0.527 -0.428 -0.083 -0.402 0.027 0.033 0.070 -0.790
## Ag:EF:RcAfA 0.050 0.267 -0.060 -0.398 -0.093 0.000 -0.006 -0.001 -0.077
## SxMn:EF:RAA -0.055 0.024 0.232 0.388 0.081 -0.013 0.026 0.014 -0.063
## 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.345
## SexMen:EF -0.263 -0.064
## SxMn:RcAfrA -0.123 -0.422 0.271
## EF:RacAfrAm -0.344 -0.056 0.336 0.086
## Age:SxMn:EF -0.496 -0.194 0.493 0.132 -0.031
## Ag:SxMn:RAA -0.175 -0.667 0.027 0.568 -0.013 0.198
## Ag:EF:RcAfA -0.648 -0.189 -0.006 0.046 0.546 -0.044 0.035
## SxMn:EF:RAA 0.032 -0.027 -0.651 -0.150 -0.599 -0.035 0.087 -0.059
Single Leg Stand (Mean) - Model 4 (Four-Way Interactions)
mm4 = lmer(PhyPerfSingleMN~(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: PhyPerfSingleMN ~ (Age + Sex + EF + Race)^4 + BMI + PovStat + WRATtotal + (Age | HNDid)
## Data: DanThesisFinal
##
## AIC BIC logLik deviance
## 8214 8337 -4084 8168
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 6.0076 2.451
## Age 0.0296 0.172 1.00
## Residual 6.0017 2.450
## Number of obs: 1569, groups: HNDid, 1112
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 2.88e+01 7.57e-01 1.33e+03 38.03 < 2e-16
## Age -1.43e-01 2.59e-02 1.16e+03 -5.53 4.0e-08
## SexMen 4.01e-01 3.60e-01 8.61e+02 1.11 0.266
## EF 1.42e-01 9.28e-02 1.32e+03 1.53 0.126
## RaceAfrAm 9.17e-01 3.17e-01 8.76e+02 2.89 0.004
## BMI -8.41e-02 1.23e-02 1.33e+03 -6.82 1.3e-11
## PovStatBelow -3.07e-01 1.87e-01 1.36e+03 -1.64 0.100
## WRATtotal 3.29e-02 1.35e-02 1.38e+03 2.43 0.015
## Age:SexMen 8.81e-04 3.76e-02 1.16e+03 0.02 0.981
## Age:EF 8.65e-03 9.17e-03 1.16e+03 0.94 0.346
## Age:RaceAfrAm 6.06e-02 3.18e-02 1.10e+03 1.90 0.057
## SexMen:EF 7.96e-03 1.35e-01 1.33e+03 0.06 0.953
## SexMen:RaceAfrAm -2.65e-01 4.58e-01 8.45e+02 -0.58 0.563
## EF:RaceAfrAm 8.99e-02 1.25e-01 1.31e+03 0.72 0.471
## Age:SexMen:EF 7.36e-03 1.37e-02 1.25e+03 0.54 0.591
## Age:SexMen:RaceAfrAm 1.59e-03 4.70e-02 1.13e+03 0.03 0.973
## Age:EF:RaceAfrAm 9.60e-04 1.21e-02 1.09e+03 0.08 0.937
## SexMen:EF:RaceAfrAm -1.91e-01 1.77e-01 1.34e+03 -1.08 0.282
## Age:SexMen:EF:RaceAfrAm -4.14e-03 1.77e-02 1.20e+03 -0.23 0.815
##
## Correlation of Fixed Effects:
## (Intr) Age SexMen EF RcAfrA BMI PvSttB WRATtt Ag:SxM
## Age 0.178
## SexMen -0.245 -0.426
## EF -0.034 -0.258 0.265
## RaceAfrAm -0.327 -0.486 0.534 0.271
## BMI -0.470 0.030 0.006 0.081 -0.065
## PovStatBelw -0.266 0.028 0.027 0.052 -0.021 0.082
## WRATtotal -0.811 0.011 0.021 -0.169 0.133 -0.012 0.179
## Age:SexMen -0.082 -0.689 0.558 0.186 0.329 -0.030 -0.013 -0.053
## Age:EF -0.086 -0.468 0.183 0.672 0.207 0.017 -0.019 -0.010 0.323
## Age:RcAfrAm -0.141 -0.812 0.347 0.215 0.613 -0.013 0.002 -0.022 0.560
## SexMen:EF 0.129 0.177 -0.350 -0.663 -0.205 -0.064 -0.020 -0.012 -0.170
## SxMn:RcAfrA 0.144 0.338 -0.785 -0.201 -0.683 0.086 -0.001 -0.012 -0.441
## EF:RacAfrAm 0.086 0.192 -0.198 -0.726 -0.155 -0.048 -0.019 0.040 -0.134
## Age:SxMn:EF 0.041 0.313 -0.181 -0.458 -0.133 -0.026 -0.004 0.038 -0.362
## Ag:SxMn:RAA 0.041 0.551 -0.445 -0.154 -0.409 0.014 0.036 0.077 -0.801
## Ag:EF:RcAfA 0.063 0.354 -0.139 -0.515 -0.121 -0.039 0.007 0.026 -0.245
## SxMn:EF:RAA -0.067 -0.135 0.266 0.507 0.112 0.026 0.009 -0.015 0.131
## A:SM:EF:RAA -0.039 -0.242 0.140 0.357 0.078 0.057 -0.017 -0.040 0.279
## 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.380
## SexMen:EF -0.461 -0.144
## SxMn:RcAfrA -0.143 -0.427 0.269
## EF:RacAfrAm -0.499 -0.128 0.493 0.109
## Age:SxMn:EF -0.670 -0.255 0.668 0.140 0.337
## Ag:SxMn:RAA -0.259 -0.679 0.135 0.569 0.089 0.291
## Ag:EF:RcAfA -0.759 -0.255 0.351 0.084 0.667 0.510 0.174
## SxMn:EF:RAA 0.351 0.090 -0.759 -0.162 -0.702 -0.509 -0.073 -0.468
## A:SM:EF:RAA 0.519 0.174 -0.518 -0.073 -0.457 -0.774 -0.218 -0.684
## 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.636
Single Leg Stand (Mean) - Model 5 (Retaining significant higher-order terms)
mm5 = lmer(PhyPerfSingleMN~Age + Sex + EF + Race + BMI + PovStat + WRATtotal + (Age|HNDid) + Age*Race + Age*EF ,data=DanThesisFinal,REML=F)
summary(mm5)
## Linear mixed model fit by maximum likelihood ['merModLmerTest']
## Formula: PhyPerfSingleMN ~ Age + Sex + EF + Race + BMI + PovStat + WRATtotal + (Age | HNDid) + Age * Race + Age * EF
## Data: DanThesisFinal
##
## AIC BIC logLik deviance
## 8201 8276 -4087 8173
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 6.0345 2.457
## Age 0.0292 0.171 1.00
## Residual 6.0267 2.455
## Number of obs: 1569, groups: HNDid, 1112
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 29.0309 0.7403 1340.0000 39.22 < 2e-16
## Age -0.1418 0.0178 1163.0000 -7.98 3.6e-15
## SexMen 0.2128 0.1746 1343.0000 1.22 0.22312
## EF 0.1409 0.0461 1441.0000 3.05 0.00230
## RaceAfrAm 0.7772 0.2287 838.0000 3.40 0.00071
## BMI -0.0857 0.0121 1327.0000 -7.06 2.7e-12
## PovStatBelow -0.3060 0.1868 1368.0000 -1.64 0.10162
## WRATtotal 0.0307 0.0135 1388.0000 2.28 0.02259
## Age:RaceAfrAm 0.0619 0.0226 1124.0000 2.74 0.00619
## Age:EF 0.0107 0.0043 1157.0000 2.48 0.01331
##
## Correlation of Fixed Effects:
## (Intr) Age SexMen EF RcAfrA BMI PvSttB WRATtt Ag:RAA
## Age 0.158
## SexMen -0.181 -0.051
## EF 0.184 -0.089 -0.018
## RaceAfrAm -0.312 -0.456 -0.029 0.157
## BMI -0.488 0.013 0.147 0.033 -0.021
## PovStatBelw -0.267 0.032 0.031 0.057 -0.033 0.086
## WRATtotal -0.827 -0.031 0.007 -0.313 0.183 -0.008 0.177
## Age:RcAfrAm -0.146 -0.787 0.017 0.108 0.590 0.006 0.033 0.036
## Age:EF -0.052 -0.262 0.033 0.570 0.137 -0.009 -0.053 0.036 0.221