Changeset 20054 in main


Ignore:
Timestamp:
03/18/20 09:04:45 (13 days ago)
Author:
Paul Leo
Message:

Promoting modifcations to Birth Average and Percentage backend files to production

File:
1 edited

Legend:

Unmodified
Added
Removed
  • adopters/mt/trunk/src/main/backend_qModules/birth23/MT_Average_data_frame.def

    r20026 r20054  
    225225 *              2a. If count > 20                                                                                                               ;
    226226 *                      Report count and rates with confidence intervals                                        ;
    227  *              2b. If count > 20                                                                                                               ;
     227 *              2b. If count < 20                                                                                                               ;
    228228 *                      Suppress rates and counts                                                                                       ;
    229229 * ZW-s program uses ".A" to identify cells for suppression.                            ;
     
    235235 * code at the end of this file to work.                                                                                ;
    236236 *******************************************************************************;
    237  
    238  **************** ASK LOIS ABOUT STDERR COMPUTATIONS ***************************;
    239237
    240238        data tmp;
     
    244242
    245243                if popcount>=300 then do;
    246 
    247244                        *   Suppress cells if absolute count <5 ;
    248245                       
    249246                        if 1<=count<5 then do;
     247                                n=.A;
    250248                                mean=.A;
    251                                 LL=.A;
    252                                 UL=.A;
     249                                LL=put('**',8.0);
     250                                UL=put('**',8.0);
    253251                                redflag=put('Suppressed', $15.);
    254                                 n=.A;
    255                         end;
    256 
    257                                 *   do not compute rates based on fewer than 20 events, but provide counts;
     252                        end;
     253
     254                        *   do not compute rates based on fewer than 20 events, but provide counts;
    258255                       
    259256                        if 5<=count<20 then do;
    260257                                mean=.A;
    261                                 LL=.A;
    262                                 UL=.A;
     258                                LL=put('**',8.0);
     259                                UL=put('**',8.0);
    263260                                redflag=put('Suppressed Avg', $16.);
    264261                        end;
    265262
    266                         * for events >= 20 use binomial exact confidence intervals;
    267                        
    268                         if count >= 20, then do;
    269                                 LL=put(t1, 8.2);
    270                                 LL=(1 – betainv (.975, popcount – count – 1, count))*100);
    271                                 UL=((1 – betainv (.025, popcount – count, count + 1))*100);
    272                                 LL=compress(LL);
    273                                 UL=compress(UL);
     263                        * for events between 20 and 100, use T distribution to calculate confidence intervals;
     264                        * note LL and UL are calculated in step 4 above.;
     265
     266                        if 20<=count<=100 then do;
     267                       
    274268                                redflag=put('-', $15.);
    275                         end;
    276                        
     269
     270                        end;
     271
     272                        * and normal approximation for more than 100 events. ;
     273                        * which is also computed in step 4 above;
     274
     275                        if count>100 then do;
     276
     277                                redflag=put('-', $15.);
     278                        end;
    277279                end;
    278280               
    279281                if popcount<300 then do;
    280282
    281                         *  suppress if popcount <300 and count < 20 events.             ;
     283                        *  suppress if popcount <300 and count < 20 events;
    282284               
    283285                        if 1<=count<20 then do;
     286                                n=.A;
    284287                                mean=.A;
    285                                 LL=.A;
    286                                 UL=.A;
     288                                LL=put('**',8.0);
     289                                UL=put('**',8.0);
    287290                                redflag=put('Suppressed', $15.);
    288                                 n=.A;
    289                         end;
    290                        
    291                        
    292                         * report rate and n if popcount <300 and count >= 20 events.    ;
    293 
    294                         * for events >= 20 use binomial exact confidence intervals;
    295                        
    296                         if count >= 20, then do;
    297                                 LL=put(t1, 8.2);
    298                                 LL=(1 – betainv (.975, popcount – count – 1, count))*100);
    299                                 UL=((1 – betainv (.025, popcount – count, count + 1))*100);
    300                                 LL=compress(LL);
    301                                 UL=compress(UL);
     291                        end;
     292                       
     293                       
     294                        * report rate and n if popcount <300 and count >= 20 events;
     295
     296                        if 20<=count<=100 then do;
     297                        * for events between 20 and 100, use T distribution to calculate confidence intervals;
     298                        * note LL and UL are calculated in step 4 above.;
     299                       
     300                                redflag=put('-', $15.);
     301
     302                        end;
     303
     304                        *  and normal approximation for more than 100 events. ;
     305                        * which is also computed in step 4 above;
     306                       
     307                        if count>100 then do;
     308
    302309                                redflag=put('-', $15.);
    303310                        end;
     
    306313                if count=0 then do;
    307314                        mean=.A;
    308                         LL=.A;
    309                         UL=.A;
     315                        LL=put('**',8.0);
     316                        UL=put('**',8.0);
    310317                end;
    311318
    312319                if count<0 then do;
    313                         mean=0;
    314320                        mean=.A;
    315                         LL=.A;
    316                         UL=.A;
     321                        LL=put('**',8.0);
     322                        UL=put('**',8.0);
    317323                        redflag=put('No Variance', $15.);  *no variance;
    318324                end;
     
    347353 ********************************************************************;
    348354--------BoNdArY--------
    349 1 suppressed_variables 4
     3551 suppressed_variables 2
    350356mean
    351357n
    352 LL
    353 UL
    354 --------BoNdArY--------
     358--------BoNdArY--------
Note: See TracChangeset for help on using the changeset viewer.