Dan's Thesis Analyses - Grip Strength (Left)
load(file="/Users/meganwilliams/Desktop/UMBC 2016/Dan's Thesis/DanThesisFinal.rdata")
library(lmerTest)
Left Grip Strength - Model 1 (Main Effects)
mm1 = lmer(PhyPerfLeft~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: PhyPerfLeft ~ Age + Sex + EF + Race + BMI + PovStat + WRATtotal + DomHand + (Age | HNDid)
## Data: DanThesisFinal
##
## AIC BIC logLik deviance
## 15764 15845 -7868 15736
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 39.8137 6.3098
## Age 0.0073 0.0854 -0.51
## Residual 30.6232 5.5338
## Number of obs: 2261, groups: HNDid, 1485
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 23.8385 1.8737 1850.0000 12.72 < 2e-16
## Age -0.2133 0.0214 806.0000 -9.99 < 2e-16
## SexMen 16.8739 0.4237 1426.0000 39.83 < 2e-16
## EF 0.3456 0.0844 2245.0000 4.09 4.4e-05
## RaceAfrAm 2.4589 0.4282 1405.0000 5.74 1.1e-08
## BMI 0.0755 0.0262 1673.0000 2.88 0.0040
## PovStatBelow -1.4035 0.4362 1429.0000 -3.22 0.0013
## WRATtotal 0.0734 0.0299 1552.0000 2.45 0.0143
## DomHandLeft 0.4508 1.0608 2198.0000 0.42 0.6709
## DomHandRight -1.8496 0.9339 2105.0000 -1.98 0.0478
##
## Correlation of Fixed Effects:
## (Intr) Age SexMen EF RcAfrA BMI PvSttB WRATtt DmHndL
## Age -0.030
## SexMen -0.226 0.008
## EF 0.220 0.071 -0.011
## RaceAfrAm -0.250 0.026 -0.009 0.087
## BMI -0.459 -0.033 0.196 0.033 0.003
## PovStatBelw -0.237 0.084 0.058 0.065 -0.039 0.053
## WRATtotal -0.712 0.022 0.022 -0.384 0.185 -0.007 0.174
## DomHandLeft -0.439 0.007 0.034 0.005 -0.010 0.045 0.019 0.006
## DomHandRght -0.484 -0.009 0.070 0.009 -0.017 0.030 0.016 -0.018 0.819
Left Grip Strength - Model 2 (Two-Way Interactions)
mm2 = lmer(PhyPerfLeft~(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: PhyPerfLeft ~ (Age + Sex + EF + Race)^2 + BMI + PovStat + WRATtotal + DomHand + (Age | HNDid)
## Data: DanThesisFinal
##
## AIC BIC logLik deviance
## 15769 15884 -7865 15729
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 39.3523 6.273
## Age 0.0102 0.101 -0.47
## Residual 30.6053 5.532
## Number of obs: 2261, groups: HNDid, 1485
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 24.0542 1.8901 1844.0000 12.73 < 2e-16
## Age -0.2255 0.0384 795.0000 -5.87 6.4e-09
## SexMen 16.8350 0.6619 1435.0000 25.43 < 2e-16
## EF 0.3310 0.1356 2165.0000 2.44 0.0147
## RaceAfrAm 2.3120 0.5604 1441.0000 4.13 3.9e-05
## BMI 0.0745 0.0264 1669.0000 2.82 0.0049
## PovStatBelow -1.3937 0.4358 1423.0000 -3.20 0.0014
## WRATtotal 0.0696 0.0300 1547.0000 2.32 0.0204
## DomHandLeft 0.4805 1.0609 2195.0000 0.45 0.6506
## DomHandRight -1.8564 0.9348 2099.0000 -1.99 0.0472
## Age:SexMen -0.0480 0.0437 828.0000 -1.10 0.2731
## Age:EF -0.0013 0.0080 1083.0000 -0.16 0.8714
## Age:RaceAfrAm 0.0532 0.0445 808.0000 1.20 0.2320
## SexMen:EF 0.2525 0.1571 2186.0000 1.61 0.1081
## SexMen:RaceAfrAm 0.1473 0.8539 1414.0000 0.17 0.8631
## EF:RaceAfrAm -0.1741 0.1554 2177.0000 -1.12 0.2627
##
## Correlation of Fixed Effects:
## (Intr) Age SexMen EF RcAfrA BMI PvSttB WRATtt DmHndL
## Age -0.018
## SexMen -0.231 0.057
## EF 0.092 0.087 0.081
## RaceAfrAm -0.251 0.129 0.468 0.145
## BMI -0.454 -0.049 0.084 0.087 -0.042
## PovStatBelw -0.236 0.014 0.039 0.056 -0.030 0.055
## WRATtotal -0.709 -0.012 0.038 -0.262 0.153 -0.006 0.174
## DomHandLeft -0.437 0.011 0.030 0.004 -0.009 0.047 0.018 0.007
## DomHandRght -0.483 -0.001 0.055 -0.007 -0.016 0.033 0.015 -0.014 0.819
## Age:SexMen 0.044 -0.429 -0.145 -0.017 -0.009 -0.024 0.026 0.009 -0.044
## Age:EF 0.007 -0.187 0.022 -0.097 -0.037 0.009 -0.028 0.000 0.031
## Age:RcAfrAm -0.023 -0.704 0.008 -0.077 -0.162 0.056 0.032 0.026 0.016
## SexMen:EF 0.115 -0.013 -0.191 -0.461 -0.128 -0.097 -0.014 -0.046 -0.018
## SxMn:RcAfrA 0.103 -0.018 -0.758 -0.073 -0.626 0.065 -0.006 -0.031 -0.003
## EF:RacAfrAm -0.020 -0.068 0.008 -0.612 -0.056 -0.034 -0.011 0.073 0.008
## DmHndR Ag:SxM Age:EF Ag:RAA SxM:EF SM:RAA
## Age
## SexMen
## EF
## RaceAfrAm
## BMI
## PovStatBelw
## WRATtotal
## DomHandLeft
## DomHandRght
## Age:SexMen -0.049
## Age:EF 0.010 -0.011
## Age:RcAfrAm 0.023 -0.014 0.200
## SexMen:EF -0.031 0.091 0.026 -0.012
## SxMn:RcAfrA -0.005 0.030 -0.003 0.012 0.176
## EF:RacAfrAm 0.041 -0.008 0.031 0.086 -0.037 -0.022
Left Grip Strength - Model 3 (Three-Way Interactions)
mm3 = lmer(PhyPerfLeft~(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: PhyPerfLeft ~ (Age + Sex + EF + Race)^3 + BMI + PovStat + WRATtotal + DomHand + (Age | HNDid)
## Data: DanThesisFinal
##
## AIC BIC logLik deviance
## 15775 15912 -7863 15727
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 39.21732 6.2624
## Age 0.00917 0.0957 -0.46
## Residual 30.63569 5.5350
## Number of obs: 2261, groups: HNDid, 1485
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 2.42e+01 1.89e+00 1.84e+03 12.79 < 2e-16
## Age -2.36e-01 4.49e-02 8.16e+02 -5.26 1.9e-07
## SexMen 1.69e+01 6.85e-01 1.40e+03 24.60 < 2e-16
## EF 3.07e-01 1.52e-01 2.05e+03 2.02 0.0440
## RaceAfrAm 2.26e+00 5.66e-01 1.38e+03 4.00 6.6e-05
## BMI 7.22e-02 2.65e-02 1.66e+03 2.73 0.0065
## PovStatBelow -1.40e+00 4.36e-01 1.42e+03 -3.21 0.0014
## WRATtotal 6.94e-02 3.00e-02 1.55e+03 2.32 0.0207
## DomHandLeft 4.55e-01 1.06e+00 2.20e+03 0.43 0.6684
## DomHandRight -1.88e+00 9.36e-01 2.10e+03 -2.00 0.0452
## Age:SexMen -4.85e-02 7.11e-02 8.62e+02 -0.68 0.4953
## Age:EF 6.92e-03 1.34e-02 1.12e+03 0.52 0.6060
## Age:RaceAfrAm 6.13e-02 5.73e-02 7.89e+02 1.07 0.2848
## SexMen:EF 2.73e-01 2.31e-01 2.10e+03 1.18 0.2380
## SexMen:RaceAfrAm 1.10e-01 8.76e-01 1.39e+03 0.13 0.9004
## EF:RaceAfrAm -1.29e-01 2.05e-01 2.12e+03 -0.63 0.5288
## Age:SexMen:EF 1.11e-02 1.63e-02 1.11e+03 0.68 0.4967
## Age:SexMen:RaceAfrAm -3.88e-04 9.11e-02 8.30e+02 0.00 0.9966
## Age:EF:RaceAfrAm -2.36e-02 1.60e-02 1.10e+03 -1.47 0.1410
## SexMen:EF:RaceAfrAm -4.69e-02 3.13e-01 2.18e+03 -0.15 0.8810
##
## Correlation of Fixed Effects:
## (Intr) Age SexMen EF RcAfrA BMI PvSttB WRATtt DmHndL
## Age -0.035
## SexMen -0.230 0.133
## EF 0.066 0.102 0.159
## RaceAfrAm -0.254 0.178 0.478 0.159
## BMI -0.455 -0.045 0.075 0.081 -0.042
## PovStatBelw -0.236 0.018 0.041 0.054 -0.027 0.055
## WRATtotal -0.708 -0.009 0.040 -0.225 0.152 -0.006 0.174
## DomHandLeft -0.437 0.031 0.027 -0.013 -0.005 0.048 0.017 0.006
## DomHandRght -0.483 0.010 0.050 -0.018 -0.014 0.034 0.014 -0.015 0.819
## Age:SexMen 0.047 -0.611 -0.233 -0.053 -0.106 0.006 0.005 0.002 -0.054
## Age:EF 0.043 -0.285 -0.033 -0.152 -0.078 -0.022 -0.009 0.002 -0.002
## Age:RcAfrAm 0.001 -0.788 -0.104 -0.081 -0.205 0.054 0.016 0.019 -0.013
## SexMen:EF 0.098 -0.057 -0.276 -0.610 -0.134 -0.064 -0.018 -0.044 0.017
## SxMn:RcAfrA 0.106 -0.106 -0.768 -0.111 -0.631 0.071 -0.010 -0.033 -0.005
## EF:RacAfrAm 0.004 -0.072 -0.115 -0.709 -0.082 -0.031 -0.016 0.044 0.032
## Age:SxMn:EF -0.012 0.139 0.102 0.070 0.034 -0.017 -0.012 0.000 0.009
## Ag:SxMn:RAA -0.023 0.483 0.183 0.043 0.125 -0.027 0.015 0.004 0.034
## Ag:EF:RcAfA -0.051 0.180 -0.012 0.099 0.059 0.056 -0.004 -0.003 0.024
## SxMn:EF:RAA -0.024 0.039 0.194 0.447 0.059 -0.001 0.012 0.018 -0.043
## 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.114
## Age:RcAfrAm 0.008 0.485 0.236
## SexMen:EF 0.002 0.131 0.081 0.045
## SxMn:RcAfrA -0.002 0.177 0.007 0.128 0.205
## EF:RacAfrAm 0.051 0.033 0.098 0.069 0.454 0.048
## Age:SxMn:EF 0.019 -0.204 -0.483 -0.139 -0.082 -0.052 -0.025
## Ag:SxMn:RAA 0.009 -0.787 -0.103 -0.621 -0.103 -0.201 -0.036 0.189
## Ag:EF:RcAfA 0.027 -0.013 -0.637 -0.106 -0.035 0.031 -0.102 -0.008 -0.005
## SxMn:EF:RAA -0.034 -0.088 -0.041 -0.040 -0.730 -0.112 -0.647 0.024 0.098
## 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.028
Left Grip Strength - Model 4 (Four-Way Interactions)
mm4 = lmer(PhyPerfLeft~(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: PhyPerfLeft ~ (Age + Sex + EF + Race)^4 + BMI + PovStat + WRATtotal + DomHand + (Age | HNDid)
## Data: DanThesisFinal
##
## AIC BIC logLik deviance
## 15777 15920 -7863 15727
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 39.23226 6.2636
## Age 0.00931 0.0965 -0.46
## Residual 30.62353 5.5339
## Number of obs: 2261, groups: HNDid, 1485
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 2.42e+01 1.89e+00 1.84e+03 12.79 < 2e-16
## Age -2.39e-01 4.54e-02 8.31e+02 -5.26 1.8e-07
## SexMen 1.68e+01 6.90e-01 1.38e+03 24.40 < 2e-16
## EF 3.02e-01 1.53e-01 1.98e+03 1.97 0.0486
## RaceAfrAm 2.25e+00 5.66e-01 1.38e+03 3.98 7.4e-05
## BMI 7.21e-02 2.65e-02 1.66e+03 2.72 0.0066
## PovStatBelow -1.40e+00 4.36e-01 1.42e+03 -3.22 0.0013
## WRATtotal 6.95e-02 3.00e-02 1.55e+03 2.32 0.0204
## DomHandLeft 4.40e-01 1.06e+00 2.20e+03 0.41 0.6793
## DomHandRight -1.88e+00 9.36e-01 2.10e+03 -2.01 0.0446
## Age:SexMen -4.19e-02 7.29e-02 8.66e+02 -0.57 0.5659
## Age:EF 9.86e-03 1.51e-02 1.12e+03 0.65 0.5145
## Age:RaceAfrAm 6.35e-02 5.75e-02 7.96e+02 1.10 0.2697
## SexMen:EF 2.81e-01 2.32e-01 2.01e+03 1.21 0.2250
## SexMen:RaceAfrAm 1.52e-01 8.81e-01 1.39e+03 0.17 0.8631
## EF:RaceAfrAm -1.22e-01 2.06e-01 2.04e+03 -0.59 0.5536
## Age:SexMen:EF 3.61e-03 2.40e-02 1.18e+03 0.15 0.8803
## Age:SexMen:RaceAfrAm -6.08e-03 9.21e-02 8.38e+02 -0.07 0.9474
## Age:EF:RaceAfrAm -2.91e-02 2.07e-02 1.04e+03 -1.41 0.1603
## SexMen:EF:RaceAfrAm -6.02e-02 3.15e-01 2.09e+03 -0.19 0.8484
## Age:SexMen:EF:RaceAfrAm 1.39e-02 3.27e-02 1.12e+03 0.42 0.6717
Left Grip Strength - Model 5 (Retaining significant higher-order terms)
mm5 = lmer(PhyPerfLeft~Age + Sex + EF + Race + BMI + PovStat + WRATtotal + DomHand + (Age|HNDid), data=DanThesisFinal,REML=F)
summary(mm5)
## Linear mixed model fit by maximum likelihood ['merModLmerTest']
## Formula: PhyPerfLeft ~ Age + Sex + EF + Race + BMI + PovStat + WRATtotal + DomHand + (Age | HNDid)
## Data: DanThesisFinal
##
## AIC BIC logLik deviance
## 15764 15845 -7868 15736
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 39.8137 6.3098
## Age 0.0073 0.0854 -0.51
## Residual 30.6232 5.5338
## Number of obs: 2261, groups: HNDid, 1485
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 23.8385 1.8737 1850.0000 12.72 < 2e-16
## Age -0.2133 0.0214 806.0000 -9.99 < 2e-16
## SexMen 16.8739 0.4237 1426.0000 39.83 < 2e-16
## EF 0.3456 0.0844 2245.0000 4.09 4.4e-05
## RaceAfrAm 2.4589 0.4282 1405.0000 5.74 1.1e-08
## BMI 0.0755 0.0262 1673.0000 2.88 0.0040
## PovStatBelow -1.4035 0.4362 1429.0000 -3.22 0.0013
## WRATtotal 0.0734 0.0299 1552.0000 2.45 0.0143
## DomHandLeft 0.4508 1.0608 2198.0000 0.42 0.6709
## DomHandRight -1.8496 0.9339 2105.0000 -1.98 0.0478
##
## Correlation of Fixed Effects:
## (Intr) Age SexMen EF RcAfrA BMI PvSttB WRATtt DmHndL
## Age -0.030
## SexMen -0.226 0.008
## EF 0.220 0.071 -0.011
## RaceAfrAm -0.250 0.026 -0.009 0.087
## BMI -0.459 -0.033 0.196 0.033 0.003
## PovStatBelw -0.237 0.084 0.058 0.065 -0.039 0.053
## WRATtotal -0.712 0.022 0.022 -0.384 0.185 -0.007 0.174
## DomHandLeft -0.439 0.007 0.034 0.005 -0.010 0.045 0.019 0.006
## DomHandRght -0.484 -0.009 0.070 0.009 -0.017 0.030 0.016 -0.018 0.819