1
次のコードを使用して、RSMパッケージのpersp関数を使用して、Rのx-y平面に輪郭投影を伴う3D応答曲面を描きました。rsmを使用してperspの後に線を描くR
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)パースビューのX-Y平面上の例えば、Y = 5X + 8追加のラインを描きますか?
(2)x-y平面上にzの散布図を投影しますか?