Main Effect Plots for Cannabis and BP

library(effects)

## Loading required package: colorspace
## 
## Attaching package: 'effects'
## 
## The following object is masked from 'package:car':
## 
##     Prestige

load(file ="/Users/meganwilliams/Desktop/Research/Cannabis and Cardiovascular Disease/exdata.rdata")

lm1 = lm(PhysMeanSBPsitimp~MedHxMarijStatusCurrent + Sex + Race + Age0 + HseHldEducation + PovStat + MedHxCigaretteStatusNFC + MedHxAlcStatusNFC + MedHxCokeStatusNFC + MedHxOpiateStatusNFC + CESimp + DMMclusterdich + CVDclusterdich + MedHxMedsBP + PhysBMI,exdata)

lm2 = lm(PhysMeanPPsitimp~MedHxMarijStatusCurrent + Sex + Race + Age0 + HseHldEducation + PovStat + MedHxCigaretteStatusNFC + MedHxAlcStatusNFC + MedHxCokeStatusNFC + MedHxOpiateStatusNFC + CESimp + DMMclusterdich + CVDclusterdich + MedHxMedsBP + PhysBMI,exdata)

lm3 = lm(PhysMeanDBPsitimp~MedHxMarijStatusCurrent + Sex + Race + Age0 + HseHldEducation + PovStat + MedHxCigaretteStatusNFC + MedHxAlcStatusNFC + MedHxCokeStatusNFC + MedHxOpiateStatusNFC + CESimp + DMMclusterdich + CVDclusterdich + MedHxMedsBP + PhysBMI,exdata)

Systolic Blood Pressure

Main Effect Plot for Marijuana Use

plot(Effect(c("MedHxMarijStatusCurrent"),lm1))

plot of chunk unnamed-chunk-2

Plot for Marijuana Use by Age

plot(Effect(c("MedHxMarijStatusCurrent","Age0"),lm1))

plot of chunk unnamed-chunk-3

Plot for Marijuana Use by Race

plot(Effect(c("MedHxMarijStatusCurrent","Race"),lm1))

plot of chunk unnamed-chunk-4

Plot for Marijuana Use by Gender

plot(Effect(c("MedHxMarijStatusCurrent","Sex"),lm1))

plot of chunk unnamed-chunk-5

Pulse Pressure

Main Effect Plot for Marijuana Use

plot(Effect(c("MedHxMarijStatusCurrent"),lm2))

plot of chunk unnamed-chunk-6

Plot for Marijuana Use by Age

plot(Effect(c("MedHxMarijStatusCurrent","Age0"),lm2))

plot of chunk unnamed-chunk-7

Plot for Marijuana Use by Race

plot(Effect(c("MedHxMarijStatusCurrent","Race"),lm2))

plot of chunk unnamed-chunk-8

Plot for Marijuana Use by Gender

plot(Effect(c("MedHxMarijStatusCurrent","Sex"),lm2))

plot of chunk unnamed-chunk-9

Diastolic Blood Pressure

Main Effect Plot for Marijuana Use

plot(Effect(c("MedHxMarijStatusCurrent"),lm3))

plot of chunk unnamed-chunk-10