opt_rotation: Finds an Optimal Rotation Vector

Description Usage Arguments Value Examples

Description

This function will produce an optimal rotation vector for some design. It will produce Resolution IV if possible and if not will find a minimal aberration. If multiple minimal aberration designs exist it will randomly choose one and ignore any "desirability" between factors.

Usage

1
opt_rotation(design)

Arguments

design

This needs to be a coded design matrix using (0,1,2) where each row corresponds to one treatment.

Value

It will produce a list

Rotation Vector

The first element of the list will be the optimum rotation vector and a small message whether Resoluiton IV is obtainable or if this is minimal abberration

psuedo-design matrix

This is the fake design matrix which is the coded effects of length 3 that are aliased with the intercept in the original design. This will not neccessarly be a regular fractional factorial design (hence psuedo)

Examples

1
2
3
4
5
6
7
x <- c("a", "b", "c", "d", "e", "f")
l <- matrix(c(1,2,0,1,1,1, 0,1,1,2,0,0, 1,1,1,2,2,0), nrow = 3, byrow = TRUE)
trees <- c(0,0,0)
m <- TripleFrac:::part(x, l, trees)[[1]]
head(m)
what_frac(m)
opt_rotation(m)

vinny-paris/TripleFrac documentation built on May 10, 2019, 8:27 a.m.