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