R 分组矩阵乘法
我在做两个矩阵的乘法,比如beta和Xs。 但我的数据有一个组变量,我按组进行了回归。 下面是我正在尝试做的一个例子R 分组矩阵乘法,r,matrix,multiplication,R,Matrix,Multiplication,我在做两个矩阵的乘法,比如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) colname
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我们可以循环通过'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