马尔可夫模型转移概率数组行和检查异常问题咨询

问题背景 #

我正在构建决策分析模型（马尔可夫模型），需要验证一个85行、85列、35切片的转移概率数组每行和是否等于1。

原代码及异常现象 #

使用以下R代码进行检查：

# 设置标记跟踪所有切片的行和是否均为1
all_slices_have_rowsums_equal_to_1 <- FALSE

# 遍历a_P_tunnels中的所有切片
for (i in 1:35) {
  # 计算当前切片的行和
  row_sums <- rowSums(a_P_tunnels [ , , i])
  
  # 检查是否存在行和不等于1的情况
  if (any(row_sums != 1)) {
    # 输出切片编号及不等于1的行和
    cat("Slice", i, "has row sums", row_sums[row_sums != 1], "\n")
  }
}

执行后异常输出了所有35个切片的行和均为1的信息，但代码逻辑是仅当行和不等于1时才输出，这一现象的原因是什么？

异常原因 #

核心问题是浮点数精度误差：在R的数值计算中，浮点数运算会产生微小的精度偏差，比如实际行和可能是0.999999999999999或1.000000000000001这类近似1的值。原代码用row_sums != 1做严格相等判断，会把这些因计算误差导致的近似值误判为不等于1，进而触发输出逻辑。

修正后可行代码 #

解决思路是引入精度容差，判断行和与1的差值是否在可接受范围内，修正后的代码如下：

# 遍历a_P_tunnels中的所有切片
for (i in 1:35) {
  # 计算当前切片的行和
  row_sums <- rowSums(a_P_tunnels [ , , i])
  
  # 检查是否存在行和与1的差异超出5位小数容差的情况
  if (!isTRUE(all.equal(row_sums, 1, tolerance = 1e-5))) {
    # 计算每个行和与1的绝对差值
    differences <- abs(row_sums - 1)
    if (any(differences > 1e-5)) {
      # 输出切片编号及不符合要求的行和（保留4位小数）
      cat("Slice", i, "has row sums", format(round(row_sums[differences > 1e-5], 4), nsmall = 20), "\n")
    }
  }
}

这段代码通过all.equal()函数设置容差tolerance = 1e-5（允许行和与1的差值不超过0.00001），同时通过计算绝对差值二次确认，避免了浮点数精度误差导致的误判，仅在行和与1的差异确实超出合理范围时才输出提示。

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