2

行列のすべての行の外積を計算する効率的な方法を見つけるのに助けが必要です。

これまで私は試しました

X <- matrix(1:100, ncol=4)
lapply(1:i, function(i) tcrossprod(X[i,]))

「現実の」(別名「大きな」)マトリックスに適用すると、かなり遅くて面倒なので、ゾッとします。

結果の型が配列かリストかは問題ではありません。

何か案は?plyr友達も遅くなります。

4

1 に答える 1

1

以下のコードで目的の出力が得られますか? もしそうなら、それはあなたのデータで十分に速いですか?

library(data.table)

# Example matrix
m <- matrix(1:100, ncol=4)

# Convert to data.table
X <- as.data.table(m)

# Add row numbers
X[,row:=1:nrow(X)]

# Reshape to long format
X.long<-data.table(reshape(X, direction="long", varying=list(rev(rev(names(X))[-1])), v.names="value", idvar="row", timevar="col", times=1:(ncol(X)-1)),key="row")

# Join the long data.table to itself, perform the products, and then reshape back to wide format
reshape(X.long[X.long,list(col=i.col, col.old=col, prod=value*i.value)], direction="wide", idvar=c("row","col"), timevar="col.old")

##     row col prod.1 prod.2 prod.3 prod.4
##  1:   1   1      1     26     51     76
##  2:   1   2     26    676   1326   1976
##  3:   1   3     51   1326   2601   3876
##  4:   1   4     76   1976   3876   5776
##  5:   2   1      4     54    104    154
##  6:   2   2     54    729   1404   2079
##  7:   2   3    104   1404   2704   4004
##  8:   2   4    154   2079   4004   5929
##  ... etc
于 2013-02-15T05:12:24.143 に答える