************************************************************************
* 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 csz1fmt 1='Infant/Child'
	            2='Adult'
	            3='Large'
	            4='Thigh'
                ;
  VALUE csz2fmt 1,2='Infant/Child'
	            3='Adult'
	            4='Large'
	            5='Thigh'
                ;
  VALUE selfmt  1='Age ge 20'
                2='Age lt 20'; 
	
data analysis_data;
set nh.analysis_data;

if ridstatr = 2;  ***examined  ;

/*variable created for study subpopulation of interest:used on domain statement for SAS Survey procedure*/
if ridageyr >=20 then sel=1;
else sel=2;

LABEL                           
       riagendr = 'Gender'  
	   sel = 'Subpop: 1=age ge 20 2=age lt 20'
          ;                                                       
                     
run;

title ' '; 








proc surveyfreq 
data=analysis_data;

  strata sdmvstra;

  cluster sdmvpsu;

  weight wtmec4yr;

table sel*riagendr*bpacsz/col row 
	  nostd nowt wchisq wllchisq chisq chisq1; 
format riagendr sexfmt. bpacsz csz2fmt. sel selfmt.; 
run;