如何在R中计算矩阵行均值/行方差并求解其标准差？

首先修正你代码里的小问题：要生成100×100的矩阵，rgamma()的第一个参数得是100*100=10000，不然只会生成100个元素，矩阵会变成100行1列，不符合你的需求。

下面提供两种实现方式：

方式一：向量化操作（R推荐的高效写法，无需手动遍历） #

这种写法利用R的内置函数批量处理，比循环效率更高：

# 生成100×100的gamma分布矩阵
data <- matrix(rgamma(100*100, 4.623371, 2.757291), nrow = 100)

# 计算每行均值
row_means <- rowMeans(data)
# 计算每行方差（apply函数对每行调用var方法）
row_vars <- apply(data, 1, var)
# 得到每行均值/方差的结果向量x
x <- row_means / row_vars
# 计算x的标准差
x_sd <- sd(x)

# 输出结果
print(x_sd)

方式二：手动遍历（for循环实现） #

如果是为了理解遍历逻辑，用循环也能完成需求：

# 生成正确的100×100矩阵
data <- matrix(rgamma(100*100, 4.623371, 2.757291), nrow = 100)

# 初始化空向量存储每行的计算结果
x <- numeric(nrow(data))

# 遍历每一行
for (i in 1:nrow(data)) {
  row_data <- data[i, ]  # 取出第i行数据
  mean_val <- mean(row_data)  # 计算该行均值
  var_val <- var(row_data)  # 计算该行方差
  x[i] <- mean_val / var_val  # 存入结果向量
}

# 计算x的标准差
x_sd <- sd(x)

print(x_sd)

说明：R中优先推荐向量化操作，尤其是数据规模较大时，执行速度远快于手动循环；循环写法更适合用来理解逐行处理的逻辑。

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