代码之家 › 专栏 › 技术社区 › wan

按组的矩阵乘法

multiplication matrix r

wan · 技术社区 · 10 年前

我在做两个矩阵的乘法,比如beta和Xs。但我的数据有一个组变量,我按组进行了回归。这是我要做的一个例子。

matx <- matrix(c(1,1,2,1,3,4,2,5,6,2,7,8,3,9,10,3,11,12),nrow=6,ncol=3,byrow=T)
colnames(matx) <- c("group","x1","x2")
b <- matrix(c(1,3,4,2,5,6,3,0,0),nrow=3,ncol=3,byrow=T)
colnames(b) <- c("group","b1","b2")

我想做的是乘法 matx 和 b 通过 group 变量组1的x1和x2应乘以组1的b1和b2。是否有任何代码可以简化数据的分组?

1 回复 | 直到 10 年前

akrun 10 年前

我们可以循环通过 unique “group”的值,同时子集矩阵和do %*%

sapply(unique(matx[,'group']), function(i) 
          matx[matx[,1]==i,-1]%*% b[b[,1]==i, -1])
#     [,1] [,2] [,3]
#[1,]   11   61    0
#[2,]   25   83    0

推荐文章

Jay · JavaScript中的乘法表-2个表

8 年前

Ryan · 使用16位DSP的33*33位

8 年前

Fizi · 熊猫groupby:获取每个群体的产品价值

8 年前

user2079550 · 如何快速执行类似于乘法的矩阵运算?

8 年前

Siebe Bosch · 为什么这个乘法会导致OverflowException?

8 年前

Micael Salomon · 根据特定列的值乘以值

8 年前

juhist · 为什么gcc或clang on ARM不使用“使用乘法除以不变整数”的技巧?

8 年前

John · 在Java中,如何将两个一维数组的对应值相乘?

8 年前

Michael26 · 乘法公式字符串-java

8 年前

Serge Kashlik · 如何在SQL中计算并添加新列

8 年前