# Plots a response surface of a polynomial equation of second degree with interactive controls Plots an RSA object, or a response surface with specified parameters, with interactive controls for coefficients. ```r 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! ## Examples ```r # 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) ``` ## See Also `plotRSA`, `RSA`

demoRSA function