1                              The SAS System   23:24 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\P0201.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};
16         BETA=J(2,1,0);
17         
18         HOLD=J(1,5,.);
19         HOLDNM={"PI1" "PI2" "N1" "N2" "POWER"};
20         INHOLD={.08  .04    50   50 ,
21                 .08  .04   100  100 ,
22                 .08  .04   200  200 ,
23                 .08  .04   400  400 ,
24         
25                 .05  .02    50   50 ,
26                 .05  .02   100  100 ,
27                 .05  .02   200  200 ,
28                 .05  .02   400  400  };
29         NCASES=NROW(INHOLD);
30         
31         ALPHA={.05};
32         REPN={1};
33           DO ICASE=1 TO NCASES;
34         
35           PI1=INHOLD[ICASE,1];
2                              The SAS System   23:24 Friday, August 5, 1994

36           PI2=INHOLD[ICASE,2];
37         
38           BETA[1,1]=PI1;
39           BETA[2,1]=PI2;
40         
41           SIGMA1=PI1#(1.-PI1);
42           SIGMA2=PI2#(1.-PI2);
43           SIGMA=(SIGMA1+SIGMA2)/2;
44         
45           N1=INHOLD[ICASE,3];
46           N2=INHOLD[ICASE,4];
47           ESSENCEX = BLOCK( J(N1,1,1) , J(N2,1,1)  );
48         
49           RUN POWER;
50         
51           POWCASES=NROW(_HOLDPOW);
52           TEMP=J(POWCASES,1,1) @ INHOLD[ICASE,*] ;
52                                                     *DUPLICATE INFO;
53           HOLD=HOLD // ( TEMP || _HOLDPOW[*,6] ) ;
54         *GRAB OTHER COLMNS AS NEEDED, BUT 1) ADD VAR NAMES IN HOLDNM,
55                                           2) DEFINE HOLD WITH CORRECT # 
COL'S;
56           END;
57         HOLD=HOLD[2:NROW(HOLD),*];
58         CREATE ONE VAR HOLDNM;
59         APPEND FROM HOLD;
60         CLOSE ONE;
61         FINISH;
NOTE: Module SHELL defined.
61                 * SHELL;
62         RUN  SHELL;
NOTE: The data set WORK.ONE has 8 observations and 5 variables.
63         
Exiting IML.
NOTE: Storage library LOAD01.POWSTORE closed.
NOTE: The PROCEDURE IML used 8.63 seconds.


64         PROC PRINT DATA=ONE NOOBS;
65         BY PI1 NOTSORTED;
66         
NOTE: The PROCEDURE PRINT printed page 1.
NOTE: The PROCEDURE PRINT used 0.96 seconds.


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