Dan's Thesis Analyses - Chair Stands
load(file="/Users/meganwilliams/Desktop/UMBC 2016/Dan's Thesis/DanThesisFinal.rdata")
library(lmerTest)
Chair Stands - Model 1 (Main Effects)
mm1 = lmer(PhyPerfStands10~Age + Sex + EF + Race + BMI + PovStat + WRATtotal + (Age|HNDid), data=DanThesisFinal,REML=F)
summary(mm1)
## Linear mixed model fit by maximum likelihood ['merModLmerTest']
## Formula: PhyPerfStands10 ~ Age + Sex + EF + Race + BMI + PovStat + WRATtotal + (Age | HNDid)
## Data: DanThesisFinal
##
## AIC BIC logLik deviance
## 14110 14178 -7043 14086
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 23.03798 4.7998
## Age 0.00776 0.0881 1.00
## Residual 36.71407 6.0592
## Number of obs: 2047, groups: HNDid, 1372
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 34.4258 1.5689 1364.0000 21.94 < 2e-16
## Age 0.2131 0.0202 1548.0000 10.54 < 2e-16
## SexMen -1.5129 0.3889 1250.0000 -3.89 0.00011
## EF -0.1864 0.0827 1767.0000 -2.25 0.02432
## RaceAfrAm -0.2556 0.4008 1256.0000 -0.64 0.52377
## BMI 0.0299 0.0258 1353.0000 1.16 0.24680
## PovStatBelow 1.6077 0.4082 1252.0000 3.94 8.7e-05
## WRATtotal -0.0260 0.0286 1394.0000 -0.91 0.36391
##
## Correlation of Fixed Effects:
## (Intr) Age SexMen EF RcAfrA BMI PvSttB
## Age -0.002
## SexMen -0.211 -0.005
## EF 0.262 0.102 -0.020
## RaceAfrAm -0.279 0.030 -0.016 0.115
## BMI -0.514 -0.015 0.186 0.040 -0.015
## PovStatBelw -0.276 0.102 0.028 0.062 -0.048 0.078
## WRATtotal -0.824 0.010 0.020 -0.400 0.180 -0.004 0.188
Chair Stands - Model 2 (Two-Way Interactions)
mm2 = lmer(PhyPerfStands10~(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: PhyPerfStands10 ~ (Age + Sex + EF + Race)^2 + BMI + PovStat + WRATtotal + (Age | HNDid)
## Data: DanThesisFinal
##
## AIC BIC logLik deviance
## 14112 14214 -7038 14076
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 22.76922 4.7717
## Age 0.00676 0.0822 1.00
## Residual 36.68273 6.0566
## Number of obs: 2047, groups: HNDid, 1372
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 3.43e+01 1.59e+00 1.38e+03 21.66 < 2e-16
## Age 1.91e-01 3.67e-02 1.54e+03 5.20 2.3e-07
## SexMen -1.15e+00 6.13e-01 1.29e+03 -1.87 0.062
## EF -9.07e-02 1.33e-01 1.73e+03 -0.68 0.496
## RaceAfrAm -2.09e-02 5.28e-01 1.32e+03 -0.04 0.968
## BMI 2.70e-02 2.61e-02 1.36e+03 1.03 0.302
## PovStatBelow 1.64e+00 4.08e-01 1.26e+03 4.01 6.3e-05
## WRATtotal -2.78e-02 2.87e-02 1.40e+03 -0.97 0.333
## Age:SexMen 1.11e-01 4.11e-02 1.54e+03 2.70 0.007
## Age:EF -5.10e-03 7.86e-03 1.87e+03 -0.65 0.516
## Age:RaceAfrAm -3.23e-02 4.21e-02 1.52e+03 -0.77 0.443
## SexMen:EF 8.29e-03 1.52e-01 1.67e+03 0.05 0.956
## SexMen:RaceAfrAm -4.43e-01 7.99e-01 1.28e+03 -0.55 0.579
## EF:RaceAfrAm -1.74e-01 1.51e-01 1.66e+03 -1.16 0.247
##
## Correlation of Fixed Effects:
## (Intr) Age SexMen EF RcAfrA BMI PvSttB WRATtt Ag:SxM
## Age 0.035
## SexMen -0.223 -0.030
## EF 0.098 0.083 0.108
## RaceAfrAm -0.286 -0.023 0.492 0.179
## BMI -0.503 -0.029 0.045 0.087 -0.068
## PovStatBelw -0.277 0.034 0.030 0.050 -0.023 0.078
## WRATtotal -0.819 -0.016 0.048 -0.267 0.160 -0.005 0.187
## Age:SexMen 0.007 -0.432 0.023 -0.027 -0.016 -0.025 0.015 -0.002
## Age:EF -0.001 -0.219 0.032 0.074 0.012 0.013 -0.047 0.004 0.002
## Age:RcAfrAm -0.053 -0.706 0.017 -0.064 0.057 0.044 0.026 0.029 -0.011
## SexMen:EF 0.126 -0.020 -0.224 -0.467 -0.161 -0.087 -0.011 -0.055 0.117
## SxMn:RcAfrA 0.108 -0.008 -0.768 -0.090 -0.643 0.103 -0.014 -0.044 0.037
## EF:RacAfrAm 0.003 -0.088 0.002 -0.609 -0.069 -0.033 -0.017 0.075 -0.008
## Age:EF Ag:RAA SxM:EF SM:RAA
## Age
## SexMen
## EF
## RaceAfrAm
## BMI
## PovStatBelw
## WRATtotal
## Age:SexMen
## Age:EF
## Age:RcAfrAm 0.211
## SexMen:EF 0.001 -0.016
## SxMn:RcAfrA -0.010 -0.006 0.207
## EF:RacAfrAm 0.023 0.111 -0.030 -0.028
Chair Stands - Model 3 (Three-Way Interactions)
mm3 = lmer(PhyPerfStands10~(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: PhyPerfStands10 ~ (Age + Sex + EF + Race)^3 + BMI + PovStat + WRATtotal + (Age | HNDid)
## Data: DanThesisFinal
##
## AIC BIC logLik deviance
## 14118 14242 -7037 14074
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 22.65756 4.7600
## Age 0.00682 0.0826 1.00
## Residual 36.70136 6.0582
## Number of obs: 2047, groups: HNDid, 1372
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 3.43e+01 1.59e+00 1.38e+03 21.57 < 2e-16
## Age 1.74e-01 4.39e-02 1.59e+03 3.96 7.8e-05
## SexMen -1.08e+00 6.31e-01 1.33e+03 -1.71 0.087
## EF -3.29e-02 1.50e-01 1.75e+03 -0.22 0.826
## RaceAfrAm 2.56e-02 5.30e-01 1.33e+03 0.05 0.962
## BMI 2.77e-02 2.61e-02 1.36e+03 1.06 0.289
## PovStatBelow 1.65e+00 4.08e-01 1.26e+03 4.05 5.5e-05
## WRATtotal -2.79e-02 2.87e-02 1.40e+03 -0.97 0.331
## Age:SexMen 1.47e-01 6.76e-02 1.52e+03 2.17 0.030
## Age:EF 3.09e-03 1.36e-02 1.83e+03 0.23 0.821
## Age:RaceAfrAm -8.15e-03 5.53e-02 1.57e+03 -0.15 0.883
## SexMen:EF -1.31e-01 2.24e-01 1.74e+03 -0.58 0.559
## SexMen:RaceAfrAm -5.40e-01 8.09e-01 1.34e+03 -0.67 0.504
## EF:RaceAfrAm -2.60e-01 2.01e-01 1.81e+03 -1.29 0.197
## Age:SexMen:EF -2.07e-02 1.59e-02 1.87e+03 -1.30 0.194
## Age:SexMen:RaceAfrAm -4.95e-02 8.59e-02 1.53e+03 -0.58 0.565
## Age:EF:RaceAfrAm 9.86e-04 1.58e-02 1.87e+03 0.06 0.950
## SexMen:EF:RaceAfrAm 2.02e-01 3.03e-01 1.66e+03 0.66 0.506
##
## Correlation of Fixed Effects:
## (Intr) Age SexMen EF RcAfrA BMI PvSttB WRATtt Ag:SxM
## Age 0.014
## SexMen -0.231 -0.035
## EF 0.057 0.044 0.196
## RaceAfrAm -0.289 -0.032 0.494 0.193
## BMI -0.503 -0.024 0.045 0.080 -0.065
## PovStatBelw -0.278 0.035 0.034 0.057 -0.021 0.078
## WRATtotal -0.819 0.006 0.052 -0.225 0.160 -0.005 0.188
## Age:SexMen 0.021 -0.625 0.018 -0.028 0.020 -0.007 -0.001 -0.030
## Age:EF 0.021 -0.363 0.008 0.105 -0.013 -0.009 -0.025 -0.008 0.173
## Age:RcAfrAm -0.026 -0.798 0.028 -0.029 0.066 0.037 0.012 -0.001 0.502
## SexMen:EF 0.137 -0.044 -0.305 -0.612 -0.166 -0.068 -0.031 -0.058 0.127
## SxMn:RcAfrA 0.114 0.032 -0.772 -0.144 -0.645 0.100 -0.017 -0.044 -0.024
## EF:RacAfrAm 0.045 -0.044 -0.154 -0.706 -0.101 -0.028 -0.030 0.039 0.038
## Age:SxMn:EF -0.005 0.179 0.058 -0.057 -0.007 -0.009 -0.014 0.008 -0.217
## Ag:SxMn:RAA -0.019 0.495 -0.013 0.019 -0.041 -0.012 0.014 0.036 -0.790
## Ag:EF:RcAfA -0.028 0.236 -0.030 -0.068 0.037 0.032 0.007 0.007 -0.059
## SxMn:EF:RAA -0.071 0.050 0.228 0.449 0.079 0.011 0.029 0.029 -0.113
## 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.300
## SexMen:EF -0.026 0.030
## SxMn:RcAfrA -0.029 -0.051 0.233
## EF:RacAfrAm -0.048 0.059 0.454 0.075
## Age:SxMn:EF -0.495 -0.170 0.105 -0.002 -0.021
## Ag:SxMn:RAA -0.144 -0.631 -0.098 0.061 -0.046 0.186
## Ag:EF:RcAfA -0.654 -0.159 -0.023 0.036 0.116 0.010 0.026
## SxMn:EF:RAA -0.025 -0.053 -0.727 -0.135 -0.651 0.012 0.127 0.010
Chair Stands - Model 4 (Four-Way Interactions)
mm4 = lmer(PhyPerfStands10~(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: PhyPerfStands10 ~ (Age + Sex + EF + Race)^4 + BMI + PovStat + WRATtotal + (Age | HNDid)
## Data: DanThesisFinal
##
## AIC BIC logLik deviance
## 14120 14249 -7037 14074
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 22.65760 4.7600
## Age 0.00702 0.0838 1.00
## Residual 36.68631 6.0569
## Number of obs: 2047, groups: HNDid, 1372
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 3.43e+01 1.59e+00 1.38e+03 21.58 < 2e-16
## Age 1.79e-01 4.48e-02 1.61e+03 4.01 6.4e-05
## SexMen -1.06e+00 6.33e-01 1.33e+03 -1.67 0.095
## EF -3.82e-02 1.50e-01 1.77e+03 -0.25 0.799
## RaceAfrAm 3.58e-02 5.30e-01 1.33e+03 0.07 0.946
## BMI 2.75e-02 2.61e-02 1.36e+03 1.05 0.292
## PovStatBelow 1.66e+00 4.08e-01 1.26e+03 4.07 5.1e-05
## WRATtotal -2.80e-02 2.87e-02 1.40e+03 -0.97 0.330
## Age:SexMen 1.36e-01 6.96e-02 1.54e+03 1.96 0.050
## Age:EF -1.56e-03 1.56e-02 1.81e+03 -0.10 0.920
## Age:RaceAfrAm -1.27e-02 5.58e-02 1.59e+03 -0.23 0.820
## SexMen:EF -1.15e-01 2.26e-01 1.82e+03 -0.51 0.612
## SexMen:RaceAfrAm -5.86e-01 8.12e-01 1.34e+03 -0.72 0.471
## EF:RaceAfrAm -2.48e-01 2.02e-01 1.84e+03 -1.22 0.221
## Age:SexMen:EF -9.64e-03 2.39e-02 1.88e+03 -0.40 0.687
## Age:SexMen:RaceAfrAm -4.01e-02 8.72e-02 1.56e+03 -0.46 0.646
## Age:EF:RaceAfrAm 9.29e-03 2.08e-02 1.86e+03 0.45 0.654
## SexMen:EF:RaceAfrAm 1.70e-01 3.07e-01 1.85e+03 0.55 0.580
## Age:SexMen:EF:RaceAfrAm -1.99e-02 3.21e-02 1.87e+03 -0.62 0.536
##
## Correlation of Fixed Effects:
## (Intr) Age SexMen EF RcAfrA BMI PvSttB WRATtt Ag:SxM
## Age 0.015
## SexMen -0.231 -0.023
## EF 0.057 0.032 0.192
## RaceAfrAm -0.289 -0.027 0.495 0.191
## BMI -0.503 -0.026 0.044 0.081 -0.065
## PovStatBelw -0.278 0.040 0.035 0.055 -0.020 0.078
## WRATtotal -0.819 0.005 0.052 -0.225 0.160 -0.005 0.188
## Age:SexMen 0.020 -0.642 0.004 -0.014 0.013 -0.004 -0.007 -0.029
## Age:EF 0.018 -0.408 -0.024 0.121 -0.027 -0.002 -0.035 -0.006 0.262
## Age:RcAfrAm -0.026 -0.801 0.020 -0.021 0.063 0.038 0.008 -0.001 0.515
## SexMen:EF 0.136 -0.019 -0.294 -0.614 -0.161 -0.069 -0.027 -0.058 0.094
## SxMn:RcAfrA 0.114 0.014 -0.773 -0.138 -0.645 0.101 -0.019 -0.044 -0.004
## EF:RacAfrAm 0.045 -0.023 -0.147 -0.707 -0.097 -0.029 -0.027 0.038 0.014
## Age:SxMn:EF -0.002 0.265 0.086 -0.081 0.019 -0.015 0.011 0.004 -0.317
## Ag:SxMn:RAA -0.018 0.513 -0.003 0.008 -0.036 -0.014 0.018 0.035 -0.797
## Ag:EF:RcAfA -0.020 0.305 0.019 -0.090 0.048 0.016 0.023 0.004 -0.196
## SxMn:EF:RAA -0.070 0.016 0.214 0.452 0.073 0.013 0.024 0.029 -0.069
## A:SM:EF:RAA -0.001 -0.198 -0.064 0.057 -0.031 0.012 -0.027 0.003 0.236
## 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.326
## SexMen:EF -0.081 0.014
## SxMn:RcAfrA 0.019 -0.039 0.220
## EF:RacAfrAm -0.090 0.044 0.460 0.066
## Age:SxMn:EF -0.650 -0.212 0.159 -0.069 0.060
## Ag:SxMn:RAA -0.209 -0.639 -0.075 0.046 -0.028 0.253
## Ag:EF:RcAfA -0.750 -0.207 0.060 -0.032 0.153 0.487 0.133
## SxMn:EF:RAA 0.059 -0.030 -0.732 -0.118 -0.655 -0.116 0.094 -0.101
## A:SM:EF:RAA 0.485 0.135 -0.119 0.091 -0.099 -0.746 -0.175 -0.647
## 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.167
Chair Stands - Model 5 (Retaining significant higher-order terms)
mm5 = lmer(PhyPerfStands10~Age + Sex + EF + Race + BMI + PovStat + WRATtotal + (Age|HNDid) + Age*Sex, data=DanThesisFinal,REML=F)
summary(mm5)
## Linear mixed model fit by maximum likelihood ['merModLmerTest']
## Formula: PhyPerfStands10 ~ Age + Sex + EF + Race + BMI + PovStat + WRATtotal + (Age | HNDid) + Age * Sex
## Data: DanThesisFinal
##
## AIC BIC logLik deviance
## 14105 14178 -7039 14079
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 22.789 4.7738
## Age 0.007 0.0837 1.00
## Residual 36.724 6.0600
## Number of obs: 2047, groups: HNDid, 1372
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 34.3769 1.5657 1369.0000 21.96 < 2e-16
## Age 0.1691 0.0258 1529.0000 6.55 7.7e-11
## SexMen -1.4029 0.3899 1310.0000 -3.60 0.00033
## EF -0.1778 0.0826 1772.0000 -2.15 0.03159
## RaceAfrAm -0.2405 0.4000 1260.0000 -0.60 0.54783
## BMI 0.0288 0.0258 1358.0000 1.12 0.26403
## PovStatBelow 1.6245 0.4075 1258.0000 3.99 7.1e-05
## WRATtotal -0.0255 0.0286 1396.0000 -0.89 0.37194
## Age:SexMen 0.1105 0.0408 1528.0000 2.71 0.00687
##
## Correlation of Fixed Effects:
## (Intr) Age SexMen EF RcAfrA BMI PvSttB WRATtt
## Age 0.004
## SexMen -0.211 -0.064
## EF 0.261 0.054 -0.016
## RaceAfrAm -0.279 0.013 -0.014 0.115
## BMI -0.514 -0.001 0.184 0.039 -0.015
## PovStatBelw -0.276 0.069 0.030 0.062 -0.048 0.078
## WRATtotal -0.824 0.004 0.021 -0.399 0.180 -0.004 0.188
## Age:SexMen -0.010 -0.624 0.096 0.040 0.015 -0.017 0.017 0.006