qrp_report
Public
Source
xxxxxxxxxx keep runid group: order day lag_episodes_atrisk %if "&curve" = "CIF" %then %do; &includevars. &eventvar. %end; %else %do; &curve._: %end;;***************************************************************************************************** PROGRAM OVERVIEW****************************************************************************************************** PROGRAM: figure_survivalcurves_createdata.sas * Created (mm/dd/yyyy): 07/14/2021**--------------------------------------------------------------------------------------------------* PURPOSE: This macro computes Kaplan-Meier (KM), Cumulative Distribution Function (CDF), and* Cumulative Incidence Function (CIF) estimates* * Program inputs: * - MSOC dataset containing 1 row per day and columns containing counts of episodes censored* * Program outputs: * * * PARAMETERS: * - dataset: aggregate dataset from %aggregate_report_tables* - rename: syntax to rename variables on input dataset* - curve: KM, CDF or CIF* - whereclause: where statement to restrict &dataset* - dayvar: variable name for censor day variable* - includegroups: Groups to include in report* - includevars: For KM and CDF curves, censor reasons to include in final dataset. * For CIF curves, censor reason to use as competing risk* - eventvar: For CIF curves, censor reason to use as event of interest* - transposedata: Y/N/T whether all groups are in 1 plot (Y=Yes, N=No, T=only execute transpose)* - discardnegativetimegroups: list of groups to remove negative ttswitch (Relevant for T6 only)* - figure: standard figure # from FIGUREFILE* * Programming Notes: * - Curve = CDF will produce a 1-CDF plot* - Requires group names to be unique across runIDs* - Can specify transposedata = T if prior dataset exists with stacked data - used to produce * 2nd plot with the same data **--------------------------------------------------------------------------------------------------* CONTACT INFO: * Sentinel Coordinating Center* info@sentinelsystem.org****************************************************************************************************;%macro figure_survivalcurves_createdata(dataset=, rename=, curve=, whereclause=, dayvar=, includegroups=, includevars=, eventvar=, transposedata=, discardnegativetimegroups=, figure=); %put =====> MACRO CALLED: figure_survivalcurves_createdata; /*--------------------------------------------------------------------------------------------*/ /* Confirm input dataset exists */ /*--------------------------------------------------------------------------------------------*/ %isdata(dataset=&dataset.); %if %eval(&nobs.<1) %then %do; %put WARNING: SENTINEL: plot &figure. requested, however dataset &dataset. does not exist; %goto skipfigure; %end; /*--------------------------------------------------------------------------------------------*/ /* If transposedata T then do no execute aggregation and computation */ /*--------------------------------------------------------------------------------------------*/ %if &transposedata. ne T %then %do; /*--------------------------------------------------------------------------------------------*/ /* Select rows and columns and aggregate data */ /*--------------------------------------------------------------------------------------------*/ proc means data=&dataset.(&rename. where=(&whereclause.)) nway noprint; var episodes &includevars. %if "&curve" = "CIF" %then %do; &eventvar. %end;; class runid group &dayvar. / missing; output out=_survivaldata(drop=_: where=(missing(day)=0) rename=&dayvar.=day) sum=; run; %isdata(dataset=_survivaldata); %if %eval(&nobs.<1) %then %do; %put WARNING: SENTINEL: plot &figure. requested, however no observations available to compute curve; %goto skipfigure; %end; /*--------------------------------------------------------------------------------------------*/ /* Remove negative time for Type 6 */ /*--------------------------------------------------------------------------------------------*/ %let minday = 0;