function [p,stats]=quantreg(x,y,tau,order,Nboot)
% Quantile Regression
%
% USAGE: [p,stats]=quantreg(x,y,tau[,order,nboot]);
%
% INPUTS:
% x,y: data that is fitted. (x and y should be columns)
% Note: that if x is a matrix with several columns then multiple
% linear regression is used and the "order" argument is not used.
% tau: quantile used in regression.
% order: polynomial order. (default=1)
% (negative orders are interpreted as zero intercept)
% nboot: number of bootstrap surrogates used in statistical inference.(default=200)
%
% stats is a structure with the following fields:
% .pse: standard error on p. (not independent)
% .pboot: the bootstrapped polynomial coefficients.
% .yfitci: 95% confidence interval on "polyval