LIBNAME NH "C:\NHANES\DATA";

*-------------------------------------------------------------------------;
* Use the PROC FORMAT procedure to assign text labels to the numeric      ;
* values of user-defined formats                                          ;
*-------------------------------------------------------------------------;

proc format;
    value YESNOS 
        1='Used any calcium supp'
        2='No supp use'; 
    value GENDER 
        1='Male'
        2='Female';
run;

*--------------------------------------------------------------------------;
* Create a dataset from the permanent dataset libname.demoadv             ;
*--------------------------------------------------------------------------;

data demoadv;
  set nh.demoadv;
  format riagendr gender. anycalsup yesnos.;
run;

*--------------------------------------------------------------------------;
* Sort the data by strata and PSU                                          ;
*--------------------------------------------------------------------------;

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

*--------------------------------------------------------------------------;
* Use proc descript to generate proportions                                ;
*--------------------------------------------------------------------------;

proc descript data=demoadv design=wr ;
  subpopn ridageyr >= 50;
  nest sdmvstra sdmvpsu; 
  weight wtint2yr;
  subgroup riagendr;
  levels 2;
  var anycalsup; 
  catlevel 1;
  table riagendr;
  print nsum="Sample Size" percent="Percent" sepercent="SE" /
    nohead notime style=NCHS nsumfmt=f8.0 percentfmt=f8.4 sepercentfmt=f8.4;
  rtitle "Prevalence of SPs age 50 and older who report calcium supplement use";
run;