Unique contribution of Marijuana Use to each model (semi-partial correlations) ========================================================

##Make Marijuana Use Numeric cannabisHyp$MarijuanaUseNum <- ifelse(cannabisHyp$MedHxMarijStatusCurrent == "Current.Use",1,0)

##Make Sex Numeric cannabisHyp$SexNum <- ifelse(cannabisHyp$Sex == "Men",1,0)

##Make Race Numeric cannabisHyp$RaceNum <- ifelse(cannabisHyp$Race == "White",1,0)

##Make PovStat Numeric cannabisHyp$PovNum <- ifelse(cannabisHyp$PovStat == "Above",1,0) save(cannabisHyp,file="~/Desktop/Research/Cannabis and Cardiovascular Disease/cannabisHyp.rdata")

load(file="~/Desktop/Research/Cannabis and Cardiovascular Disease/cannabisHyp.rdata")
library(ppcor)

Systolic Blood Pressure

spcor.test(cannabisHyp$MarijuanaUseNum,cannabisHyp$PhysMeanSBPsitimp,cannabisHyp[,c("SexNum","RaceNum","Age0","HseHldEducation","PovNum","MedHxCigaretteStatusNFC","MedHxAlcStatusNFC","MedHxCokeStatusNFC","MedHxOpiateStatusNFC","CESimp", "DMMclusterdich","CVDclusterdich","MedHxMedsBP","PhysBMI")])

##   estimate   p.value statistic    n gp  Method
## 1  0.06638 0.0005307     3.465 2728 14 pearson

Diastolic Blood Pressure

spcor.test(cannabisHyp$MarijuanaUseNum,cannabisHyp$PhysMeanDBPsitimp,cannabisHyp[,c("SexNum","RaceNum","Age0","HseHldEducation","PovNum","MedHxCigaretteStatusNFC","MedHxAlcStatusNFC","MedHxCokeStatusNFC","MedHxOpiateStatusNFC","CESimp", "DMMclusterdich","CVDclusterdich","MedHxMedsBP","PhysBMI")])

##   estimate p.value statistic    n gp  Method
## 1   0.0243  0.2056     1.266 2728 14 pearson

Pulse Pressure

spcor.test(cannabisHyp$MarijuanaUseNum,cannabisHyp$PhysMeanPPsitimp,cannabisHyp[,c("SexNum","RaceNum","Age0","HseHldEducation","PovNum","MedHxCigaretteStatusNFC","MedHxAlcStatusNFC","MedHxCokeStatusNFC","MedHxOpiateStatusNFC","CESimp", "DMMclusterdich","CVDclusterdich","MedHxMedsBP","PhysBMI")])

##   estimate   p.value statistic    n gp  Method
## 1  0.07048 0.0002337      3.68 2728 14 pearson