clear
capture log close
set more 1
capture log using NHANES_analysis.log, replace

****************************************************
*     Data Analysis  - Sample Stata Program        *
*       For CSP Summer Seminar                     *
****************************************************


** Analyses for Patient Factors Associated with dichotomous dependent variable: numbness in past 3 months

use NHANES_transformed

** Make a global macro commands
global demog "age female race"
global social "educ3 income"
global dmhealth "duration bmxbmi meds3 healthst"
global xiall "age female i.race i.educ3 i.income duration bmxbmi i.meds3 i.healthst" //need separate global when using "i."


**For Tables 1 & 2: Descriptive Statistics

sum a1c age duration bmxbmi
for var a1c_3 female race educ3 income meds3 healthst: tab X


**again taking into account dependent var

bysort numb: sum a1c age duration bmxbmi
for var a1c_3 female race educ3 income meds3 healthst: tab X numb, chi ro co




** For Table 3:  Results for DepVar = numb (Dichotomous DepVar)
** Chi-squares for DepVar = numb
tab numb meds3, co ro ch
tab numb a1c_3, co ro ch

** Mean A1cs 
tab numb,sum(a1c) 


** Get Odds Ratios
* Bivariate Comparisons
logistic numb a1c