mixDesign {r6qualitytools} | R Documentation |
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
)
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. |
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
|
randomize |
Logical value. If 'TRUE' the RunOrder of the mixture design will be randomized (default). |
seed |
Nmerical value giving the input for set.seed. |
Value
The function mixDesig()
returns an object of class mixDesig()
.
Note
In this version the creation of (augmented) lattice, centroid mixture designs is fully supported. Getters and Setter methods for the mixDesign object exist just as for objects of class facDesign
(i.e. factorial designs).
The creation of constrained component proportions is partially supported but don't rely on it. Visualization (i.e. ternary plots) for some of these designs can be done with the help of the wirePlot3
and contourPlot3
function.
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()