Marijuana and Blood Pressure Analyses 10-04-16
Systolic Blood Pressure
model1 = lm(PhysMeanSBPsitimp~MedHxMarijStatusCurrent + Sex + Race + Age0 + HseHldEducation + PovStat + MedHxCigaretteStatusNFC + MedHxAlcStatusNFC + MedHxCokeStatusNFC + MedHxOpiateStatusNFC + CESimp + DMMclusterdich + CVDclusterdich + MedHxMedsBP + PhysBMI + MedHxMarijStatusCurrent*HTN, data = exdata)
summary(model1)
##
## Call:
## lm(formula = PhysMeanSBPsitimp ~ MedHxMarijStatusCurrent + Sex +
## Race + Age0 + HseHldEducation + PovStat + MedHxCigaretteStatusNFC +
## MedHxAlcStatusNFC + MedHxCokeStatusNFC + MedHxOpiateStatusNFC +
## CESimp + DMMclusterdich + CVDclusterdich + MedHxMedsBP +
## PhysBMI + MedHxMarijStatusCurrent * HTN, data = exdata)
##
## Residuals:
## Min 1Q Median 3Q Max
## -44.521574 -9.674918 -0.303762 9.242150 66.562431
##
## Coefficients:
## Estimate Std. Error t value
## (Intercept) 85.141690793 3.090930055 27.54565
## MedHxMarijStatusCurrent1 1.250305548 1.163846647 1.07429
## SexMen 2.582222614 0.601573284 4.29245
## RaceAfrAm 2.161980923 0.594147944 3.63879
## Age0 0.367989213 0.033974397 10.83137
## HseHldEducation -0.139043872 0.100519778 -1.38325
## PovStatBelow 1.138546578 0.611855684 1.86081
## MedHxCigaretteStatusNFC -0.315287970 0.389005401 -0.81050
## MedHxAlcStatusNFC 0.703616579 0.372089843 1.89099
## MedHxCokeStatusNFC 0.475756138 0.612296863 0.77700
## MedHxOpiateStatusNFC -2.006321995 0.734817559 -2.73037
## CESimp 0.006559648 0.026440309 0.24809
## DMMclusterdich1 -0.685567482 0.679248470 -1.00930
## CVDclusterdich1 -2.369450531 1.294792146 -1.82999
## MedHxMedsBP -10.438176223 0.972703586 -10.73110
## PhysBMI 0.390620450 0.039781785 9.81908
## HTNYes 20.585931135 0.957518972 21.49924
## MedHxMarijStatusCurrent1:HTNYes 3.129540557 1.752626381 1.78563
## Pr(>|t|)
## (Intercept) < 2.22e-16
## MedHxMarijStatusCurrent1 0.28279001
## SexMen 1.8294e-05
## RaceAfrAm 0.00027906
## Age0 < 2.22e-16
## HseHldEducation 0.16670297
## PovStatBelow 0.06287982
## MedHxCigaretteStatusNFC 0.41772566
## MedHxAlcStatusNFC 0.05873325
## MedHxCokeStatusNFC 0.43722544
## MedHxOpiateStatusNFC 0.00636748
## CESimp 0.80408156
## DMMclusterdich1 0.31291984
## CVDclusterdich1 0.06736238
## MedHxMedsBP < 2.22e-16
## PhysBMI < 2.22e-16
## HTNYes < 2.22e-16
## MedHxMarijStatusCurrent1:HTNYes 0.07427159
##
## Residual standard error: 14.74233 on 2699 degrees of freedom
## (85 observations deleted due to missingness)
## Multiple R-squared: 0.3179475, Adjusted R-squared: 0.3136515
## F-statistic: 74.0102 on 17 and 2699 DF, p-value: < 2.2204e-16
Pulse Pressure
model2 = lm(PhysMeanPPsitimp~MedHxMarijStatusCurrent + Sex + Race + Age0 + HseHldEducation + PovStat + MedHxCigaretteStatusNFC + MedHxAlcStatusNFC + MedHxCokeStatusNFC + MedHxOpiateStatusNFC + CESimp + DMMclusterdich + CVDclusterdich + MedHxMedsBP + PhysBMI + MedHxMarijStatusCurrent*HTN, data = exdata)
summary(model2)
##
## Call:
## lm(formula = PhysMeanPPsitimp ~ MedHxMarijStatusCurrent + Sex +
## Race + Age0 + HseHldEducation + PovStat + MedHxCigaretteStatusNFC +
## MedHxAlcStatusNFC + MedHxCokeStatusNFC + MedHxOpiateStatusNFC +
## CESimp + DMMclusterdich + CVDclusterdich + MedHxMedsBP +
## PhysBMI + MedHxMarijStatusCurrent * HTN, data = exdata)
##
## Residuals:
## Min 1Q Median 3Q Max
## -49.393732 -7.413964 -1.036415 6.370976 49.990828
##
## Coefficients:
## Estimate Std. Error t value
## (Intercept) 17.919929446 2.346457551 7.63701
## MedHxMarijStatusCurrent1 2.218134921 0.883525898 2.51055
## SexMen -0.547994391 0.456680077 -1.19995
## RaceAfrAm 1.782507885 0.451043182 3.95197
## Age0 0.403807853 0.025791422 15.65667
## HseHldEducation -0.053284176 0.076308874 -0.69827
## PovStatBelow 0.108153026 0.464485887 0.23284
## MedHxCigaretteStatusNFC 0.368347653 0.295310681 1.24732
## MedHxAlcStatusNFC 0.364754856 0.282469356 1.29131
## MedHxCokeStatusNFC -0.644836117 0.464820805 -1.38728
## MedHxOpiateStatusNFC -1.051254732 0.557831519 -1.88454
## CESimp 0.004291475 0.020071973 0.21380
## DMMclusterdich1 1.157480009 0.515646641 2.24472
## CVDclusterdich1 0.371809954 0.982932242 0.37827
## MedHxMedsBP -5.637601703 0.738421004 -7.63467
## PhysBMI 0.248360826 0.030200059 8.22385
## HTNYes 9.713759598 0.726893712 13.36338
## MedHxMarijStatusCurrent1:HTNYes 0.193625012 1.330493843 0.14553
## Pr(>|t|)
## (Intercept) 3.0654e-14
## MedHxMarijStatusCurrent1 0.012112
## SexMen 0.230263
## RaceAfrAm 7.9500e-05
## Age0 < 2.22e-16
## HseHldEducation 0.485069
## PovStatBelow 0.815900
## MedHxCigaretteStatusNFC 0.212387
## MedHxAlcStatusNFC 0.196708
## MedHxCokeStatusNFC 0.165471
## MedHxOpiateStatusNFC 0.059600
## CESimp 0.830716
## DMMclusterdich1 0.024867
## CVDclusterdich1 0.705263
## MedHxMedsBP 3.1204e-14
## PhysBMI 3.0265e-16
## HTNYes < 2.22e-16
## MedHxMarijStatusCurrent1:HTNYes 0.884304
##
## Residual standard error: 11.19154 on 2699 degrees of freedom
## (85 observations deleted due to missingness)
## Multiple R-squared: 0.2479523, Adjusted R-squared: 0.2432155
## F-statistic: 52.34519 on 17 and 2699 DF, p-value: < 2.2204e-16
Diastolic Blood Pressure
model3 = lm(PhysMeanDBPsitimp~MedHxMarijStatusCurrent + Sex + Race + Age0 + HseHldEducation + PovStat + MedHxCigaretteStatusNFC + MedHxAlcStatusNFC + MedHxCokeStatusNFC + MedHxOpiateStatusNFC + CESimp + DMMclusterdich + CVDclusterdich + MedHxMedsBP + PhysBMI + MedHxMarijStatusCurrent*HTN, data = exdata)
summary(model3)
##
## Call:
## lm(formula = PhysMeanDBPsitimp ~ MedHxMarijStatusCurrent + Sex +
## Race + Age0 + HseHldEducation + PovStat + MedHxCigaretteStatusNFC +
## MedHxAlcStatusNFC + MedHxCokeStatusNFC + MedHxOpiateStatusNFC +
## CESimp + DMMclusterdich + CVDclusterdich + MedHxMedsBP +
## PhysBMI + MedHxMarijStatusCurrent * HTN, data = exdata)
##
## Residuals:
## Min 1Q Median 3Q Max
## -32.678437 -6.910081 0.115360 7.005594 36.093898
##
## Coefficients:
## Estimate Std. Error t value
## (Intercept) 67.221761347 2.064743081 32.55696
## MedHxMarijStatusCurrent1 -0.967829373 0.777450239 -1.24488
## SexMen 3.130217005 0.401851305 7.78949
## RaceAfrAm 0.379473038 0.396891173 0.95611
## Age0 -0.035818640 0.022694917 -1.57827
## HseHldEducation -0.085759695 0.067147270 -1.27719
## PovStatBelow 1.030393551 0.408719954 2.52103
## MedHxCigaretteStatusNFC -0.683635624 0.259855835 -2.63083
## MedHxAlcStatusNFC 0.338861723 0.248556232 1.36332
## MedHxCokeStatusNFC 1.120592255 0.409014662 2.73974
## MedHxOpiateStatusNFC -0.955067263 0.490858558 -1.94571
## CESimp 0.002268173 0.017662142 0.12842
## DMMclusterdich1 -1.843047491 0.453738374 -4.06192
## CVDclusterdich1 -2.741260485 0.864921910 -3.16937
## MedHxMedsBP -4.800574521 0.649766563 -7.38815
## PhysBMI 0.142259624 0.026574256 5.35329
## HTNYes 10.872171537 0.639623232 16.99777
## MedHxMarijStatusCurrent1:HTNYes 2.935915545 1.170755446 2.50771
## Pr(>|t|)
## (Intercept) < 2.22e-16
## MedHxMarijStatusCurrent1 0.2132853
## SexMen 9.5191e-15
## RaceAfrAm 0.3391005
## Age0 0.1146215
## HseHldEducation 0.2016457
## PovStatBelow 0.0117585
## MedHxCigaretteStatusNFC 0.0085662
## MedHxAlcStatusNFC 0.1728953
## MedHxCokeStatusNFC 0.0061893
## MedHxOpiateStatusNFC 0.0517935
## CESimp 0.8978262
## DMMclusterdich1 5.0058e-05
## CVDclusterdich1 0.0015448
## MedHxMedsBP 1.9743e-13
## PhysBMI 9.3634e-08
## HTNYes < 2.22e-16
## MedHxMarijStatusCurrent1:HTNYes 0.0122099
##
## Residual standard error: 9.847887 on 2699 degrees of freedom
## (85 observations deleted due to missingness)
## Multiple R-squared: 0.1849221, Adjusted R-squared: 0.1797883
## F-statistic: 36.02001 on 17 and 2699 DF, p-value: < 2.2204e-16
bargraph.CI(exdata$HTN,exdata$PhysMeanDBPsitimp,group=exdata$MedHxMarijStatusCurrent, split=F, col=NULL, angle=NULL, density=NULL,
lc=TRUE, uc=TRUE, legend=T, ncol=1, leg.lab=NULL, x.leg=NULL, y.leg=NULL,
cex.leg=1, bty="n", bg="white", err.width=if(length(levels(as.factor(exdata$MedHxmMarijStatusCurrent)))>10)
0 else .1, err.col="black", err.lty=1, fun = function(x) mean(x, na.rm=TRUE),
xlab="Marijuana Use",ylab="Diastolic Blood Pressure", ylim=c(35,115), xpd=FALSE, cex.axis=1,
yaxt="s", xaxt="s", data=NULL, subset=NULL,main="Marijuana Use by HTN Status")
box()
axis(4,labels=F)
