Crude_Associations.Rmd

---
title: "Crude Associations"
author: "Bryon Langford, Matthew Hoctor"
date: "7/20/2021"
output:
  html_document:
    number_sections: no
    theme: lumen
    toc: yes
    toc_float:
      collapsed: yes
      smooth_scroll: no
  pdf_document:
    toc: yes
  word_document:
    toc: yes
---

```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = TRUE)

if (Sys.info()["sysname"] == "Windows") 
  knitr::opts_chunk$set(engine.path = list(stata = "C:/Program Files/Stata17/StataBE-64")) 
if ((Sys.info()["sysname"] == "Darwin"))
  knitr::opts_chunk$set(engine.path = list(stata = "/Applications/Stata/StataIC"))
```

```{r libraries, include=FALSE}
library(Statamarkdown)
```

# Associations between MJ exposure and covariates

## Gender

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): ta gndr MJ,col percent
```

## Education

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): ta EDUC_cat MJ,col percent
```

## Race

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): ta race_eth MJ,col percent
```

## Survey Year

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): ta sddsrvyr MJ,col percent
svy, subpop(if include==1): ta sddsrvyr MJ,row percent
svy: ta sddsrvyr MJ,row percent
```

## Income

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): mean indfmpir,over(MJ)coeflegend
# svy, subpop(if include==1): mean indfmpir,over(MJ)

test _b[c.indfmpir@0bn.MJ] = _b[c.indfmpir@1.MJ] = _b[c.indfmpir@2.MJ] = _b[c.indfmpir@3.MJ] = _b[c.indfmpir@4.MJ]
```

## Age

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): mean ridageyr,over(MJ)coeflegend
# svy, subpop(if include==1): mean ridageyr,over(MJ)

test _b[c.ridageyr@0bn.MJ] = _b[c.ridageyr@1.MJ] = _b[c.ridageyr@2.MJ] = _b[c.ridageyr@3.MJ] = _b[c.ridageyr@4.MJ]
```

## Smoking

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): ta SMK_cat MJ,col percent
```

## Alcohol

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): ta AL_cat MJ,col percent
```

## BMI

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): mean bmxbmi,over(MJ)coeflegend
# svy, subpop(if include==1): mean bmxbmi,over(MJ)

test _b[c.bmxbmi@0bn.MJ] = _b[c.bmxbmi@1.MJ] = _b[c.bmxbmi@2.MJ] = _b[c.bmxbmi@3.MJ] = _b[c.bmxbmi@4.MJ]
```

## Healthy Eating Index

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): mean hei2015,over(MJ)coeflegend
# svy, subpop(if include==1): mean hei2015,over(MJ)

test _b[c.hei2015@0bn.MJ] = _b[c.hei2015@1.MJ] = _b[c.hei2015@2.MJ] = _b[c.hei2015@3.MJ] = _b[c.hei2015@4.MJ]
```

# Associations between outcome (HTN category) & covariates

## MJ Category

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): ta MJ BP_cat,col percent
```

## Gender

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): ta gndr BP_cat,col percent
```

## Education

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): ta EDUC_cat BP_cat,col percent
```

## Race

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): ta race_eth BP_cat,col percent
```

## Survey Year

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): ta sddsrvyr BP_cat,col percent
svy, subpop(if include==1): ta sddsrvyr BP_cat,row percent
svy: ta sddsrvyr BP_cat,row percent
```

## Income

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): mean indfmpir,over(BP_cat)coeflegend
# svy, subpop(if include==1): mean indfmpir,over(BP_cat)

test _b[c.indfmpir@0bn.BP_cat] = _b[c.indfmpir@1.BP_cat] = _b[c.indfmpir@2.BP_cat] = _b[c.indfmpir@3.BP_cat]
```

## Age

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): mean ridageyr,over(BP_cat)coeflegend
# svy, subpop(if include==1): mean ridageyr,over(BP_cat)

test _b[c.ridageyr@0bn.BP_cat] = _b[c.ridageyr@1.BP_cat] = _b[c.ridageyr@2.BP_cat] = _b[c.ridageyr@3.BP_cat]
```

## Smoking

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): ta SMK_cat BP_cat,col percent
```

## Alcohol

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): ta AL_cat BP_cat,col percent
```

## BMI

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): mean bmxbmi,over(BP_cat)coeflegend
# svy, subpop(if include==1): mean bmxbmi,over(BP_cat)

test _b[c.bmxbmi@0bn.BP_cat] = _b[c.bmxbmi@1.BP_cat] = _b[c.bmxbmi@2.BP_cat] = _b[c.bmxbmi@3.BP_cat]
```

## Healthy Eating Index

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): mean hei2015,over(BP_cat)coeflegend
# svy, subpop(if include==1): mean hei2015,over(BP_cat)

test _b[c.hei2015@0bn.BP_cat] = _b[c.hei2015@1.BP_cat] = _b[c.hei2015@2.BP_cat] = _b[c.hei2015@3.BP_cat]
```

# Crude coefficients and ORs

## Coefficients

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

xi: svy,subpop(if include==1): logit BP_cat i.MJ
```

## ORs

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

xi: svy,subpop(if include==1): logit BP_cat i.MJ,or
```

# Correlations among potential confounders

## Gender

### Education

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): ta gndr EDUC_cat,row percent
```

### Race

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): ta gndr race_eth,row percent
```

### Survey Year

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): ta gndr sddsrvyr,row percent
```



### Smoking

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): ta gndr SMK_cat,row percent
```

### Alcohol

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): ta gndr AL_cat,row percent
```

### BMI

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): mean bmxbmi,over(gndr)coeflegend

test _b[c.bmxbmi@0bn.gndr] = _b[c.bmxbmi@1.gndr]
```

### HEI

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): mean hei2015,over(gndr)coeflegend

test _b[c.hei2015@0bn.gndr] = _b[c.hei2015@1.gndr]
```

### Income

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): mean indfmpir,over(gndr)coeflegend

test _b[c.indfmpir@0bn.gndr] = _b[c.indfmpir@1.gndr]
```

### Age

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): mean ridageyr,over(gndr)coeflegend

test _b[c.ridageyr@0bn.gndr] = _b[c.ridageyr@1.gndr]
```

## Education

### Race

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): ta EDUC_cat race_eth,row percent
```

### Survey Year

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): ta EDUC_cat sddsrvyr,row percent
```



### Smoking

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): ta EDUC_cat SMK_cat,row percent
```

### Alcohol

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): ta EDUC_cat AL_cat,row percent
```

### BMI

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): mean bmxbmi,over(EDUC_cat)coeflegend

test _b[c.bmxbmi@1bn.EDUC_cat] = _b[c.bmxbmi@2.EDUC_cat] = _b[c.bmxbmi@3.EDUC_cat] = _b[c.bmxbmi@4.EDUC_cat] = _b[c.bmxbmi@5.EDUC_cat]
```

### HEI

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): mean hei2015,over(EDUC_cat)coeflegend

test _b[c.hei2015@1bn.EDUC_cat] = _b[c.hei2015@2.EDUC_cat] = _b[c.hei2015@3.EDUC_cat] = _b[c.hei2015@4.EDUC_cat] = _b[c.hei2015@5.EDUC_cat]
```

### Income

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): mean indfmpir,over(EDUC_cat)coeflegend

test _b[c.indfmpir@1bn.EDUC_cat] = _b[c.indfmpir@2.EDUC_cat] = _b[c.indfmpir@3.EDUC_cat] = _b[c.indfmpir@4.EDUC_cat] = _b[c.indfmpir@5.EDUC_cat]
```

### Age

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): mean ridageyr,over(EDUC_cat)coeflegend

test _b[c.ridageyr@1bn.EDUC_cat] = _b[c.ridageyr@2.EDUC_cat] = _b[c.ridageyr@3.EDUC_cat] = _b[c.ridageyr@4.EDUC_cat] = _b[c.ridageyr@5.EDUC_cat]
```

## Race

### Survey Year

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): ta race_eth sddsrvyr,row percent
```

### Smoking

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): ta race_eth SMK_cat,row percent
```

### Alcohol

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): ta race_eth AL_cat,row percent
```

### BMI

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): mean bmxbmi,over(race_eth)coeflegend

test _b[c.bmxbmi@1bn.race_eth] = _b[c.bmxbmi@2.race_eth] = _b[c.bmxbmi@3.race_eth] = _b[c.bmxbmi@4.race_eth]
```

### HEI

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): mean hei2015,over(race_eth)coeflegend

test _b[c.hei2015@1bn.race_eth] = _b[c.hei2015@2.race_eth] = _b[c.hei2015@3.race_eth] = _b[c.hei2015@4.race_eth]
```

### Income

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): mean indfmpir,over(race_eth)coeflegend

test _b[c.indfmpir@1bn.race_eth] = _b[c.indfmpir@2.race_eth] = _b[c.indfmpir@3.race_eth] = _b[c.indfmpir@4.race_eth]
```

### Age

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): mean ridageyr,over(race_eth)coeflegend

test _b[c.ridageyr@1bn.race_eth] = _b[c.ridageyr@2.race_eth] = _b[c.ridageyr@3.race_eth] = _b[c.ridageyr@4.race_eth]
```

## Smoking

### Survey Year

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): ta SMK_cat sddsrvyr,row percent
```

### Alcohol

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): ta SMK_cat AL_cat,row percent
```

### BMI

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): mean bmxbmi,over(SMK_cat)coeflegend

test _b[c.bmxbmi@1bn.SMK_cat] = _b[c.bmxbmi@2.SMK_cat] = _b[c.bmxbmi@3.SMK_cat] = _b[c.bmxbmi@4.SMK_cat]
```

### HEI

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): mean hei2015,over(SMK_cat)coeflegend

test _b[c.hei2015@1bn.SMK_cat] = _b[c.hei2015@2.SMK_cat] = _b[c.hei2015@3.SMK_cat] = _b[c.hei2015@4.SMK_cat]
```

### Income

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): mean indfmpir,over(SMK_cat)coeflegend

test _b[c.indfmpir@1bn.SMK_cat] = _b[c.indfmpir@2.SMK_cat] = _b[c.indfmpir@3.SMK_cat] = _b[c.indfmpir@4.SMK_cat]
```

### Age

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): mean ridageyr,over(SMK_cat)coeflegend

test _b[c.ridageyr@1bn.SMK_cat] = _b[c.ridageyr@2.SMK_cat] = _b[c.ridageyr@3.SMK_cat] = _b[c.ridageyr@4.SMK_cat]
```

## Alcohol

### Survey Year

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): ta AL_cat sddsrvyr,row percent
```

### BMI

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): mean bmxbmi,over(AL_cat)coeflegend

test _b[c.bmxbmi@0bn.AL_cat] = _b[c.bmxbmi@1bn.AL_cat] = _b[c.bmxbmi@2.AL_cat]
```

### HEI

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): mean hei2015,over(AL_cat)coeflegend

test _b[c.hei2015@0bn.AL_cat] = _b[c.hei2015@1bn.AL_cat] = _b[c.hei2015@2.AL_cat]
```

### Income

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): mean indfmpir,over(AL_cat)coeflegend

test _b[c.indfmpir@0bn.AL_cat] = _b[c.indfmpir@1bn.AL_cat] = _b[c.indfmpir@2.AL_cat]
```

### Age

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): mean ridageyr,over(AL_cat)coeflegend

test _b[c.ridageyr@0bn.AL_cat] = _b[c.ridageyr@1bn.AL_cat] = _b[c.ridageyr@2.AL_cat]
```

## Survey Year

### BMI

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): mean bmxbmi,over(sddsrvyr)coeflegend

test _b[c.bmxbmi@4bn.sddsrvyr] = _b[c.bmxbmi@5.sddsrvyr] = _b[c.bmxbmi@6.sddsrvyr] = _b[c.bmxbmi@7.sddsrvyr] = _b[c.bmxbmi@8.sddsrvyr] = _b[c.bmxbmi@9.sddsrvyr] = _b[c.bmxbmi@10.sddsrvyr]
```

### HEI

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): mean hei2015,over(sddsrvyr)coeflegend

test _b[c.hei2015@4bn.sddsrvyr] = _b[c.hei2015@5.sddsrvyr] = _b[c.hei2015@6.sddsrvyr] = _b[c.hei2015@7.sddsrvyr] = _b[c.hei2015@8.sddsrvyr] = _b[c.hei2015@9.sddsrvyr] = _b[c.hei2015@10.sddsrvyr]
```

### Income

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): mean indfmpir,over(sddsrvyr)coeflegend

test _b[c.indfmpir@4bn.sddsrvyr] = _b[c.indfmpir@5.sddsrvyr] = _b[c.indfmpir@6.sddsrvyr] = _b[c.indfmpir@7.sddsrvyr] = _b[c.indfmpir@8.sddsrvyr] = _b[c.indfmpir@9.sddsrvyr] = _b[c.indfmpir@10.sddsrvyr]
```

### Age

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

svy, subpop(if include==1): mean ridageyr,over(sddsrvyr)coeflegend

test _b[c.ridageyr@4bn.sddsrvyr] = _b[c.ridageyr@5.sddsrvyr] = _b[c.ridageyr@6.sddsrvyr] = _b[c.ridageyr@7.sddsrvyr] = _b[c.ridageyr@8.sddsrvyr] = _b[c.ridageyr@9.sddsrvyr] = _b[c.ridageyr@10.sddsrvyr]
```

## continuous varaibles

```{stata}
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

pwcorr indfmpir bmxbmi hei2015 ridageyr if include==1 [w=wtmec12yr]
```