Changeset 20072 in main


Ignore:
Timestamp:
03/19/20 16:21:49 (2 weeks ago)
Author:
Paul Leo
Message:

finally corrected the confidence intervals for percentages

File:
1 edited

Legend:

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

    r20053 r20072  
    225225                       
    226226                        if count >= 20 then do;
    227                                 diffcountL=popcount-count-1;
    228                                 diffcountU=popcount-count;
    229                                 countplus1=count+1;
    230                         /************* Keeping lines for steps to obtain  LL and UL    *************
    231                                 Lowerquant=betainv(.975,diffcountL,count);
    232                                 Upperquant=betainv(.025,diffcountU,countplus1);
    233                                 Lowerval=(1-Lowerquant)*100;
    234                                 Upperval=(1-Upperquant)*100;
    235                                 LL=rateper-Lowerval;
    236                                 UL=rateper+Upperval;
    237                         ***************************************************************************/
    238                                 LL=rateper-(1-(betainv(.975,diffcountL,count)))*100;
    239                                 UL=rateper+(1-(betainv(.025,diffcountU,countplus1)))*100;
     227                                LL=(1-betainv(.975,totalnum-count+1,count))*100;
     228                                UL=(1-betainv(.025,totalnum-count,count+1))*100;
    240229                                LL=compress(LL);
    241230                                UL=compress(UL);
     
    271260                       
    272261                        if count >= 20 then do;
    273                                 diffcountL=popcount-count-1;
    274                                 diffcountU=popcount-count;
    275                                 countplus1=count+1;
    276                         /************* Keeping lines for steps to obtain  LL and UL    *************
    277                                 Lowerquant=betainv(.975,diffcountL,count);
    278                                 Upperquant=betainv(.025,diffcountU,countplus1);
    279                                 Lowerval=(1-Lowerquant)*100;
    280                                 Upperval=(1-Upperquant)*100;
    281                                 LL=rateper-Lowerval;
    282                                 UL=rateper+Upperval;
    283                         ***************************************************************************/
    284                                 LL=rateper-(1-(betainv(.975,diffcountL,count)))*100;
    285                                 UL=rateper+(1-(betainv(.025,diffcountU,countplus1)))*100;
     262
     263                                LL=(1-betainv(.975,totalnum-count+1,count))*100;
     264                                UL=(1-betainv(.025,totalnum-count,count+1))*100;
    286265                                LL=compress(LL);
    287266                                UL=compress(UL);
Note: See TracChangeset for help on using the changeset viewer.