NOTE: Copyright (c) 1999-2001 by SAS Institute Inc., Cary, NC, USA.
NOTE: SAS (r) Proprietary Software Release 8.2 (TS2M0)
      Licensed to GEORGE WASHINGTON UNIVERSITY-CAMPUSWIDE-T/R, Site 0001180006.
NOTE: This session is executing on the WIN_ME  platform.
 
 
 
NOTE: SAS initialization used:
      real time           15.04 seconds
 
1    ****  ST210A_F SAS;
2    filename sa_flst 'c:\sasjobs\stat210\jl2003\sa_flst.pdf';
3    ods printer pdf file=sa_flst;
NOTE: Writing ODS PRINTER output to DISK destination "c:\sasjobs\stat210\jl2003\sa_flst.pdf",
      printer "PDF".
4
5    TITLE1 'ANALYSIS OF MULTIPLE 2 X 2 TABLES';
6
7    %MACRO FREQANAL;
8    options nonotes;
9    DATA ONE; SET ONE;
10   STRATUM = '        ';STRATUM = K;
11
12   PROC MEANS SUM NOPRINT;
13      VAR D1 N1 D2 N2;
14      OUTPUT OUT=ONESUM
15      SUM= D1 N1 D2 N2;
16   DATA ONESUM; SET ONESUM; STRATUM = 'MARGINAL';
17
18   DATA TWO; SET ONE ONESUM;
19
20   *** ESTIMATORS WITHIN STRATA;
21
22   A=D1; B=D2; C=N1-D1; D=N2-D2;
23   N=N1+N2; DT=D1+D2;  P=DT/N;
24   M=1/(1/N1 + 1/N2);
25   P1=D1/N1;P2=D2/N2; RR=P1/P2; LRR=LOG(RR);
26   RD=P1-P2;  OR=A*D/(B*C); LOR=LOG(OR);
27
28   VLOR = (1/(N1*P1*(1-P1))) + (1/(N2*P2*(1-P2)));
29   VLRR = ((1-P1)/(N1*P1)) + ((1-P2)/(N2*P2));
30   VRD = (P1*(1-P1)/N1) + (P2*(1-P2)/N2);
31   U95RD = RD + 1.96*SQRT(VRD);
32   L95RD = RD - 1.96*SQRT(VRD);
33   U95LRR = LRR + 1.96*SQRT(VLRR);
34   L95LRR = LRR - 1.96*SQRT(VLRR);
35   U95RR = EXP(U95LRR);
36   L95RR = EXP(L95LRR);
37   U95LOR = LOR + 1.96*SQRT(VLOR);
38   L95LOR = LOR - 1.96*SQRT(VLOR);
39   U95OR = EXP(U95LOR);
40   L95OR = EXP(L95LOR);
41
42   PROC PRINT; VAR STRATUM A B C D N1 N2 P1 P2 P;
43   PROC PRINT; VAR STRATUM
44    RD VRD L95RD U95RD;
45   PROC PRINT; VAR STRATUM
46    RR LRR VLRR
47    L95LRR U95LRR L95RR U95RR;
48   PROC PRINT; VAR STRATUM
49    OR LOR VLOR
50    L95LOR U95LOR L95OR U95OR
51    ;
52
53   ************ Generate data set for use with proc FREQ **************;
54
55   DATA FOUR; SET ONE;
56   KEEP I J K F;
57   D3=N1-D1;
58   D4=N2-D2;
59   I = 1; J = 1; F =D1; OUTPUT;
60   I = 2; J = 1; F =D2; OUTPUT;
61   I = 1; J = 2; F =D3; OUTPUT;
62   I = 2; J = 2; F =D4; OUTPUT;
63
64   %LABELIT;
65   options notes;
66   PROC FREQ; TABLE I*J / CHISQ CMH NOCOL NOPERCENT;
67   WEIGHT F;
68   format i group. j response.;
69   TITLE3 'ESTIMATES AND TESTS WITHIN COMBINED SAMPLE -- '
70   'MARGINAL, UNADJUSTED ANALYSES FOR GROUP X RESPONSE (I X J)';
71
72   PROC FREQ; TABLE K*J / CHISQ NOCOL NOPERCENT;
73   WEIGHT F;
74   format j response. k stratum.;
75   TITLE3 'ESTIMATES AND TESTS WITHIN COMBINED SAMPLE -- '
76   'MARGINAL, UNADJUSTED ANALYSES FOR STRATUM X RESPONSE (K X J)';
77
78   PROC FREQ; TABLE K*I / CHISQ NOCOL NOPERCENT;
79   WEIGHT F;
80   format i group. k stratum.;
81   TITLE3 'ESTIMATES AND TESTS WITHIN COMBINED SAMPLE -- '
82   'MARGINAL, UNADJUSTED ANALYSES FOR STRATUM X GROUP (K X I)';
83
84   PROC FREQ; TABLE I*J / CHISQ CMH NOCOL NOPERCENT; BY K;
85   WEIGHT F;
86   format i group. j response. k stratum.;
87   TITLE3 'ESTIMATES AND TESTS WITHIN STRATA --'
88   ' CONDITIONAL NESTED ANALYSES OF GROUP X RESPONSE (I X J) WITHIN STRATA (K)';
 
89
90   PROC FREQ; TABLE K*I*J / CMH NOPRINT;
91   WEIGHT F;
92   format i group. j response. k stratum.;
93   TITLE3 'PARTIAL, ADJUSTED ESTIMATES AND TESTS OF GROUP X RESPONSE COMBINED OVER STRATA';
 
94   %MEND;
95
96   DATA ONE;
97   INPUT K N1 D1 N2 D2;
98   CARDS;
 
NOTE: The data set WORK.ONE has 3 observations and 5 variables.
NOTE: DATA statement used:
      real time           0.71 seconds
 
 
102  *****;
103  TITLE2 'COUNTRY A DATA';
104  %MACRO LABELIT;
105  label i='group (independent variable)'
106        j='response (dependent variable)'
107        k='stratum (confounding variable)';
108  proc format;
109    value group 1='drug' 2='placebo';
110    value stratum 1='acid dependent' 2='drug dependent' 3='intermediate';
111    value response 1='ulcer healed' 2='not healed';
112  %MEND;
113  %FREQANAL;
NOTE: PROCEDURE FORMAT used:
      real time           0.22 seconds
 
 
 
NOTE: There were 12 observations read from the data set WORK.FOUR.
NOTE: PROCEDURE FREQ used:
      real time           0.76 seconds
 
 
 
NOTE: There were 12 observations read from the data set WORK.FOUR.
NOTE: PROCEDURE FREQ used:
      real time           0.10 seconds
 
 
 
NOTE: There were 12 observations read from the data set WORK.FOUR.
NOTE: PROCEDURE FREQ used:
      real time           0.05 seconds
 
 
 
NOTE: There were 12 observations read from the data set WORK.FOUR.
NOTE: PROCEDURE FREQ used:
      real time           0.54 seconds
 
 
114
115  /****************************************************************;
116  DATA ONE;
117  INPUT K N1 D1 N2 D2;
118  CARDS;
119   1 18  8 15 4
120   2 27 21 18 9
121  *****;
122  %MACRO LABELIT;
123  label i='group (independent variable)'
124        j='response (dependent variable)'
125        k='stratum (confounding variable)';
126  proc format;
127    value group 1='ST depress' 2='none';
128    value stratum 1='female' 2='male';
129    value response 1='CAD yes' 2='no';
130  %mend;
131  TITLE2 'Koch, Imrey et al. data from Stokes, Davis, Koch, 1995';
132  title3 'Stratum 1=f, 2=M. Group 1= ECG ST segment depression, 2=none';
133  %FREQANAL;
134  */;
135
136  ods printer close;
NOTE: ODS PRINTER printed 17 pages to c:\sasjobs\stat210\jl2003\sa_flst.pdf.
137  RUN;
 
NOTE: There were 12 observations read from the data set WORK.FOUR.
NOTE: PROCEDURE FREQ used:
      real time           0.93 seconds
 
 
NOTE: SCL source line.
138  
     -
     180
 
ERROR 180-322: Statement is not valid or it is used out of proper order.