conf.design (version 2.0.0)

conf.set: Find all confounded effects

Description

Find minimal complete sets of confounded effects from a defining set for symmetric confounded factorial designs. Useful for checking if a low order interaction will be unintentionally confounded with blocks. As in the usual convention, only effects whose leading factor has an index of one are listed.

All factors must have the same (prime) number of levels.

Usage

conf.set(G, p)

Arguments

G

Matrix whose rows define the effects to be confounded with blocks, in the same way as for conf.design.

p

Number of levels for each factor. Must be a prime number.

Value

A matrix like G with a minimal set of confounded with blocks defined in the rows.

Side Effects

None

Details

The function constructs all linear functions of the rows of G (over GF(p)), and removes those rows whose leading non-zero component is not one.

See Also

conf.design

Examples

### If A B^2 C and B C D are confounded with blocks, then so are A C^2 D
### and A B D^2.

G <- rbind(c(A = 1, B = 2, C = 1, D = 0),
           c(A = 0, B = 1, C = 1, D = 1))
conf.set(G, 3)
###      A B C D
### [1,] 1 2 1 0
### [2,] 0 1 1 1
### [3,] 1 0 2 1
### [4,] 1 1 0 2

### Only three-factor interactions are confounded, so the design is
### presumably useful.

as.matrix(replications( ~ .^2, conf.design(G, 3)))
###          [,1]
### Blocks      9
### A          27
### B          27
### C          27
### D          27
### Blocks:A    3
### Blocks:B    3
### Blocks:C    3
### Blocks:D    3
### A:B         9
### A:C         9
### A:D         9
### B:C         9
### B:D         9
### C:D         9