如何在R语言中求解单个不等式方程？以Gerrodette不等式为例

核心思路 #

你不需要用quadprog或limsolve这类多元方程工具——这是一元二次不等式求解问题，用R基础包就能搞定，还能轻松替换参数。

步骤1：参数转化与不等式推导 #

原始不等式：
r² * n³ > 12*cv²(zα/2 + zB)²

整理为关于zB的标准二次不等式形式：

0 < zB² + 2zα/2zB + (zα/2² - (r²n³)/(12cv²))

代入你提供的参数（r=0.02，n=40，cv=0.67，zα/2=1.645）后，得到：
0 < zB² + 3.3zB - 2.04
同时限定zB的范围为-3.4 < zB < 3.4

步骤2：R实现方案 #

方法1：直接求解当前参数的zB范围 #

利用R基础函数计算二次函数的根，再结合给定区间筛选有效解：

# 定义二次函数
f <- function(z) z^2 + 3.3*z - 2.04

# 计算两个实根
root1 <- (-3.3 - sqrt(3.3^2 + 4*2.04))/2  # 负根≈-3.832
root2 <- (-3.3 + sqrt(3.3^2 + 4*2.04))/2  # 正根≈0.532

# 结合zB的限定范围筛选有效区间
lower_bound <- -3.4
upper_bound <- 3.4

# 开口向上的二次函数，>0的区间为z<root1或z>root2，筛选在限定范围内的部分
valid_intervals <- list()
if (root2 < upper_bound) {
  valid_intervals <- c(valid_intervals, list(c(root2, upper_bound)))
}

# 输出结果
cat("满足条件的zB范围：\n")
for (interval in valid_intervals) {
  cat(sprintf("(%.4f, %.2f)\n", interval[1], interval[2]))
}

运行后得到有效区间：(0.5323, 3.40)

方法2：可替换参数的通用函数 #

如果需要频繁调整r、n、cv等参数，封装一个通用函数更高效：

solve_gerrodette_zB <- function(r, n, cv, z_alpha2, zB_min = -3.4, zB_max = 3.4) {
  # 计算二次不等式的系数
  a <- 1
  b <- 2 * z_alpha2
  c <- z_alpha2^2 - (r^2 * n^3) / (12 * cv^2)
  
  # 判别式判断是否有实根
  discriminant <- b^2 - 4*a*c
  if (discriminant <= 0) {
    return(list(valid_zB_intervals = c(zB_min, zB_max)))
  }
  
  # 计算排序后的根
  root1 <- (-b - sqrt(discriminant)) / (2*a)
  root2 <- (-b + sqrt(discriminant)) / (2*a)
  if (root1 > root2) {
    root1 <- root2 -> root2
  }
  
  # 筛选有效区间
  valid_intervals <- list()
  if (root1 > zB_min) {
    valid_intervals <- c(valid_intervals, list(c(zB_min, root1)))
  }
  if (root2 < zB_max) {
    valid_intervals <- c(valid_intervals, list(c(root2, zB_max)))
  }
  
  return(list(
    quadratic_coefficients = list(a=a, b=b, c=c),
    roots = c(root1=round(root1,4), root2=round(root2,4)),
    valid_zB_intervals = valid_intervals
  ))
}

# 测试你的参数
result <- solve_gerrodette_zB(r=0.02, n=40, cv=0.67, z_alpha2=1.645)
print(result$valid_zB_intervals)

# 更换参数示例
result_new <- solve_gerrodette_zB(r=0.03, n=35, cv=0.7, z_alpha2=1.96)
print(result_new$valid_zB_intervals)

为什么不用quadprog/limsolve？ #

这类工具针对的是多元线性/二次规划或联立方程组，你的问题是简单的一元二次不等式，用基础R的数学计算就能高效解决，无需额外依赖。

内容的提问来源于stack exchange，提问作者megsryder