矩阵列运算

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

比方说有一个n*p矩阵，每一列分别是A1,A2,....,Ap, 想要做成一个新的p*(p-1)/2列矩阵分别是(A1*A2, A1*A3,.....,A(p-1)*Ap)
是否可以不用循环写出？

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

1. p <- ncol(matr)
   
2. a <- lapply(2:p, fun <- function(x){matr[,c(x:p)]})
   
3. b <- lapply(1:(p-1), fun <- function(x){matr[,x]*a[[x]]})
   
4. result <- do.call("cbind", b)

复制代码