如何在R语言中统计下三角矩阵的元素数量

嘿，针对你遇到的这个问题，我来给你梳理几种在R里的实现方式，刚好匹配你给出的示例场景：

首先先明确：你示例里的矩阵是严格下三角矩阵（对角线及以上全为0，非零元素都在对角线下方），这类矩阵的严格下三角区域元素总数有通用的数学公式，当然也可以用R的内置函数来灵活计算，分两种常见情况说明：

1. 计算严格下三角区域的总元素数（不管元素是否为0） #

如果你只是想知道严格下三角区域一共有多少个元素（比如你的5x5矩阵就是10个），可以先获取矩阵的阶数，再用经典的求和公式计算：

# 先还原你的矩阵
Matrix <- c(0, 4, 3, 1, 2, 0, 0, 3, 2, 1, 0, 0, 0, 2, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0)
Matrix <- matrix(Matrix, 5, 5)

# 获取矩阵的行数/列数（因为是方阵，nrow和ncol结果一致）
n <- nrow(Matrix)
# 严格下三角元素数公式：1+2+...+(n-1) = n*(n-1)/2
lower_tri_total <- n * (n - 1) / 2
print(lower_tri_total)
# 输出结果：10

这个公式是通用的，任何n阶方阵的严格下三角区域元素数都可以这么算。

2. 统计严格下三角区域中的非零元素数 #

如果你的矩阵里严格下三角区域可能混有0，而你只想统计其中的非零元素（就像你的示例那样），可以用R内置的lower.tri()函数来快速筛选：

# lower.tri()会返回和原矩阵同维度的逻辑矩阵，严格下三角位置为TRUE
# 用它索引原矩阵，再筛选非零元素，最后sum统计数量
non_zero_count <- sum(Matrix[lower.tri(Matrix)] != 0)
print(non_zero_count)
# 输出结果：10，和你的示例一致

原理很简单：Matrix[lower.tri(Matrix)]会提取出所有严格下三角的元素，然后!= 0会把非零元素转为TRUE（数值1），零转为FALSE（数值0），sum()就会直接统计出非零元素的数量。

补充：如果需要包含对角线的下三角区域 #

要是你需要统计的是包含对角线的下三角区域元素数，那公式换成n*(n+1)/2，筛选时给lower.tri()加个diag=TRUE参数就行：

# 包含对角线的下三角总元素数
lower_tri_with_diag_total <- n * (n + 1) / 2
print(lower_tri_with_diag_total)
# 输出：15

# 包含对角线的下三角非零元素数
non_zero_with_diag <- sum(Matrix[lower.tri(Matrix, diag = TRUE)] != 0)
print(non_zero_with_diag)
# 输出：10（因为你的示例里对角线全是0）

这样就能根据你的实际需求灵活选择啦～

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