Changeset 21542 in main


Ignore:
Timestamp:
09/15/20 16:58:12 (5 days ago)
Author:
Paul Leo
Message:

Maryland BRFSS backend .def and .SAS programs
Changed SurveyAARateDist8.sas so that if wgt_percent > .95 it will use RSE=(StdErr?/wgt_percent) instead of RSE=(StdErr/?(1-wgt_percent).
NOTE: need to do this for all AARate.sas programs for MD
Changed SurveyCrudeRate?.sas to use Denom <50 instead of Denom < 30 which is what CDC suggests
NOTE: need to do this for all CrudeRate?.sas programs for MD
Changed RateCrude?.def to include the correct number of suppressed variables which is actually 5, not 4

Location:
adopters/md/trunk/src/main/backend_qModules/brfss23
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • adopters/md/trunk/src/main/backend_qModules/brfss23/RateCrude.def

    r18175 r21542  
    2828--------BoNdArY--------
    2929--------BoNdArY--------
    30 1 suppressed_variables 4
     301 suppressed_variables 5
    3131percent
    3232wgtN
  • adopters/md/trunk/src/main/backend_qModules/brfss23/SurveyAARateDist8.sas

    r19541 r21542  
    525525* redflag is the statistical stability indicator, based on                ;
    526526* the relative standard error (RSE, or coefficient of  variation.         ;
    527 *  now using wgt_percent  which is mean*weight to calculate RSE           ;
     527* now using wgt_percent  which is mean*weight to calculate RSE            ;
    528528**************************************************************************;
    529529data Stability;
     
    531531   
    532532        redflag=put('Stable', $14.);
    533        
    534   if 0<wgt_percent<.50 then RSE=(StdErr/wgt_percent);
    535   if .50<=wgt_percent<1 then RSE=(StdErr/(1-wgt_percent));
     533 
     534  RSE=(StdErr/wgt_percent);
     535  /* if 0<wgt_percent<.50 then RSE=(StdErr/wgt_percent);  */
     536  if .50<=wgt_percent<.95 then RSE=(StdErr/(1-wgt_percent));
     537
     538 
     539/*      RSE=(StdErr/wgt_percent); */
    536540 
    537541   if ((wgt_percent = 0) or (wgt_percent = 1)) then RSE=0; /* No Variance */
    538542   if RSE = . then RSE = 0; * is this reasonable?;
    539    ******************/
    540543   
    541544  if RSE >.3 then redflag=put('Unstable', $14.);
    542545  if (wgt_percent = 0) OR (wgt_percent = 1) then redflag=put('No Variance', $14.);
    543546 
    544   proc print data=Stability noobs;  title2 'Stability labels prior to setting totals to '.' and rounding and multiply by 100,';
     547  if &varname. = '.' then delete;
     548 
     549  proc print data=Stability noobs;  title2 'Stability labels prior to setting totals to '.' and rounding and multiply by 100, deleting cross1 and cross2 totals';
    545550run;
    546551
     
    560565   LL=100*LL;
    561566   UL=100*UL;
     567
     568/* ask Lois about the following  */
     569
     570        if UL > 100 then UL = 100;
    562571
    563572   run;
     
    595604   
    596605   /* if &varname in (%spvar2%); *This is the value for the indicator dimension passed in by the URL.;*/
    597         if &varname = '%spvar2%';
     606        if &varname in (%spvar2% .);
    598607        if (0<SampleN < 50) OR  (RSE >.3) then do;
    599608                 ageadj_percent = .A;
  • adopters/md/trunk/src/main/backend_qModules/brfss23/SurveyCrudeRate.sas

    r21330 r21542  
    284284  ************************************************************************;
    285285  if &varname in (%spvar2% .);
    286   if (0<Denom<30)  OR  (CV >.3) then do;        /* MD cell supression Rule  Denom (number of folks who answered either y or n)*/
     286  if (0<Denom<50)  OR  (CV >.3) then do;        /* MD cell supression Rule  Denom (number of folks who answered either y or n)*/
    287287        wgtN  = .A;
    288288        percent = .A;
Note: See TracChangeset for help on using the changeset viewer.