qrp_report

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

*                                         SENTINEL PROGRAM

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

* PROGRAM:  aggregate_tree.sas 

* Created (mm/dd/yy): 04/28/2021

*

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

* PURPOSE: Create analytic datasets that can be used as inputs to TreeScan software

* 

* Program inputs:       

*    input.&tree_table.

*    Type 3 Analysis:

*       msoc.[RUNID]_t#_treeads_[TREEANALYSISID]_[LEVELID]_[LEVELNUM]_[PERIODID].sas7bdat 

*       msoc.[RUNID]_t3_tree_wkdays_[PERIODID]_agg.sas7bdat

*    Type 2 and Type 4:

*       msoc.[RUNID]_t#_treeads_[TREEANALYSISID]_[LEVELID]_[LEVELNUM]_[PERIODID].sas7bdat

*

*  PARAMETERS:  

*          

*  Programming Notes: 

*

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

* CONTACT INFO: 

*  Sentinel Coordinating Center

*  info@sentinelsystem.org

*

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

​

%macro aggregate_tree();

​

  %do n = 1 %to &numrunid.;

   %let runid = %scan(&runidlist., &n.); 

   %if %index(&reporttype., TREE) > 0 %then %do;

     %let type = %substr(&reporttype,5,1);

   %end;

   %else %do;

     %let type = %substr(&reporttype.,2,1);

   %end;

   /***********************************************************************************************

    Identify EOI and REF for type 2 and 4

   ***********************************************************************************************/

   *Get EOI/REF for type 2/4 analysis;

    %if %eval(&type. ne 3) %then %do;

        proc sql noprint;

            create table comparison as

            select x.analysisgrp, y.eoi, y.ref

            from infolder.&&&runid._psmatchfile. (keep=analysisgrp psestimategrp) as x

            inner join infolder.&&&runid._psestimationfile. as y

            on x.psestimategrp = y.psestimategrp;

        quit;

    %end;

  /***********************************************************************************************

    Determine if there is a t3treewkdays dataset created

   ***********************************************************************************************/

   %let run_t3wk = N;

   %if %sysfunc(exist(agg_t3_tree_wkdays_&periodid)) =1 %then %do;

      proc sql;

       select case when count(levelid) > 0 then 'Y'

              else 'N' end into: run_t3wk

       from infolder.&&&runid._userstrata. (where = (lowcase(tableid) = 't3treewkdays'));

      quit;