求助：R中反向计算边际税的方法及相互依赖变量b、c求解

一、反向计算边际税 #

R没有内置的反向边际税计算函数，但可以自己编写循环或自定义函数实现。核心逻辑是依据税率阶梯规则，从税后收入倒推税前收入：

• 先明确各税率区间的起征点、适用税率及速算扣除数
• 通过循环从最高税率区间开始判断，逐步扣除对应区间的税额，最终得到税前收入

以下是两档税率的反推示例（假设：收入≤50000税率10%，超过部分20%）：

reverse_marginal_tax <- function(after_tax) {
  # 定义税率区间：(下限, 上限, 税率, 速算扣除数)
  tax_brackets <- list(
    c(0, 50000, 0.1, 0),
    c(50000, Inf, 0.2, 5000)
  )
  
  pre_tax <- after_tax
  for (bracket in tax_brackets) {
    lower <- bracket[1]
    upper <- bracket[2]
    rate <- bracket[3]
    deduction <- bracket[4]
    
    # 计算当前区间对应的税后收入范围
    bracket_after_min <- lower * (1 - rate) + deduction
    bracket_after_max <- if (is.infinite(upper)) Inf else upper * (1 - rate) + deduction
    
    if (after_tax >= bracket_after_min && after_tax <= bracket_after_max) {
      pre_tax <- (after_tax - deduction) / (1 - rate)
      break
    }
  }
  return(pre_tax)
}

# 测试
reverse_marginal_tax(45000)  # 对应税前50000
reverse_marginal_tax(85000)  # 对应税前100000

二、求解相互依赖的变量b和c #

你的问题本质是解线性方程组，先手动推导逻辑：
已知条件：

a = 20000
b = 0.1 * c
c = a + b

将b代入c的公式：
c = 20000 + 0.1c → 0.9c = 20000 → c = 20000 / 0.9 ≈ 22222.22
b = 0.1 * c ≈ 2222.22

用R代码实现有两种方式：

# 方式一：直接计算
a <- 20000
c <- a / 0.9
b <- 0.1 * c

cat("b =", round(b, 2), "\nc =", round(c, 2))

# 方式二：用solve函数解线性方程组
# 转化为方程组：b - 0.1c = 0；c - b = a
A <- matrix(c(1, -0.1, -1, 1), nrow=2)
B <- c(0, a)
solution <- solve(A, B)
b <- solution[1]
c <- solution[2]

cat("\n用solve求解结果：\nb =", round(b, 2), "\nc =", round(c, 2))

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