0

最近我进行了一个二元混合实验，目前我正在尝试使用 persp 函数绘制模型。然而，由于混合实验的性质，两个混合分量的值不应超过 1（即 max(x+y)=1）。然而，在 persp 函数中加入这种累积约束时，我遇到了困难。

下面的代码显示了一个没有累积约束的绘图示例：

x<-c(0,0.25,0.50,0.75,1)
y<-c(1,0.75,0.50,0.25,0)
z<-rnorm(5,40,1)
xyz<-bind_cols(x,y,z)
lmxyz<-lm(z~x+y,data=xyz)
persp(lmxyz,~x+y,image=TRUE,theta=-50)

没有累积约束的绘图

在 persp 函数中使用“bounds”参数会产生以下错误：

代码：

persp(lmxyz,~x+y,image=TRUE,bounds = list((x+y)<1),theta=-50)

结果：

Error in contour.lm(x, form, at, bounds, zlim, xlabs, atpos = atpos, plot.it = FALSE) : 
  'bounds' must be a NAMED list of bounds for each variable

问题：有没有办法在 persp 函数中加入累积约束？或者是否有更可行的方法来可视化这些类型的模型（包括 rsm 模型），包括这种类型的约束？

4

0 回答 0