Changeset 21817 in main


Ignore:
Timestamp:
10/16/20 18:03:06 (2 weeks 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

Location:
adopters/hi/branches/2.3/src/main/backend_qModules/brfss
Files:
2 edited

Legend:

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

    r21797 r21817  
    5151                from tmp, yrnum
    5252        quit;
    53    
    54  data tmp;
     53      proc print data=tmpyrwgt (obs=10);
     54            title1 ' '; title2 'tmpyrwgt';
     55          run;
     56
     57data tmp;
    5558     set tmpyrwgt;
    5659         if upcase(%cross1%) ^=YEAR
     
    172175data tmp4;
    173176  set tmp2 tmp1;
    174   f=log(mean)-log(1-mean); 
    175   s=stderr/(mean*(1-mean));
    176   Lf=f-1.96*s; 
    177   Uf=f+1.96*s; 
    178   percent=mean*100;
    179   lower=(exp(Lf)/(1+exp(Lf)))*100;   
    180   upper=(exp(Uf)/(1+exp(Uf)))*100;
     177        if mean >0 then do;
     178                f=log(mean)-log(1-mean); 
     179                s=stderr/(mean*(1-mean));
     180                Lf=f-1.96*s; 
     181                Uf=f+1.96*s; 
     182                percent=mean*100;
     183          lower=(exp(Lf)/(1+exp(Lf)))*100;   
     184          upper=(exp(Uf)/(1+exp(Uf)))*100;
     185  end;
    181186  drop f s Lf Uf VarName;
    182187  proc sort data=tmp4; by %cross1%
  • adopters/hi/branches/2.3/src/main/backend_qModules/brfss/SurveyCrudeRateNoVarLevel.sas

    r21802 r21817  
    5050                from tmp, yrnum
    5151        quit;
    52    
     52        /*
     53      proc print data=tmpyrwgt (obs=10);
     54        title1 ' '; title2 'tmpyrwgt';
     55    */   
     56           run;
     57
    5358 data tmp;
    5459     set tmpyrwgt;
Note: See TracChangeset for help on using the changeset viewer.