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")
Это построит плоскость регрессии с точками, но я бы хотел построить это без очков.Пропуск x, y, z при вызове функции scatter3D приводит к ошибке.