************************************************************************
* Program:  C:\NHANES\hypothesistesting_chisquare_SUDAAN.sas           *
* Proposal: Perform Chi-Square Test Using SUDAAN                       * 
************************************************************************;
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'
                ;
	
data analysis_data;
set nh.analysis_data;

if ridstatr = 2;  ***examined  ;

run;

title ' '; 



proc sort data=analysis_data;
  by sdmvstra sdmvpsu;

proc crosstab 
data=analysis_data design=wr;

  nest sdmvstra sdmvpsu;

  weight wtmec4yr;

  subpopn ridageyr >= 20;

  recode bpacsz = (1 3 4 5);

  class riagendr bpacsz/NoFreq;

  table riagendr*bpacsz;

  print nsum rowper colper/stest =all;

  test chisq;

rformat riagendr sexfmt.;
rformat bpacsz csz1fmt.;
rtitle "Chi-square test for blood pressure cuff size: NHANES 1999-2002";
run;