demoRSA | R Documentation |
Plots an RSA object, or a response surface with specified parameters, with interactive controls for coefficients.
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, ... )
x |
Either an RSA object (returned by the |
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 |
|
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) |
No details so far. Just play around with the interface!
plotRSA
, RSA
# 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)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.