如何執行的矩陣計算（例如協方差矩陣），用於不同等級因子可變的

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)