Ignore:
Timestamp:
10/13/20 13:10:40 (7 weeks ago)
Author:
LoisHaggard_NM
Message:

NMBIS Backend QModules

File:
1 edited

Legend:

Unmodified
Added
Removed
  • adopters/nm/trunk/src/main/backend_qModules/3.0/_EPHT/childlead/CrudeRateNM_AnnZERO.def

    r19939 r21775  
    139139  if number=. then n=0; else n=number;
    140140  rate=n/popnum;
    141   rateper=rate*100;
     141  PctPersons=rate*100;
    142142  stderr=sqrt(rate*(1-rate)/popnum)*100;
    143143  if n=0 then stderr=sqrt((3/popnum)*(1-(3/popnum))/popnum)*100;
    144   t1=(rateper-(1.96*stderr));
     144  t1=(PctPersons-(1.96*stderr));
    145145  if (t1<0) then t1=0;
    146146  LL=put(t1, 8.2);
    147   UL=put((rateper+(1.96*stderr)), 8.2);
     147  UL=put((PctPersons+(1.96*stderr)), 8.2);
    148148  LL=compress(LL);
    149149  UL=compress(UL);
     
    159159/*************************************************/
    160160
    161   rse=(stderr/rateper);
    162   if rse>.3 then redflag=1; *statistical stability – 1 red flag;
    163   if rse>.5 then redflag=2; *statistical stability - 2 red flags;
    164   if stderr=. then redflag=2; *no variance, n=1;
    165 keep %cross1% %cross2% n popnum rate rateper LL UL rse redflag;
     161        if n>0 then do;
     162                rse=(stderr/PctPersons);
     163                redflag=put('', $12.);
     164                if rse>.3 then redflag=put('Unstable', $12.);
     165                if rse>.5 then redflag=put('VeryUnstable', $12.);
     166                if stderr=. then redflag=put('Unstable', $12.);
     167        end;
     168        if n<=0 then redflag=put('Unstable', $12.);
     169
     170keep %cross1% %cross2% n popnum rate PctPersons LL UL rse redflag;
    166171proc print data=tmp; title2 'TMP AGAIN';
    167172run;
     
    176181/*********************************************
    177182
     183
    178184data tmp;
    179   set tmp;
    180   if ((0<n<4) and (1<=popnum<5000)) then do;
    181   n = .A;
    182   rate = .A;
    183 end;
     185        set tmp;
     186        if ((0<n<4) and (.<popcount<5000)) then do;
     187                rateper = .A;
     188                n = .A;
     189                popcount = .A;
     190                LL = put('**', 8.0);
     191                UL = put('**', 8.0);
     192                redflag = put('Suppressed', $12.);
     193        end;
     194if popcount=. then redflag=put('', $12.);
     195
     196
    184197*/
    185198proc print data=tmp; title2 'FINAL'; run;
     
    188201
    189202--------BoNdArY--------
    190 f out_variable rateper
     203f out_variable PctPersons
    191204# definition for output file
    192205f xml_out_map_file XMLRateNumerDenomLCLUCL.map
     
    196209--------BoNdArY--------
    197210f out_detail lbl_not_used__see_xml_out_map_file
    198   rateper 15.3
     211  PctPersons 15.3
    199212  n 15.0
    200213  popnum 15.0
     
    2062191 suppressed_variables 2
    207220n
    208 rateper
    209 --------BoNdArY--------
     221PctPersons
     222--------BoNdArY--------
Note: See TracChangeset for help on using the changeset viewer.