Lifetime IPV and Hypertension Regression Models (Updated 3/1/15)

Systolic Blood Pressure

m1sBP = lmer(adjSBP~Age * LifeIPV * Sex  * Race + BMIcat + PovStat + (Age|HNDid) + (1|subclass),LifeIPVbp)
(st = step(m1sBP))

Final Model by backwards elimination

m1sBP = lmer(adjSBP ~Age + LifeIPV + Sex + Race + BMIcat + Age:LifeIPV + Age:Sex + LifeIPV:Sex + Age:LifeIPV:Sex + (Age|HNDid) + (1|subclass),data = LifeIPVbp)
summary(m1sBP)
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: adjSBP ~ Age + LifeIPV + Sex + Race + BMIcat + Age:LifeIPV +      Age:Sex + LifeIPV:Sex + Age:LifeIPV:Sex + (Age | HNDid) +      (1 | subclass) 
##    Data: LifeIPVbp 
## 
## REML criterion at convergence: 5533 
## 
## Random effects:
##  Groups   Name        Variance Std.Dev. Corr
##  HNDid    (Intercept) 165.906  12.880       
##           Age           0.031   0.176   1.00
##  subclass (Intercept)   0.000   0.000       
##  Residual             133.399  11.550       
## Number of obs: 666, groups: HNDid, 334; subclass, 222
## 
## Fixed effects:
##                     Estimate Std. Error      df t value Pr(>|t|)
## (Intercept)          113.413      2.272 386.000   49.93  < 2e-16
## Age                    0.795      0.149 477.000    5.34  1.4e-07
## LifeIPV1              -0.860      2.217 280.000   -0.39   0.6984
## SexMen                -2.414      2.095 278.000   -1.15   0.2502
## RaceAfrAm              4.542      1.793 314.000    2.53   0.0118
## BMIcatObese            8.975      1.663 604.000    5.40  9.8e-08
## Age:LifeIPV1          -0.211      0.224 477.000   -0.94   0.3472
## Age:SexMen            -0.389      0.220 489.000   -1.77   0.0772
## LifeIPV1:SexMen       15.249      4.823 288.000    3.16   0.0017
## Age:LifeIPV1:SexMen    1.370      0.513 502.000    2.67   0.0078
## 
## Correlation of Fixed Effects:
##             (Intr) Age    LfIPV1 SexMen RcAfrA BMIctO Ag:LIPV1 Ag:SxM
## Age          0.201                                                   
## LifeIPV1    -0.499 -0.212                                            
## SexMen      -0.451 -0.220  0.453                                     
## RaceAfrAm   -0.519 -0.010  0.194 -0.051                              
## BMIcatObese -0.511  0.015 -0.071  0.048 -0.089                       
## Age:LifIPV1 -0.122 -0.664  0.289  0.143  0.032 -0.056                
## Age:SexMen  -0.115 -0.676  0.135  0.287 -0.034 -0.007  0.448         
## LfIPV1:SxMn  0.163  0.097 -0.447 -0.429 -0.010  0.071 -0.133   -0.124
## Ag:LIPV1:SM  0.001  0.290 -0.118 -0.122  0.044  0.060 -0.437   -0.430
##             LIPV1:
## Age               
## LifeIPV1          
## SexMen            
## RaceAfrAm         
## BMIcatObese       
## Age:LifIPV1       
## Age:SexMen        
## LfIPV1:SxMn       
## Ag:LIPV1:SM  0.414

Age * LifeIPV * Sex Interaction Plot (sBP)
plot of chunk unnamed-chunk-5

Diastolic Blood Pressure

m1dBP = lmer(adjDBP~Age * LifeIPV * Sex * Race + PovStat + BMIcat + (Age|HNDid) + (1|subclass),LifeIPVbp)
(st = step(m1dBP))

Final Model by backwards elimination

m1dBP = lmer(adjDBP~Age + LifeIPV + Sex + Race + PovStat + BMIcat + Age:LifeIPV + Age:Sex + LifeIPV:Sex + Age:LifeIPV:Sex + (Age|HNDid) + (1|subclass), LifeIPVbp)
summary(m1dBP)
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: adjDBP ~ Age + LifeIPV + Sex + Race + PovStat + BMIcat + Age:LifeIPV +      Age:Sex + LifeIPV:Sex + Age:LifeIPV:Sex + (Age | HNDid) +      (1 | subclass) 
##    Data: LifeIPVbp 
## 
## REML criterion at convergence: 4887 
## 
## Random effects:
##  Groups   Name        Variance Std.Dev. Corr 
##  HNDid    (Intercept) 3.38e+01 5.8133        
##           Age         7.28e-04 0.0270   -1.00
##  subclass (Intercept) 5.72e-04 0.0239        
##  Residual             6.38e+01 7.9866        
## Number of obs: 666, groups: HNDid, 334; subclass, 222
## 
## Fixed effects:
##                     Estimate Std. Error       df t value Pr(>|t|)
## (Intercept)          64.0079     1.2794 359.0000   50.03  < 2e-16
## Age                  -0.0262     0.0835 416.0000   -0.31   0.7535
## LifeIPV1              1.1155     1.1695 306.0000    0.95   0.3409
## SexMen                2.3086     1.1020 304.0000    2.09   0.0370
## RaceAfrAm             3.2575     0.9821 321.0000    3.32   0.0010
## PovStatBelow          2.8070     0.9415 320.0000    2.98   0.0031
## BMIcatObese           5.3263     0.9808 536.0000    5.43  8.5e-08
## Age:LifeIPV1         -0.0370     0.1254 413.0000   -0.30   0.7680
## Age:SexMen            0.1214     0.1235 421.0000    0.98   0.3261
## LifeIPV1:SexMen       4.4797     2.5417 302.0000    1.76   0.0790
## Age:LifeIPV1:SexMen   0.6703     0.2987 428.0000    2.24   0.0253
## 
## Correlation of Fixed Effects:
##             (Intr) Age    LfIPV1 SexMen RcAfrA PvSttB BMIctO Ag:LIPV1
## Age          0.036                                                   
## LifeIPV1    -0.435 -0.094                                            
## SexMen      -0.438 -0.081  0.439                                     
## RaceAfrAm   -0.476 -0.001  0.204 -0.057                              
## PovStatBelw -0.251  0.100 -0.092  0.068 -0.042                       
## BMIcatObese -0.531  0.025 -0.083  0.055 -0.109  0.022                
## Age:LifIPV1 -0.009 -0.667  0.108  0.050  0.024 -0.070 -0.062         
## Age:SexMen  -0.011 -0.679  0.062  0.097 -0.016 -0.089 -0.014  0.452  
## LfIPV1:SxMn  0.158  0.033 -0.435 -0.430 -0.001 -0.070  0.074 -0.043  
## Ag:LIPV1:SM -0.036  0.277 -0.032 -0.043  0.053 -0.009  0.052 -0.418  
##             Ag:SxM LIPV1:
## Age                      
## LifeIPV1                 
## SexMen                   
## RaceAfrAm                
## PovStatBelw              
## BMIcatObese              
## Age:LifIPV1              
## Age:SexMen               
## LfIPV1:SxMn -0.039       
## Ag:LIPV1:SM -0.411  0.282

Age * LifeIPV * Sex Interaction Plot (dBP)
plot of chunk unnamed-chunk-9