Dan's Thesis Analyses - Grip Strength (Right)
load(file="/Users/meganwilliams/Desktop/UMBC 2016/Dan's Thesis/DanThesisFinal.rdata")
library(lmerTest)
Right Grip Strength - Model 1 (Main Effects)
mm1 = lmer(PhyPerfRight~Age + Sex + EF + Race + BMI + PovStat + WRATtotal + DomHand + (Age|HNDid), data=DanThesisFinal,REML=F)
summary(mm1)
## Linear mixed model fit by maximum likelihood ['merModLmerTest']
## Formula: PhyPerfRight ~ Age + Sex + EF + Race + BMI + PovStat + WRATtotal + DomHand + (Age | HNDid)
## Data: DanThesisFinal
##
## AIC BIC logLik deviance
## 15652 15732 -7812 15624
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 34.18190 5.847
## Age 0.00231 0.048 -1.00
## Residual 29.81064 5.460
## Number of obs: 2273, groups: HNDid, 1489
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 22.2124 1.7686 1881.0000 12.56 < 2e-16
## Age -0.1857 0.0201 1796.0000 -9.23 < 2e-16
## SexMen 15.7058 0.3985 1444.0000 39.41 < 2e-16
## EF 0.3585 0.0805 2247.0000 4.45 8.8e-06
## RaceAfrAm 2.1611 0.4029 1425.0000 5.36 9.5e-08
## BMI 0.0853 0.0247 1679.0000 3.45 0.00057
## PovStatBelow -1.3955 0.4098 1446.0000 -3.41 0.00068
## WRATtotal 0.0569 0.0282 1577.0000 2.02 0.04376
## DomHandLeft -0.8075 1.0139 2238.0000 -0.80 0.42588
## DomHandRight 0.3894 0.8916 2168.0000 0.44 0.66236
##
## Correlation of Fixed Effects:
## (Intr) Age SexMen EF RcAfrA BMI PvSttB WRATtt DmHndL
## Age -0.031
## SexMen -0.224 0.010
## EF 0.219 0.071 -0.011
## RaceAfrAm -0.253 0.029 -0.007 0.091
## BMI -0.458 -0.031 0.193 0.034 0.008
## PovStatBelw -0.235 0.082 0.057 0.064 -0.036 0.050
## WRATtotal -0.708 0.020 0.022 -0.387 0.185 -0.007 0.172
## DomHandLeft -0.444 0.006 0.033 0.008 -0.008 0.045 0.020 0.005
## DomHandRght -0.488 -0.009 0.069 0.013 -0.015 0.030 0.018 -0.020 0.820
Right Grip Strength - Model 2 (Two-Way Interactions)
mm2 = lmer(PhyPerfRight~(Age + Sex + EF + Race)^2 + BMI + PovStat + WRATtotal + DomHand + (Age|HNDid), data=DanThesisFinal,REML=F)
summary(mm2)
## Linear mixed model fit by maximum likelihood ['merModLmerTest']
## Formula: PhyPerfRight ~ (Age + Sex + EF + Race)^2 + BMI + PovStat + WRATtotal + DomHand + (Age | HNDid)
## Data: DanThesisFinal
##
## AIC BIC logLik deviance
## 15654 15769 -7807 15614
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 33.68114 5.8035
## Age 0.00241 0.0491 -1.00
## Residual 29.89537 5.4677
## Number of obs: 2273, groups: HNDid, 1489
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 2.23e+01 1.78e+00 1.88e+03 12.54 < 2e-16
## Age -2.22e-01 3.60e-02 1.79e+03 -6.18 7.9e-10
## SexMen 1.55e+01 6.22e-01 1.45e+03 24.97 < 2e-16
## EF 2.88e-01 1.29e-01 2.15e+03 2.23 0.02564
## RaceAfrAm 1.83e+00 5.27e-01 1.45e+03 3.47 0.00054
## BMI 8.58e-02 2.49e-02 1.68e+03 3.45 0.00058
## PovStatBelow -1.36e+00 4.09e-01 1.44e+03 -3.33 0.00088
## WRATtotal 5.59e-02 2.82e-02 1.58e+03 1.98 0.04743
## DomHandLeft -7.88e-01 1.01e+00 2.24e+03 -0.78 0.43681
## DomHandRight 4.34e-01 8.93e-01 2.17e+03 0.49 0.62679
## Age:SexMen -4.09e-02 4.11e-02 1.81e+03 -1.00 0.31975
## Age:EF -7.76e-03 7.57e-03 2.21e+03 -1.03 0.30498
## Age:RaceAfrAm 9.11e-02 4.18e-02 1.80e+03 2.18 0.02936
## SexMen:EF 1.90e-01 1.49e-01 2.18e+03 1.27 0.20306
## SexMen:RaceAfrAm 3.41e-01 8.03e-01 1.44e+03 0.43 0.67077
## EF:RaceAfrAm -1.31e-02 1.48e-01 2.17e+03 -0.09 0.92911
##
## Correlation of Fixed Effects:
## (Intr) Age SexMen EF RcAfrA BMI PvSttB WRATtt DmHndL
## Age -0.022
## SexMen -0.230 0.064
## EF 0.091 0.088 0.082
## RaceAfrAm -0.253 0.140 0.469 0.149
## BMI -0.453 -0.049 0.082 0.087 -0.039
## PovStatBelw -0.234 0.010 0.038 0.056 -0.027 0.052
## WRATtotal -0.705 -0.011 0.038 -0.263 0.153 -0.007 0.172
## DomHandLeft -0.443 0.011 0.031 0.005 -0.006 0.047 0.019 0.006
## DomHandRght -0.489 0.000 0.056 -0.005 -0.013 0.033 0.016 -0.016 0.821
## Age:SexMen 0.045 -0.426 -0.152 -0.016 -0.012 -0.022 0.031 0.010 -0.046
## Age:EF 0.010 -0.191 0.024 -0.100 -0.039 0.011 -0.029 -0.003 0.026
## Age:RcAfrAm -0.019 -0.703 0.005 -0.078 -0.172 0.055 0.034 0.024 0.015
## SexMen:EF 0.116 -0.013 -0.195 -0.462 -0.130 -0.099 -0.014 -0.044 -0.019
## SxMn:RcAfrA 0.104 -0.022 -0.758 -0.074 -0.625 0.065 -0.006 -0.032 -0.004
## EF:RacAfrAm -0.020 -0.068 0.008 -0.610 -0.057 -0.033 -0.012 0.070 0.011
## DmHndR Ag:SxM Age:EF Ag:RAA SxM:EF SM:RAA
## Age
## SexMen
## EF
## RaceAfrAm
## BMI
## PovStatBelw
## WRATtotal
## DomHandLeft
## DomHandRght
## Age:SexMen -0.050
## Age:EF 0.006 -0.011
## Age:RcAfrAm 0.023 -0.016 0.204
## SexMen:EF -0.033 0.092 0.029 -0.013
## SxMn:RcAfrA -0.006 0.032 -0.004 0.014 0.179
## EF:RacAfrAm 0.044 -0.009 0.030 0.087 -0.036 -0.023
Right Grip Strength - Model 3 (Three-Way Interactions)
mm3 = lmer(PhyPerfRight~(Age + Sex + EF + Race)^3 + BMI + PovStat + WRATtotal + DomHand + (Age|HNDid), data=DanThesisFinal,REML=F)
summary(mm3)
## Linear mixed model fit by maximum likelihood ['merModLmerTest']
## Formula: PhyPerfRight ~ (Age + Sex + EF + Race)^3 + BMI + PovStat + WRATtotal + DomHand + (Age | HNDid)
## Data: DanThesisFinal
##
## AIC BIC logLik deviance
## 15655 15793 -7804 15607
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 33.36550 5.776
## Age 0.00212 0.046 -1.00
## Residual 29.90691 5.469
## Number of obs: 2273, groups: HNDid, 1489
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 22.5613 1.7824 1880.0000 12.66 < 2e-16
## Age -0.2301 0.0420 1815.0000 -5.48 5e-08
## SexMen 15.6206 0.6453 1392.0000 24.21 < 2e-16
## EF 0.2490 0.1448 2015.0000 1.72 0.08553
## RaceAfrAm 1.7910 0.5321 1375.0000 3.37 0.00078
## BMI 0.0823 0.0249 1671.0000 3.31 0.00096
## PovStatBelow -1.3529 0.4081 1442.0000 -3.32 0.00094
## WRATtotal 0.0558 0.0281 1575.0000 1.98 0.04760
## DomHandLeft -0.8230 1.0142 2243.0000 -0.81 0.41720
## DomHandRight 0.3801 0.8926 2176.0000 0.43 0.67029
## Age:SexMen -0.0700 0.0668 1825.0000 -1.05 0.29505
## Age:EF 0.0155 0.0127 2185.0000 1.23 0.22048
## Age:RaceAfrAm 0.0849 0.0537 1798.0000 1.58 0.11382
## SexMen:EF 0.2024 0.2190 2101.0000 0.92 0.35545
## SexMen:RaceAfrAm 0.1665 0.8241 1381.0000 0.20 0.83992
## EF:RaceAfrAm 0.0352 0.1947 2110.0000 0.18 0.85675
## Age:SexMen:EF -0.0082 0.0154 2202.0000 -0.53 0.59466
## Age:SexMen:RaceAfrAm 0.0562 0.0856 1800.0000 0.66 0.51127
## Age:EF:RaceAfrAm -0.0377 0.0152 2201.0000 -2.48 0.01305
## SexMen:EF:RaceAfrAm 0.0049 0.2978 2169.0000 0.02 0.98688
##
## Correlation of Fixed Effects:
## (Intr) Age SexMen EF RcAfrA BMI PvSttB WRATtt DmHndL
## Age -0.039
## SexMen -0.229 0.143
## EF 0.065 0.104 0.161
## RaceAfrAm -0.257 0.191 0.480 0.163
## BMI -0.453 -0.046 0.071 0.081 -0.039
## PovStatBelw -0.234 0.013 0.041 0.055 -0.024 0.052
## WRATtotal -0.704 -0.008 0.040 -0.226 0.152 -0.007 0.172
## DomHandLeft -0.442 0.030 0.028 -0.013 -0.003 0.047 0.018 0.005
## DomHandRght -0.488 0.010 0.050 -0.017 -0.012 0.033 0.016 -0.017 0.821
## Age:SexMen 0.047 -0.609 -0.245 -0.054 -0.114 0.009 0.010 0.003 -0.054
## Age:EF 0.045 -0.293 -0.034 -0.160 -0.081 -0.018 -0.010 -0.001 -0.005
## Age:RcAfrAm 0.005 -0.787 -0.111 -0.083 -0.217 0.055 0.020 0.017 -0.012
## SexMen:EF 0.099 -0.059 -0.282 -0.611 -0.137 -0.066 -0.019 -0.043 0.018
## SxMn:RcAfrA 0.106 -0.114 -0.768 -0.113 -0.631 0.071 -0.010 -0.033 -0.005
## EF:RacAfrAm 0.005 -0.073 -0.116 -0.707 -0.083 -0.030 -0.017 0.041 0.037
## Age:SxMn:EF -0.011 0.143 0.105 0.074 0.035 -0.020 -0.012 -0.003 0.012
## Ag:SxMn:RAA -0.023 0.482 0.192 0.044 0.132 -0.029 0.013 0.003 0.032
## Ag:EF:RcAfA -0.052 0.185 -0.013 0.104 0.061 0.055 -0.004 0.001 0.022
## SxMn:EF:RAA -0.025 0.040 0.198 0.446 0.059 0.000 0.014 0.018 -0.045
## DmHndR Ag:SxM Age:EF Ag:RAA SxM:EF SM:RAA EF:RAA A:SM:E A:SM:R
## Age
## SexMen
## EF
## RaceAfrAm
## BMI
## PovStatBelw
## WRATtotal
## DomHandLeft
## DomHandRght
## Age:SexMen -0.038
## Age:EF -0.021 0.119
## Age:RcAfrAm 0.008 0.483 0.243
## SexMen:EF 0.001 0.134 0.088 0.047
## SxMn:RcAfrA -0.004 0.187 0.008 0.135 0.209
## EF:RacAfrAm 0.054 0.033 0.103 0.069 0.455 0.048
## Age:SxMn:EF 0.022 -0.208 -0.483 -0.143 -0.083 -0.055 -0.027
## Ag:SxMn:RAA 0.008 -0.787 -0.106 -0.619 -0.105 -0.211 -0.035 0.192
## Ag:EF:RcAfA 0.023 -0.015 -0.637 -0.109 -0.040 0.032 -0.107 -0.008 -0.004
## SxMn:EF:RAA -0.036 -0.089 -0.044 -0.040 -0.729 -0.114 -0.647 0.022 0.100
## A:EF:R
## Age
## SexMen
## EF
## RaceAfrAm
## BMI
## PovStatBelw
## WRATtotal
## DomHandLeft
## DomHandRght
## Age:SexMen
## Age:EF
## Age:RcAfrAm
## SexMen:EF
## SxMn:RcAfrA
## EF:RacAfrAm
## Age:SxMn:EF
## Ag:SxMn:RAA
## Ag:EF:RcAfA
## SxMn:EF:RAA 0.031
Right Grip Strength - Model 4 (Four-Way Interactions)
mm4 = lmer(PhyPerfRight~(Age + Sex + EF + Race)^4 + BMI + PovStat + WRATtotal + DomHand + (Age|HNDid), data=DanThesisFinal,REML=F)
summary(mm4)
## Linear mixed model fit by maximum likelihood ['merModLmerTest']
## Formula: PhyPerfRight ~ (Age + Sex + EF + Race)^4 + BMI + PovStat + WRATtotal + DomHand + (Age | HNDid)
## Data: DanThesisFinal
##
## AIC BIC logLik deviance
## 15657 15800 -7803 15607
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 33.34594 5.775
## Age 0.00203 0.045 -1.00
## Residual 29.90384 5.468
## Number of obs: 2273, groups: HNDid, 1489
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 2.25e+01 1.78e+00 1.88e+03 12.64 < 2e-16
## Age -2.26e-01 4.25e-02 1.83e+03 -5.32 1.2e-07
## SexMen 1.57e+01 6.49e-01 1.38e+03 24.13 < 2e-16
## EF 2.56e-01 1.45e-01 1.94e+03 1.77 0.07772
## RaceAfrAm 1.81e+00 5.33e-01 1.37e+03 3.39 0.00072
## BMI 8.24e-02 2.49e-02 1.67e+03 3.31 0.00094
## PovStatBelow -1.35e+00 4.08e-01 1.44e+03 -3.30 0.00100
## WRATtotal 5.56e-02 2.81e-02 1.58e+03 1.98 0.04839
## DomHandLeft -8.03e-01 1.01e+00 2.24e+03 -0.79 0.42887
## DomHandRight 3.85e-01 8.93e-01 2.18e+03 0.43 0.66600
## Age:SexMen -7.87e-02 6.85e-02 1.86e+03 -1.15 0.25070
## Age:EF 1.17e-02 1.43e-02 2.17e+03 0.82 0.41103
## Age:RaceAfrAm 8.20e-02 5.39e-02 1.81e+03 1.52 0.12830
## SexMen:EF 1.91e-01 2.20e-01 2.00e+03 0.87 0.38480
## SexMen:RaceAfrAm 1.11e-01 8.29e-01 1.38e+03 0.13 0.89377
## EF:RaceAfrAm 2.53e-02 1.95e-01 2.02e+03 0.13 0.89716
## Age:SexMen:EF 1.41e-03 2.26e-02 2.20e+03 0.06 0.95048
## Age:SexMen:RaceAfrAm 6.37e-02 8.65e-02 1.83e+03 0.74 0.46160
## Age:EF:RaceAfrAm -3.06e-02 1.96e-02 2.20e+03 -1.56 0.11871
## SexMen:EF:RaceAfrAm 2.25e-02 2.99e-01 2.08e+03 0.08 0.94019
## Age:SexMen:EF:RaceAfrAm -1.79e-02 3.09e-02 2.20e+03 -0.58 0.56393
Right Grip Strength - Model 5 (Retaining significant higher-order terms)
mm5 = lmer(PhyPerfRight~Age + Sex + EF + Race + BMI + PovStat + WRATtotal + DomHand + (Age|HNDid) + Age*Race + Age*EF + EF*Race + Age*EF*Race ,data=DanThesisFinal,REML=F)
summary(mm5)
## Linear mixed model fit by maximum likelihood ['merModLmerTest']
## Formula: PhyPerfRight ~ Age + Sex + EF + Race + BMI + PovStat + WRATtotal + DomHand + (Age | HNDid) + Age * Race + Age * EF + EF * Race + Age * EF * Race
## Data: DanThesisFinal
##
## AIC BIC logLik deviance
## 15647 15750 -7806 15611
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 33.66616 5.8023
## Age 0.00214 0.0462 -1.00
## Residual 29.81225 5.4601
## Number of obs: 2273, groups: HNDid, 1489
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 22.3426 1.7648 1879.0000 12.66 < 2e-16
## Age -0.2571 0.0334 1850.0000 -7.71 2.1e-14
## SexMen 15.7034 0.3971 1443.0000 39.54 < 2e-16
## EF 0.3376 0.1146 2068.0000 2.94 0.00327
## RaceAfrAm 1.8785 0.4134 1366.0000 4.54 6.0e-06
## BMI 0.0852 0.0247 1670.0000 3.45 0.00057
## PovStatBelow -1.3378 0.4086 1443.0000 -3.27 0.00108
## WRATtotal 0.0580 0.0281 1576.0000 2.06 0.03956
## DomHandLeft -0.8661 1.0119 2241.0000 -0.86 0.39213
## DomHandRight 0.3775 0.8906 2173.0000 0.42 0.67169
## Age:RaceAfrAm 0.1059 0.0422 1824.0000 2.51 0.01212
## Age:EF 0.0120 0.0111 2202.0000 1.08 0.27986
## EF:RaceAfrAm 0.0370 0.1484 2064.0000 0.25 0.80341
## Age:EF:RaceAfrAm -0.0376 0.0151 2203.0000 -2.48 0.01318
##
## Correlation of Fixed Effects:
## (Intr) Age SexMen EF RcAfrA BMI PvSttB WRATtt DmHndL
## Age -0.015
## SexMen -0.223 -0.009
## EF 0.157 0.135 0.003
## RaceAfrAm -0.245 0.198 -0.009 0.144
## BMI -0.458 -0.048 0.194 0.051 0.005
## PovStatBelw -0.235 0.024 0.057 0.055 -0.040 0.052
## WRATtotal -0.706 -0.007 0.020 -0.319 0.169 -0.009 0.171
## DomHandLeft -0.445 -0.004 0.034 -0.001 -0.009 0.046 0.020 0.005
## DomHandRght -0.489 -0.017 0.069 -0.019 -0.019 0.030 0.017 -0.017 0.820
## Age:RcAfrAm -0.012 -0.790 0.018 -0.107 -0.217 0.045 0.035 0.024 0.011
## Age:EF 0.042 -0.308 0.020 -0.138 -0.111 -0.025 -0.017 -0.001 0.001
## EF:RacAfrAm -0.008 -0.102 -0.019 -0.710 -0.104 -0.041 -0.012 0.068 0.008
## Ag:EF:RcAfA -0.051 0.224 0.004 0.101 0.105 0.049 -0.004 0.000 0.022
## DmHndR Ag:RAA Age:EF EF:RAA
## Age
## SexMen
## EF
## RaceAfrAm
## BMI
## PovStatBelw
## WRATtotal
## DomHandLeft
## DomHandRght
## Age:RcAfrAm 0.018
## Age:EF -0.012 0.243
## EF:RacAfrAm 0.040 0.102 0.105
## Ag:EF:RcAfA 0.022 -0.144 -0.731 -0.114