r - 在 R 中使用 rsm 在 persp 之后画一条线

Question

我使用以下代码使用 RSM 包的 persp 函数在 R 中的 xy 平面上绘制具有轮廓投影的 3d 响应曲面。

x <- seq(-3,3,by=0.25) 
y <- seq(-3,3,by=0.25) 
d <- expand.grid(x=x,y=y)
z <- c(data=NA,1089)
b0 = 5.628; b1 = 0; b2 = 0; b3 = -.1 ; b4 =.1; b5 = -.1
k=1
for (i in 1:25) {
for (j in 1:25) {
z[k]=b0+b1*x[i]+b2*y[j]+b3*x[i]*x[i]+b4*x[i]*y[j]+ b5*y[j]*y[j]
k=k+1
}  }
library(rsm)
data.lm <- lm(z~poly(x,y,degree=2),data=d)
persp(data.lm,x~y, zlim=c(0,max(z)),contour=list(z="bottom",col="colors"),theta=-55,phi=25)

画完这张图后，有没有办法可以：

(1) 在 persp xy 平面上画一条额外的线，例如 y=5x+8？

(2) 在 xy 平面上投影 z 的散点图？

Answer 1

该行可以添加trans3d：

res1 <- persp(data.lm,x~y, zlim=c(0,max(z)),contour=list(z="bottom",col="colors"),theta=-55,phi=25)    
xy <- matrix(c((-3-8)/5,-3,(3-8)/5,3),ncol=2,byrow = T)
lines(trans3d(xy[,2], xy[,1], 0, pmat = res1$`y ~ x`$transf), col = 3)