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 + PhysBMI + MedHxMarijStatusCurrent*HTN, data = exdata)
summary(model1)
##
## Call:
## lm(formula = PhysMeanSBPsitimp ~ MedHxMarijStatusCurrent + Sex +
## Race + Age0 + HseHldEducation + PovStat + MedHxCigaretteStatusNFC +
## MedHxAlcStatusNFC + MedHxCokeStatusNFC + MedHxOpiateStatusNFC +
## CESimp + DMMclusterdich + CVDclusterdich + PhysBMI + MedHxMarijStatusCurrent *
## HTN, data = exdata)
##
## Residuals:
## Min 1Q Median 3Q Max
## -44.646919 -10.142703 -0.426330 9.732826 65.110968
##
## Coefficients:
## Estimate Std. Error t value
## (Intercept) 88.312533082 3.141144059 28.11477
## MedHxMarijStatusCurrent1 1.233779430 1.188194712 1.03836
## SexMen 2.850856837 0.613626957 4.64591
## RaceAfrAm 1.922727119 0.606151032 3.17203
## Age0 0.341313250 0.034592213 9.86677
## HseHldEducation -0.145598533 0.102620881 -1.41880
## PovStatBelow 1.294006233 0.624481344 2.07213
## MedHxCigaretteStatusNFC -0.431439526 0.396990110 -1.08678
## MedHxAlcStatusNFC 0.422536571 0.378932043 1.11507
## MedHxCokeStatusNFC 0.728724368 0.624643415 1.16662
## MedHxOpiateStatusNFC -2.183716271 0.750000989 -2.91162
## CESimp 0.001285249 0.026988809 0.04762
## DMMclusterdich1 -1.749630195 0.686030494 -2.55037
## CVDclusterdich1 -4.233284454 1.309933889 -3.23168
## PhysBMI 0.357116468 0.040488797 8.82013
## HTNYes 13.456871826 0.703991324 19.11511
## MedHxMarijStatusCurrent1:HTNYes 4.755759164 1.782592223 2.66789
## Pr(>|t|)
## (Intercept) < 2.22e-16
## MedHxMarijStatusCurrent1 0.2991933
## SexMen 3.5481e-06
## RaceAfrAm 0.0015308
## Age0 < 2.22e-16
## HseHldEducation 0.1560727
## PovStatBelow 0.0383481
## MedHxCigaretteStatusNFC 0.2772326
## MedHxAlcStatusNFC 0.2649187
## MedHxCokeStatusNFC 0.2434650
## MedHxOpiateStatusNFC 0.0036251
## CESimp 0.9620214
## DMMclusterdich1 0.0108157
## CVDclusterdich1 0.0012455
## PhysBMI < 2.22e-16
## HTNYes < 2.22e-16
## MedHxMarijStatusCurrent1:HTNYes 0.0076785
##
## Residual standard error: 15.05076 on 2700 degrees of freedom
## (85 observations deleted due to missingness)
## Multiple R-squared: 0.2888468, Adjusted R-squared: 0.2846326
## F-statistic: 68.54065 on 16 and 2700 DF, p-value: < 2.2204e-16
bargraph.CI(exdata$HTN,exdata$PhysMeanSBPsitimp,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="Systolic Blood Pressure", ylim=c(75,200), xpd=FALSE, cex.axis=1,
yaxt="s", xaxt="s", data=NULL, subset=NULL,main="Marijuana Use by HTN Status")
box()
axis(4,labels=F)
Pulse Pressure
model2 = lm(PhysMeanPPsitimp~MedHxMarijStatusCurrent + Sex + Race + Age0 + HseHldEducation + PovStat + MedHxCigaretteStatusNFC + MedHxAlcStatusNFC + MedHxCokeStatusNFC + MedHxOpiateStatusNFC + CESimp + DMMclusterdich + CVDclusterdich + PhysBMI + MedHxMarijStatusCurrent*HTN, data = exdata)
summary(model2)
##
## Call:
## lm(formula = PhysMeanPPsitimp ~ MedHxMarijStatusCurrent + Sex +
## Race + Age0 + HseHldEducation + PovStat + MedHxCigaretteStatusNFC +
## MedHxAlcStatusNFC + MedHxCokeStatusNFC + MedHxOpiateStatusNFC +
## CESimp + DMMclusterdich + CVDclusterdich + PhysBMI + MedHxMarijStatusCurrent *
## HTN, data = exdata)
##
## Residuals:
## Min 1Q Median 3Q Max
## -45.784295 -7.600849 -1.004901 6.492415 51.533973
##
## Coefficients:
## Estimate Std. Error t value
## (Intercept) 19.632483967 2.360360524 8.31758
## MedHxMarijStatusCurrent1 2.209209255 0.892849179 2.47434
## SexMen -0.402906521 0.461099783 -0.87379
## RaceAfrAm 1.653288218 0.455482124 3.62975
## Age0 0.389400312 0.025993744 14.98054
## HseHldEducation -0.056824313 0.077112756 -0.73690
## PovStatBelow 0.192115933 0.469256133 0.40941
## MedHxCigaretteStatusNFC 0.305614835 0.298311624 1.02448
## MedHxAlcStatusNFC 0.212945086 0.284742189 0.74785
## MedHxCokeStatusNFC -0.508209364 0.469377918 -1.08273
## MedHxOpiateStatusNFC -1.147064407 0.563575785 -2.03533
## CESimp 0.001442801 0.020280292 0.07114
## DMMclusterdich1 0.582785576 0.515506219 1.13051
## CVDclusterdich1 -0.634836521 0.984328061 -0.64494
## PhysBMI 0.230265509 0.030424634 7.56839
## HTNYes 5.863393784 0.529002586 11.08387
## MedHxMarijStatusCurrent1:HTNYes 1.071936759 1.339499314 0.80025
## Pr(>|t|)
## (Intercept) < 2.22e-16
## MedHxMarijStatusCurrent1 0.01340942
## SexMen 0.38230781
## RaceAfrAm 0.00028896
## Age0 < 2.22e-16
## HseHldEducation 0.46124780
## PovStatBelow 0.68227475
## MedHxCigaretteStatusNFC 0.30569953
## MedHxAlcStatusNFC 0.45461441
## MedHxCokeStatusNFC 0.27902515
## MedHxOpiateStatusNFC 0.04191485
## CESimp 0.94328920
## DMMclusterdich1 0.25836128
## CVDclusterdich1 0.51901824
## PhysBMI 5.152e-14
## HTNYes < 2.22e-16
## MedHxMarijStatusCurrent1:HTNYes 0.42363532
##
## Residual standard error: 11.30964 on 2700 degrees of freedom
## (85 observations deleted due to missingness)
## Multiple R-squared: 0.2317109, Adjusted R-squared: 0.2271581
## F-statistic: 50.89389 on 16 and 2700 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 + PhysBMI + MedHxMarijStatusCurrent*HTN, data = exdata)
summary(model3)
##
## Call:
## lm(formula = PhysMeanDBPsitimp ~ MedHxMarijStatusCurrent + Sex +
## Race + Age0 + HseHldEducation + PovStat + MedHxCigaretteStatusNFC +
## MedHxAlcStatusNFC + MedHxCokeStatusNFC + MedHxOpiateStatusNFC +
## CESimp + DMMclusterdich + CVDclusterdich + PhysBMI + MedHxMarijStatusCurrent *
## HTN, data = exdata)
##
## Residuals:
## Min 1Q Median 3Q Max
## -33.997012 -7.026092 0.162221 7.001323 35.871242
##
## Coefficients:
## Estimate Std. Error t value
## (Intercept) 68.6800491149 2.0755815723 33.08954
## MedHxMarijStatusCurrent1 -0.9754298257 0.7851263751 -1.24239
## SexMen 3.2537633581 0.4054678104 8.02471
## RaceAfrAm 0.2694389008 0.4005279249 0.67271
## Age0 -0.0480870614 0.0228575826 -2.10377
## HseHldEducation -0.0887742202 0.0678090548 -1.30918
## PovStatBelow 1.1018902998 0.4126400910 2.67034
## MedHxCigaretteStatusNFC -0.7370543614 0.2623201424 -2.80975
## MedHxAlcStatusNFC 0.2095914848 0.2503878683 0.83707
## MedHxCokeStatusNFC 1.2369337319 0.4127471832 2.99683
## MedHxOpiateStatusNFC -1.0366518638 0.4955800190 -2.09180
## CESimp -0.0001575522 0.0178334624 -0.00883
## DMMclusterdich1 -2.3324157715 0.4533100759 -5.14530
## CVDclusterdich1 -3.5984479334 0.8655682739 -4.15732
## PhysBMI 0.1268509584 0.0267538833 4.74140
## HTNYes 7.5934780429 0.4651780983 16.32381
## MedHxMarijStatusCurrent1:HTNYes 3.6838224054 1.1778878967 3.12748
## Pr(>|t|)
## (Intercept) < 2.22e-16
## MedHxMarijStatusCurrent1 0.2142022
## SexMen 1.5025e-15
## RaceAfrAm 0.5011898
## Age0 0.0354909
## HseHldEducation 0.1905850
## PovStatBelow 0.0076228
## MedHxCigaretteStatusNFC 0.0049937
## MedHxAlcStatusNFC 0.4026288
## MedHxCokeStatusNFC 0.0027529
## MedHxOpiateStatusNFC 0.0365500
## CESimp 0.9929517
## DMMclusterdich1 2.8626e-07
## CVDclusterdich1 3.3203e-05
## PhysBMI 2.2326e-06
## HTNYes < 2.22e-16
## MedHxMarijStatusCurrent1:HTNYes 0.0017819
##
## Residual standard error: 9.945128 on 2700 degrees of freedom
## (85 observations deleted due to missingness)
## Multiple R-squared: 0.1684379, Adjusted R-squared: 0.1635102
## F-statistic: 34.18134 on 16 and 2700 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)
