Существует ли функция R, которая разрешает линейную модель второго порядка?

Question

Я - новичок в R и программировании, и я стараюсь выполнить, вероятно, простую задачу. Я сделал код, который создает второй порядок модели, и я хочу ввести переменные в этой модели и найти «значение Y»

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

Я сделал это до сих пор:

modFOI <- rsm(Rendimento~FO(x1,x2,x3,x4)+TWI(x1,x2,x3,x4)+PQ(x1,x2,x3,x4),data=CR) # com interações
summary(modFOI)
print(modFOI)

После этого я нашел модель SO, но теперь я хочу создать переменные, такие как x1, x2, x3, и ввести их в модель и найти Y. Я также хотел бы найти оптимальный Y

Answer 1

Оптимальное значение отображается как постоянная точка на выходе summary(modFOI). Вы также можете запустить steepest(modFOI), чтобы увидеть трассу оценочных значений по пути наискорейшего подъема.

Для прогнозирования создайте фрейм данных с требуемыми наборами значений x. Например,

testdat <- data.frame(x1 = -1:1, x2 = 0, x3 = 0, x4 = 1)

Затем используйте функцию predict() с newdata:

predict(modFOI, newdata = testdat)

Answer 2

Самый простой способ создать многочлен (2-й порядок), о котором я могу подумать, следующий:

DF <- data.frame(x = runif(10,0,1),
                 y = runif(10,0,1) )
mod <- lm(DF$y ~ DF$x + I(DF$x^2))

predict(mod, new.data=data.frame(x=c(1,2,3,4,5)))

Примечание.при использовании predict new.data должен быть в формате data.frame, а переменная должна иметь то же имя, что и переменная в модели (здесь x)

Надеюсь, это поможет