用r写一个函数，input为多个x和y值，通过比较huber函数和y=x^2函数的值，输出huber函数的delta的范围

这是一个简单的例子代码，仅供参考：

huber_delta_range <- function(x, y) {
  # 计算每个数据点的 huber 函数和 y=x^2 函数的值
  huber_vals <- huber_loss(x, y, delta=1)
  y_vals <- y - x^2
  
  # 计算 huber 函数和 y=x^2 函数的值的差
  delta_vals <- huber_vals - y_vals
  
  # 输出 delta 的范围
  delta_min <- min(delta_vals)
  delta_max <- max(delta_vals)
  cat("Huber函数的delta的范围为:", delta_min, "到", delta_max, "\n")
}

# 计算 huber 函数的值
huber_loss <- function(x, y, delta) {
  abs_diff <- abs(y - x^2)
  is_large <- abs_diff > delta
  large_diff <- delta * (abs_diff - 0.5 * delta)
  small_diff <- 0.5 * abs_diff^2
  huber_vals <- is_large * large_diff + (1 - is_large) * small_diff
  return(huber_vals)
}

# 例子数据
x <- c(1, 2, 3, 4, 5)
y <- c(1.2, 4, 9.1, 16.5, 26.8)

# 输出 delta 的范围
huber_delta_range(x, y)

这个函数的输入是 x 和 y，输出是 Huber 函数的 delta 的范围。函数中的 huber_loss 函数计算每个数据点的 Huber 函数的值，然后计算 Huber 函数和 y=x^2 函数的值的差，最后输出 delta 的范围。