Function performing the first bootstrap procedure to yield the parameter estimates
boot1(y, phi1, arp, nbs, x, allb, method, scores)the response variable
the Durbin two-stage estimate of the autoregressive parameter rho
the order of autoregressive errors
the bootstrap size
the original design matrix (including intercept), without centering
all the Durbin two-stage estimates of the regression coefficients
If "OLS", uses the ordinary least square; If "RANK", uses the rank-based fit
Default is Wilcoxon scores
An estimate of the bias is returned