cqr.cd | R Documentation |
Composite quantile regression (cqr) find the estimated coefficient which minimize the absolute error for various quantile level. The algorithm base on greedy coordinate descent and Edgeworth's for ordinary l_1 regression.
cqr.cd(X,y,tau,beta,maxit,toler)
X |
the design matrix |
y |
response variable |
tau |
vector of quantile level |
beta |
initial value of estimate coefficient (default naive guess by least square estimation) |
maxit |
maxim iteration (default 200) |
toler |
the tolerance critical for stop the algorithm (default 1e-3) |
a list
structure is with components
beta |
the vector of estimated coefficient |
b |
intercept |
cqr.cd(x,y,tau) work properly only if the least square estimation is good.
Wu, T.T. and Lange, K. (2008). Coordinate Descent Algorithms for Lasso Penalized Regression. Annals of Applied Statistics, 2, No 1, 224–244.
Hui Zou and Ming Yuan(2008). Composite Quantile Regression and the Oracle Model Selection Theory, The Annals of Statistics, 36, Number 3, Page 1108–1126.
set.seed(1) n=100 p=2 a=rnorm(n*p, mean = 1, sd =1) x=matrix(a,n,p) beta=rnorm(p,1,1) beta=matrix(beta,p,1) y=x%*%beta-matrix(rnorm(n,0.1,1),n,1) tau=1:5/6 # x is 1000*10 matrix, y is 1000*1 vector, beta is 10*1 vector cqr.cd(x,y,tau)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.