qrp_report

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

*                                           PROGRAM OVERVIEW

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

*

* PROGRAM: figure_survivalcurves_output.sas  

* Created (mm/dd/yyyy): 07/28/2021

*

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

* PURPOSE: This macro includes a proc sgplot to produce Kaplan-Meier, Cumulative Distribution

*          Function (CDF), and Cumulative Incidence Function (CIF) curves with an at-risk table

*                                        

*  Program inputs:                                                                                   

*   - Dataset(s) computed in figure_survivalcurves_createdata.sas (L1 plots) or 

*     l2_effect_estimate_km_createdata.sas (L2 plots)

* 

*  Program outputs: 

*   - Dataset(s) to output/repdata for each figure

* 

* 

*  PARAMETERS:  

*   

*            

*  Programming Notes:         

*  Utility macro %output_survivalcurves created to execute proc sgplot for each km/cdf plot 

*

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

* CONTACT INFO: 

*  Sentinel Coordinating Center

*  info@sentinelsystem.org

*

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

​

%macro figure_survivalcurves_output();

​

    %put =====> MACRO CALLED: figure_survivalcurves_output;

​

    %macro output_survivalcurves(dataset=,

                         where=,

                         figtitle=,

                         figfn=,

                         xaxislabel=,

                         yaxislabel=,

                         figure=,

                         font=,

                         analysis=,

                         analysisgrp=,

                         monitoringperiod=,

                         eoi=,

                         ref=,

                         eoilabel=,

                         reflabel=,

                         kmrefpop=);

​

        /* Obtain x and y axis values */

        %let num_fn = 0;

        %let kmxtickmarks = ;

        %let kmytickmarks = ;

​

        data _null_;

            set figurefile(where=(figure="&figure."));

            call symputx('atrisktable', includeatrisktable);

        run;

​

        /*assign axes*/

        %figure_axes(data=&dataset., figure=&figure., figuresub=overall, where=&where., xtickmarks=kmxtickmarks, xvar=day, ytickmarks=kmytickmarks, yvar=);

​

        %tableletter(); 

        %isdata(dataset=repdata.Figure&figurenum.&tableletter.);

        %if %eval(&nobs=0) %then %do;

        data repdata.Figure&figurenum.&tableletter.;

        set &dataset(where=(&where));

        /* Only selects days for corresponding tickmarks */

        %if &atrisktable = Y %then %do;

        if day in (&kmxtickmarks) then xaxisatrisk=day;

        %end;

        /* Add columns for Sentinel Views */

        %if &reporttype = T2L2 and (&figure=F3 or &figure=F4 or &figure=F5) %then %do;

        length monitoringperiod 3 analysis $13 analysisgrp eoi ref $40;

        %if &analysis = Adjusted %then %do;

            %if &figure = F4 and &pscsfile. = psmatchfile %then %let analysis=Conditional;

            %else %if &figure = F5 and &pscsfile. = psmatchfile %then %let analysis=Unconditional;

            %else %if &pscsfile. = stratificationfile | &pscsfile. = iptwfile %then %let analysis=Weighted; 

        %end;

        analysis="&analysis";

        analysisgrp="&analysisgrp";

        monitoringperiod=&monitoringperiod;

        eoi="&eoi";

        ref="&ref";

        eoilabel="&eoilabel";

        reflabel="&reflabel";

        %end;