1                              The SAS System   23:57 Friday, August 5, 1994

NOTE: Copyright(c) 1989 by SAS Institute Inc., Cary, NC USA. 
NOTE: SAS (r) Proprietary Software Release 6.08  TS405
      Licensed to UNIVERSITY OF NORTH CAROLINA AT CHAPEL HILL, Site 
      0003944001.



1          TITLE1 "\KEM\POWER\TWOBYTWO\P0202.PGM";
2          TITLE2 "POOLED VAR T INDEP BINOMIAL, 2X2 TABLE";
3          LIBNAME LOAD01 "C:\PROGRAMS\POWER\STORED\" ;
NOTE: Libref LOAD01 was successfully assigned as follows: 
      Engine:        V608 
      Physical Name: C:\PROGRAMS\POWER\STORED
4          
5          PROC IML WORKSIZE=200 SYMSIZE=400;
Worksize    = 204800
Symbol size = 409600
IML Ready
6          RESET STORAGE=LOAD01.POWSTORE;
NOTE: New storage library = LOAD01.POWSTORE
7          LOAD  MODULE = _ALL_           ;
8          
9          START SHELL;
9                        * JUST TO ALLOW DO LOOP AROUND LINK;
10         
11         ROUND={3};
12         U = I(1);
13         C = {1 -1};
14         OPT_ON={COLLAPSE NOPRINT };
15         OPT_OFF = {ALPHA WARN CASE RHOSCAL SIGSCAL BETASCAL};
16         BETA=J(2,1,0);
17         
18         NSET={50}#(1:20)`;
19         NSET={25}//NSET;
20         PISET={.08 .04 ,
21                .05 .02 } ;
22         RN=NROW(NSET);
23         INHOLD= ( PISET @ J(RN,1,1) ) || ( J(NROW(PISET),2,1) @ NSET ) ;
24         
25         NCASES=NROW(INHOLD);
26         
27         ALPHA={.05};
28         REPN={1};
29           DO ICASE=1 TO NCASES;
30         
31           PI1=INHOLD[ICASE,1];
32           PI2=INHOLD[ICASE,2];
33         
34           BETA[1,1]=PI1;
35           BETA[2,1]=PI2;
2                              The SAS System   23:57 Friday, August 5, 1994

36         
37           SIGMA1=PI1#(1.-PI1);
38           SIGMA2=PI2#(1.-PI2);
39           SIGMA=(SIGMA1+SIGMA2)/2;
40         
41           N1=INHOLD[ICASE,3];
42           N2=INHOLD[ICASE,4];
43           ESSENCEX = BLOCK( J(N1,1,1) , J(N2,1,1)  );
44         
45           RUN POWER;
46         
47           POWCASES=NROW(_HOLDPOW);
48           TEMP=J(POWCASES,1,1) @ INHOLD[ICASE,*] ;
48                                                     *DUPLICATE INFO;
49           HOLD=HOLD // ( TEMP || _HOLDPOW ) ;
50         
51           END;
52         HOLDNM={"PI1" "PI2" "N1" "N2" "TOTAL_N" "POWER"};
53         CREATE ONE VAR HOLDNM;
54         APPEND FROM HOLD;
55         CLOSE ONE;
56         FINISH;
NOTE: Module SHELL defined.
56                 * SHELL;
57         RUN  SHELL;
NOTE: The data set WORK.ONE has 42 observations and 6 variables.
58         
Exiting IML.
NOTE: Storage library LOAD01.POWSTORE closed.
NOTE: 12 workspace compresses.
NOTE: The PROCEDURE IML used 15.67 seconds.


59         PROC PRINT DATA=ONE NOOBS;
60         BY PI1 PI2 NOTSORTED;
61         PAGEBY PI1;
62         LABEL PI1     ="Control rate, PI1"
63               PI2     ="Treated rate, PI2" ;
64         

NOTE: The PROCEDURE PRINT printed pages 1-2.
NOTE: The PROCEDURE PRINT used 1.09 seconds.


65         PROC PLOT DATA=ONE;
66         BY PI1 PI2 NOTSORTED;
67         PLOT POWER*TOTAL_N="/" / VPOS=35 VZERO VAXIS=0 TO 1 BY .1
68                                  HPOS=65 HAXIS= 0 TO 2000 BY 250  ;
69         LABEL PI1     ="Control rate, PI1"
70               PI2     ="Treated rate, PI2" ;
3                              The SAS System   23:57 Friday, August 5, 1994

71         
NOTE: The PROCEDURE PLOT printed pages 3-4.
NOTE: The PROCEDURE PLOT used 0.61 seconds.


NOTE: SAS Institute Inc., SAS Campus Drive, Cary, NC USA 27513-2414