Function intended for multiple response quantile regression
called from boot.rq
for wild bootstrap option.
rqs.fit(x, y, tau=0.5, tol = 0.0001)
the design matrix an n by p matrix.
the response variable as a n by m matrix
the quantile desired, if tau lies outside (0,1)
tolerance parameter for Barrodale and Roberts exterior point method.