aux_1 <- numeric(length(y))
for(i in seq_along(y)){
aux_1[i] <- optim(0, objetivo, k = y[i], model = model, method = "Brent",
lower = min(x), upper = max(x))$par
}
# Tomando os resultados mais próximos:
result <- round(aux_1,0)
result
[1] 100 51 50 27 125 25 76 100 124 25 99 124 125 76 75 50 50 74
[19] 125 100 25 74 100 77 26 102 25 74 49 98 75 76 125 125 75 25
extract <- function(x,y){
y <- sort(y, decreasing = T)
for(i in seq_along(y)){
x[x/y[i] < 1.25 & x/y[i] >= 0.98] = y[i]}
return(x)
}
result <- extract(result,v)
result
[1] 100 50 50 25 125 25 75 100 125 25 100 125 125 75 75 50 50 75
[19] 125 100 25 75 100 75 25 100 25 75 2550 100 75 75 125 125 75 25
EDIT_2 (15/06/18, 20:20): Refiz o vetor de valores para ficar mais similar ao do data.frame apresentado como exemplo e acrescentei as boas sugestões do comentário do Rui Barradas abaixo. Ainda, como essa resposta carecia de aproximação para ficar apenas com valores do vetor (25, 50, 75, 100, 125), criei a função 'extract'. No entanto, é uma função específica para os valores desse exemplo, pode se ajustar os para casos mais gerais.