Systolic Blood Pressure: Model 1
ms1 = lmer(SBP~Age + Sex + (Age|HNDid),bp)
(st=step(ms1))
Systolic Blood Pressure: Model 1: Summary
summary(ms1)
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: SBP ~ Age + Sex + (Age | HNDid)
## Data: bp
##
## REML criterion at convergence: 36816
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 6.2497 2.500
## Age 0.0262 0.162 1.00
## Residual 164.9819 12.845
## Number of obs: 4389, groups: HNDid, 2390
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 9.32e+01 1.54e+00 2.23e+03 60.47 <2e-16
## Age 5.58e-01 3.03e-02 2.80e+03 18.43 <2e-16
## SexMen 1.06e-02 5.93e-01 2.15e+03 0.02 0.99
##
## Correlation of Fixed Effects:
## (Intr) Age
## Age -0.968
## SexMen -0.158 -0.004
Systolic Blood Pressure: Model 1: Plots
plot(st)
Diastolic Blood Pressure: Model 1
md1 = lmer(DBP~Age + Sex + (Age|HNDid),bp)
(st=step(md1))
Diastolic Blood Pressure: Model 1:Summary
summary(md1)
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: DBP ~ Age + Sex + (Age | HNDid)
## Data: bp
##
## REML criterion at convergence: 32619
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 78.6899 8.8707
## Age 0.0027 0.0519 -1.00
## Residual 67.6361 8.2241
## Number of obs: 4383, groups: HNDid, 2389
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 72.3075 0.9960 2205.0000 72.60 <2e-16
## Age -0.0250 0.0188 2858.0000 -1.33 0.18
## SexMen 2.0823 0.3624 2308.0000 5.75 1e-08
##
## Correlation of Fixed Effects:
## (Intr) Age
## Age -0.972
## SexMen -0.160 0.007
Diastolic Blood Pressure: Model 1: Plots
plot(st)
Mean Arterial Pressure: Model 1
mm1 = lmer(MAP~Age + Sex + (Age|HNDid),bp)
(st=step(mm1))
Mean Arterial Pressure: Model 1: Summary
summary(mm1)
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: MAP ~ Age + Sex + (Age | HNDid)
## Data: bp
##
## REML criterion at convergence: 33557
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 5.63e+01 7.50381
## Age 4.77e-05 0.00691 -1.00
## Residual 8.18e+01 9.04553
## Number of obs: 4383, groups: HNDid, 2389
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 79.076 1.100 2237.000 71.89 <2e-16
## Age 0.174 0.021 2881.000 8.30 <2e-16
## SexMen 1.278 0.408 2311.000 3.13 0.0018
##
## Correlation of Fixed Effects:
## (Intr) Age
## Age -0.971
## SexMen -0.160 0.004
Mean Arterial Pressure: Model 1: Plots
plot(st)
Systolic Blood Pressure: Model 2
ms2 = lmer(SBP~Age + Sex + Age50 + (Age|HNDid),bp)
(st=step(ms2))
Systolic Blood Pressure: Model 2: Summary
summary(ms2)
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: SBP ~ Age + Sex + Age50 + (Age | HNDid)
## Data: bp
##
## REML criterion at convergence: 36808
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 5.7839 2.405
## Age 0.0267 0.163 1.00
## Residual 164.8294 12.839
## Number of obs: 4389, groups: HNDid, 2390
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 9.68e+01 2.11e+00 3.53e+03 45.88 <2e-16
## Age 4.66e-01 4.81e-02 4.17e+03 9.69 <2e-16
## SexMen 9.87e-03 5.92e-01 2.16e+03 0.02 0.987
## Age50 2.08e+00 8.40e-01 3.72e+03 2.48 0.013
##
## Correlation of Fixed Effects:
## (Intr) Age SexMen
## Age -0.976
## SexMen -0.117 -0.001
## Age50 0.684 -0.777 -0.003
Systolic Blood Pressure: Model 2: Plots
plot(st)
Diastolic Blood Pressure: Model 2
md2 = lmer(DBP~Age + Sex + Age50 + (Age|HNDid),bp)
(st=step(md2))
Diastolic Blood Pressure: Model 2:Summary
summary(md2)
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: DBP ~ Age + Sex + Age50 + (Age | HNDid)
## Data: bp
##
## REML criterion at convergence: 32612
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 81.57427 9.032
## Age 0.00302 0.055 -1.00
## Residual 67.46996 8.214
## Number of obs: 4383, groups: HNDid, 2389
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 74.5849 1.3474 3946.0000 55.35 < 2e-16
## Age -0.0837 0.0301 4352.0000 -2.78 0.0054
## SexMen 2.0742 0.3623 2308.0000 5.72 1.2e-08
## Age50 1.3100 0.5301 3864.0000 2.47 0.0135
##
## Correlation of Fixed Effects:
## (Intr) Age SexMen
## Age -0.974
## SexMen -0.125 0.012
## Age50 0.673 -0.781 -0.009
Diastolic Blood Pressure: Model 2: Plots
plot(st)
Mean Arterial Pressure: Model 2
mm2 = lmer(MAP~Age + Sex + Age50 + (Age|HNDid),bp)
(st=step(mm2))
Mean Arterial Pressure: Model 2: Summary
summary(mm2)
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: MAP ~ Age + Sex + Age50 + (Age | HNDid)
## Data: bp
##
## REML criterion at convergence: 33550
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 5.76e+01 7.58759
## Age 7.25e-05 0.00851 -1.00
## Residual 8.17e+01 9.03622
## Number of obs: 4383, groups: HNDid, 2389
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 8.17e+01 1.49e+00 3.86e+03 54.73 <2e-16
## Age 1.07e-01 3.35e-02 4.37e+03 3.19 0.0014
## SexMen 1.27e+00 4.08e-01 2.31e+03 3.11 0.0019
## Age50 1.51e+00 5.88e-01 3.81e+03 2.57 0.0102
##
## Correlation of Fixed Effects:
## (Intr) Age SexMen
## Age -0.975
## SexMen -0.123 0.008
## Age50 0.676 -0.780 -0.007
Mean Arterial Pressure: Model 2: Plots
plot(st)
Systolic Blood Pressure: Model 3
ms3 = lmer(SBP~Age + Sex + Age50 + Race + PovStat + Race*PovStat + (Age|HNDid),bp)
(st=step(ms3))
Systolic Blood Pressure: Model 3: Summary
summary(ms3)
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: SBP ~ Age + Sex + Age50 + Race + PovStat + Race * PovStat + (Age | HNDid)
## Data: bp
##
## REML criterion at convergence: 36758
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 6.2597 2.502
## Age 0.0249 0.158 1.00
## Residual 165.0642 12.848
## Number of obs: 4389, groups: HNDid, 2390
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 93.8302 2.1604 3478.0000 43.43 < 2e-16
## Age 0.4704 0.0479 4160.0000 9.83 < 2e-16
## SexMen 0.0606 0.5875 2142.0000 0.10 0.918
## Age50 2.0707 0.8375 3730.0000 2.47 0.013
## RaceAfrAm 3.8032 0.7515 2134.0000 5.06 4.5e-07
## PovStatBelow 1.0076 1.0021 2098.0000 1.01 0.315
## RaceAfrAm:PovStatBelow -0.0571 1.2489 2122.0000 -0.05 0.964
##
## Correlation of Fixed Effects:
## (Intr) Age SexMen Age50 RcAfrA PvSttB
## Age -0.952
## SexMen -0.122 0.001
## Age50 0.667 -0.778 -0.003
## RaceAfrAm -0.186 -0.001 0.002 -0.007
## PovStatBelw -0.156 0.009 0.044 -0.006 0.407
## RcAfrAm:PSB 0.106 0.007 -0.012 0.009 -0.602 -0.801
Systolic Blood Pressure: Model 3: Plots
plot(st)
Diastolic Blood Pressure: Model 3
md3 = lmer(DBP~Age + Sex + Age50 + Race + PovStat + Race*PovStat + (Age|HNDid),bp)
(st=step(md3))
Diastolic Blood Pressure: Model 3:Summary
summary(md3)
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: DBP ~ Age + Sex + Age50 + Race + PovStat + Race * PovStat + (Age | HNDid)
## Data: bp
##
## REML criterion at convergence: 32569
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 81.96018 9.0532
## Age 0.00331 0.0575 -1.00
## Residual 67.51401 8.2167
## Number of obs: 4383, groups: HNDid, 2389
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 72.7675 1.3785 3883.0000 52.79 < 2e-16
## Age -0.0792 0.0299 4346.0000 -2.65 0.0082
## SexMen 2.1381 0.3597 2299.0000 5.94 3.2e-09
## Age50 1.3081 0.5286 3882.0000 2.47 0.0134
## RaceAfrAm 1.9175 0.4576 2249.0000 4.19 2.9e-05
## PovStatBelow 0.8541 0.6132 2249.0000 1.39 0.1638
## RaceAfrAm:PovStatBelow 0.2361 0.7652 2284.0000 0.31 0.7577
##
## Correlation of Fixed Effects:
## (Intr) Age SexMen Age50 RcAfrA PvSttB
## Age -0.952
## SexMen -0.130 0.014
## Age50 0.657 -0.782 -0.010
## RaceAfrAm -0.179 0.000 0.005 -0.007
## PovStatBelw -0.157 0.017 0.044 -0.008 0.405
## RcAfrAm:PSB 0.102 0.005 -0.013 0.010 -0.598 -0.800
Diastolic Blood Pressure: Model 3: Plots
plot(st)
Mean Arterial Pressure: Model 3
mm3 = lmer(MAP~Age + Sex + Age50 + Race + PovStat + Race*PovStat + (Age|HNDid),bp)
(st=step(mm3))
Mean Arterial Pressure: Model 3: Summary
summary(mm3)
## Linear mixed model fit by REML ['merModLmerTest']
## Formula: MAP ~ Age + Sex + Age50 + Race + PovStat + Race * PovStat + (Age | HNDid)
## Data: bp
##
## REML criterion at convergence: 33496
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## HNDid (Intercept) 5.92e+01 7.6967
## Age 1.85e-04 0.0136 -1.00
## Residual 8.17e+01 9.0415
## Number of obs: 4383, groups: HNDid, 2389
##
## Fixed effects:
## Estimate Std. Error df t value Pr(>|t|)
## (Intercept) 7.95e+01 1.53e+00 3.81e+03 52.06 < 2e-16
## Age 1.11e-01 3.33e-02 4.36e+03 3.34 0.00085
## SexMen 1.33e+00 4.04e-01 2.30e+03 3.29 0.00101
## Age50 1.51e+00 5.86e-01 3.84e+03 2.57 0.01015
## RaceAfrAm 2.59e+00 5.15e-01 2.26e+03 5.03 5.2e-07
## PovStatBelow 9.76e-01 6.89e-01 2.25e+03 1.42 0.15670
## RaceAfrAm:PovStatBelow 7.13e-02 8.60e-01 2.28e+03 0.08 0.93388
##
## Correlation of Fixed Effects:
## (Intr) Age SexMen Age50 RcAfrA PvSttB
## Age -0.952
## SexMen -0.128 0.010
## Age50 0.660 -0.781 -0.008
## RaceAfrAm -0.182 -0.001 0.004 -0.007
## PovStatBelw -0.157 0.015 0.044 -0.008 0.406
## RcAfrAm:PSB 0.103 0.005 -0.013 0.009 -0.599 -0.800
Mean Arterial Pressure: Model 3: Plots
plot(st)
