qualityTools (version 1.31.1)

wirePlot3: function to create a ternary plot (3D wire plot) for an object of class mixDesign

Description

This function creates a ternary plot for mixture designs (i.e. object of class mixDesign)

Usage

wirePlot3(x, y, z, response, data = NULL, xlim, ylim, zlim, main, xlab, ylab, zlab, border, sub, form = "linear", phi, theta, ticktype, col = 1, col.text, cex.axis, axes = TRUE, steps, factors, fun, plot)

Arguments

x

factor 1 of the mixDesign object

y

factor 2 of the mixDesign object

z

factor 3 of the mixDesign object

response

the response of the mixDesign object

data

the mixDesign object

xlim

...

ylim

...

zlim

...

main

...

xlab

...

ylab

...

zlab

...

border

...

sub

...

form

"linear", "quadratic", "fullCubic", "specialCubic", or formula such as y ~ A + B + C

phi

...

theta

...

ticktype

...

col

a color or color.palette

col.text

...

cex.axis

...

axes

TRUE/FALSE

steps

resolution of the plot (number of rows for the square matrix)

factors

list of factors for categorizing with setting in case there are more than 3 factors (not yet implemented)

fun

...

plot

if FALSE no plot is being produced

Value

• a matrix containing NA's and numerics

References

Cornell (2002). Experiments with Mixtures. Third Edition. Wiley. Myers; Montgomery; Anderson-Cook (2009) Response Surface Methodology - Wiley

See Also

contourPlot3

Examples

#yarn elongation example p.564 Response Surface Methodology
mdo = mixDesign(3,2, center = FALSE, axial = FALSE, randomize = FALSE, replicates  = c(1,1,2,3))
names(mdo) = c("polyethylene", "polystyrene", "polypropylene")
units(mdo) = "percent"
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)
response(mdo) = elongation

dev.new(14,14); par(mfrow = c(2,2))
wirePlot3(A, B, C, elongation, data = mdo, form = "linear")
wirePlot3(A, B, C, elongation, data = mdo, form = "quadratic", col = 2)
wirePlot3(A, B, C, elongation, data = mdo, form = "elongation ~ I(A^2) - B:A + I(C^2)", col = 3)
wirePlot3(A, B, C, elongation, data = mdo, form = "quadratic", col = colorRampPalette(c("yellow", "white", "red")))