我在绘制交互和使用预测函数时遇到问题。这是我的最终模型：
`mymodel <- glm(y ~ SITE + TREAT*HET, family= quasibinomial, data = PRESU)`.
我想绘制Treat（分类，2级）和Het（连续）之间的互动。继Crawley第16章（或多或少，但示例中没有其他变量）之后，我尝试了这段代码：
x <-seq(0,1,length=100)
v <- rep("treat1", length(x))
y <- predict(mymodel, list(TREAT=factor(v),HET=x), type="response")
plot(HET[TREAT=="treat1"], fitted(mymodel)[TREAT=="treat1"],col="black")
lines(x,y,col="Black",lwd=1.5,lty=3)
t <- rep("treat2", length(x))
s <- predict(mymodel, list(TREAT=factor(t),HET=x), type="response")
points(HET[TREAT=="treat2"],fitted(mymodel)[TREAT=="treat2"],col="red")
lines(x,s,col="red",lwd=1.5)
我在预测中得到了这个输出
eval（expr，envir，enclos）出错：找不到对象'SITE'。
我无法理解如何从模型中绘制我的交互。如果我删除SITE代码有效，但我想将其保留在最终模型中。我希望我能清楚地解释一下。感谢谁能提供帮助。
数据集
SITE TREAT1 HET Tot Prey
1 treat1 0.3 30 9
1 treat2 0.5 30 16
2 treat1 0.8 30 19
2 treat2 0.76 30 25
3 treat2 0.97 30 1
3 treat1 0.75 30 23
4 treat1 0.67 30 18
4 treat2 0.89 30 29
我对每个网站都有更多的重复，但是我们说这与我的数据集类似。响应变量是y =（猎物，tot-prey）