2006 HCSDB Adult Sampling Report - Tricare

COM = &POPSIZE*SQRT(&POPSIZE/(&POPSIZE-1));NUM = COM/DSUM&ITE;RUN;PROC MEANS DATA=&DAT NOPRINT;VAR NUM DEN COM;BY &DOMAIN;OUTPUT OUT=DSIZEA SUM=NUMS DENS COMS;RUN;DATA &ODAT;MERGE &DAT DSIZEA;BY &DOMAIN;ND=(&SSQUARE*NUMS**2)/(&V0+&SSQUARE*DENS);NHO=ND*COM/COMS;DROP ND NUM DEN COM NUMS DENS COMS;RUN;%MEND OPTALLO;/*----------------------------------------------------------------TO RETREIVE THE NUMBER OF OBSERVATIONS IN A SAS DATA SET------------------------------------------------------------------*/%MACRO NUMOBS(DSN);%GLOBAL NUM;/* THIS MACRO CONTAINS THE NUMBER OF OBS IN THE DATA*/DATA _NULL_;IF 0 THEN SET &DSN NOBS=COUNT;CALL SYMPUT(’NUM’,LEFT(PUT(COUNT,8.)));STOP;RUN;%MEND NUMOBS;/*----------------------------------------------------------------ITERATE UNTIL THE REMAINING DOMAINS HAVE NHO GREATER THANTHE PREVIOUS SAMPLE SIZES-----------------------------------------------------------------*/%MACRO ITERATE;%OPTALLO(STE,DOM&ITE,POPSIZE,VSTAR,OSTAT);DATA FIN&I STE;SET OSTAT;IF NHF < NHO THEN FIN = FIN +1;IF FIN=&I then output FIN&I;IF FIN = &I + 1 then output STE;RUN;%VAR(FIN&I,DOM&ITE,POPSIZE,NHF,SUMMARY);DATA STE;MERGE STE (IN=A) SUMMARY ;BY DOM&ITE;IF A;IF VSUM=. THEN VSUM=0;****SHOULD EXIST!!!;VSTAR= VSTAR - VSUM/DSUM&ITE**2;DROP VSUM;RUN;%MEND ITERATE;/*--------------------------------------------------------------------MAIN PART OF THE PROGRAM: ’ITE’ INDICATES THE LEVEL OF DOMAINS---------------------------------------------------------------------*/%MACRO MPART(ITE);PROC SORT data=indata;BY DOM&ITE;RUN;F-2