demoRSA {RSA}

Plots a response surface of a polynomial equation of second degree with interactive controls

Description

Plots an RSA object, or a response surface with specified parameters, with interactive controls for coefficients.

Usage

demoRSA(
  x = NULL,
  y = 0,
  x2 = 0,
  y2 = 0,
  xy = 0,
  w = 0,
  wx = 0,
  wy = 0,
  x3 = 0,
  xy2 = 0,
  x2y = 0,
  y3 = 0,
  b0 = 0,
  type = "3d",
  zlim = c(-2, 2),
  xlim = c(-2, 2),
  ylim = c(-2, 2),
  xlab = NULL,
  ylab = NULL,
  zlab = NULL,
  points = TRUE,
  model = "full",
  project = c("PA1", "PA2"),
  extended = FALSE,
  ...
)

Arguments

x

Either an RSA object (returned by the RSA function), or the coefficient for the X predictor

y

Y coefficient

x2

X^2 coefficient

y2

Y^2 coefficient

xy

XY interaction coefficient

w

W coefficient (for (un)constrained absolute difference model)

wx

WX coefficient (for (un)constrained absolute difference model)

wy

WY coefficient (for (un)constrained absolute difference model)

x3

X^3 coefficient

xy2

XY^2 coefficient

x2y

X^2Y coefficient

y3

Y^3 coefficient

b0

Intercept

type

3d for 3d surface plot, contour for 2d contour plot. Shortcuts (i.e., first letter of string) are sufficient; be careful: "contour" is very slow at the moment

zlim

Limits of the z axis

xlim

Limits of the x axis

ylim

Limits of the y axis

xlab

Label of the x axis

ylab

Label of the y axis

zlab

Label of the z axis

points

A list of parameters which define the appearance of the raw scatter points: show = TRUE: Should the original data points be overplotted? value="raw": Plot the original z value, "predicted": plot the predicted z value. jitter=0: Amount of jitter for the raw data points. cex = .5: multiplication factor for point size. See ?plotRSA for details.

model

If x is an RSA object: from which model should the response surface be computed?

project

Which features should be projected on the floor? See ?plotRSA for details.

extended

Show additional controls (not implemented yet)

...

Other parameters passed through to plot.RSA (e.g., xlab, ylab, zlab, cex, legend)

Details

No details so far. Just play around with the interface!

See Also

plotRSA, RSA

Examples

# Plot response surfaces from known parameters
# example of Edwards (2002), Figure 3
## Not run: 
demoRSA(x=.314, y=-.118, x2=-.145, y2=-.102, xy=.299, b0=5.628, type="3d")
demoRSA(x=.314, y=-.118, x2=-.145, y2=-.102, xy=.299, b0=5.628, legend=FALSE, type="c")

## End(Not run)

# Plot response surface from an RSA object
## Not run: 
set.seed(0xBEEF)
n <- 300
err <- 2
x <- rnorm(n, 0, 5)
y <- rnorm(n, 0, 5)
df <- data.frame(x, y)
df <- within(df, {
	diff <- x-y
	absdiff <- abs(x-y)
	SD <- (x-y)^2
	z.diff <- diff + rnorm(n, 0, err)
	z.abs <- absdiff + rnorm(n, 0, err)
	z.sq <- SD + rnorm(n, 0, err)
	z.add <- diff + 0.4*x + rnorm(n, 0, err)
	z.complex <- 0.4*x + - 0.2*x*y + + 0.1*x^2 - 0.03*y^2 + rnorm(n, 0, err)
})

r1 <- RSA(z.sq~x*y, df)
demoRSA(r1)
demoRSA(r1, points=TRUE, model="SQD")

## End(Not run)