
 z <- .Fortran("rq1", as.integer(m0fn), as.integer(p), 
        as.integer(mofn + 5), as.integer(p + 2), 
        as.double(x[s[,i],]), as.double(y[s[,i]]), as.double(tau), as.double(tol), 
        flag = as.integer(1), coef = double(p), resid = double(mofn), 
        integer(mofn), double((mofn + 5) * (p + 2)), double(m), as.integer(1), 
        sol = double((p + 2)), dsol = double(m), lsol = as.integer(1), 
        PACKAGE = "quantreg")

