Question

如何在R（Matrix包）中使用sparseMatrix命令构建

下面的矩阵

[1,] 1.002 0.210 0.002  .      .      .       2.943   0.051
[2,] 0.210 1.002 0.210  .      .      .       7.515   2.943
[3,] 0.002 0.210 1.002  .      .      .       0.843   7.515
[4,] .     .     .     16.003  3.354  0.031  18.691   1.122
[5,] .     .     .      3.354 16.003  3.354  13.675  18.691
[6,] .     .     .      0.031  3.354 16.003   0.440  13.675
[7,] 2.943 7.515 0.843 18.691 13.675  0.440 109.002  22.848
[8,] 0.051 2.943 7.515  1.122 18.691 13.675  22.848 109.002

我还提供了dput命令，可以将矩阵复制到你的控制台

new("dsCMatrix"
    , i = c(0L, 0L, 1L, 0L, 1L, 2L, 3L, 3L, 4L, 3L, 4L, 5L, 0L, 1L, 2L, 
3L, 4L, 5L, 6L, 0L, 1L, 2L, 3L, 4L, 5L, 6L, 7L)
    , p = c(0L, 1L, 3L, 6L, 7L, 9L, 12L, 19L, 27L)
    , Dim = c(8L, 8L)
    , Dimnames = list(NULL, NULL)
    , x = c(1.002, 0.21, 1.002, 0.002, 0.21, 1.002, 16.003, 3.354, 16.003, 
0.031, 3.354, 16.003, 2.943, 7.515, 0.843, 18.691, 13.675, 0.44, 
109.002, 0.051, 2.943, 7.515, 1.122, 18.691, 13.675, 22.848, 
109.002)
    , uplo = "U"
    , factors = list()
)

Answer 1

你有一个CsparseMatrix并且他们没有列索引，但很容易强制到拥有它们的TsparseMatrix。 i和j参数（内部x值中存在的值）基于0而不是1，所以如果你想用它们解决问题，你需要加1。

 CSM <- new("dsCMatrix"
     , i = c(0L, 0L, 1L, 0L, 1L, 2L, 3L, 3L, 4L, 3L, 4L, 5L, 0L, 1L, 2L, 
 3L, 4L, 5L, 6L, 0L, 1L, 2L, 3L, 4L, 5L, 6L, 7L)
     , p = c(0L, 1L, 3L, 6L, 7L, 9L, 12L, 19L, 27L)
     , Dim = c(8L, 8L)
     , Dimnames = list(NULL, NULL)
     , x = c(1.002, 0.21, 1.002, 0.002, 0.21, 1.002, 16.003, 3.354, 16.003, 
 0.031, 3.354, 16.003, 2.943, 7.515, 0.843, 18.691, 13.675, 0.44, 
 109.002, 0.051, 2.943, 7.515, 1.122, 18.691, 13.675, 22.848, 
 109.002)
     , uplo = "U"
     , factors = list()
 )
 TSM <-as(CSM, "TsparseMatrix")
 TSM@i
# [1] 0 0 1 0 1 2 3 3 4 3 4 5 0 1 2 3 4 5 6 0 1 2 3 4 5 6 7
 TSM@j
# [1] 0 1 1 2 2 2 3 4 4 5 5 5 6 6 6 6 6 6 6 7 7 7 7 7 7 7 7
 TSM@j +1
# [1] 1 2 2 3 3 3 4 5 5 6 6 6 7 7 7 7 7 7 7 8 8 8 8 8 8 8 8
 TSM@i +1
# [1] 1 1 2 1 2 3 4 4 5 4 5 6 1 2 3 4 5 6 7 1 2 3 4 5 6 7 8

R中的sparseMatrix命令

1 个答案: