RGL Embedded Least Squares / ANOVA Plot

[coatless]

Interactive vs. stagnate plots requires use of rgl

To use with rmarkdown, might need to enable hook_webgl()? Need fallback for pdf version. (JS Error)

open3d()
x = rnorm(100)
y = rnorm(100)
z = 0.2*x - 0.3*y + rnorm(100, sd = 0.3)
fit = lm(z ~ x + y)
plot3d(x, y, z, type = "s", col = "red", size = 1)
coefs = coef(fit)
a = coefs["x"]
b = coefs["y"]
c = -1
d = coefs["(Intercept)"]
planes3d(a, b, c, d, alpha = 0.5)

Better version in car::scatter3d:

http://www.inside-r.org/packages/cran/car/docs/identify3d

See:

http://www.ats.ucla.edu/stat/r/dae/intreg.htm

For anova plot