Changeset 19761 in main
 Timestamp:
 01/28/20 09:27:18 (3 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

adopters/nj/trunk/src/main/ibisq/qModules2.2/vdrs/AgeRate.def
r19409 r19761 4 4 #f data_where %cross1% %cross2%^=99 5 5 # pop_where %cross2%=1&and&15<=age<=44 6 f small_num 20 7 f small_pop 20 8 6 9 ######################################### 7 10 BoNdArY 8 11 1 script 9 data tmp; 10 set tmp; 11 x=1; 12 agepop=agegroup; 13 run; 14 12 15 13 proc summary data=tmp; 16 14 var x; … … 30 28 proc summary data=poptmp; 31 29 var popcount; 32 class age grp %popcross1% %popcross2%;30 class agepop %popcross1% %popcross2%; 33 31 output out=pop sum=popnum; 34 32 run; … … 36 34 data pop; 37 35 set pop; 38 if age grp=. then agegrp=1;36 if agepop=. then agepop=1; 39 37 ?popcross1? if %popcross1%=. then %cross1%=1; 40 38 ?popcross2? if %popcross2%=. then %cross2%=1; … … 46 44 set pop.morstd; 47 45 if (stdyr=2000); 48 stdpop=stdpop/100; /*add to 1*/49 46 run; 50 47 proc sort data=stdpop; 51 48 by agepop; 52 49 run; 53 proc print data=stdpop;54 title 'stdpop';55 run;56 50 57 51 proc sql; … … 59 53 select tmp.*, pop.*, stdpop.* 60 54 from tmp, pop, stdpop 61 where tmp.agepop=pop.age grp55 where tmp.agepop=pop.agepop 62 56 and tmp.agepop=stdpop.agepop 63 57 ?popcross1? and tmp.%cross1%=pop.%popcross1% … … 90 84 select number.*, pop.* 91 85 from number, pop 92 where pop.age grp=186 where pop.agepop=1 93 87 ?popcross1? and number.%cross1%=pop.%popcross1% 94 88 ?popcross2? and number.%cross2%=pop.%popcross2% … … 105 99 ?cross2? if %cross2%=. then delete; 106 100 107 *rate=rate/100;101 rate=rate/100; 108 102 stdev=1.96*sqrt(rate*(1rate)/popnum); 109 103 t1=(ratestdev)*100000; … … 111 105 r1=put(t1, 8.2); 112 106 r2=put((rate+stdev)*100000, 8.2); 113 rate=rate*100000; 114 cv=100*sqrt((100000/popnum)/rate); 115 107 rate=rate*100000; 116 108 r1=compress(r1); 117 109 r2=compress(r2); … … 138 130 n 15.0 139 131 popnum 15.0 140 r1 15.2 141 r2 15.2 142 cv 7.4 132 r1 15. 133 r2 15. 143 134 BoNdArY 144 BoNdArY145 1 suppressed_variables 2146 rate147 n148 BoNdArY149
Note: See TracChangeset
for help on using the changeset viewer.