/*fluorosis.sas*/
/*June 2, 2005*/
/*Laurie Barker*/
/*Program to recreate estimates for NHANES 1999-2002 MMWR*/
/*Fluorosis and Dean's Index*/

/*Reads NH4MMWR */
/*Fluorosis data not collected in NHANES 1988-1994*/

/*Variables used in this section */
/*
  RIDAGEYR   Age in Years  
  RIAGENDR   Gender  
  RIDRETH2   Race/Ethnicity  
  PIRGRP3    Poverty in 3 groups  
  WTMEC4YR   Weights for 1999-2002 MEC examined  
  SDDSRVYR   Survey cycle  
  SDMVSTRA   Strata  
  SDMVPSU    PSU  

  ASDIST5B  ="Dist 5B 6 to 39 - age and sex, 4 age groups - overall"
  ASDIST5D  ="Dist 5D 6 to 39 - age and sex, single years of age - overall"
  OHDDFIX ="Deans Fluorosis Index - cannot be assessed set to . "
  OHDDFIaX ="Deans Fluorosis Index (mod and severe combined) - cannot be assessed set to . "
  OHDCFIX ="Community Fluorosis Index - cannot be assessed set to . "
*/

proc format;

value sex
      1="Male"
      2="Female"
;
value phase
      1="1999-2000"
      2="2001-2002"
;
value age
      1="6-11"
	  2="12-15"
	  3="16-19"
	  4="20-39"
;
value race
      1="NH White"
      2="NH Black"
      3="Mexican American"
	  4="Other/including Multi-racial"
	  5="Other Hispanic"
;
value pir
      1="<100"
	  2="100-199"
	  3="200+"
;
value yesno
      1="yes"
	  2="No"
;
run;

/*NHANES 1999-2002*/
LIBNAME NHANESAB "K:/DOH/EVERYONE/NHANES/NHANES19992002/AnalysisFiles";
proc sort data=NHANESAB.nh4mmwr;
by SDMVSTRA SDMVPSU;
run;

/* Table E1. Prevalence and severity of enamel fluorosis among persons aged 6 to 39 years,
             by selected characteristics -- United States 1999-2002 
             one column for each category
             no, questionable, very mild, mild, moderate, severe*/
proc descript data=NHANESAB.nh4mmwr design=wr;/*E1. fluorosis prevalence [No]  - not age adjusted*/
weight WTMEC4YR;
nest SDMVSTRA SDMVPSU/strlev=1 psulev=2 missunit;
var OHDDFISX;
catlevel 1;
subpopn ridageyr >=6 and ridageyr<=39;
subgroup  RIDRETH2  ASDIST5B  RIAGENDR PIRGRP3;
level         5         4          2        3 ;
tables    ASDIST5B RIAGENDR RIDRETH2 PIRGRP3 
          ASDIST5B*RIAGENDR ASDIST5B*RIDRETH2 ASDIST5B*PIRGRP3;
rformat RIDRETH2 race. ;
rformat riagendr sex. ;
rformat ASDIST5B age. ;
rformat PIRGRP3 pir. ;

print nsum wsum total percent sepercent / style =nchs totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
output nsum wsum total percent sepercent / 
       filename=NHANESAB.N4E1Au REPLACE filetype=SAS 
       totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
rtitle "E1. N IV phase I + II not age/sex adjusted (1 yr) fluorosis prevalence [No] among ages 6 to 39 (N4E1Au)";
run;
data NHANESAB.N4E1Au;
set  NHANESAB.N4E1Au;
format tablename $8.;
tablename="N4E1Au   ";
run;
proc descript data=NHANESAB.nh4mmwr design=wr;/*E1. fluorosis prevalence [No]   - age adjusted*/
weight WTMEC4YR;
nest SDMVSTRA SDMVPSU/strlev=1 psulev=2 missunit;
var OHDDFISX;
catlevel 1;
subpopn ridageyr >=6 and ridageyr<=39;
subgroup  RIDRETH2  ASDIST5B  RIAGENDR PIRGRP3 ASDIST5F;
level         5         4          2        3    34;
tables      RIAGENDR 
          ASDIST5B*RIAGENDR;
Stdvar ASDIST5F;    
Stdwgt	0.029554 0.030007 0.029449 0.031862 0.031908 0.030384
	    0.029728 0.029592 0.029675 0.029811 0.029071 0.030263
		0.029698 0.030656 0.030459 0.028400 0.026824 0.026318
		0.025722 0.025805 0.025104 0.025911 0.026032 0.030829
        0.030165 0.021196 0.028468 0.028973 0.030829 0.032345
        0.032813 0.033710 0.032360 0.036079
;
rformat RIDRETH2 race. ;
rformat riagendr sex. ;
rformat ASDIST5B age. ;
rformat PIRGRP3 pir. ;

print nsum wsum total percent sepercent / style =nchs totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
output nsum wsum total percent sepercent / 
       filename=NHANESAB.N4E1Aaa REPLACE filetype=SAS 
       totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2  ;
rtitle "E1. N IV phase I + II age adjusted (1 yr) fluorosis prevalence [No] among ages 6 to 39 (N4E1Aaa)";
run;
data NHANESAB.N4E1Aaa;
set  NHANESAB.N4E1Aaa;
format tablename $8.;
tablename="N4E1Aaa  ";
run;
proc descript data=NHANESAB.nh4mmwr design=wr;/*E1. fluorosis prevalence [No]   - age/sex adjusted*/
weight WTMEC4YR;
nest SDMVSTRA SDMVPSU/strlev=1 psulev=2 missunit;
var OHDDFISX;
catlevel 1;
subpopn ridageyr >=6 and ridageyr<=39;
subgroup  RIDRETH2  ASDIST5B  RIAGENDR PIRGRP3 ASDIST5D;
level         5         4          2        3    68;
tables    ASDIST5B  RIDRETH2 PIRGRP3 
          ASDIST5B*RIDRETH2 ASDIST5B*PIRGRP3;
Stdvar ASDIST5D;    
Stdwgt 0.015018 0.014307 0.015250 0.014539 0.014943 0.014284 0.016208 0.015415
       0.016246 0.015422 0.015445 0.014711 0.015093 0.014419 0.015025 0.014351
       0.015100 0.014351 0.015168 0.014419 0.014831 0.014022 0.015467 0.014569
       0.015070 0.014404 0.015489 0.014943 0.015400 0.014831 0.014351 0.013835
       0.013513 0.013109 0.013199 0.012929 0.012847 0.012682 0.012862 0.012749
       0.012427 0.012487 0.012832 0.012884 0.012862 0.012974 0.015212 0.015385
       0.014868 0.015070 0.014172 0.014351 0.014037 0.014224 0.014232 0.014524
       0.015265 0.015332 0.015984 0.016118 0.016193 0.016365 0.016642 0.016815
       0.015969 0.016156 0.017878 0.017926 	
; 
rformat RIDRETH2 race. ;
rformat riagendr sex. ;
rformat ASDIST5B age. ;
rformat PIRGRP3 pir. ;

print nsum wsum total percent sepercent / style =nchs totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
output nsum wsum total percent sepercent/ 
       filename=NHANESAB.N4E1Aasa REPLACE filetype=SAS 
       totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
rtitle "E1. N IV phase I + II age/sex adjusted (1 yr) fluorosis prevalence [No] among ages 6 to 39 (N4E1Aasa)";
run;
data NHANESAB.N4E1Aasa;
set  NHANESAB.N4E1Aasa;
format tablename $8.;
tablename="N4E1Aasa ";
run;

proc descript data=NHANESAB.nh4mmwr design=wr;/*E1. fluorosis prevalence [Questionable]  - not age adjusted*/
weight WTMEC4YR;
nest SDMVSTRA SDMVPSU/strlev=1 psulev=2 missunit;
var OHDDFISX;
catlevel 2;
subpopn ridageyr >=6 and ridageyr<=39;
subgroup  RIDRETH2  ASDIST5B  RIAGENDR PIRGRP3;
level         5         4          2        3 ;
tables    ASDIST5B RIAGENDR RIDRETH2 PIRGRP3 
          ASDIST5B*RIAGENDR ASDIST5B*RIDRETH2 ASDIST5B*PIRGRP3;
rformat RIDRETH2 race. ;
rformat riagendr sex. ;
rformat ASDIST5B age. ;
rformat PIRGRP3 pir. ;

print nsum wsum total percent sepercent / style =nchs totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
output nsum wsum total percent sepercent / 
       filename=NHANESAB.N4E1Bu REPLACE filetype=SAS 
       totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
rtitle "E1. N IV phase I + II not age/sex adjusted (1 yr) fluorosis prevalence [Questionable] among ages 6 to 39 (N4E1Bu)";
run;
data NHANESAB.N4E1Bu;
set  NHANESAB.N4E1Bu;
format tablename $8.;
tablename="N4E1Bu   ";
run;
proc descript data=NHANESAB.nh4mmwr design=wr;/*E1. fluorosis prevalence [Questionable]   - age adjusted*/
weight WTMEC4YR;
nest SDMVSTRA SDMVPSU/strlev=1 psulev=2 missunit;
var OHDDFISX;
catlevel 2;
subpopn ridageyr >=6 and ridageyr<=39;
subgroup  RIDRETH2  ASDIST5B  RIAGENDR PIRGRP3 ASDIST5F;
level         5         4          2        3    34;
tables      RIAGENDR 
          ASDIST5B*RIAGENDR;
Stdvar ASDIST5F;    
Stdwgt	0.029554 0.030007 0.029449 0.031862 0.031908 0.030384
	    0.029728 0.029592 0.029675 0.029811 0.029071 0.030263
		0.029698 0.030656 0.030459 0.028400 0.026824 0.026318
		0.025722 0.025805 0.025104 0.025911 0.026032 0.030829
        0.030165 0.021196 0.028468 0.028973 0.030829 0.032345
        0.032813 0.033710 0.032360 0.036079
;
rformat RIDRETH2 race. ;
rformat riagendr sex. ;
rformat ASDIST5B age. ;
rformat PIRGRP3 pir. ;

print nsum wsum total percent sepercent / style =nchs totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
output nsum wsum total percent sepercent / 
       filename=NHANESAB.N4E1Baa REPLACE filetype=SAS 
       totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2  ;
rtitle "E1. N IV phase I + II age adjusted (1 yr) fluorosis prevalence [Questionable] among ages 6 to 39 (N4E1Baa)";
run;
data NHANESAB.N4E1Baa;
set  NHANESAB.N4E1Baa;
format tablename $8.;
tablename="N4E1Baa  ";
run;
proc descript data=NHANESAB.nh4mmwr design=wr;/*E1. fluorosis prevalence [Questionable]   - age/sex adjusted*/
weight WTMEC4YR;
nest SDMVSTRA SDMVPSU/strlev=1 psulev=2 missunit;
var OHDDFISX;
catlevel 2;
subpopn ridageyr >=6 and ridageyr<=39;
subgroup  RIDRETH2  ASDIST5B  RIAGENDR PIRGRP3 ASDIST5D;
level         5         4          2        3    68;
tables    ASDIST5B  RIDRETH2 PIRGRP3 
          ASDIST5B*RIDRETH2 ASDIST5B*PIRGRP3;
Stdvar ASDIST5D;    
Stdwgt 0.015018 0.014307 0.015250 0.014539 0.014943 0.014284 0.016208 0.015415
       0.016246 0.015422 0.015445 0.014711 0.015093 0.014419 0.015025 0.014351
       0.015100 0.014351 0.015168 0.014419 0.014831 0.014022 0.015467 0.014569
       0.015070 0.014404 0.015489 0.014943 0.015400 0.014831 0.014351 0.013835
       0.013513 0.013109 0.013199 0.012929 0.012847 0.012682 0.012862 0.012749
       0.012427 0.012487 0.012832 0.012884 0.012862 0.012974 0.015212 0.015385
       0.014868 0.015070 0.014172 0.014351 0.014037 0.014224 0.014232 0.014524
       0.015265 0.015332 0.015984 0.016118 0.016193 0.016365 0.016642 0.016815
       0.015969 0.016156 0.017878 0.017926 	
; 
rformat RIDRETH2 race. ;
rformat riagendr sex. ;
rformat ASDIST5B age. ;
rformat PIRGRP3 pir. ;

print nsum wsum total percent sepercent / style =nchs totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
output nsum wsum total percent sepercent/ 
       filename=NHANESAB.N4E1Basa REPLACE filetype=SAS 
       totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
rtitle "E1. N IV phase I + II age/sex adjusted (1 yr) fluorosis prevalence [Questionable] among ages 6 to 39 (N4E1Basa)";
run;
data NHANESAB.N4E1Basa;
set  NHANESAB.N4E1Basa;
format tablename $8.;
tablename="N4E1Basa ";
run;

proc descript data=NHANESAB.nh4mmwr design=wr;/*E1. fluorosis prevalence [Very Mild] - not age adjusted*/
weight WTMEC4YR;
nest SDMVSTRA SDMVPSU/strlev=1 psulev=2 missunit;
var OHDDFISX;
catlevel 3;
subpopn ridageyr >=6 and ridageyr<=39;
subgroup  RIDRETH2  ASDIST5B  RIAGENDR PIRGRP3;
level         5         4          2        3 ;
tables    ASDIST5B RIAGENDR RIDRETH2 PIRGRP3 
          ASDIST5B*RIAGENDR ASDIST5B*RIDRETH2 ASDIST5B*PIRGRP3;
rformat RIDRETH2 race. ;
rformat riagendr sex. ;
rformat ASDIST5B age. ;
rformat PIRGRP3 pir. ;

print nsum wsum total percent sepercent / style =nchs totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
output nsum wsum total percent sepercent / 
       filename=NHANESAB.N4E1Cu REPLACE filetype=SAS 
       totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
rtitle "E1. N IV phase I + II not age/sex adjusted (1 yr) fluorosis prevalence [Very Mild] among ages 6 to 39 (N4E1Cu)";
run;
data NHANESAB.N4E1Cu;
set  NHANESAB.N4E1Cu;
format tablename $8.;
tablename="N4E1Cu   ";
run;
proc descript data=NHANESAB.nh4mmwr design=wr;/*E1. fluorosis prevalence [Very Mild] - age adjusted*/
weight WTMEC4YR;
nest SDMVSTRA SDMVPSU/strlev=1 psulev=2 missunit;
var OHDDFISX;
catlevel 3;
subpopn ridageyr >=6 and ridageyr<=39;
subgroup  RIDRETH2  ASDIST5B  RIAGENDR PIRGRP3 ASDIST5F;
level         5         4          2        3    34;
tables      RIAGENDR 
          ASDIST5B*RIAGENDR;
Stdvar ASDIST5F;    
Stdwgt	0.029554 0.030007 0.029449 0.031862 0.031908 0.030384
	    0.029728 0.029592 0.029675 0.029811 0.029071 0.030263
		0.029698 0.030656 0.030459 0.028400 0.026824 0.026318
		0.025722 0.025805 0.025104 0.025911 0.026032 0.030829
        0.030165 0.021196 0.028468 0.028973 0.030829 0.032345
        0.032813 0.033710 0.032360 0.036079
;
rformat RIDRETH2 race. ;
rformat riagendr sex. ;
rformat ASDIST5B age. ;
rformat PIRGRP3 pir. ;

print nsum wsum total percent sepercent / style =nchs totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
output nsum wsum total percent sepercent / 
       filename=NHANESAB.N4E1Caa REPLACE filetype=SAS 
       totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2  ;
rtitle "E1. N IV phase I + II age adjusted (1 yr) fluorosis prevalence [Very Mild] among ages 6 to 39 (N4E1Caa)";
run;
data NHANESAB.N4E1Caa;
set  NHANESAB.N4E1Caa;
format tablename $8.;
tablename="N4E1Caa  ";
run;
proc descript data=NHANESAB.nh4mmwr design=wr;/*E1. fluorosis prevalence [Very Mild] - age/sex adjusted*/
weight WTMEC4YR;
nest SDMVSTRA SDMVPSU/strlev=1 psulev=2 missunit;
var OHDDFISX;
catlevel 3;
subpopn ridageyr >=6 and ridageyr<=39;
subgroup  RIDRETH2  ASDIST5B  RIAGENDR PIRGRP3 ASDIST5D;
level         5         4          2        3    68;
tables    ASDIST5B  RIDRETH2 PIRGRP3 
          ASDIST5B*RIDRETH2 ASDIST5B*PIRGRP3;
Stdvar ASDIST5D;    
Stdwgt 0.015018 0.014307 0.015250 0.014539 0.014943 0.014284 0.016208 0.015415
       0.016246 0.015422 0.015445 0.014711 0.015093 0.014419 0.015025 0.014351
       0.015100 0.014351 0.015168 0.014419 0.014831 0.014022 0.015467 0.014569
       0.015070 0.014404 0.015489 0.014943 0.015400 0.014831 0.014351 0.013835
       0.013513 0.013109 0.013199 0.012929 0.012847 0.012682 0.012862 0.012749
       0.012427 0.012487 0.012832 0.012884 0.012862 0.012974 0.015212 0.015385
       0.014868 0.015070 0.014172 0.014351 0.014037 0.014224 0.014232 0.014524
       0.015265 0.015332 0.015984 0.016118 0.016193 0.016365 0.016642 0.016815
       0.015969 0.016156 0.017878 0.017926 	
; 
rformat RIDRETH2 race. ;
rformat riagendr sex. ;
rformat ASDIST5B age. ;
rformat PIRGRP3 pir. ;

print nsum wsum total percent sepercent / style =nchs totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
output nsum wsum total percent sepercent/ 
       filename=NHANESAB.N4E1Casa REPLACE filetype=SAS 
       totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
rtitle "E1. N IV phase I + II age/sex adjusted (1 yr) fluorosis prevalence [Very Mild] among ages 6 to 39 (N4E1Casa)";
run;
data NHANESAB.N4E1Casa;
set  NHANESAB.N4E1Casa;
format tablename $8.;
tablename="N4E1Casa ";
run;

proc descript data=NHANESAB.nh4mmwr design=wr;/*E1. fluorosis prevalence [Mild] - not age adjusted*/
weight WTMEC4YR;
nest SDMVSTRA SDMVPSU/strlev=1 psulev=2 missunit;
var OHDDFISX;
catlevel 4;
subpopn ridageyr >=6 and ridageyr<=39;
subgroup  RIDRETH2  ASDIST5B  RIAGENDR PIRGRP3;
level         5         4          2        3 ;
tables    ASDIST5B RIAGENDR RIDRETH2 PIRGRP3 
          ASDIST5B*RIAGENDR ASDIST5B*RIDRETH2 ASDIST5B*PIRGRP3;
rformat RIDRETH2 race. ;
rformat riagendr sex. ;
rformat ASDIST5B age. ;
rformat PIRGRP3 pir. ;

print nsum wsum total percent sepercent / style =nchs totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
output nsum wsum total percent sepercent / 
       filename=NHANESAB.N4E1Du REPLACE filetype=SAS 
       totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
rtitle "E1. N IV phase I + II not age/sex adjusted (1 yr) fluorosis prevalence [Mild] among ages 6 to 39 (N4E1Du)";
run;
data NHANESAB.N4E1Du;
set  NHANESAB.N4E1Du;
format tablename $8.;
tablename="N4E1Du   ";
run;
proc descript data=NHANESAB.nh4mmwr design=wr;/*E1. fluorosis prevalence [Mild] - age adjusted*/
weight WTMEC4YR;
nest SDMVSTRA SDMVPSU/strlev=1 psulev=2 missunit;
var OHDDFISX;
catlevel 4;
subpopn ridageyr >=6 and ridageyr<=39;
subgroup  RIDRETH2  ASDIST5B  RIAGENDR PIRGRP3 ASDIST5F;
level         5         4          2        3    34;
tables      RIAGENDR 
          ASDIST5B*RIAGENDR;
Stdvar ASDIST5F;    
Stdwgt	0.029554 0.030007 0.029449 0.031862 0.031908 0.030384
	    0.029728 0.029592 0.029675 0.029811 0.029071 0.030263
		0.029698 0.030656 0.030459 0.028400 0.026824 0.026318
		0.025722 0.025805 0.025104 0.025911 0.026032 0.030829
        0.030165 0.021196 0.028468 0.028973 0.030829 0.032345
        0.032813 0.033710 0.032360 0.036079
;
rformat RIDRETH2 race. ;
rformat riagendr sex. ;
rformat ASDIST5B age. ;
rformat PIRGRP3 pir. ;

print nsum wsum total percent sepercent / style =nchs totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
output nsum wsum total percent sepercent / 
       filename=NHANESAB.N4E1Daa REPLACE filetype=SAS 
       totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2  ;
rtitle "E1. N IV phase I + II age adjusted (1 yr) fluorosis prevalence [Mild] among ages 6 to 39 (N4E1Daa)";
run;
data NHANESAB.N4E1Daa;
set  NHANESAB.N4E1Daa;
format tablename $8.;
tablename="N4E1Daa  ";
run;
proc descript data=NHANESAB.nh4mmwr design=wr;/*E1. fluorosis prevalence [Mild] - age/sex adjusted*/
weight WTMEC4YR;
nest SDMVSTRA SDMVPSU/strlev=1 psulev=2 missunit;
var OHDDFISX;
catlevel 4;
subpopn ridageyr >=6 and ridageyr<=39;
subgroup  RIDRETH2  ASDIST5B  RIAGENDR PIRGRP3 ASDIST5D;
level         5         4          2        3    68;
tables    ASDIST5B  RIDRETH2 PIRGRP3 
          ASDIST5B*RIDRETH2 ASDIST5B*PIRGRP3;
Stdvar ASDIST5D;    
Stdwgt 0.015018 0.014307 0.015250 0.014539 0.014943 0.014284 0.016208 0.015415
       0.016246 0.015422 0.015445 0.014711 0.015093 0.014419 0.015025 0.014351
       0.015100 0.014351 0.015168 0.014419 0.014831 0.014022 0.015467 0.014569
       0.015070 0.014404 0.015489 0.014943 0.015400 0.014831 0.014351 0.013835
       0.013513 0.013109 0.013199 0.012929 0.012847 0.012682 0.012862 0.012749
       0.012427 0.012487 0.012832 0.012884 0.012862 0.012974 0.015212 0.015385
       0.014868 0.015070 0.014172 0.014351 0.014037 0.014224 0.014232 0.014524
       0.015265 0.015332 0.015984 0.016118 0.016193 0.016365 0.016642 0.016815
       0.015969 0.016156 0.017878 0.017926 	
; 
rformat RIDRETH2 race. ;
rformat riagendr sex. ;
rformat ASDIST5B age. ;
rformat PIRGRP3 pir. ;

print nsum wsum total percent sepercent / style =nchs totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
output nsum wsum total percent sepercent/ 
       filename=NHANESAB.N4E1Dasa REPLACE filetype=SAS 
       totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
rtitle "E1. N IV phase I + II age/sex adjusted (1 yr) fluorosis prevalence [Mild] among ages 6 to 39 (N4E1Dasa)";
run;
data NHANESAB.N4E1Dasa;
set  NHANESAB.N4E1Dasa;
format tablename $8.;
tablename="N4E1Dasa ";
run;


proc descript data=NHANESAB.nh4mmwr design=wr;/*E1. fluorosis prevalence [Moderate] - not age adjusted*/
weight WTMEC4YR;
nest SDMVSTRA SDMVPSU/strlev=1 psulev=2 missunit;
var OHDDFISX;
catlevel 5;
subpopn ridageyr >=6 and ridageyr<=39;
subgroup  RIDRETH2  ASDIST5B  RIAGENDR PIRGRP3;
level         5         4          2        3 ;
tables    ASDIST5B RIAGENDR RIDRETH2 PIRGRP3 
          ASDIST5B*RIAGENDR ASDIST5B*RIDRETH2 ASDIST5B*PIRGRP3;
rformat RIDRETH2 race. ;
rformat riagendr sex. ;
rformat ASDIST5B age. ;
rformat PIRGRP3 pir. ;

print nsum wsum total percent sepercent / style =nchs totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
output nsum wsum total percent sepercent / 
       filename=NHANESAB.N4E1Eu REPLACE filetype=SAS 
       totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
rtitle "E1. N IV phase I + II not age/sex adjusted (1 yr) fluorosis prevalence [Moderate] among ages 6 to 39 (N4E1Eu)";
run;
data NHANESAB.N4E1Eu;
set  NHANESAB.N4E1Eu;
format tablename $8.;
tablename="N4E1Eu   ";
run;
proc descript data=NHANESAB.nh4mmwr design=wr;/*E1. fluorosis prevalence [Moderate] - age adjusted*/
weight WTMEC4YR;
nest SDMVSTRA SDMVPSU/strlev=1 psulev=2 missunit;
var OHDDFISX;
catlevel 5;
subpopn ridageyr >=6 and ridageyr<=39;
subgroup  RIDRETH2  ASDIST5B  RIAGENDR PIRGRP3 ASDIST5F;
level         5         4          2        3    34;
tables      RIAGENDR 
          ASDIST5B*RIAGENDR;
Stdvar ASDIST5F;    
Stdwgt	0.029554 0.030007 0.029449 0.031862 0.031908 0.030384
	    0.029728 0.029592 0.029675 0.029811 0.029071 0.030263
		0.029698 0.030656 0.030459 0.028400 0.026824 0.026318
		0.025722 0.025805 0.025104 0.025911 0.026032 0.030829
        0.030165 0.021196 0.028468 0.028973 0.030829 0.032345
        0.032813 0.033710 0.032360 0.036079
;
rformat RIDRETH2 race. ;
rformat riagendr sex. ;
rformat ASDIST5B age. ;
rformat PIRGRP3 pir. ;

print nsum wsum total percent sepercent / style =nchs totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
output nsum wsum total percent sepercent / 
       filename=NHANESAB.N4E1Eaa REPLACE filetype=SAS 
       totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2  ;
rtitle "E1. N IV phase I + II age adjusted (1 yr) fluorosis prevalence [Moderate] among ages 6 to 39 (N4E1Eaa)";
run;
data NHANESAB.N4E1Eaa;
set  NHANESAB.N4E1Eaa;
format tablename $8.;
tablename="N4E1Eaa  ";
run;
proc descript data=NHANESAB.nh4mmwr design=wr;/*E1. fluorosis prevalence [Moderate] - age/sex adjusted*/
weight WTMEC4YR;
nest SDMVSTRA SDMVPSU/strlev=1 psulev=2 missunit;
var OHDDFISX;
catlevel 5;
subpopn ridageyr >=6 and ridageyr<=39;
subgroup  RIDRETH2  ASDIST5B  RIAGENDR PIRGRP3 ASDIST5D;
level         5         4          2        3    68;
tables    ASDIST5B  RIDRETH2 PIRGRP3 
          ASDIST5B*RIDRETH2 ASDIST5B*PIRGRP3;
Stdvar ASDIST5D;    
Stdwgt 0.015018 0.014307 0.015250 0.014539 0.014943 0.014284 0.016208 0.015415
       0.016246 0.015422 0.015445 0.014711 0.015093 0.014419 0.015025 0.014351
       0.015100 0.014351 0.015168 0.014419 0.014831 0.014022 0.015467 0.014569
       0.015070 0.014404 0.015489 0.014943 0.015400 0.014831 0.014351 0.013835
       0.013513 0.013109 0.013199 0.012929 0.012847 0.012682 0.012862 0.012749
       0.012427 0.012487 0.012832 0.012884 0.012862 0.012974 0.015212 0.015385
       0.014868 0.015070 0.014172 0.014351 0.014037 0.014224 0.014232 0.014524
       0.015265 0.015332 0.015984 0.016118 0.016193 0.016365 0.016642 0.016815
       0.015969 0.016156 0.017878 0.017926 	
; 
rformat RIDRETH2 race. ;
rformat riagendr sex. ;
rformat ASDIST5B age. ;
rformat PIRGRP3 pir. ;

print nsum wsum total percent sepercent / style =nchs totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
output nsum wsum total percent sepercent/ 
       filename=NHANESAB.N4E1Easa REPLACE filetype=SAS 
       totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
rtitle "E1. N IV phase I + II age/sex adjusted (1 yr) fluorosis prevalence [Moderate] among ages 6 to 39 (N4E1Easa)";
run;
data NHANESAB.N4E1Easa;
set  NHANESAB.N4E1Easa;
format tablename $8.;
tablename="N4E1Easa ";
run;


proc descript data=NHANESAB.nh4mmwr design=wr;/*E1. fluorosis prevalence [Severe] - not age adjusted*/
weight WTMEC4YR;
nest SDMVSTRA SDMVPSU/strlev=1 psulev=2 missunit;
var OHDDFISX;
catlevel 6;
subpopn ridageyr >=6 and ridageyr<=39;
subgroup  RIDRETH2  ASDIST5B  RIAGENDR PIRGRP3;
level         5         4          2        3 ;
tables    ASDIST5B RIAGENDR RIDRETH2 PIRGRP3 
          ASDIST5B*RIAGENDR ASDIST5B*RIDRETH2 ASDIST5B*PIRGRP3;
rformat RIDRETH2 race. ;
rformat riagendr sex. ;
rformat ASDIST5B age. ;
rformat PIRGRP3 pir. ;

print nsum wsum total percent sepercent / style =nchs totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
output nsum wsum total percent sepercent / 
       filename=NHANESAB.N4E1Fu REPLACE filetype=SAS 
       totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
rtitle "E1. N IV phase I + II not age/sex adjusted (1 yr) fluorosis prevalence [Severe] among ages 6 to 39 (N4E1Fu)";
run;
data NHANESAB.N4E1Fu;
set  NHANESAB.N4E1Fu;
format tablename $8.;
tablename="N4E1Fu   ";
run;
proc descript data=NHANESAB.nh4mmwr design=wr;/*E1. fluorosis prevalence [Severe] - age adjusted*/
weight WTMEC4YR;
nest SDMVSTRA SDMVPSU/strlev=1 psulev=2 missunit;
var OHDDFISX;
catlevel 6;
subpopn ridageyr >=6 and ridageyr<=39;
subgroup  RIDRETH2  ASDIST5B  RIAGENDR PIRGRP3 ASDIST5F;
level         5         4          2        3    34;
tables      RIAGENDR 
          ASDIST5B*RIAGENDR;
Stdvar ASDIST5F;    
Stdwgt	0.029554 0.030007 0.029449 0.031862 0.031908 0.030384
	    0.029728 0.029592 0.029675 0.029811 0.029071 0.030263
		0.029698 0.030656 0.030459 0.028400 0.026824 0.026318
		0.025722 0.025805 0.025104 0.025911 0.026032 0.030829
        0.030165 0.021196 0.028468 0.028973 0.030829 0.032345
        0.032813 0.033710 0.032360 0.036079
;
rformat RIDRETH2 race. ;
rformat riagendr sex. ;
rformat ASDIST5B age. ;
rformat PIRGRP3 pir. ;

print nsum wsum total percent sepercent / style =nchs totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
output nsum wsum total percent sepercent / 
       filename=NHANESAB.N4E1Faa REPLACE filetype=SAS 
       totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2  ;
rtitle "E1. N IV phase I + II age adjusted (1 yr) fluorosis prevalence [Severe] among ages 6 to 39 (N4E1Faa)";
run;
data NHANESAB.N4E1Faa;
set  NHANESAB.N4E1Faa;
format tablename $8.;
tablename="N4E1Faa  ";
run;
proc descript data=NHANESAB.nh4mmwr design=wr;/*E1. fluorosis prevalence [Severe] - age/sex adjusted*/
weight WTMEC4YR;
nest SDMVSTRA SDMVPSU/strlev=1 psulev=2 missunit;
var OHDDFISX;
catlevel 6;
subpopn ridageyr >=6 and ridageyr<=39;
subgroup  RIDRETH2  ASDIST5B  RIAGENDR PIRGRP3 ASDIST5D;
level         5         4          2        3    68;
tables    ASDIST5B  RIDRETH2 PIRGRP3 
          ASDIST5B*RIDRETH2 ASDIST5B*PIRGRP3;
Stdvar ASDIST5D;    
Stdwgt 0.015018 0.014307 0.015250 0.014539 0.014943 0.014284 0.016208 0.015415
       0.016246 0.015422 0.015445 0.014711 0.015093 0.014419 0.015025 0.014351
       0.015100 0.014351 0.015168 0.014419 0.014831 0.014022 0.015467 0.014569
       0.015070 0.014404 0.015489 0.014943 0.015400 0.014831 0.014351 0.013835
       0.013513 0.013109 0.013199 0.012929 0.012847 0.012682 0.012862 0.012749
       0.012427 0.012487 0.012832 0.012884 0.012862 0.012974 0.015212 0.015385
       0.014868 0.015070 0.014172 0.014351 0.014037 0.014224 0.014232 0.014524
       0.015265 0.015332 0.015984 0.016118 0.016193 0.016365 0.016642 0.016815
       0.015969 0.016156 0.017878 0.017926 	
; 
rformat RIDRETH2 race. ;
rformat riagendr sex. ;
rformat ASDIST5B age. ;
rformat PIRGRP3 pir. ;

print nsum wsum total percent sepercent / style =nchs totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
output nsum wsum total percent sepercent/ 
       filename=NHANESAB.N4E1Fasa REPLACE filetype=SAS 
       totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
rtitle "E1. N IV phase I + II age/sex adjusted (1 yr) fluorosis prevalence [Severe] among ages 6 to 39 (N4E1Fasa)";
run;
data NHANESAB.N4E1Fasa;
set  NHANESAB.N4E1Fasa;
format tablename $8.;
tablename="N4E1Fasa ";
run;

/* Table E1b. Prevalence and severity of enamel fluorosis among persons aged 6 to 39 years,
             by selected characteristics -- United States 1999-2002 
             one column for each category - moderate and severe combined
             plus weighed dean's index
             no, questionable, very mild, mild, moderate/severe Weighted Dean Index*/

proc descript data=NHANESAB.nh4mmwr design=wr;/*E1b. fluorosis prevalence [No]  - not age adjusted*/
weight WTMEC4YR;
nest SDMVSTRA SDMVPSU/strlev=1 psulev=2 missunit;
var OHDDFIaSX;
catlevel 1;
subpopn ridageyr >=6 and ridageyr<=39;
subgroup  RIDRETH2  ASDIST5B  RIAGENDR PIRGRP3;
level         5         4          2        3 ;
tables    ASDIST5B RIAGENDR RIDRETH2 PIRGRP3 
          ASDIST5B*RIAGENDR ASDIST5B*RIDRETH2 ASDIST5B*PIRGRP3;
rformat RIDRETH2 race. ;
rformat riagendr sex. ;
rformat ASDIST5B age. ;
rformat PIRGRP3 pir. ;

print nsum wsum total percent sepercent / style =nchs totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
output nsum wsum total percent sepercent / 
       filename=NHANESAB.N4E2Au REPLACE filetype=SAS 
       totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
rtitle "E1b. N IV phase I + II not age/sex adjusted (1 yr) fluorosis prevalence [No] among ages 6 to 39 (N4E2Au)";
run;
data NHANESAB.N4E2Au;
set  NHANESAB.N4E2Au;
format tablename $8.;
tablename="N4E2Au   ";
run;
proc descript data=NHANESAB.nh4mmwr design=wr;/*E1b. fluorosis prevalence [No]   - age adjusted*/
weight WTMEC4YR;
nest SDMVSTRA SDMVPSU/strlev=1 psulev=2 missunit;
var OHDDFIaSX;
catlevel 1;
subpopn ridageyr >=6 and ridageyr<=39;
subgroup  RIDRETH2  ASDIST5B  RIAGENDR PIRGRP3 ASDIST5F;
level         5         4          2        3    34;
tables      RIAGENDR 
          ASDIST5B*RIAGENDR;
Stdvar ASDIST5F;    
Stdwgt	0.029554 0.030007 0.029449 0.031862 0.031908 0.030384
	    0.029728 0.029592 0.029675 0.029811 0.029071 0.030263
		0.029698 0.030656 0.030459 0.028400 0.026824 0.026318
		0.025722 0.025805 0.025104 0.025911 0.026032 0.030829
        0.030165 0.021196 0.028468 0.028973 0.030829 0.032345
        0.032813 0.033710 0.032360 0.036079
;
rformat RIDRETH2 race. ;
rformat riagendr sex. ;
rformat ASDIST5B age. ;
rformat PIRGRP3 pir. ;

print nsum wsum total percent sepercent / style =nchs totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
output nsum wsum total percent sepercent / 
       filename=NHANESAB.N4E2Aaa REPLACE filetype=SAS 
       totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2  ;
rtitle "E1b. N IV phase I + II age adjusted (1 yr) fluorosis prevalence [No] among ages 6 to 39 (N4E2Aaa)";
run;
data NHANESAB.N4E2Aaa;
set  NHANESAB.N4E2Aaa;
format tablename $8.;
tablename="N4E2Aaa  ";
run;
proc descript data=NHANESAB.nh4mmwr design=wr;/*E1b. fluorosis prevalence [No]   - age/sex adjusted*/
weight WTMEC4YR;
nest SDMVSTRA SDMVPSU/strlev=1 psulev=2 missunit;
var OHDDFIaSX;
catlevel 1;
subpopn ridageyr >=6 and ridageyr<=39;
subgroup  RIDRETH2  ASDIST5B  RIAGENDR PIRGRP3 ASDIST5D;
level         5         4          2        3    68;
tables    ASDIST5B  RIDRETH2 PIRGRP3 
          ASDIST5B*RIDRETH2 ASDIST5B*PIRGRP3;
Stdvar ASDIST5D;    
Stdwgt 0.015018 0.014307 0.015250 0.014539 0.014943 0.014284 0.016208 0.015415
       0.016246 0.015422 0.015445 0.014711 0.015093 0.014419 0.015025 0.014351
       0.015100 0.014351 0.015168 0.014419 0.014831 0.014022 0.015467 0.014569
       0.015070 0.014404 0.015489 0.014943 0.015400 0.014831 0.014351 0.013835
       0.013513 0.013109 0.013199 0.012929 0.012847 0.012682 0.012862 0.012749
       0.012427 0.012487 0.012832 0.012884 0.012862 0.012974 0.015212 0.015385
       0.014868 0.015070 0.014172 0.014351 0.014037 0.014224 0.014232 0.014524
       0.015265 0.015332 0.015984 0.016118 0.016193 0.016365 0.016642 0.016815
       0.015969 0.016156 0.017878 0.017926 	
; 
rformat RIDRETH2 race. ;
rformat riagendr sex. ;
rformat ASDIST5B age. ;
rformat PIRGRP3 pir. ;

print nsum wsum total percent sepercent / style =nchs totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
output nsum wsum total percent sepercent/ 
       filename=NHANESAB.N4E2Aasa REPLACE filetype=SAS 
       totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
rtitle "E1b. N IV phase I + II age/sex adjusted (1 yr) fluorosis prevalence [No] among ages 6 to 39 (N4E2Aasa)";
run;
data NHANESAB.N4E2Aasa;
set  NHANESAB.N4E2Aasa;
format tablename $8.;
tablename="N4E2Aasa ";
run;

proc descript data=NHANESAB.nh4mmwr design=wr;/*E1b. fluorosis prevalence [Questionable]  - not age adjusted*/
weight WTMEC4YR;
nest SDMVSTRA SDMVPSU/strlev=1 psulev=2 missunit;
var OHDDFIaSX;
catlevel 2;
subpopn ridageyr >=6 and ridageyr<=39;
subgroup  RIDRETH2  ASDIST5B  RIAGENDR PIRGRP3;
level         5         4          2        3 ;
tables    ASDIST5B RIAGENDR RIDRETH2 PIRGRP3 
          ASDIST5B*RIAGENDR ASDIST5B*RIDRETH2 ASDIST5B*PIRGRP3;
rformat RIDRETH2 race. ;
rformat riagendr sex. ;
rformat ASDIST5B age. ;
rformat PIRGRP3 pir. ;

print nsum wsum total percent sepercent / style =nchs totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
output nsum wsum total percent sepercent / 
       filename=NHANESAB.N4E2Bu REPLACE filetype=SAS 
       totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
rtitle "E1b. N IV phase I + II not age/sex adjusted (1 yr) fluorosis prevalence [Questionable] among ages 6 to 39 (N4E2Bu)";
run;
data NHANESAB.N4E2Bu;
set  NHANESAB.N4E2Bu;
format tablename $8.;
tablename="N4E2Bu   ";
run;
proc descript data=NHANESAB.nh4mmwr design=wr;/*E1b. fluorosis prevalence [Questionable]   - age adjusted*/
weight WTMEC4YR;
nest SDMVSTRA SDMVPSU/strlev=1 psulev=2 missunit;
var OHDDFIaSX;
catlevel 2;
subpopn ridageyr >=6 and ridageyr<=39;
subgroup  RIDRETH2  ASDIST5B  RIAGENDR PIRGRP3 ASDIST5F;
level         5         4          2        3    34;
tables      RIAGENDR 
          ASDIST5B*RIAGENDR;
Stdvar ASDIST5F;    
Stdwgt	0.029554 0.030007 0.029449 0.031862 0.031908 0.030384
	    0.029728 0.029592 0.029675 0.029811 0.029071 0.030263
		0.029698 0.030656 0.030459 0.028400 0.026824 0.026318
		0.025722 0.025805 0.025104 0.025911 0.026032 0.030829
        0.030165 0.021196 0.028468 0.028973 0.030829 0.032345
        0.032813 0.033710 0.032360 0.036079
;
rformat RIDRETH2 race. ;
rformat riagendr sex. ;
rformat ASDIST5B age. ;
rformat PIRGRP3 pir. ;

print nsum wsum total percent sepercent / style =nchs totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
output nsum wsum total percent sepercent / 
       filename=NHANESAB.N4E2Baa REPLACE filetype=SAS 
       totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2  ;
rtitle "E1b. N IV phase I + II age adjusted (1 yr) fluorosis prevalence [Questionable] among ages 6 to 39 (N4E2Baa)";
run;
data NHANESAB.N4E2Baa;
set  NHANESAB.N4E2Baa;
format tablename $8.;
tablename="N4E2Baa  ";
run;
proc descript data=NHANESAB.nh4mmwr design=wr;/*E1b. fluorosis prevalence [Questionable]   - age/sex adjusted*/
weight WTMEC4YR;
nest SDMVSTRA SDMVPSU/strlev=1 psulev=2 missunit;
var OHDDFIaSX;
catlevel 2;
subpopn ridageyr >=6 and ridageyr<=39;
subgroup  RIDRETH2  ASDIST5B  RIAGENDR PIRGRP3 ASDIST5D;
level         5         4          2        3    68;
tables    ASDIST5B  RIDRETH2 PIRGRP3 
          ASDIST5B*RIDRETH2 ASDIST5B*PIRGRP3;
Stdvar ASDIST5D;    
Stdwgt 0.015018 0.014307 0.015250 0.014539 0.014943 0.014284 0.016208 0.015415
       0.016246 0.015422 0.015445 0.014711 0.015093 0.014419 0.015025 0.014351
       0.015100 0.014351 0.015168 0.014419 0.014831 0.014022 0.015467 0.014569
       0.015070 0.014404 0.015489 0.014943 0.015400 0.014831 0.014351 0.013835
       0.013513 0.013109 0.013199 0.012929 0.012847 0.012682 0.012862 0.012749
       0.012427 0.012487 0.012832 0.012884 0.012862 0.012974 0.015212 0.015385
       0.014868 0.015070 0.014172 0.014351 0.014037 0.014224 0.014232 0.014524
       0.015265 0.015332 0.015984 0.016118 0.016193 0.016365 0.016642 0.016815
       0.015969 0.016156 0.017878 0.017926 	
; 
rformat RIDRETH2 race. ;
rformat riagendr sex. ;
rformat ASDIST5B age. ;
rformat PIRGRP3 pir. ;

print nsum wsum total percent sepercent / style =nchs totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
output nsum wsum total percent sepercent/ 
       filename=NHANESAB.N4E2Basa REPLACE filetype=SAS 
       totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
rtitle "E1b. N IV phase I + II age/sex adjusted (1 yr) fluorosis prevalence [Questionable] among ages 6 to 39 (N4E2Basa)";
run;
data NHANESAB.N4E2Basa;
set  NHANESAB.N4E2Basa;
format tablename $8.;
tablename="N4E2Basa ";
run;

proc descript data=NHANESAB.nh4mmwr design=wr;/*E1b. fluorosis prevalence [Very Mild] - not age adjusted*/
weight WTMEC4YR;
nest SDMVSTRA SDMVPSU/strlev=1 psulev=2 missunit;
var OHDDFIaSX;
catlevel 3;
subpopn ridageyr >=6 and ridageyr<=39;
subgroup  RIDRETH2  ASDIST5B  RIAGENDR PIRGRP3;
level         5         4          2        3 ;
tables    ASDIST5B RIAGENDR RIDRETH2 PIRGRP3 
          ASDIST5B*RIAGENDR ASDIST5B*RIDRETH2 ASDIST5B*PIRGRP3;
rformat RIDRETH2 race. ;
rformat riagendr sex. ;
rformat ASDIST5B age. ;
rformat PIRGRP3 pir. ;

print nsum wsum total percent sepercent / style =nchs totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
output nsum wsum total percent sepercent / 
       filename=NHANESAB.N4E2Cu REPLACE filetype=SAS 
       totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
rtitle "E1b. N IV phase I + II not age/sex adjusted (1 yr) fluorosis prevalence [Very Mild] among ages 6 to 39 (N4E2Cu)";
run;
data NHANESAB.N4E2Cu;
set  NHANESAB.N4E2Cu;
format tablename $8.;
tablename="N4E2Cu   ";
run;
proc descript data=NHANESAB.nh4mmwr design=wr;/*E1b. fluorosis prevalence [Very Mild] - age adjusted*/
weight WTMEC4YR;
nest SDMVSTRA SDMVPSU/strlev=1 psulev=2 missunit;
var OHDDFIaSX;
catlevel 3;
subpopn ridageyr >=6 and ridageyr<=39;
subgroup  RIDRETH2  ASDIST5B  RIAGENDR PIRGRP3 ASDIST5F;
level         5         4          2        3    34;
tables      RIAGENDR 
          ASDIST5B*RIAGENDR;
Stdvar ASDIST5F;    
Stdwgt	0.029554 0.030007 0.029449 0.031862 0.031908 0.030384
	    0.029728 0.029592 0.029675 0.029811 0.029071 0.030263
		0.029698 0.030656 0.030459 0.028400 0.026824 0.026318
		0.025722 0.025805 0.025104 0.025911 0.026032 0.030829
        0.030165 0.021196 0.028468 0.028973 0.030829 0.032345
        0.032813 0.033710 0.032360 0.036079
;
rformat RIDRETH2 race. ;
rformat riagendr sex. ;
rformat ASDIST5B age. ;
rformat PIRGRP3 pir. ;

print nsum wsum total percent sepercent / style =nchs totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
output nsum wsum total percent sepercent / 
       filename=NHANESAB.N4E2Caa REPLACE filetype=SAS 
       totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2  ;
rtitle "E1b. N IV phase I + II age adjusted (1 yr) fluorosis prevalence [Very Mild] among ages 6 to 39 (N4E2Caa)";
run;
data NHANESAB.N4E2Caa;
set  NHANESAB.N4E2Caa;
format tablename $8.;
tablename="N4E2Caa  ";
run;
proc descript data=NHANESAB.nh4mmwr design=wr;/*E1b. fluorosis prevalence [Very Mild]  - age/sex adjusted*/
weight WTMEC4YR;
nest SDMVSTRA SDMVPSU/strlev=1 psulev=2 missunit;
var OHDDFIaSX;
catlevel 3;
subpopn ridageyr >=6 and ridageyr<=39;
subgroup  RIDRETH2  ASDIST5B  RIAGENDR PIRGRP3 ASDIST5D;
level         5         4          2        3    68;
tables    ASDIST5B  RIDRETH2 PIRGRP3 
          ASDIST5B*RIDRETH2 ASDIST5B*PIRGRP3;
Stdvar ASDIST5D;    
Stdwgt 0.015018 0.014307 0.015250 0.014539 0.014943 0.014284 0.016208 0.015415
       0.016246 0.015422 0.015445 0.014711 0.015093 0.014419 0.015025 0.014351
       0.015100 0.014351 0.015168 0.014419 0.014831 0.014022 0.015467 0.014569
       0.015070 0.014404 0.015489 0.014943 0.015400 0.014831 0.014351 0.013835
       0.013513 0.013109 0.013199 0.012929 0.012847 0.012682 0.012862 0.012749
       0.012427 0.012487 0.012832 0.012884 0.012862 0.012974 0.015212 0.015385
       0.014868 0.015070 0.014172 0.014351 0.014037 0.014224 0.014232 0.014524
       0.015265 0.015332 0.015984 0.016118 0.016193 0.016365 0.016642 0.016815
       0.015969 0.016156 0.017878 0.017926 	
; 
rformat RIDRETH2 race. ;
rformat riagendr sex. ;
rformat ASDIST5B age. ;
rformat PIRGRP3 pir. ;

print nsum wsum total percent sepercent / style =nchs totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
output nsum wsum total percent sepercent/ 
       filename=NHANESAB.N4E2Casa REPLACE filetype=SAS 
       totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
rtitle "E1b. N IV phase I + II age/sex adjusted (1 yr) fluorosis prevalence [Very Mild] among ages 6 to 39 (N4E2Casa)";
run;
data NHANESAB.N4E2Casa;
set  NHANESAB.N4E2Casa;
format tablename $8.;
tablename="N4E2Casa ";
run;

proc descript data=NHANESAB.nh4mmwr design=wr;/*E1b. fluorosis prevalence [Mild] - not age adjusted*/
weight WTMEC4YR;
nest SDMVSTRA SDMVPSU/strlev=1 psulev=2 missunit;
var OHDDFIaSX;
catlevel 4;
subpopn ridageyr >=6 and ridageyr<=39;
subgroup  RIDRETH2  ASDIST5B  RIAGENDR PIRGRP3;
level         5         4          2        3 ;
tables    ASDIST5B RIAGENDR RIDRETH2 PIRGRP3 
          ASDIST5B*RIAGENDR ASDIST5B*RIDRETH2 ASDIST5B*PIRGRP3;
rformat RIDRETH2 race. ;
rformat riagendr sex. ;
rformat ASDIST5B age. ;
rformat PIRGRP3 pir. ;

print nsum wsum total percent sepercent / style =nchs totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
output nsum wsum total percent sepercent / 
       filename=NHANESAB.N4E2Du REPLACE filetype=SAS 
       totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
rtitle "E1b. N IV phase I + II not age/sex adjusted (1 yr) fluorosis prevalence [Mild] among ages 6 to 39 (N4E2Du)";
run;
data NHANESAB.N4E2Du;
set  NHANESAB.N4E2Du;
format tablename $8.;
tablename="N4E2Du   ";
run;
proc descript data=NHANESAB.nh4mmwr design=wr;/*E1b. fluorosis prevalence [Mild] - age adjusted*/
weight WTMEC4YR;
nest SDMVSTRA SDMVPSU/strlev=1 psulev=2 missunit;
var OHDDFIaSX;
catlevel 4;
subpopn ridageyr >=6 and ridageyr<=39;
subgroup  RIDRETH2  ASDIST5B  RIAGENDR PIRGRP3 ASDIST5F;
level         5         4          2        3    34;
tables      RIAGENDR 
          ASDIST5B*RIAGENDR;
Stdvar ASDIST5F;    
Stdwgt	0.029554 0.030007 0.029449 0.031862 0.031908 0.030384
	    0.029728 0.029592 0.029675 0.029811 0.029071 0.030263
		0.029698 0.030656 0.030459 0.028400 0.026824 0.026318
		0.025722 0.025805 0.025104 0.025911 0.026032 0.030829
        0.030165 0.021196 0.028468 0.028973 0.030829 0.032345
        0.032813 0.033710 0.032360 0.036079
;
rformat RIDRETH2 race. ;
rformat riagendr sex. ;
rformat ASDIST5B age. ;
rformat PIRGRP3 pir. ;

print nsum wsum total percent sepercent / style =nchs totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
output nsum wsum total percent sepercent / 
       filename=NHANESAB.N4E2Daa REPLACE filetype=SAS 
       totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2  ;
rtitle "E1b. N IV phase I + II age adjusted (1 yr) fluorosis prevalence [Mild] among ages 6 to 39 (N4E2Daa)";
run;
data NHANESAB.N4E2Daa;
set  NHANESAB.N4E2Daa;
format tablename $8.;
tablename="N4E2Daa  ";
run;
proc descript data=NHANESAB.nh4mmwr design=wr;/*E1b. fluorosis prevalence [Mild] - age/sex adjusted*/
weight WTMEC4YR;
nest SDMVSTRA SDMVPSU/strlev=1 psulev=2 missunit;
var OHDDFIaSX;
catlevel 4;
subpopn ridageyr >=6 and ridageyr<=39;
subgroup  RIDRETH2  ASDIST5B  RIAGENDR PIRGRP3 ASDIST5D;
level         5         4          2        3    68;
tables    ASDIST5B  RIDRETH2 PIRGRP3 
          ASDIST5B*RIDRETH2 ASDIST5B*PIRGRP3;
Stdvar ASDIST5D;    
Stdwgt 0.015018 0.014307 0.015250 0.014539 0.014943 0.014284 0.016208 0.015415
       0.016246 0.015422 0.015445 0.014711 0.015093 0.014419 0.015025 0.014351
       0.015100 0.014351 0.015168 0.014419 0.014831 0.014022 0.015467 0.014569
       0.015070 0.014404 0.015489 0.014943 0.015400 0.014831 0.014351 0.013835
       0.013513 0.013109 0.013199 0.012929 0.012847 0.012682 0.012862 0.012749
       0.012427 0.012487 0.012832 0.012884 0.012862 0.012974 0.015212 0.015385
       0.014868 0.015070 0.014172 0.014351 0.014037 0.014224 0.014232 0.014524
       0.015265 0.015332 0.015984 0.016118 0.016193 0.016365 0.016642 0.016815
       0.015969 0.016156 0.017878 0.017926 	
; 
rformat RIDRETH2 race. ;
rformat riagendr sex. ;
rformat ASDIST5B age. ;
rformat PIRGRP3 pir. ;

print nsum wsum total percent sepercent / style =nchs totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
output nsum wsum total percent sepercent/ 
       filename=NHANESAB.N4E2Dasa REPLACE filetype=SAS 
       totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
rtitle "E1b. N IV phase I + II age/sex adjusted (1 yr) fluorosis prevalence [Mild] among ages 6 to 39 (N4E2Dasa)";
run;
data NHANESAB.N4E2Dasa;
set  NHANESAB.N4E2Dasa;
format tablename $8.;
tablename="N4E2Dasa ";
run;


proc descript data=NHANESAB.nh4mmwr design=wr;/*E1b. fluorosis prevalence [Mod/Severe] - not age adjusted*/
weight WTMEC4YR;
nest SDMVSTRA SDMVPSU/strlev=1 psulev=2 missunit;
var OHDDFIaSX;
catlevel 5;
subpopn ridageyr >=6 and ridageyr<=39;
subgroup  RIDRETH2  ASDIST5B  RIAGENDR PIRGRP3;
level         5         4          2        3 ;
tables    ASDIST5B RIAGENDR RIDRETH2 PIRGRP3 
          ASDIST5B*RIAGENDR ASDIST5B*RIDRETH2 ASDIST5B*PIRGRP3;
rformat RIDRETH2 race. ;
rformat riagendr sex. ;
rformat ASDIST5B age. ;
rformat PIRGRP3 pir. ;

print nsum wsum total percent sepercent / style =nchs totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
output nsum wsum total percent sepercent / 
       filename=NHANESAB.N4E2Eu REPLACE filetype=SAS 
       totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
rtitle "E1b. N IV phase I + II not age/sex adjusted (1 yr) fluorosis prevalence [Mod/Severe] among ages 6 to 39 (N4E2Eu)";
run;
data NHANESAB.N4E2Eu;
set  NHANESAB.N4E2Eu;
format tablename $8.;
tablename="N4E2Eu   ";
run;
proc descript data=NHANESAB.nh4mmwr design=wr;/*E1b. fluorosis prevalence [Mod/Severe] - age adjusted*/
weight WTMEC4YR;
nest SDMVSTRA SDMVPSU/strlev=1 psulev=2 missunit;
var OHDDFIaSX;
catlevel 5;
subpopn ridageyr >=6 and ridageyr<=39;
subgroup  RIDRETH2  ASDIST5B  RIAGENDR PIRGRP3 ASDIST5F;
level         5         4          2        3    34;
tables      RIAGENDR 
          ASDIST5B*RIAGENDR;
Stdvar ASDIST5F;    
Stdwgt	0.029554 0.030007 0.029449 0.031862 0.031908 0.030384
	    0.029728 0.029592 0.029675 0.029811 0.029071 0.030263
		0.029698 0.030656 0.030459 0.028400 0.026824 0.026318
		0.025722 0.025805 0.025104 0.025911 0.026032 0.030829
        0.030165 0.021196 0.028468 0.028973 0.030829 0.032345
        0.032813 0.033710 0.032360 0.036079
;
rformat RIDRETH2 race. ;
rformat riagendr sex. ;
rformat ASDIST5B age. ;
rformat PIRGRP3 pir. ;

print nsum wsum total percent sepercent / style =nchs totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
output nsum wsum total percent sepercent / 
       filename=NHANESAB.N4E2Eaa REPLACE filetype=SAS 
       totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2  ;
rtitle "E1b. N IV phase I + II age adjusted (1 yr) fluorosis prevalence [Mod/Severe] among ages 6 to 39 (N4E2Eaa)";
run;
data NHANESAB.N4E2Eaa;
set  NHANESAB.N4E2Eaa;
format tablename $8.;
tablename="N4E2Eaa  ";
run;
proc descript data=NHANESAB.nh4mmwr design=wr;/*E1b. fluorosis prevalence [Mod/Severe] - age/sex adjusted*/
weight WTMEC4YR;
nest SDMVSTRA SDMVPSU/strlev=1 psulev=2 missunit;
var OHDDFIaSX;
catlevel 5;
subpopn ridageyr >=6 and ridageyr<=39;
subgroup  RIDRETH2  ASDIST5B  RIAGENDR PIRGRP3 ASDIST5D;
level         5         4          2        3    68;
tables    ASDIST5B  RIDRETH2 PIRGRP3 
          ASDIST5B*RIDRETH2 ASDIST5B*PIRGRP3;
Stdvar ASDIST5D;    
Stdwgt 0.015018 0.014307 0.015250 0.014539 0.014943 0.014284 0.016208 0.015415
       0.016246 0.015422 0.015445 0.014711 0.015093 0.014419 0.015025 0.014351
       0.015100 0.014351 0.015168 0.014419 0.014831 0.014022 0.015467 0.014569
       0.015070 0.014404 0.015489 0.014943 0.015400 0.014831 0.014351 0.013835
       0.013513 0.013109 0.013199 0.012929 0.012847 0.012682 0.012862 0.012749
       0.012427 0.012487 0.012832 0.012884 0.012862 0.012974 0.015212 0.015385
       0.014868 0.015070 0.014172 0.014351 0.014037 0.014224 0.014232 0.014524
       0.015265 0.015332 0.015984 0.016118 0.016193 0.016365 0.016642 0.016815
       0.015969 0.016156 0.017878 0.017926 	
; 
rformat RIDRETH2 race. ;
rformat riagendr sex. ;
rformat ASDIST5B age. ;
rformat PIRGRP3 pir. ;

print nsum wsum total percent sepercent / style =nchs totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
output nsum wsum total percent sepercent/ 
       filename=NHANESAB.N4E2Easa REPLACE filetype=SAS 
       totalfmt=F13.1 wsumfmt=f15.1 percentfmt=f6.2 sepercentfmt=f6.2 ;
rtitle "E1b. N IV phase I + II age/sex adjusted (1 yr) fluorosis prevalence [Mod/Severe] among ages 6 to 39 (N4E2Easa)";
run;
data NHANESAB.N4E2Easa;
set  NHANESAB.N4E2Easa;
format tablename $8.;
tablename="N4E2Easa ";
run;


proc descript data=NHANESAB.nh4mmwr design=wr;/*E1b. Mean Dean Score  - not age adjusted*/
weight WTMEC4YR;
nest SDMVSTRA SDMVPSU/strlev=1 psulev=2 missunit;
var OHDMDSX;

subpopn ridageyr >=6 and ridageyr<=39;
subgroup  RIDRETH2  ASDIST5B  RIAGENDR PIRGRP3;
level         5         4          2        3 ;
tables    ASDIST5B RIAGENDR RIDRETH2 PIRGRP3 
          ASDIST5B*RIAGENDR ASDIST5B*RIDRETH2 ASDIST5B*PIRGRP3;
rformat RIDRETH2 race. ;
rformat riagendr sex. ;
rformat ASDIST5B age. ;
rformat PIRGRP3 pir. ;

print nsum wsum total mean semean / style =nchs totalfmt=F13.1 wsumfmt=f15.1 meanfmt=f6.2 semeanfmt=f6.2 ;
output nsum wsum total mean semean / 
       filename=NHANESAB.N4E2Fu REPLACE filetype=SAS 
       totalfmt=F13.1 wsumfmt=f15.1 meanfmt=f6.2 semeanfmt=f6.2 ;
rtitle "E1b. N IV phase I + II not age/sex adjusted (1 yr) Mean Dean Score among ages 6 to 39 (N4E2Fu)";
run;
data NHANESAB.N4E2Fu;
set  NHANESAB.N4E2Fu;
format tablename $8.;
tablename="N4E2Fu   ";
run;
proc descript data=NHANESAB.nh4mmwr design=wr;/*E1b. Mean Dean Score   - age adjusted*/
weight WTMEC4YR;
nest SDMVSTRA SDMVPSU/strlev=1 psulev=2 missunit;
var OHDMDSX;

subpopn ridageyr >=6 and ridageyr<=39;
subgroup  RIDRETH2  ASDIST5B  RIAGENDR PIRGRP3 ASDIST5F;
level         5         4          2        3    34;
tables      RIAGENDR 
          ASDIST5B*RIAGENDR;
Stdvar ASDIST5F;    
Stdwgt	0.029554 0.030007 0.029449 0.031862 0.031908 0.030384
	    0.029728 0.029592 0.029675 0.029811 0.029071 0.030263
		0.029698 0.030656 0.030459 0.028400 0.026824 0.026318
		0.025722 0.025805 0.025104 0.025911 0.026032 0.030829
        0.030165 0.021196 0.028468 0.028973 0.030829 0.032345
        0.032813 0.033710 0.032360 0.036079
;
rformat RIDRETH2 race. ;
rformat riagendr sex. ;
rformat ASDIST5B age. ;
rformat PIRGRP3 pir. ;

print nsum wsum total mean semean / style =nchs totalfmt=F13.1 wsumfmt=f15.1 meanfmt=f6.2 semeanfmt=f6.2 ;
output nsum wsum total mean semean / 
       filename=NHANESAB.N4E2Faa REPLACE filetype=SAS 
       totalfmt=F13.1 wsumfmt=f15.1 meanfmt=f6.2 semeanfmt=f6.2  ;
rtitle "E1b. N IV phase I + II age adjusted (1 yr) Mean Dean Score among ages 6 to 39 (N4E2Faa)";
run;
data NHANESAB.N4E2Faa;
set  NHANESAB.N4E2Faa;
format tablename $8.;
tablename="N4E2Faa  ";
run;
proc descript data=NHANESAB.nh4mmwr design=wr;/*E1b. Mean Weighted Dean's   - age/sex adjusted*/
weight WTMEC4YR;
nest SDMVSTRA SDMVPSU/strlev=1 psulev=2 missunit;
var OHDMDSX;

subpopn ridageyr >=6 and ridageyr<=39;
subgroup  RIDRETH2  ASDIST5B  RIAGENDR PIRGRP3 ASDIST5D;
level         5         4          2        3    68;
tables    ASDIST5B  RIDRETH2 PIRGRP3 
          ASDIST5B*RIDRETH2 ASDIST5B*PIRGRP3;
Stdvar ASDIST5D;    
Stdwgt 0.015018 0.014307 0.015250 0.014539 0.014943 0.014284 0.016208 0.015415
       0.016246 0.015422 0.015445 0.014711 0.015093 0.014419 0.015025 0.014351
       0.015100 0.014351 0.015168 0.014419 0.014831 0.014022 0.015467 0.014569
       0.015070 0.014404 0.015489 0.014943 0.015400 0.014831 0.014351 0.013835
       0.013513 0.013109 0.013199 0.012929 0.012847 0.012682 0.012862 0.012749
       0.012427 0.012487 0.012832 0.012884 0.012862 0.012974 0.015212 0.015385
       0.014868 0.015070 0.014172 0.014351 0.014037 0.014224 0.014232 0.014524
       0.015265 0.015332 0.015984 0.016118 0.016193 0.016365 0.016642 0.016815
       0.015969 0.016156 0.017878 0.017926 	
; 
rformat RIDRETH2 race. ;
rformat riagendr sex. ;
rformat ASDIST5B age. ;
rformat PIRGRP3 pir. ;

print nsum wsum total mean semean / style =nchs totalfmt=F13.1 wsumfmt=f15.1 meanfmt=f6.2 semeanfmt=f6.2 ;
output nsum wsum total mean semean/ 
       filename=NHANESAB.N4E2Fasa REPLACE filetype=SAS 
       totalfmt=F13.1 wsumfmt=f15.1 meanfmt=f6.2 semeanfmt=f6.2 ;
rtitle "E1b. N IV phase I + II age/sex adjusted (1 yr) Mean Dean Score among ages 6 to 39 (N4E2Fasa)";
run;
data NHANESAB.N4E2Fasa;
set  NHANESAB.N4E2Fasa;
format tablename $8.;
tablename="N4E2Fasa ";
run;

/*Compile estimates into one data set for export to MSExcel file (.xls)*/
data NHANESAB.N4Fsis;
 set
  NHANESAB.N4E1Au
  NHANESAB.N4E1Aaa
  NHANESAB.N4E1Aasa
  NHANESAB.N4E1Bu
  NHANESAB.N4E1Baa
  NHANESAB.N4E1Basa
  NHANESAB.N4E1Cu
  NHANESAB.N4E1Caa
  NHANESAB.N4E1Casa
  NHANESAB.N4E1Du
  NHANESAB.N4E1Daa
  NHANESAB.N4E1Dasa
  NHANESAB.N4E1Eu
  NHANESAB.N4E1Eaa
  NHANESAB.N4E1Easa
  NHANESAB.N4E1Fu
  NHANESAB.N4E1Faa
  NHANESAB.N4E1Fasa
  NHANESAB.N4E2Au
  NHANESAB.N4E2Aaa
  NHANESAB.N4E2Aasa
  NHANESAB.N4E2Bu
  NHANESAB.N4E2Baa
  NHANESAB.N4E2Basa
  NHANESAB.N4E2Cu
  NHANESAB.N4E2Caa
  NHANESAB.N4E2Casa
  NHANESAB.N4E2Du
  NHANESAB.N4E2Daa
  NHANESAB.N4E2Dasa
  NHANESAB.N4E2Eu
  NHANESAB.N4E2Eaa
  NHANESAB.N4E2Easa
  NHANESAB.N4E2Fu
  NHANESAB.N4E2Faa
  NHANESAB.N4E2Fasa

;

if       ASDIST5B=0 then AgeGrp="Total 6 to 19      ";
else if  ASDIST5B=1 then AgeGrp="Age 6 to 11        ";
else if  ASDIST5B=2 then AgeGrp="Age 12 to 15       ";
else if  ASDIST5B=3 then AgeGrp="Age 16 to 19       ";
else if  ASDIST5B=4 then AgeGrp="Age 20 to 39       ";

if      RIDRETH2=0 then Category="Total Race/Ethn   ";
else if RIDRETH2=1 then Category="NH White          ";
else if RIDRETH2=2 then Category="NH Black          ";
else if RIDRETH2=3 then Category="Mexican American  ";
else if RIDRETH2=4 then Category="Other/MultiRace   ";
else if RIDRETH2=5 then Category="Other Hispanic    ";

/*
if      RACE=0 then Category="Total Race DMARACER   ";
else if RACE=1 then Category="White                 ";
else if RACE=2 then Category="Black                 ";
else if RACE=3 then Category="Other                 ";
else if RACE=4 then Category="MexAm unknown race    ";

if		HISPANIC=0 then Category="Total HISPANIC    ";
else if HISPANIC=1 then Category="Hispanic          ";
else if HISPANIC=2 then Category="Not Hispanic      ";

if 		DMAETHNR=0 then Category="Total DMAETHNR    ";
else if DMAETHNR=1 then Category="Mex. Am DMAETHNR  ";
else if DMAETHNR=2 then Category="OthH. DMAETHNR    ";
else if DMAETHNR=3 then Category="NotH. DMAETHNR    ";
*/

if      RIAGENDR=0 then Category="Total Gender      ";
else if RIAGENDR=1 then Category="Male              ";
else if RIAGENDR=2 then Category="Female            ";

if      PIRGRP3=0 then  Category="Total PIR         ";
else if PIRGRP3=1 then  Category="<100              ";
else if PIRGRP3=2 then  Category="100-199           ";
else if PIRGRP3=3 then  Category=">=200             ";

/* if      DMDEDUx=0 then Category="Total Education   ";
else if DMDEDUx=1 then Category="<HS               ";
else if DMDEDUx=2 then Category="HS/GED            ";
else if DMDEDUx=3 then Category=">HS               ";
*/

/* if      HFHEDUx=0 then Category="Total HH Education";
else if HFHEDUx=1 then Category="<HS               ";
else if HFHEDUx=2 then Category="HS/GED            ";
else if HFHEDUx=3 then Category=">HS               ";
*/

/* if 		DISABLEX=0 then Category="Total DisableX   ";
else if DISABLEX=1 then Category="With Disabilities";
else if DISABLEX=2 then Category="No Disabilities  ";
*/

if percent>0 then do;
 RSE=sepercent/percent; /*Relative Standard Error*/
end;

if mean>0 then do;
 RSE=semean/mean; /*Relative Standard Error*/
end;

keep
 tablename AgeGrp Category nsum rse percent sepercent mean semean
;

run;
PROC EXPORT DATA= NHANESAB.N4Fsis 
            OUTFILE= "K:\DOH\EVERYONE\NHANES\MMWR19992002\N4Fsis.xls" 
            DBMS=EXCEL2000 REPLACE;
RUN;

/*fluorosis.sas*/