************************************************************************
* Program:  C:\NHANES\.sas                                          *
* Proposal:                                                          * 
************************************************************************;
LIBNAME NH "C:\NHANES\DATA";
OPTIONS NODATE NOCENTER;
options ls=72;

proc format;                                       
                                                   
  VALUE sexfmt  1 = 'Male'                          
                2 = 'Female'                        
                ;                                  
  

  VALUE racefmt  1 = 'NH White'
                 2 = 'NH Black'
			     3 = 'Mexican American'
			     4 = 'Other'
			      ;
                       
  VALUE agefmt                             
                1 = '20-39'                               
                2 = '40-59'                               
				3 = '60+'       
				;         

 run;   
data ANALYSIS_DATA;
set NH.ANALYSIS_DATA; 

if ridstatr = 2;  ***examined  ;

if ridageyr >=20 then sel=1;
else sel=2;
 
age = .;
if 20 LE ridageyr LE 39 then age=1;
if 40 LE ridageyr LE 59 then age=2;   
if ridageyr GE 60 then age=3;  

race=.;
if ridreth1=3 then race=1;
if ridreth1=4 then race=2; 
if ridreth1=1 then race=3;
if ridreth1=2 or ridreth1=5 then race=4;

LABEL                           
       age  = 'AGE GROUP'    
	   race = 'Race Ethnicity'	
       riagendr = 'Gender'  
          ;                                                       
                                                                                                              
run;

proc sort data=analysis_data;
  by sdmvstra sdmvpsu; 
run;

proc descript 
data=analysis_data design=wr;

  nest sdmvstra sdmvpsu;

  weight wtmec4yr;

  subpopn ridageyr >= 20;

  class riagendr/NoFREQ;

  var bpxsar;

  print nsum mean semean/style=nchs;

rformat riagendr sexfmt.;
rtitle "Mean systolic blood pressure: NHANES 1999-2002"; 
run;

proc sort data=analysis_data;
  by sdmvstra sdmvpsu; 
run;

proc descript 
data=analysis_data design=wr;

  nest sdmvstra sdmvpsu;

  weight wtmec4yr;

  subpopn ridageyr >= 20;

  class riagendr/NoFREQ;

  var bpxsar;

  contrast riagendr = (1 -1)/name = "Males vs. Females";

  print nsum t_mean p_mean/style=nchs;

rformat riagendr sexfmt.;
rtitle "Significance test for difference between mean systolic blood pressure for males and females";
rtitle2  "NHANES 1999-2002"; 
run;