R KNex Koenker-Ng 稀疏模型矩阵和响应向量示例

说明

Koenker 和 Ng 在示例中使用的模型矩阵 mm 和相应的响应向量 y。矩阵 mm 是一个稀疏矩阵，有 1850 行和 712 列，但只有 8758 个非零条目。它是一个"dgCMatrix" 对象。向量y就是长度为1850的numeric。

用法

data(KNex)

例子

data(KNex, package = "Matrix")
class(KNex$mm)
dim(KNex$mm)
image(KNex$mm)
str(KNex)

system.time( # a fraction of a second
  sparse.sol <- with(KNex, solve(crossprod(mm), crossprod(mm, y))))

head(round(sparse.sol,3))

## Compare with QR-based solution ("more accurate, but slightly slower"):
system.time(
  sp.sol2 <- with(KNex, qr.coef(qr(mm), y) ))

all.equal(sparse.sol, sp.sol2, tolerance = 1e-13) # TRUE

参考

Roger Koenker and Pin Ng (2003). SparseM: A sparse matrix package for R; J. of Statistical Software, 8 (6), doi:10.18637/jss.v008.i06