Контурный график для квадратичного уравнения в R

Question

Я пытаюсь воссоздать контурный график, показанный на рисунке ниже, используя R. Это контурный график постоянной потери сиропа в зависимости от скорости и давления.

Данные, используемые для этого примера:

speed  = c(100,100,100,100,100,100,120,120,120,120,120,120,140,140,140,140,140,140) 
pressure = c(10,10,15,15,20,20,10,10,15,15,20,20,10,10,15,15,20,20) 
syrup_loss = c(-35,-25,110,75,4,5,-45,-60,-10,30,-40,-30,-40,15,80,54,31,36) 

с этими данными квадратного уравнения была создана:

model <- lm(syrup_loss~speed + pressure + speed^2 + pressure^2 + speed*pressure) 
summary(model) 

Что дает:

г = 1217,30556 - 31,25625 * x + 86,01667 * y + 0,1291 * x^2 - 2,87333 * y^2 + 0,02875 * x * у

Я попытался создать контур участка, используя этот код, но он не дает приемлемого результата:

x = seq(100, 140, len=100) 
y = seq(10, 20, len=100) 
z = outer(x, 1217.30556 - 31.25625*x + 86.01667*y + 0.1291*x^2 - 
      2.87333*y^2 + 0.02875*x*y) 
contour(x, y, z, nlev=12) 

Answer 1

В функции outer, то Y аргумент отсутствует. Код

z <- outer(x, y, function(x, y) 1217.30556 - 31.25625*x + 86.01667*y + 0.1291*x^2 - 2.87333*y^2 + 0.02875*x*y) 

исправлена ​​проблема.