r6qualitytools (version 1.0.1)

mixDesign: mixDesign: Mixture Designs

Description

Function to generate simplex lattice and simplex centroid mixture designs with optional center points and axial points.

Usage

mixDesign(
  p,
  n = 3,
  type = "lattice",
  center = TRUE,
  axial = FALSE,
  delta,
  replicates = 1,
  lower,
  total = 1,
  randomize,
  seed = 1234
)

Value

The function mixDesig() returns an object of class mixDesig().

Arguments

p

Numerical value giving the amount of factors.

n

Numerical value specifying the degree (ignored if type = `centroid`).

type

Character string giving the type of design. type can be `lattice` or `centroid` (referencing to the first source under the section references]. By default type is set to `lattice`.

center

Logical value specifying whether (optional) center points will be added. By default `center` is set to `TRUE`.

axial

Logical value specifying whether (optional) axial points will be added. By default `axial` is set to `FALSE`.

delta

Numerical value giving the delta (see references) for axial runs. No default setting.

replicates

Vector with the number of replicates for the different design points i.e. c(center = 1, axial = 1, pureBlend = 1, BinaryBlend = 1, p-3 blend, p-2 blend, p-1 blend). By default `replicates` is set to `1`.

lower

Numeric vector of lower-bound constraints on the component proportions (i.e. must be given in percent).

total

Numeric vector with

• [1] the percentage of the mixture made up by the q - components (e.g. q = 3 and x1 + x2 + x3 = 0.8, total = 0.8 with 0.2 for the other factors being held constant)

• [2] overall total in corresponding units (e.g. 200ml for the overall mixture)

randomize

Logical value. If `TRUE` the RunOrder of the mixture design will be randomized (default).

seed

Nmerical value giving the input for set.seed.

See Also

mixDesign.c, facDesign.c, facDesign, fracDesign, rsmDesign, wirePlot3, contourPlot3.

Examples

# Example usage of mixDesign
mdo <- mixDesign(3, 2, center = FALSE, axial = FALSE, randomize = FALSE, replicates = c(1, 1, 2, 3))

mdo$names(c("polyethylene", "polystyrene", "polypropylene"))
elongation <- c(11.0, 12.4, 15.0, 14.8, 16.1, 17.7,
                16.4, 16.6, 8.8, 10.0, 10.0, 9.7,
                11.8, 16.8, 16.0)
mdo$.response(elongation)

mdo$units()
mdo$summary()