Marijuana Use and Puse Wave Velocity Analyses
Marijuana use and PULSE WAVE VELOCITY - no interactions
##
## Call:
## lm(formula = PWV ~ MarijStatus + Sex + Race + Age0 + HseHldEducation +
## PovStat + MedHxCigaretteStatusNFC + MedHxAlcStatusNFC + MedHxCokeStatusNFC +
## MedHxOpiateStatusNFC + CESimp + DMMclusterdich + CVDclusterdich +
## MedHxMedsBP + PhysBMI, data = exdata)
##
## Residuals:
## Pulse wave velocity
## Min 1Q Median 3Q Max
## -7.19 -1.11 -0.26 0.72 55.12
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 2.76393 0.49333 5.60 2.4e-08
## MarijStatus 0.14208 0.15786 0.90 0.3682
## SexMen 0.31347 0.10500 2.99 0.0029
## RaceAfrAm 0.31665 0.10392 3.05 0.0023
## Age0 0.07331 0.00592 12.39 < 2e-16
## HseHldEducation -0.05382 0.01774 -3.03 0.0024
## PovStatBelow -0.26615 0.10721 -2.48 0.0131
## MedHxCigaretteStatusNFC -0.01194 0.06904 -0.17 0.8627
## MedHxAlcStatusNFC 0.16775 0.06570 2.55 0.0107
## MedHxCokeStatusNFC -0.19148 0.10602 -1.81 0.0710
## MedHxOpiateStatusNFC -0.16869 0.12630 -1.34 0.1818
## CESimp 0.00101 0.00469 0.22 0.8292
## DMMclusterdich 0.63588 0.11961 5.32 1.2e-07
## CVDclusterdich 0.03853 0.23225 0.17 0.8683
## MedHxMedsBP 0.31728 0.12070 2.63 0.0086
## PhysBMI 0.06059 0.00710 8.54 < 2e-16
##
## Residual standard error: 2.38 on 2286 degrees of freedom
## (63 observations deleted due to missingness)
## Multiple R-squared: 0.168, Adjusted R-squared: 0.162
## F-statistic: 30.8 on 15 and 2286 DF, p-value: <2e-16
Marijuana use and PULSE WAVE VELOCITY - (w/ Marijuana Use * Sex interaction)
summary(model2)
##
## Call:
## lm(formula = PWV ~ MarijStatus + Sex + Race + Age0 + HseHldEducation +
## PovStat + MedHxCigaretteStatusNFC + MedHxAlcStatusNFC + MedHxCokeStatusNFC +
## MedHxOpiateStatusNFC + CESimp + DMMclusterdich + CVDclusterdich +
## MedHxMedsBP + PhysBMI + MarijStatus * Sex, data = exdata)
##
## Residuals:
## Pulse wave velocity
## Min 1Q Median 3Q Max
## -7.10 -1.11 -0.26 0.72 55.10
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 2.761808 0.493355 5.60 2.4e-08
## MarijStatus 0.289969 0.229617 1.26 0.2068
## SexMen 0.349055 0.112404 3.11 0.0019
## RaceAfrAm 0.316331 0.103923 3.04 0.0024
## Age0 0.073240 0.005918 12.38 < 2e-16
## HseHldEducation -0.054493 0.017756 -3.07 0.0022
## PovStatBelow -0.267809 0.107230 -2.50 0.0126
## MedHxCigaretteStatusNFC -0.012926 0.069057 -0.19 0.8515
## MedHxAlcStatusNFC 0.165848 0.065741 2.52 0.0117
## MedHxCokeStatusNFC -0.194163 0.106064 -1.83 0.0673
## MedHxOpiateStatusNFC -0.164684 0.126382 -1.30 0.1927
## CESimp 0.000825 0.004697 0.18 0.8606
## DMMclusterdich 0.634420 0.119626 5.30 1.2e-07
## CVDclusterdich 0.034495 0.232305 0.15 0.8820
## MedHxMedsBP 0.318082 0.120708 2.64 0.0085
## PhysBMI 0.060805 0.007100 8.56 < 2e-16
## MarijStatus:SexMen -0.264529 0.298249 -0.89 0.3752
##
## Residual standard error: 2.38 on 2285 degrees of freedom
## (63 observations deleted due to missingness)
## Multiple R-squared: 0.168, Adjusted R-squared: 0.162
## F-statistic: 28.9 on 16 and 2285 DF, p-value: <2e-16
Marijuana use and PULSE WAVE VELOCITY - (w/ Marijuana Use * Race interaction)
summary(model3)
##
## Call:
## lm(formula = PWV ~ MarijStatus + Sex + Race + Age0 + HseHldEducation +
## PovStat + MedHxCigaretteStatusNFC + MedHxAlcStatusNFC + MedHxCokeStatusNFC +
## MedHxOpiateStatusNFC + CESimp + DMMclusterdich + CVDclusterdich +
## MedHxMedsBP + PhysBMI + MarijStatus * Race, data = exdata)
##
## Residuals:
## Pulse wave velocity
## Min 1Q Median 3Q Max
## -7.13 -1.12 -0.26 0.73 55.11
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 2.74662 0.49387 5.56 3.0e-08
## MarijStatus 0.30806 0.26552 1.16 0.2461
## SexMen 0.31387 0.10501 2.99 0.0028
## RaceAfrAm 0.34489 0.11009 3.13 0.0018
## Age0 0.07334 0.00592 12.39 < 2e-16
## HseHldEducation -0.05386 0.01774 -3.04 0.0024
## PovStatBelow -0.26478 0.10723 -2.47 0.0136
## MedHxCigaretteStatusNFC -0.01188 0.06905 -0.17 0.8634
## MedHxAlcStatusNFC 0.16800 0.06571 2.56 0.0106
## MedHxCokeStatusNFC -0.19476 0.10611 -1.84 0.0666
## MedHxOpiateStatusNFC -0.16627 0.12634 -1.32 0.1883
## CESimp 0.00103 0.00469 0.22 0.8262
## DMMclusterdich 0.63478 0.11963 5.31 1.2e-07
## CVDclusterdich 0.03923 0.23227 0.17 0.8659
## MedHxMedsBP 0.31810 0.12071 2.64 0.0085
## PhysBMI 0.06058 0.00710 8.54 < 2e-16
## MarijStatus:RaceAfrAm -0.24501 0.31514 -0.78 0.4370
##
## Residual standard error: 2.38 on 2285 degrees of freedom
## (63 observations deleted due to missingness)
## Multiple R-squared: 0.168, Adjusted R-squared: 0.162
## F-statistic: 28.9 on 16 and 2285 DF, p-value: <2e-16
Marijuana use and PULSE WAVE VELOCITY - (w/ Marijuana Use * Age interaction)
summary(model4)
##
## Call:
## lm(formula = PWV ~ MarijStatus + Sex + Race + Age0 + HseHldEducation +
## PovStat + MedHxCigaretteStatusNFC + MedHxAlcStatusNFC + MedHxCokeStatusNFC +
## MedHxOpiateStatusNFC + CESimp + DMMclusterdich + CVDclusterdich +
## MedHxMedsBP + PhysBMI + MarijStatus * Age0, data = exdata)
##
## Residuals:
## Pulse wave velocity
## Min 1Q Median 3Q Max
## -7.06 -1.11 -0.27 0.72 55.12
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 2.662438 0.503604 5.29 1.4e-07
## MarijStatus 0.912628 0.784582 1.16 0.2449
## SexMen 0.309778 0.105062 2.95 0.0032
## RaceAfrAm 0.316315 0.103918 3.04 0.0024
## Age0 0.075330 0.006249 12.05 < 2e-16
## HseHldEducation -0.053376 0.017744 -3.01 0.0027
## PovStatBelow -0.265508 0.107211 -2.48 0.0133
## MedHxCigaretteStatusNFC -0.014565 0.069094 -0.21 0.8331
## MedHxAlcStatusNFC 0.169995 0.065741 2.59 0.0098
## MedHxCokeStatusNFC -0.179296 0.106710 -1.68 0.0931
## MedHxOpiateStatusNFC -0.166693 0.126311 -1.32 0.1871
## CESimp 0.000845 0.004695 0.18 0.8572
## DMMclusterdich 0.635987 0.119608 5.32 1.2e-07
## CVDclusterdich 0.036601 0.232257 0.16 0.8748
## MedHxMedsBP 0.315831 0.120708 2.62 0.0089
## PhysBMI 0.060597 0.007095 8.54 < 2e-16
## MarijStatus:Age0 -0.017204 0.017159 -1.00 0.3162
##
## Residual standard error: 2.38 on 2285 degrees of freedom
## (63 observations deleted due to missingness)
## Multiple R-squared: 0.168, Adjusted R-squared: 0.162
## F-statistic: 28.9 on 16 and 2285 DF, p-value: <2e-16
PWV and PP correlation coefficient
cor.test(exdata$PhysMeanPPsitimp,exdata$PWV)
##
## Pearson's product-moment correlation
##
## data: exdata$PhysMeanPPsitimp and exdata$PWV
## t = 17.64, df = 2363, p-value < 2.2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
## 0.3050 0.3763
## sample estimates:
## cor
## 0.3411