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

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+'       
				; 
 
  VALUE age1fmt  1='20-29 years' 
			     2='30-39 years' 
			     3='40-49 years' 
			     4='50-59 years'
			     5='60-69 years' 
			     6='70 years and older'; 
 run;   
data ANALYSIS_DATA;
set NH.ANALYSIS_DATA; 

if ridstatr = 2;  ***examined  ;


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; 

age1=.; 
age1=1+(ridageyr>29)+(ridageyr>39)+(ridageyr>49)+(ridageyr>59)+(ridageyr>69);

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                           
       age1  = 'AGE GROUP'    
	   race = 'Race Ethnicity'	
       lbxtr='Serum triglyceride' 
       riagendr = 'Gender' 
          ;                                                       
                                                                                                              
run;





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

proc descript geometric atlevel1=1 atlevel2=2;
nest sdmvstra sdmvpsu;
weight wtsaf4yr;
subpopn ridageyr>=20/name="Adults 20 years and older"; 
class age1 riagendr/nofreq;
var lbxtr;
table riagendr*age1;
print nsum geomean segeomean/style=NCHS geomeanfmt=f6.0 segeomeanfmt=f6.1;
output nsum geomean segeomean atlev1 atlev2/filename=tg9902 replace;
rformat age1 age1fmt.; 
rformat riagendr sexfmt.; 
RUN;

data newtg9902;
set tg9902;
df=atlev2-atlev1;
drop PROCNUM TABLENO VARIABLE _C1 _C2 ATLEV1 ATLEV2;
ll=round(geomean+tinv(.025,df)*segeomean);
ul=round(geomean+tinv(.975,df)*segeomean);
geomean=round(geomean);segeomean=round(segeomean,.1);
ciwidth=ul-ll;

proc print split='/' noobs; format age1 age1fmt. riagendr sexfmt. nsum 7.0 geomean 6.0 segeomean 6.1 df 2.0;
label ll='Lower'/'Limit' ul='Upper'/'limit' df='Degrees'/'of'/'freedom'
ciwidth='Confidence'/'interval'/'width';
title1 'Geometric mean of serum triglyceride and 95 % Confidence';
title2 'interval of adults 20 years and older:';
title2 'United States, 1999-2002'; run;