如何执行的矩阵计算（例如协方差矩阵），用于不同等级因子可变的

Question

我有以下简化问题：

temp <- matrix(rnorm(900), ncol = 3) 
lev <- as.factor(rep(c("a", "b", "c"), each = 100)) 
dfr <- data.frame(lev = lev, temp = temp) 

我想calcualte方差 - 协方差矩阵的每个的三个层次。我如何有效地做到这一点？

非常感谢花时间阅读（并希望回答）我的问题！

Answer 1

没关系。我找到了解决办法：

lapply(split(dfr, dfr$lev), function(x) var(x[,-1])) 

Answer 2

我们可以使用data.table

library(data.table) 
setDT(dfr)[, as.list(var(.SD)), lev] 

或者使用by从base R

by(dfr[-1], list(dfr[,1]), FUN = var)