如何从具有2个预测变量的线性模型生成热图

Question

我已经在2个变量x和y上拟合了模型。我想可视化表格或热图中的所有z值。尽管收到了“ x必须是数字矩阵”的错误，但我在下面的代码中尝试了此操作。有更好的方法吗？

model_1 <- glm(z ~ x + y, 
                   data = mydata2, 
                   family=binomial()
    )

x <- seq(0,100,10)
y <- seq(0,100,10)

vb <- expand.grid(x=x, y=y) 

p <- predict(model_1, type="response", newdata = vb)
heatmap(p)