0
require("quantreg")
require("SparseM")
Adata <- read.csv("A.CSV", header=FALSE)
Amat = as.matrix(Adata)
A <- as.matrix.csr(Amat)
bdata <- read.csv("b.CSV", header=FALSE)
bmat = as.matrix(bdata)
b <- as.matrix.csr(bmat)
dim(A)
# [1] 156 39
dim(b)
# [1] 156 1
rq.fit.sfn(A, b, tau = 0.5)
# Error in rq.fit.sfn(A, b, tau = 0.5) :
# Dimensions of design matrix and the response vector not compatible
有人可以請解釋爲什麼我得到上述錯誤信息? A和b具有相同的行數。R稀疏矩陣支持quantreg
從這裏看一個例子:http://astrostatistics.psu.edu/datasets/2006tutorial/html/quantreg/html/rq.fit.sfn.html該函數似乎想要在您的設計矩陣中的一列1'w大概是攔截?,試試'A < - as.matrix.csr(cbind(1,Amat) )' – Vlo