在data.table中按行计算连续列方差遇错：rowVars参数需为矩阵/向量

我明白你遇到的问题了——当你用data.table的.SD直接传给rowVars()时，会触发Error in rowVars(.SD): Argument 'x' must be a matrix or a vector.的错误，核心原因很简单：

rowVars()来自matrixStats包，它仅接受矩阵或向量作为输入，但data.table中的.SD返回的是一个data.table对象（本质是特殊的data.frame），不符合函数的参数要求。

快速解决方案 #

只需要把.SD转换成矩阵即可，用as.matrix()包裹.SD，修改后的代码如下：

x[, variance := rowVars(as.matrix(.SD)), .SDcols = varQuant]

补充说明 #

1. 处理缺失值：如果你的数据里有NA，可以加上na.rm = TRUE参数跳过缺失值计算方差：

   x[, variance := rowVars(as.matrix(.SD), na.rm = TRUE), .SDcols = varQuant]
   

2. 确保列类型正确：要保证varQuant指定的所有列都是数值型连续列，如果混入了因子、字符等非数值列，as.matrix()会把整个矩阵转换成字符型，导致rowVars()再次报错。

示例演示 #

我们用一个简单的例子验证：

library(data.table)
library(matrixStats)

# 创建示例data.table
x <- data.table(
  id = 1:5,
  temp = c(23, 25, 22, 24, 26),
  humidity = c(60, 65, 58, 62, 68),
  pressure = c(1013, 1012, 1015, 1014, 1011)
)

# 指定要计算方差的连续型列
varQuant <- c("temp", "humidity", "pressure")

# 计算行方差
x[, variance := rowVars(as.matrix(.SD)), .SDcols = varQuant]

# 查看结果
x

运行后就能得到每行三个连续列的方差值，不会再触发错误。

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