qualityTools (version 1.31.1)

contourPlot3: function to create a ternary plot (contour plot) for an object of class mixDesign

Description

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

Usage

contourPlot3(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

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

Myers; Montgomery; Anderson-Cook (2009) Response Surface Methodology - Wiley

See Also

wirePlot3

Examples

#yarn 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))
contourPlot3(A, B, C, elongation, data = mdo, form = "linear")
contourPlot3(A, B, C, elongation, data = mdo, form = "quadratic", col = 2)
contourPlot3(A, B, C, elongation, data = mdo, form = "elongation ~ I(A^2) - B:A + I(C^2)", col = 3)
contourPlot3(A, B, C, elongation, data = mdo, form = "quadratic", col = colorRampPalette(c("yellow", "white", "red")))