NOTE: Copyright (c) 2016 by SAS Institute Inc., Cary, NC, USA.
NOTE: SAS (r) Proprietary Software 9.4 (TS1M7)
      Licensed to COLORADO PREVENTION CENTER, Site 70237720.
NOTE: This session is executing on the X64_DSRV19  platform.



NOTE: Analytical products:

      SAS/STAT 15.2
      SAS/IML 15.2

NOTE: Additional host information:

 X64_DSRV19 WIN 10.0.17763  Server

NOTE: SAS initialization used:
      real time           0.73 seconds
      cpu time            0.59 seconds

1    Options MPRINT;
2    ****************************************************************;
3
4    %macro SURV_CONTNRI(DSNAME,DETAIL,EVENT,PYRS,T,PROB1,PROB2);
5    * Macro to compute NRI with any up or down (no categories);
6    *  Variables:
7    *  DSNAME = dataset name;
8    *  DETAIL = 2 for detailed printout, 1 for limited, 0 for none;
9    *  EVENT = outcome variable (coded 0,1);
10   *  PYS = follow-up time;
11   *  T = time of predicted risk (eg., 10 years) = scalar;
12   *  PROB1 = probability for model 1;
13   *  PROB2 = probability for model 2;
14   *** NOTE: data need to be bootstrapped to get correct SE;
15
16   title2 Continuous NRI for Event &EVENT at Time &T : &PROB2 vs &PROB1 ;
17   title3 "Need to Bootstrap for SE";
18
19   data nri1; set &dsname;
20   _id_=_N_;
21   * compute diffs in probs;
22   event=&EVENT;
23   pyrs=&pyrs;
24   prob1=&PROB1;
25   prob2=&PROB2;
26   diffp=prob2-prob1;
27   if prob1>. and prob2>.;
28   keep _id_ event pyrs prob1 prob2 diffp;
29   run;
30
31   data nri1; set nri1;
32   if prob2>prob1 then disc=1;
33     else if prob2=prob1 then disc=0;
34     else if prob2<prob1 then disc=-1;
35   run;
36
37   *** Control output as needed ***;
38   %if &detail=2 %then %do;
39   proc means data=nri1 n sum mean median stddev min max;
40   var prob1 prob2 diffp disc event;
41   run;
42   %end;
43
44   %if &detail=2 %then %do;
45   proc freq data=nri1;
46   tables event*disc ;
47   title4 "Crude Reclassification Ignoring Time";
48   run;
49   title4; run;
50   %end;
51
52   ods listing exclude all;
53   ods html exclude all;
54   ods output OneWayFreqs=freqs;
55   run;
56   proc freq data=nri1;
57   tables disc ;
58   run;
59   ods listing exclude none;
60   ods html exclude none;
61   run;
62
63   * Compute KM estimator for survival for those moving up or down;
64   data two;
65   _id_=0;
66   event=.; pyrs=&t;
67   * Output for up and down;
68   disc=1; bd=1; output;
69   disc=0; bd=1; output;
70   disc=-1; bd=1; output;
71   * Output for overall;
72   disc=.; bd=.; output;
73   run;
74
75   data _all_; set nri1 two;
76   run;
77
78   proc phreg data=_all_(where=(disc=1 or disc=-1)) noprint;
79      model pyrs*event(0) =  ;
80      id _id_;
81      strata disc;
82      output out=survdat survival=_survest_ ;
83   run;
84
85   proc phreg data=_all_(where=(bd ne 1)) noprint;
86      model pyrs*event(0) =  ;
87      id _id_;
88      output out=survdat0 survival=survall ;
89   run;
90
91   data base; set survdat;
92   if _id_=0;
93   * survival estimate for referent person at time &t;
94   _basesurv_=_survest_;
95   inc=1-_basesurv_;
96   keep disc _basesurv_ inc;
97   run;
98
99   data base0; set survdat0;
100  * overall survival estimate;
101  if _id_=0;
102  * survival estimate for referent person at time &t;
103  inc0=1-survall;
104  disc=.;
105  keep disc survall inc0;
106  run;
107
108  proc sort data=base; by disc;
109  proc sort data=freqs; by disc;
110
111  data calc ; merge freqs base base0; by disc;
112  keep disc _basesurv_ inc survall inc0 Frequency Percent;
113  run;
114
115  data _all_; set calc end=eof;
116  retain n_up n_down totn inc_up inc_down inc_all ecase_up ecase_down
117         econt_up econt_down 0;
118  if disc=1 then do;
119     n_up=Frequency; inc_up=inc; ecase_up=n_up*inc_up; econt_up=n_up-ecase_up;
120     totn=totn+Frequency;
121     end;
122  else if disc=0 then do;
123     n_same=Frequency;
124     totn=totn+Frequency;
125     end;
126  else if disc=-1 then do;
127     n_down=Frequency; inc_down=inc; ecase_down=n_down*inc_down; econt_down=n_down-ecase_down;
128     totn=totn+Frequency;
129     end;
130  else if disc=. then inc_all=inc0;
131
132  if eof then do;
133    ecase=totn*inc_all;
134    econt=totn-ecase;
135    %if ecase=0 or econt=0 %then %do;
136       nri=.; znri=.; p2nri=.;
137       %goto exit;
138    %end;
139    pup_case=ecase_up/ecase;
140    pdown_case=ecase_down/ecase;
141    pup_cont=econt_up/econt;
142    pdown_cont=econt_down/econt;
143    ri_case=pup_case-pdown_case;
144    ri_contr=pdown_cont-pup_cont;
145    nri=ri_case+ri_contr;
146    * corrected variance to reflect binomial;
147    evri_case=pup_case*(1-pup_case)/ecase;
148    evri_contr=pup_cont*(1-pup_cont)/econt;
149    evnri=4*(evri_case+evri_contr);
150    esenri=sqrt(evnri);
151    enricil=nri-1.96*esenri;
152    enriciu=nri+1.96*esenri;
153    ezri_case=ri_case/sqrt(evri_case);
154    ezri_contr=ri_contr/sqrt(evri_contr);
155    eznri=nri/sqrt(evnri);
156    ep2ri_case=2*(1-probnorm(abs(ezri_case)));
157    ep2ri_contr=2*(1-probnorm(abs(ezri_contr)));
158    ep2nri=2*(1-probnorm(abs(eznri)));
159    output;
160    keep totn ecase econt inc_all n_up n_same n_down inc_up inc_down
161         pup_case pdown_case pup_cont pdown_cont
162         ri_case ri_contr nri evri_case evri_contr evnri esenri enricil enriciu
163         ezri_case ezri_contr eznri ep2ri_case ep2ri_contr ep2nri
164         ;
165  end;
166
167  %if &detail>0 %then %do;
168  proc print data=_all_;
169  var    totn ecase econt inc_all n_up n_same n_down inc_up inc_down
170         pup_case pdown_case pup_cont pdown_cont
171         ri_case ri_contr nri evri_case evri_contr evnri esenri enricil enriciu
172         ezri_case ezri_contr eznri ep2ri_case ep2ri_contr ep2nri
173         ;
174  run;
175  %end;
176  title2; run;
177  * Usage:
178  * %surv_contnri(probs,1,outxy,pyrs,10,pdx,pdxy);
179  %exit: %mend SURV_CONTNRI;
180
181  ****************************************************************;
182
183
184  ****************************************************************;
185  proc import datafile="C:\Users\Julian.Chen\Downloads\Voyager WIFI\07082024\newdata1.csv"
186  out=newdata
187  dbms=csv
188  replace;
189  run;

Name "" is not a valid SAS name.
Problems were detected with provided names.  See LOG.
190   /**********************************************************************
191   *   PRODUCT:   SAS
192   *   VERSION:   9.4
193   *   CREATOR:   External File Interface
194   *   DATE:      09JUL24
195   *   DESC:      Generated SAS Datastep Code
196   *   TEMPLATE SOURCE:  (None Specified.)
197   ***********************************************************************/
198      data WORK.NEWDATA    ;
199      %let _EFIERR_ = 0; /* set the ERROR detection macro variable */
200      infile 'C:\Users\Julian.Chen\Downloads\Voyager WIFI\07082024\newdata1.csv' delimiter = ','
200! MISSOVER DSD lrecl=32767 firstobs=2 ;
201         informat VAR1 $4. ;
202         informat time best32. ;
203         informat status best32. ;
204         informat DM1 best32. ;
205         informat DM2 best32. ;
206         format VAR1 $4. ;
207         format time best12. ;
208         format status best12. ;
209         format DM1 best12. ;
210         format DM2 best12. ;
211      input
212                  VAR1  $
213                  time
214                  status
215                  DM1
216                  DM2
217      ;
218      if _ERROR_ then call symputx('_EFIERR_',1);  /* set ERROR detection macro variable */
219      run;

NOTE: The infile 'C:\Users\Julian.Chen\Downloads\Voyager WIFI\07082024\newdata1.csv' is:
      Filename=C:\Users\Julian.Chen\Downloads\Voyager WIFI\07082024\newdata1.csv,
      RECFM=V,LRECL=32767,File Size (bytes)=195596,
      Last Modified=08Jul2024:12:16:00,
      Create Time=08Jul2024:12:16:13

NOTE: 6347 records were read from the infile 'C:\Users\Julian.Chen\Downloads\Voyager
      WIFI\07082024\newdata1.csv'.
      The minimum record length was 13.
      The maximum record length was 31.
NOTE: The data set WORK.NEWDATA has 6347 observations and 5 variables.
NOTE: DATA statement used (Total process time):
      real time           0.02 seconds
      cpu time            0.00 seconds


6347 rows created in WORK.NEWDATA from C:\Users\Julian.Chen\Downloads\Voyager
WIFI\07082024\newdata1.csv.



NOTE: WORK.NEWDATA data set was successfully created.
NOTE: The data set WORK.NEWDATA has 6347 observations and 5 variables.
NOTE: PROCEDURE IMPORT used (Total process time):
      real time           0.07 seconds
      cpu time            0.04 seconds


220  proc print data=newdata(obs=10); run;
NOTE: Writing HTML Body file: sashtml.htm

NOTE: There were 10 observations read from the data set WORK.NEWDATA.
NOTE: PROCEDURE PRINT used (Total process time):
      real time           0.32 seconds
      cpu time            0.20 seconds


221  proc phreg data=newdata;
222  class DM1(ref="0");
223  model time*status(0)=DM1;
224  output out=surv1 survival=s1;
225  run;

NOTE: Convergence criterion (GCONV=1E-8) satisfied.
NOTE: The data set WORK.SURV1 has 6347 observations and 6 variables.
NOTE: PROCEDURE PHREG used (Total process time):
      real time           0.11 seconds
      cpu time            0.04 seconds


226  proc phreg data=newdata;
227  class DM1(ref="0") DM2(ref="0");
228  model time*status(0)=DM1 DM2;
229  output out=surv2 survival=s2;
230  run;

NOTE: Convergence criterion (GCONV=1E-8) satisfied.
NOTE: The data set WORK.SURV2 has 6347 observations and 6 variables.
NOTE: PROCEDURE PHREG used (Total process time):
      real time           0.10 seconds
      cpu time            0.04 seconds


MPRINT(SURV_CONTNRI):   * Macro to compute NRI with any up or down (no categories);
MPRINT(SURV_CONTNRI):   * Variables: * DSNAME = dataset name;
MPRINT(SURV_CONTNRI):   * DETAIL = 2 for detailed printout, 1 for limited, 0 for none;
MPRINT(SURV_CONTNRI):   * EVENT = outcome variable (coded 0,1);
MPRINT(SURV_CONTNRI):   * PYS = follow-up time;
231
MPRINT(SURV_CONTNRI):   * T = time of predicted risk (eg., 10 years) = scalar;
MPRINT(SURV_CONTNRI):   * PROB1 = probability for model 1;
MPRINT(SURV_CONTNRI):   * PROB2 = probability for model 2;
MPRINT(SURV_CONTNRI):   *** NOTE: data need to be bootstrapped to get correct SE;
232  %SURV_CONTNRI(newdata,2,status,time,4,surv1,surv2);
MPRINT(SURV_CONTNRI):   title2 Continuous NRI for Event status at Time 4 : surv2 vs surv1 ;
MPRINT(SURV_CONTNRI):   title3 "Need to Bootstrap for SE";
MPRINT(SURV_CONTNRI):   data nri1;
MPRINT(SURV_CONTNRI):   set newdata;
MPRINT(SURV_CONTNRI):   _id_=_N_;
MPRINT(SURV_CONTNRI):   * compute diffs in probs;
MPRINT(SURV_CONTNRI):   event=status;
MPRINT(SURV_CONTNRI):   pyrs=time;
MPRINT(SURV_CONTNRI):   prob1=surv1;
MPRINT(SURV_CONTNRI):   prob2=surv2;
MPRINT(SURV_CONTNRI):   diffp=prob2-prob1;
MPRINT(SURV_CONTNRI):   if prob1>. and prob2>.;
MPRINT(SURV_CONTNRI):   keep _id_ event pyrs prob1 prob2 diffp;
MPRINT(SURV_CONTNRI):   run;

NOTE: Variable surv1 is uninitialized.
NOTE: Variable surv2 is uninitialized.
NOTE: Missing values were generated as a result of performing an operation on missing values.
      Each place is given by: (Number of times) at (Line):(Column).
      6347 at 4:156
NOTE: There were 6347 observations read from the data set WORK.NEWDATA.
NOTE: The data set WORK.NRI1 has 0 observations and 6 variables.
NOTE: DATA statement used (Total process time):
      real time           0.01 seconds
      cpu time            0.00 seconds


MPRINT(SURV_CONTNRI):   data nri1;
MPRINT(SURV_CONTNRI):   set nri1;
MPRINT(SURV_CONTNRI):   if prob2>prob1 then disc=1;
MPRINT(SURV_CONTNRI):   else if prob2=prob1 then disc=0;
MPRINT(SURV_CONTNRI):   else if prob2<prob1 then disc=-1;
MPRINT(SURV_CONTNRI):   run;

NOTE: There were 0 observations read from the data set WORK.NRI1.
NOTE: The data set WORK.NRI1 has 0 observations and 7 variables.
NOTE: DATA statement used (Total process time):
      real time           0.01 seconds
      cpu time            0.00 seconds


MPRINT(SURV_CONTNRI):   *** Control output as needed ***;
MPRINT(SURV_CONTNRI):   proc means data=nri1 n sum mean median stddev min max;
MPRINT(SURV_CONTNRI):   var prob1 prob2 diffp disc event;
MPRINT(SURV_CONTNRI):   run;

NOTE: No observations in data set WORK.NRI1.
NOTE: PROCEDURE MEANS used (Total process time):
      real time           0.00 seconds
      cpu time            0.00 seconds


MPRINT(SURV_CONTNRI):   proc freq data=nri1;
MPRINT(SURV_CONTNRI):   tables event*disc ;
MPRINT(SURV_CONTNRI):   title4 "Crude Reclassification Ignoring Time";
MPRINT(SURV_CONTNRI):   run;

NOTE: No observations in data set WORK.NRI1.
NOTE: PROCEDURE FREQ used (Total process time):
      real time           0.00 seconds
      cpu time            0.00 seconds


MPRINT(SURV_CONTNRI):   title4;
MPRINT(SURV_CONTNRI):   run;
ERROR: The LISTING destination is not active; no select/exclude lists are available.
MPRINT(SURV_CONTNRI):   ods listing exclude all;
MPRINT(SURV_CONTNRI):   ods html exclude all;
MPRINT(SURV_CONTNRI):   ods output OneWayFreqs=freqs;
MPRINT(SURV_CONTNRI):   run;
MPRINT(SURV_CONTNRI):   proc freq data=nri1;
MPRINT(SURV_CONTNRI):   tables disc ;
MPRINT(SURV_CONTNRI):   run;

NOTE: No observations in data set WORK.NRI1.
NOTE: PROCEDURE FREQ used (Total process time):
      real time           0.00 seconds
      cpu time            0.00 seconds

WARNING: Output 'OneWayFreqs' was not created.  Make sure that the output object name, label, or path
         is spelled correctly.  Also, verify that the appropriate procedure options are used to
         produce the requested output object.  For example, verify that the NOPRINT option is not
         used.

ERROR: The LISTING destination is not active; no select/exclude lists are available.
MPRINT(SURV_CONTNRI):   ods listing exclude none;
MPRINT(SURV_CONTNRI):   ods html exclude none;
MPRINT(SURV_CONTNRI):   run;
MPRINT(SURV_CONTNRI):   * Compute KM estimator for survival for those moving up or down;
MPRINT(SURV_CONTNRI):   data two;
MPRINT(SURV_CONTNRI):   _id_=0;
MPRINT(SURV_CONTNRI):   event=.;
MPRINT(SURV_CONTNRI):   pyrs=4;
MPRINT(SURV_CONTNRI):   * Output for up and down;
MPRINT(SURV_CONTNRI):   disc=1;
MPRINT(SURV_CONTNRI):   bd=1;
MPRINT(SURV_CONTNRI):   output;
MPRINT(SURV_CONTNRI):   disc=0;
MPRINT(SURV_CONTNRI):   bd=1;
MPRINT(SURV_CONTNRI):   output;
MPRINT(SURV_CONTNRI):   disc=-1;
MPRINT(SURV_CONTNRI):   bd=1;
MPRINT(SURV_CONTNRI):   output;
MPRINT(SURV_CONTNRI):   * Output for overall;
MPRINT(SURV_CONTNRI):   disc=.;
MPRINT(SURV_CONTNRI):   bd=.;
MPRINT(SURV_CONTNRI):   output;
MPRINT(SURV_CONTNRI):   run;

NOTE: The data set WORK.TWO has 4 observations and 5 variables.
NOTE: DATA statement used (Total process time):
      real time           0.01 seconds
      cpu time            0.01 seconds


MPRINT(SURV_CONTNRI):   data _all_;
MPRINT(SURV_CONTNRI):   set nri1 two;
MPRINT(SURV_CONTNRI):   run;
NOTE: There were 0 observations read from the data set WORK.NRI1.
NOTE: There were 4 observations read from the data set WORK.TWO.
NOTE: The data set WORK._ALL_ has 4 observations and 8 variables.
NOTE: DATA statement used (Total process time):
      real time           0.01 seconds
      cpu time            0.00 seconds


MPRINT(SURV_CONTNRI):   proc phreg data=_all_(where=(disc=1 or disc=-1)) noprint;
MPRINT(SURV_CONTNRI):   model pyrs*event(0) = ;
MPRINT(SURV_CONTNRI):   id _id_;
MPRINT(SURV_CONTNRI):   strata disc;
MPRINT(SURV_CONTNRI):   output out=survdat survival=_survest_ ;
MPRINT(SURV_CONTNRI):   run;

NOTE: There are no explanatory variables in the MODEL statement.
NOTE: 1 observations were deleted in the 1st stratum due either to missing or invalid values for the
      time, censoring, frequency or explanatory variables or to invalid operations in generating the
      values for some of the explanatory variables.
NOTE: 1 observations were deleted in the 2nd stratum due either to missing or invalid values for the
      time, censoring, frequency or explanatory variables or to invalid operations in generating the
      values for some of the explanatory variables.
ERROR: There are no valid observations for the analysis.
NOTE: The SAS System stopped processing this step because of errors.
WARNING: The data set WORK.SURVDAT may be incomplete.  When this step was stopped there were 2
         observations and 9 variables.
NOTE: PROCEDURE PHREG used (Total process time):
      real time           0.01 seconds
      cpu time            0.01 seconds


MPRINT(SURV_CONTNRI):   proc phreg data=_all_(where=(bd ne 1)) noprint;
MPRINT(SURV_CONTNRI):   model pyrs*event(0) = ;
MPRINT(SURV_CONTNRI):   id _id_;
MPRINT(SURV_CONTNRI):   output out=survdat0 survival=survall ;
MPRINT(SURV_CONTNRI):   run;

NOTE: There are no explanatory variables in the MODEL statement.
NOTE: 1 observations were deleted due either to missing or invalid values for the time, censoring,
      frequency or explanatory variables or to invalid operations in generating the values for some
      of the explanatory variables.
ERROR: There are no valid observations for the analysis.
NOTE: The SAS System stopped processing this step because of errors.
WARNING: The data set WORK.SURVDAT0 may be incomplete.  When this step was stopped there were 1
         observations and 9 variables.
NOTE: PROCEDURE PHREG used (Total process time):
      real time           0.01 seconds
      cpu time            0.01 seconds



MPRINT(SURV_CONTNRI):   data base;
MPRINT(SURV_CONTNRI):   set survdat;
MPRINT(SURV_CONTNRI):   if _id_=0;
MPRINT(SURV_CONTNRI):   * survival estimate for referent person at time &t;
MPRINT(SURV_CONTNRI):   _basesurv_=_survest_;
MPRINT(SURV_CONTNRI):   inc=1-_basesurv_;
MPRINT(SURV_CONTNRI):   keep disc _basesurv_ inc;
MPRINT(SURV_CONTNRI):   run;
NOTE: Missing values were generated as a result of performing an operation on missing values.
      Each place is given by: (Number of times) at (Line):(Column).
      2 at 337:111
NOTE: There were 2 observations read from the data set WORK.SURVDAT.
NOTE: The data set WORK.BASE has 2 observations and 3 variables.
NOTE: DATA statement used (Total process time):
      real time           0.00 seconds
      cpu time            0.00 seconds


MPRINT(SURV_CONTNRI):   data base0;
MPRINT(SURV_CONTNRI):   set survdat0;
MPRINT(SURV_CONTNRI):   * overall survival estimate;
MPRINT(SURV_CONTNRI):   if _id_=0;
MPRINT(SURV_CONTNRI):   * survival estimate for referent person at time &t;
MPRINT(SURV_CONTNRI):   inc0=1-survall;
MPRINT(SURV_CONTNRI):   disc=.;
MPRINT(SURV_CONTNRI):   keep disc survall inc0;
MPRINT(SURV_CONTNRI):   run;

NOTE: Missing values were generated as a result of performing an operation on missing values.
      Each place is given by: (Number of times) at (Line):(Column).
      1 at 17:27
NOTE: There were 1 observations read from the data set WORK.SURVDAT0.
NOTE: The data set WORK.BASE0 has 1 observations and 3 variables.
NOTE: DATA statement used (Total process time):
      real time           0.00 seconds
      cpu time            0.01 seconds


MPRINT(SURV_CONTNRI):   proc sort data=base;
MPRINT(SURV_CONTNRI):   by disc;

NOTE: There were 2 observations read from the data set WORK.BASE.
NOTE: The data set WORK.BASE has 2 observations and 3 variables.
NOTE: PROCEDURE SORT used (Total process time):
      real time           0.00 seconds
      cpu time            0.01 seconds


MPRINT(SURV_CONTNRI):   proc sort data=freqs;
ERROR: File WORK.FREQS.DATA does not exist.
MPRINT(SURV_CONTNRI):   by disc;

NOTE: The SAS System stopped processing this step because of errors.
NOTE: PROCEDURE SORT used (Total process time):
      real time           0.00 seconds
      cpu time            0.00 seconds



MPRINT(SURV_CONTNRI):   data calc ;
MPRINT(SURV_CONTNRI):   merge freqs base base0;
ERROR: File WORK.FREQS.DATA does not exist.
MPRINT(SURV_CONTNRI):   by disc;
MPRINT(SURV_CONTNRI):   keep disc _basesurv_ inc survall inc0 Frequency Percent;
MPRINT(SURV_CONTNRI):   run;

WARNING: The variable Frequency in the DROP, KEEP, or RENAME list has never been referenced.
WARNING: The variable Percent in the DROP, KEEP, or RENAME list has never been referenced.
NOTE: The SAS System stopped processing this step because of errors.
WARNING: The data set WORK.CALC may be incomplete.  When this step was stopped there were 0
         observations and 5 variables.
NOTE: DATA statement used (Total process time):
      real time           0.00 seconds
      cpu time            0.01 seconds



MPRINT(SURV_CONTNRI):   data _all_;
MPRINT(SURV_CONTNRI):   set calc end=eof;
MPRINT(SURV_CONTNRI):   retain n_up n_down totn inc_up inc_down inc_all ecase_up ecase_down econt_up
econt_down 0;
MPRINT(SURV_CONTNRI):   if disc=1 then do;
MPRINT(SURV_CONTNRI):   n_up=Frequency;
MPRINT(SURV_CONTNRI):   inc_up=inc;
MPRINT(SURV_CONTNRI):   ecase_up=n_up*inc_up;
MPRINT(SURV_CONTNRI):   econt_up=n_up-ecase_up;
MPRINT(SURV_CONTNRI):   totn=totn+Frequency;
MPRINT(SURV_CONTNRI):   end;
MPRINT(SURV_CONTNRI):   else if disc=0 then do;
MPRINT(SURV_CONTNRI):   n_same=Frequency;
MPRINT(SURV_CONTNRI):   totn=totn+Frequency;
MPRINT(SURV_CONTNRI):   end;
MPRINT(SURV_CONTNRI):   else if disc=-1 then do;
MPRINT(SURV_CONTNRI):   n_down=Frequency;
MPRINT(SURV_CONTNRI):   inc_down=inc;
MPRINT(SURV_CONTNRI):   ecase_down=n_down*inc_down;
MPRINT(SURV_CONTNRI):   econt_down=n_down-ecase_down;
MPRINT(SURV_CONTNRI):   totn=totn+Frequency;
MPRINT(SURV_CONTNRI):   end;
MPRINT(SURV_CONTNRI):   else if disc=. then inc_all=inc0;
MPRINT(SURV_CONTNRI):   if eof then do;
MPRINT(SURV_CONTNRI):   ecase=totn*inc_all;
MPRINT(SURV_CONTNRI):   econt=totn-ecase;
MPRINT(SURV_CONTNRI):   pup_case=ecase_up/ecase;
MPRINT(SURV_CONTNRI):   pdown_case=ecase_down/ecase;
MPRINT(SURV_CONTNRI):   pup_cont=econt_up/econt;
MPRINT(SURV_CONTNRI):   pdown_cont=econt_down/econt;
MPRINT(SURV_CONTNRI):   ri_case=pup_case-pdown_case;
MPRINT(SURV_CONTNRI):   ri_contr=pdown_cont-pup_cont;
MPRINT(SURV_CONTNRI):   nri=ri_case+ri_contr;
MPRINT(SURV_CONTNRI):   * corrected variance to reflect binomial;
MPRINT(SURV_CONTNRI):   evri_case=pup_case*(1-pup_case)/ecase;
MPRINT(SURV_CONTNRI):   evri_contr=pup_cont*(1-pup_cont)/econt;
MPRINT(SURV_CONTNRI):   evnri=4*(evri_case+evri_contr);
MPRINT(SURV_CONTNRI):   esenri=sqrt(evnri);
MPRINT(SURV_CONTNRI):   enricil=nri-1.96*esenri;
MPRINT(SURV_CONTNRI):   enriciu=nri+1.96*esenri;
MPRINT(SURV_CONTNRI):   ezri_case=ri_case/sqrt(evri_case);
MPRINT(SURV_CONTNRI):   ezri_contr=ri_contr/sqrt(evri_contr);
MPRINT(SURV_CONTNRI):   eznri=nri/sqrt(evnri);
MPRINT(SURV_CONTNRI):   ep2ri_case=2*(1-probnorm(abs(ezri_case)));
MPRINT(SURV_CONTNRI):   ep2ri_contr=2*(1-probnorm(abs(ezri_contr)));
MPRINT(SURV_CONTNRI):   ep2nri=2*(1-probnorm(abs(eznri)));
MPRINT(SURV_CONTNRI):   output;
MPRINT(SURV_CONTNRI):   keep totn ecase econt inc_all n_up n_same n_down inc_up inc_down pup_case
pdown_case pup_cont pdown_cont ri_case ri_contr nri evri_case evri_contr evnri esenri enricil enriciu
ezri_case ezri_contr eznri ep2ri_case ep2ri_contr ep2nri ;
MPRINT(SURV_CONTNRI):   end;

NOTE: Variable Frequency is uninitialized.
NOTE: There were 0 observations read from the data set WORK.CALC.
NOTE: The data set WORK._ALL_ has 0 observations and 28 variables.
NOTE: DATA statement used (Total process time):
      real time           0.02 seconds
      cpu time            0.01 seconds


MPRINT(SURV_CONTNRI):   proc print data=_all_;
MPRINT(SURV_CONTNRI):   var totn ecase econt inc_all n_up n_same n_down inc_up inc_down pup_case
pdown_case pup_cont pdown_cont ri_case ri_contr nri evri_case evri_contr evnri esenri enricil enriciu
ezri_case ezri_contr eznri ep2ri_case ep2ri_contr ep2nri ;
MPRINT(SURV_CONTNRI):   run;

NOTE: No observations in data set WORK._ALL_.
NOTE: PROCEDURE PRINT used (Total process time):
      real time           0.00 seconds
      cpu time            0.00 seconds


MPRINT(SURV_CONTNRI):   title2;
MPRINT(SURV_CONTNRI):   run;
MPRINT(SURV_CONTNRI):   * Usage: * %surv_contnri(probs,1,outxy,pyrs,10,pdx,pdxy);