小编omu*_*mun的帖子

我需要计算:

我可以进一步将其分解为:

在RI中写了这段代码

att_num <- dim(X)[2]
A <- matrix(0, att_num, att_num)
for(i in seq(att_num)) A[,i] <- colSums(dx * X * X[,i])

但由于循环,它是可怕的慢.这一行占用了我脚本中的大部分计算时间.有没有办法可以改进这个计算？

• dx 是一个大小为1 xm的向量
• X 是一个大小为[nxm]的矩阵

例:

dx <- sample(1:100, 30, replace=T)
X <- data.frame(replicate(30,sample(0:1,100,rep=TRUE)))

att_num <- dim(X)[2]
A <- matrix(0, att_num, att_num)
for(i in seq(att_num)) A[,i] <- colSums(dx * X * X[,i])