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