************************************************************************ * Program: C:\NHANES\.sas * * Proposal: * ************************************************************************; LIBNAME NH "C:\NHANES\DATA"; OPTIONS NODATE NOCENTER; option ls=72; 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 hbpfmt 1 = 'High BP' 2 = 'No high BP' ; VALUE hbp2fmt 100 = 'High BP' 0 = 'No high BP' ; run; data ANALYSIS_DATA; merge 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; 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; *** code to define mean blood pressure measures and HBP ****; n_sbp = n(of bpxsy1-bpxsy4); n_dbp = n(of bpxdi1-bpxdi4); *Setting DBP values of 0 as missing for calculating average; array _DBP bpxdi1-bpxdi4; do over _DBP; if (_DBP = 0) then _DBP = .; end; mean_sbp = mean(of bpxsy1-bpxsy4); mean_dbp = mean(of bpxdi1-bpxdi4); if BPQ050a=1 then HBP_trt=1; else if BPQ020 in (1,2) and BPQ050a < 7 then HBP_trt=0; if n_sbp>0 and n_dbp>0 then do; if mean_sbp>=140 then SBP140=1; else SBP140=0; if mean_dbp>=90 then DBP90=1; else DBP90=0; end; if HBP_trt>=0 and SBP140>=0 and DBP90>=0 then do; if HBP_trt=1 or SBP140=1 or DBP90=1 then HBP=1; else HBP=2; end; hbpx=.; if hbp =1 then hbpx =100; if hbp =2 then hbpx =0; LABEL age = 'AGE GROUP' race = 'Race Ethnicity' riagendr = 'Gender' ; RUN; PROC SORT DATA=ANALYSIS_DATA OUT=ANALYSIS_DATA; BY sdmvstra sdmvpsu ; RUN; PROC descript data=ANALYSIS_DATA design=wr ; subpopn ridageyr >=20 ; **subset the data to reflect the subgroup to be analyzed; NEST sdmvstra sdmvpsu; weight wtmec4yr; subgroup riagendr age race; levels 2 3 4 ; var hbp; catlevel 1; table riagendr * age * race ; print nsum= "Sample Size" percent="Percent" sepercent="SE" / nohead notime style=NCHS nsumfmt=f8.0 percentfmt=f8.4 sepercentfmt=f8.4 ; rformat riagendr sexfmt.; rformat age agefmt.; rformat race racefmt.; rformat hbp hbpfmt.; rtitle "Prevalence of SPs with measured high blood pressure : NHANES 1999-2002"; run;