2016-03-14 3 views
-1

Я пытаюсь найти способ найти минимальный/максимальный из квадратичной модели. В этом случае минимум.Как найти min max от lm

x.lm <- lm(Y ~ X + I(X^2)) 

Edit: Для того, чтобы уточнить, я уже могу найти минимум у через мин (предсказать (x.lm)). Как я могу перевести это в соответствующее значение x.

+0

Пожалуйста, покажите воспроизводимый пример. – akrun

+1

Некоторые [основные математические школы] (http://www.wikihow.com/Find-the-Maximum-or-Minimum-Value-of-a-Quadratic-Function-Easily) позволяют рассчитать это из коэффициентов. – Roland

ответ

2

Проверьте это. Идея в том, что вы должны принять подогнанные значения образуют x.lm подходят

#example data 

X <- 1:100 
Y <- 1:100 + rnorm(n = 100, mean = 0, sd = 4) 

x.lm <- lm(Y ~ X + I(X^2)) 

fits <- x.lm$fitted.values #getting fits, you can take residuals, 
# and other parameters too 

# I guess you are looking for this. 
min.fit = min(fits) 
max.fit = max(fits) 

После очередного вопроса

df <- cbind(X, Y, fits) 
df <- as.data.frame(df) 
index <- which.min(df$fits) #very usefull command 
row.in.df <- df[index,] 
+0

Это дает мне самый низкий у кривой. Как насчет использования min (x.lm $ fit.values), чтобы найти соответствующее значение x? – thomasrive

+0

Я обновил свой ответ. Ofcoruse вам не нужно связывать X, Y и вписывается в один фрейм данных. Вам просто нужно получить индекс с помощью функции min, или max, и по этому индексу у вас есть соответствующий номер из X, your.desired.min <- X [who.min (подходит)] –

+0

Спасибо. Это было очень близко к ручному методу, который я выполнил, чтобы найти ответ. – thomasrive

Смежные вопросы