Changeset 21802 in main


Ignore:
Timestamp:
10/16/20 11:05:29 (5 days ago)
Author:
Paul Leo
Message:

Hawaii IBIS v2.3 BRFSS back-end .def and .sas programs
Correcting most, if not all issues with rounding and cell suppression.
Improving a comment in the code, preventing divide by zero
Still have questions though:
1) need to verify cell suppression rules
2) need to verify Statistical Stability output rules, and possibly the data_notes section describing the Statistical Stability output rules.
3) need to get Statistical Stability column working for SurveyAARateDist9plus.sas, which could be the .def file, or the module.xml

File:
1 edited

Legend:

Unmodified
Added
Removed
  • adopters/hi/branches/2.3/src/main/backend_qModules/brfss/SurveyCrudeRateNoVarLevel.sas

    r21796 r21802  
    1 /**********************************************************************;
     1***********************************************************************;
    22* Program filename: SurveyCrudeRateNoVarLevel.sas                      ;
    33*                                                                      ;
     
    5858     _llcpwt=_llcpwt/nyears;
    5959   run;
    60    
    61 /*proc print data=tmp (obs=10);
     60/*   
     61proc print data=tmp (obs=10);
    6262        title1 ' '; title2 'checktmp .'; run;
    6363*/
     
    161161data tmp4;
    162162  set tmp2 tmp1;
    163   f=log(mean)-log(1-mean); 
    164   s=stderr/(mean*(1-mean));
    165   Lf=f-1.96*s; 
    166   Uf=f+1.96*s; 
    167   percent=mean*100;
    168   lower=(exp(Lf)/(1+exp(Lf)))*100;   
    169   upper=(exp(Uf)/(1+exp(Uf)))*100;
     163        if mean >0 then do;
     164                f=log(mean)-log(1-mean); 
     165                s=stderr/(mean*(1-mean));
     166                Lf=f-1.96*s; 
     167                Uf=f+1.96*s; 
     168                percent=mean*100;
     169          lower=(exp(Lf)/(1+exp(Lf)))*100;   
     170          upper=(exp(Uf)/(1+exp(Uf)))*100;
     171  end;
    170172  drop f s Lf Uf VarName;
    171173  proc sort data=tmp4; by %cross1%
     
    229231run;
    230232
     233************************************************************;
    231234************************************************************;
    232235********** tmp7: Add RSE, redflag to dataset, clean up *****;
     
    282285data tmp;
    283286  set tmp7;
    284  
     287
    285288  if (SampleN<50) OR  (rse >.3) then do;        /* Hawaii cell supression Rule */
    286289
Note: See TracChangeset for help on using the changeset viewer.