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/qModules-2.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*(1-rate)/popnum); 109 103 t1=(rate-stdev)*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 --------BoNdArY--------145 1 suppressed_variables 2146 rate147 n148 --------BoNdArY--------149
Note: See TracChangeset
for help on using the changeset viewer.