Changeset 19761 in main


Ignore:
Timestamp:
01/28/20 09:27:18 (4 weeks ago)
Author:
Maria Baron
Message:

latest NJVDRS edits

File:
1 edited

Legend:

Unmodified
Added
Removed
  • adopters/nj/trunk/src/main/ibisq/qModules-2.2/vdrs/AgeRate.def

    r19409 r19761  
    44#f data_where %cross1% %cross2%^=99
    55# pop_where %cross2%=1&and&15<=age<=44
     6f small_num 20
     7f small_pop 20
     8
    69#########################################
    710--------BoNdArY--------
    8111 script
    9 data tmp;
    10   set tmp;
    11 x=1;
    12 agepop=agegroup;
    13 run;
    14 
     12 
    1513 proc summary data=tmp;
    1614 var x;
     
    3028 proc summary data=poptmp;
    3129 var popcount;
    32  class agegrp  %popcross1% %popcross2%;
     30 class agepop  %popcross1% %popcross2%;
    3331 output out=pop sum=popnum;
    3432 run;
     
    3634data pop;
    3735  set pop;
    38 if agegrp=. then agegrp=-1;
     36if agepop=. then agepop=-1;
    3937?popcross1? if %popcross1%=. then %cross1%=-1;
    4038?popcross2? if %popcross2%=. then %cross2%=-1;
     
    4644 set pop.morstd;
    4745 if (stdyr=2000);
    48 stdpop=stdpop/100; /*add to 1*/
    4946 run;
    5047 proc sort data=stdpop;
    5148 by agepop;
    5249 run;
    53 proc print data=stdpop;
    54 title 'stdpop';
    55 run;
    5650
    5751 proc sql;
     
    5953 select tmp.*, pop.*, stdpop.*
    6054 from tmp, pop, stdpop
    61  where tmp.agepop=pop.agegrp
     55 where tmp.agepop=pop.agepop
    6256 and tmp.agepop=stdpop.agepop
    6357 ?popcross1? and tmp.%cross1%=pop.%popcross1%
     
    9084 select number.*, pop.*
    9185 from number, pop
    92  where pop.agegrp=-1
     86 where pop.agepop=-1
    9387?popcross1? and number.%cross1%=pop.%popcross1%
    9488?popcross2? and number.%cross2%=pop.%popcross2%
     
    10599?cross2? if %cross2%=. then delete;
    106100
    107  *rate=rate/100;
     101 rate=rate/100;
    108102 stdev=1.96*sqrt(rate*(1-rate)/popnum);
    109103 t1=(rate-stdev)*100000;
     
    111105 r1=put(t1, 8.2);
    112106 r2=put((rate+stdev)*100000, 8.2);
    113 rate=rate*100000;
    114 cv=100*sqrt((100000/popnum)/rate);
    115 
     107 rate=rate*100000;
    116108 r1=compress(r1);
    117109 r2=compress(r2);
     
    138130 n 15.0
    139131 popnum 15.0
    140  r1 15.2
    141  r2 15.2
    142  cv 7.4
     132 r1 15.
     133 r2 15.
    143134--------BoNdArY--------
    144 --------BoNdArY--------
    145 1 suppressed_variables 2
    146 rate
    147 n
    148 --------BoNdArY--------
    149 
Note: See TracChangeset for help on using the changeset viewer.