qrp_report

****************************************************************************************************

*                                           PROGRAM OVERVIEW

****************************************************************************************************

*

* PROGRAM: l2_psdistribution_createdata.sas  

* Created (mm/dd/yyyy): 03/17/2021

*

*--------------------------------------------------------------------------------------------------

* PURPOSE: The macro creates both aggregated and DP specific PS distribution datasets

*                                        

*  Program inputs:                                                                                   

* 

*  Program outputs:                                                                                                                                       

*

*  PARAMETERS:                                                                       

*            

*  Programming Notes:                                                                                

*                                                                           

*

*--------------------------------------------------------------------------------------------------

* CONTACT INFO: 

*  Sentinel Coordinating Center

*  info@sentinelsystem.org

*

***************************************************************************************************;

​

%macro l2_psdistribution_createdata;

​

    %do loopcount = 1 %to &numl2comparisons.; 

        /*Initialize macro variables for the analysisgrp loop*/

        %let OutputPSDistribution = ;

        %let ratio = ;

        %let strataweight = ;

​

        /*set parameters from l2comparisonfile for this loop*/

        data _null_;

            set l2comparisonfile(where=(order=&loopcount.));

            call symputx('runid', runid);

            call symputx('analysisgrp', analysisgrp);

            call symputx('OutputPSDistribution', OutputPSDistribution);

        run;

​

        %if &OutputPSDistribution.=Y %then %do;

        /* only psmatchfile, stratificationfile, and iptwfile have histogram */

            proc sql noprint;

                select strip(file) into: psfile

                from pscs_masterinputs

                where analysisgrp = "&analysisgrp.";

            quit;

​

            %put now computing PS distributions for &analysisgrp.;

​

            %if &psfile. = psmatchfile | &psfile. = stratificationfile | &psfile. = iptwfile %then %do;

                data _null_; 

                  set pscs_masterinputs (where=(lowcase(analysisgrp)="&analysisgrp."));

                    call symputx("psestimategrp", lowcase(psestimategrp));

                         %if &psfile. = psmatchfile  %then %do;

                            call symputx('RATIO',upcase(ratio)) ;

                         %end;

                         %if &psfile. = stratificationfile %then %do;

                            call symputx("strataweight",strataweight);

                         %end;

                run; 

                data _null_; 

                  set psest_masterinputs (where=(psestimategrp="&psestimategrp."));

                  call symputx('eoi', eoi) ;

                  call symputx('ref', ref) ;    

                run;

            %end;

​

            %if &psfile. = psmatchfile | &psfile. = stratificationfile | &psfile. = iptwfile %then %do;

​

                /* Determine weights on input dataset to use for squaring */

                proc sql noprint;

                  select count(distinct(weight)) into: numweights trimmed

                  from agg_psdistribution_&periodid. (where=(runid="&runid." and lowcase(analysisgrp)="&analysisgrp."));

                  select distinct(weight) into: weight1 - :weight&numweights.

                  from agg_psdistribution_&periodid. (where=(runid="&runid." and lowcase(analysisgrp)="&analysisgrp."));

                quit;

​

                /*Create bins dataset*/

                %let start = -0.0125;

                    data bins;

                    length type $20;

                      %do bin = 1 %to 41;

                          ps_cat = &bin.;

                          bin_eoi = &start.;

                          bin_ref = &start.;

                          %let start = &start. + 0.025;