1
假设我的数据包含3个类别。我想找到每个类别的协方差矩阵,并期望返回是一个数组。如何在data.table中按类别查找协方差矩阵R
set.seed(7)
x <- rbind(cbind(rnorm(120,-1,0.1), rnorm(120,-0.5,0.1)),
cbind(rnorm(120,-0.4,0.1), rnorm(120,0,0.1)),
cbind(rnorm(120,.2,0.1), rnorm(120,0,0.1)))
lab <- c(rep(1, 120), rep(2, 120), rep(3, 120))
newx <- cbind(x, lab = lab)
s <- sapply(1:3, function(k){ var(newx[lab == k, -3]) })
dim(s) <- c(2,2,3)
回报,
, , 1
[,1] [,2]
[1,] 0.008880447 -0.001116058
[2,] -0.001116058 0.009229061
, , 2
[,1] [,2]
[1,] 0.012193536 -0.001217923
[2,] -0.001217923 0.009391710
, , 3
[,1] [,2]
[1,] 0.010752319 0.0
[2,] 0.0.008226595
我知道如何使用sapply
做到这一点,但如果我想使用data.table
,我该怎么办呢?
我曾尝试:
library(data.table)
dt <- data.table(newx)
dt[,lapply(.SD, var), by = lab]
,但如我所料不提供回报。
你想要什么输出?一个数据表是一个数据框,而不是一个三维数组。 – rawr 2014-11-03 03:23:11
@rawr我想要的输出是包含每个类别的协方差的上述数组。在这种情况下,我有3个类别,所以我期望有3个协方差矩阵。 – newbie 2014-11-03 04:16:20