back to DISC 7353 main page

mail E-mail: JLPeixoto@uh.edu

Home Page Julio L. Peixoto's home page


Replicated Blocked Design with Partial Confounding (SAS)

NOTE: This program can be run interactively by submitting after each run.

PROGRAM:

/* First create a 2 to the 3 design in 2 blocks for estimating
main effects */
proc factex;
factors A B C;
blocks nblocks=2;
model est=(A B C); /* we want to estimate A, B, C */
examine confounding aliasing;
output out=rep1 blockname=block nvals=(1 2);
run;
/* Sinc proc factex is interactive, we can proceed directly to
the second replicate, without giving proc factex again. In this
second replicate we want ABC to be estimable. */
model est=(A B C A*B*C); /* Include ABC, additional effect */
output out=rep2
blockname=block nvals=(3 4); /* Use new blocks */
run;
/* Now we estimate BC with the third replicate */
model est=(A B C A*B*C B*C);
output out=rep3 blockname=block nvals=(5 6);
run;
/* Finally, we estimate AC with the fourth replicate */
model est=(A B C A*B*C B*C A*C);
output out=rep4 blockname=block nvals=(7 8);
run;
/* Combine the four replicates: */
data combine;
set rep1 rep2 rep3 rep4;
run;
proc print data=combine;
run;

OUTPUT:

The SAS System 17:05 Sunday, April 20, 1997 9

Block Pseudo-factor Confounding Rules

[B1] = A*B*C
The SAS System 17:05 Sunday, April 20, 1997 10

Aliasing Structure

A
B
C
A*B
A*C
B*C

Block Pseudo-factor Confounding Rules

[B1] = B*C


Aliasing Structure

A
B
C
A*B
A*C
[B] = B*C
A*B*C

Block Pseudo-factor Confounding Rules

[B1] = A*C
Aliasing Structure

A
B
C
A*B
[B] = A*C
B*C
A*B*C

Block Pseudo-factor Confounding Rules

[B1] = A*B

Aliasing Structure

A
B
C
[B] = A*B
A*C
B*C
A*B*C

OBS BLOCK A B C

1 1 -1 -1 -1
2 1 -1 1 1
3 1 1 -1 1
4 1 1 1 -1
5 2 -1 -1 1
6 2 -1 1 -1
7 2 1 -1 -1
8 2 1 1 1
9 3 -1 -1 1
10 3 -1 1 -1
11 3 1 -1 1
12 3 1 1 -1
13 4 -1 -1 -1
14 4 -1 1 1
15 4 1 -1 -1
16 4 1 1 1
17 5 -1 -1 1
18 5 -1 1 1
19 5 1 -1 -1
20 5 1 1 -1
21 6 -1 -1 -1
22 6 -1 1 -1
23 6 1 -1 1
24 6 1 1 1
25 7 -1 1 -1
26 7 -1 1 1
27 7 1 -1 -1
28 7 1 -1 1
29 8 -1 -1 -1
30 8 -1 -1 1
31 8 1 1 -1
32 8 1 1 1