3D Plot

Share:

Description

Creates a wireframe diagramm for an object of class facDesign.

Usage

1
2
3
wirePlot(x, y, z, data = NULL, xlim, ylim, zlim, main, xlab, ylab, border, sub, 
         zlab, form = "fit", phi, theta, ticktype, col = 1, steps, factors, fun,
         plot)

Arguments

x

name providing the Factor A for the plot.

y

name providing the Factor B for the plot.

z

name giving the Response variable.

data

needs to be an object of class facDesign and contains the names of x,y,z.

xlim

vector giving the range of the x-axis.

ylim

vector giving the range of the y-axis.

zlim

vector giving the range of the y-axis.

main

an overall title for the plot: see title.

xlab

a title for the x axis: title.

ylab

a title for the y axis: title.

zlab

a title for the z axis: title.

border

numerical value or character string giving the color of the line drawn around the surface facets.

sub

a sub title for the plot: title.

form

a character string or a formula with the syntax “y~ x+y + x*y”.
If form is a character it has to be one out of the following:

  • “quadratic”

  • “full”

  • “interaction”

  • “linear”

  • “fit”

“fit” takes the formula from the fit in the facDesign object fdo. Quadratic or higher orders should be given as I(Variable^2).
By default form is set as “fit”.

phi

numerical value. Angle (in degree) defining the viewing direction. Phi gives the colatitude.
By default phi is set to ‘30’.

theta

numerical value. Angle (in degree) defining the viewing direction. Theta gives the azimuthal direction.
By default theta is set to ‘-30’.

ticktype

graphical parameter. Character string specifying wheter the ticks on the axes should be “simple” or “detailed” (default).

col

a predefined (1, 2, 3 or 4) or self defined colorRampPalette or color to be used (i.e. “red”).

steps

number of grid points per factor. By default steps = 25.

factors

list of 4th 5th factor with value i.e. factors = list(D = 1.2, E = -1), if nothing is specified values will be the mean of the low and the high value of the factors.

fun

function to be applied to z “desirability”.

plot

logical value. If ‘TRUE’ (default) a plot is created (most likely to disappear with time)

Details

This function can be used to display the desirability of each response by specifying fun = “desirability” or the fun = “overall” (i.e. composed) desirability of all responses. The required desirabilities can be set using desires.

Value

The function wirePlot() returns an invisible list with the following entries:

  • x - locations of grid lines for x at which the values in z are measured

  • y - locations of grid lines for y at which the values in z are measured

  • z - a matrix containing the values of z to be plotted

Note

For an example in context which shows the usage of the function wirePlot() to an object of class facDesign, please read the vignette for the package qualityTools at http://www.r-qualitytools.org/html/Improve.html.

Author(s)

Thomas Roth thomas.roth@tu-berlin.de

See Also

contourPlot
filled.contour
persp
http://www.r-qualitytools.org/html/Improve.html

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
#create a response surface design and assign random data to response y
fdo = rsmDesign(k = 3, blocks = 2)
response(fdo) = data.frame(y = rnorm(nrow(fdo)))

par(mfrow = c(3,2))

#I - display linear fit
wirePlot(A,B,y, data = fdo, form = "linear")

#II - display full fit (i.e. effect, interactions and quadratic effects
wirePlot(A,B,y, data = fdo, form = "full")

#III - display a fit specified before
fits(fdo) = lm(y ~ B + I(A^2) , data = fdo)
wirePlot(A,B,y, data = fdo, form = "fit")

#IV - display a fit given directly
wirePlot(A,B,y, data = fdo, form = "y ~ A*B + I(A^2)")

#V - display a fit using a different colorRamp
wirePlot(A,B,y, data = fdo, form = "full", col = 2)

#V - display a fit using a self defined colorRamp
myColour = colorRampPalette(c("green", "gray","blue"))
wirePlot(A,B,y, data = fdo, form = "full", col = myColour)

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.