私のコードでは、外部関数は行列を出力します。その出力をスパース行列にすることは可能ですか?この出力行列をスパース行列に変換することはできますが、最初に密行列を作成してから変換するのではなく、outer(......、sparse = T)のようなことをしたいと思います。私は非常に大きなデータを持っているからです。私のコードと短いサンプルベクトルは以下のとおりです。関連するすべての投稿を読みましたが、理解できません。
# myvec is a vector of length 100.
# mycol is a column vector of length 10.
myvec <- c( 0, 1, 1, 2, 1, 2, 2, 3, 1, 2, 1, 0, 2, 1, 2, 1, 3, 2, 2, 1, 0, 1, 2, 1, 0, 2, 1, 2, 1, 3, 2, 2, 1, 0,1,
2, 2, 3, 1, 2, 1, 0, 2, 1, 2, 1,0, 1, 1, 2, 1, 2, 2, 3, 1, 2, 1, 0, 2, 1, 2, 1, 3, 2, 2, 1, 0, 1, 2,
1, 0, 2, 1, 2, 1, 3, 2, 2, 1, 0,1, 2, 2, 3, 1, 2, 1, 0, 2, 1, 2, 1,0, 1, 1, 2, 1, 2, 2, 3)
mycol <- c(-15.07,-7.64, -15.48, -18.59, -12.35, -10.44, -11.45, -18.49, -15.24, -13.23)
res <- outer(1:length(mycol), 1:length(mycol), function(i,j) {
ifelse(myvec > 1, 0,
ifelse(myvec == 1, (mycol[j]-mycol[i])/(1-exp(-(mycol[j]-mycol[i]))),
ifelse(myvec ==0, 1, myvec)))})
また、行、列、値のトリプレットを使用してスパース行列を最初から作成しようとしました。コードは次のとおりです。ただし、エントリ(resvec)の長さが100であるため、スパース行列が作成されますが、対角へのエントリの最初の10のみが使用され、残りの非対角要素はゼロになります。しかし、私は100のエントリすべてをそれらの場所に書き込みたかったのです。
よろしくお願いします。
resvec<- as.vector(res)
matSparse <- sparseMatrix(
i = 1:length(mycol),
j = 1:length(mycol),
x = resvec,
dims = c(length(mycol),length(mycol)),
dimnames = list(rownames,colnames)
)
ps:行名と列名は同じです。
このようなもので、ゼロ以外の結果を生成する値に対してouter
のみ、何をするかをやり直していmyvec
ます。関数f
を少し単純化して、おそらくそれを関数にすることmyvec
もできます。
f <- function(i,j) {
ifelse(myvec > 1, 0,
ifelse(myvec == 1, (mycol[j]-mycol[i])/(1-exp(-(mycol[j]-mycol[i]))),
ifelse(myvec ==0, 1, myvec)))}
library(Matrix)
i <- rep(1:length(mycol), times = length(mycol))[myvec < 2]
j <- rep(1:length(mycol), each = length(mycol))[myvec < 2]
myvec <- myvec[myvec < 2]
sparseMatrix(i, j, x = f(i, j))
#> 10 x 10 sparse Matrix of class "dgCMatrix"
#>
#> [1,] 1.000000000 7.434410 1.000000000 . 2.9118151 4.675609 . 1.0000000 0.91740717 .
#> [2,] 0.004409999 1.000000 0.003087581 . 1.0000000 . 0.08629586 . . 0.02095708
#> [3,] 1.218969243 . . 0.14518433 . . . 0.1560604 . 1.00000000
#> [4,] . 10.950192 3.255184331 1.00000000 6.2521909 . . . . 5.38531586
#> [5,] 0.191815105 . 1.000000000 0.01219087 . 2.241996 . 0.0132598 0.17006693 0.62371549
#> [6,] . 2.981293 . . 0.3319956 . 0.57859728 . . .
#> [7,] . . 0.072926774 . 1.0000000 1.588597 1.00000000 . 0.08761709 0.36106528
#> [8,] . . . . 6.1532598 1.000000 7.04617334 . 1.00000000 .
#> [9,] 1.087407174 . 0.884795398 0.12180607 3.0600669 . . 0.1310995 . .
#> [10,] . 5.610957 . . . 2.972579 2.14106528 1.0000000 0.31098581 .
この記事はインターネットから収集されたものであり、転載の際にはソースを示してください。
侵害の場合は、連絡してください[email protected]
コメントを追加