R构造稀疏矩阵

我正在阅读 R 中 Matrix 包的说明。但我无法理解p函数中的参数：

sparseMatrix(i = ep, j = ep, p, x, dims, dimnames,
         symmetric = FALSE, index1 = TRUE,
         giveCsparse = TRUE, check = TRUE)

根据http://stat.ethz.ch/R-manual/R-devel/library/Matrix/html/sparseMatrix.html http://stat.ethz.ch/R-manual/R-devel/library/Matrix/html/sparseMatrix.html

p:
指针的数字（整数值）向量，每一列（或行）一个，指向 列（或行）中元素的初始（从零开始）索引。恰好是 i、j 或 p 之一 一定失踪了。

我估计p用于行或列索引的压缩表示，因为在其中包含多个元素是浪费的i or j具有相同的值来表示单个行/列。但是当我尝试提供的示例时，我仍然不明白如何p正在控制哪个元素x转到哪一行/哪一列

dn <- list(LETTERS[1:3], letters[1:5])
## pointer vectors can be used, and the (i,x) slots are sorted if necessary:
m <- sparseMatrix(i = c(3,1, 3:2, 2:1), p= c(0:2, 4,4,6), x = 1:6, dimnames = dn)

只需再往下读一点?SparseMatrix学习如何p被解释。 （特别要注意关于“扩展形式”的部分p.)

如果缺少“i”或“j”，则“p”必须是非递减整数 第一个元素为零的向量。它提供了压缩的、 或行或列索引的“指针”表示， 以缺失者为准。 ‘p’的扩展形式， ‘rep(seq_along(dp),dp)’，其中‘dp

这里有一个小功能可以帮助你see这在实践中意味着什么：

pex <- function(p) {
    dp <- diff(p)
    rep(seq_along(dp), dp)
}

## Play around with the function to discover the indices encoded by p.
pex(p = c(0,1,2,3))
# [1] 1 2 3

pex(p = c(0,0,1,2,3))
# [1] 2 3 4

pex(p = c(10,11,12,13))
# [1] 1 2 3

pex(p = c(0,0,2,5))
# [1] 2 2 3 3 3

pex(p = c(0,1,3,3,3,3,8))
# [1] 1 2 2 6 6 6 6 6