R 中没有点的回归平面的 3D 图

Question

dat <- data.frame(nitrogen = runif(50, 0, 10), temperature= rnorm(50, 10, 3))
modmat <- model.matrix(~ nitrogen * temperature, dat)
coeff <- c(1, 2, -1, 1.5)
dat$soil <- rnorm(50, mean = modmat %*% coeff, sd = 0.5)

我正在尝试绘制一个回归平面，显示 temperature 如何调节 nitrogen 和 soil 之间的关系。

library(plot3D)

x <- dat$nitrogen
y <- dat$temperature
z <- dat$soil

fit <- lm(z ~ x*y)

grid.lines = 26
x.pred <- seq(min(x), max(x), length.out = grid.lines)
y.pred <- seq(min(y), max(y), length.out = grid.lines)
xy <- expand.grid( x = x.pred, y = y.pred)
z.pred <- matrix(predict(fit, newdata = xy), 
                 nrow = grid.lines, ncol = grid.lines)

fitpoints <- predict(fit)

scatter3D(x, y, z, pch = 18, cex = 2, 
    theta = 20, phi = 20, ticktype = "detailed",
    xlab = "nitrogen", ylab = "soil", zlab = "temperature",  
    surf = list(x = x.pred, y = y.pred, z = z.pred,  
    facets = NA, fit = fitpoints), main = "dat")

这将绘制带点的回归平面，但我想绘制不带点的回归平面。调用 scatter3D() 函数时省略 x, y, z 会导致错误。

Answer 1

这是一个疯狂的想法。尝试设置 cex = 0.

scatter3D(x, y, z, pch = 18, cex = 0, 
    theta = 20, phi = 20, ticktype = "detailed",
    xlab = "nitrogen", ylab = "soil", zlab = "temperature",  
    surf = list(x = x.pred, y = y.pred, z = z.pred,  
    facets = NA, fit = fitpoints), main = "dat")

R 中没有点的回归平面的 3D 图

3D plot of regression plane without points in R

regression

r

plot3d