View source: R/4.3_Mix_designs_Functions.R
mixDesign | R Documentation |
Function to generate simplex lattice and simplex centroid mixture designs with optional center points and axial points.
mixDesign(
p,
n = 3,
type = "lattice",
center = TRUE,
axial = FALSE,
delta,
replicates = 1,
lower,
total = 1,
randomize,
seed = 1234
)
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. |
The function mixDesig()
returns an object of class mixDesig()
.
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.
mixDesign.c
, facDesign.c
, facDesign
, fracDesign
, rsmDesign
, wirePlot3
, contourPlot3
.
# 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()
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.