1## This is just from ?anova.rq (extended)
2library(quantreg)
3data(barro)
4fit0 <- rq(y.net ~  lgdp2 + fse2 + gedy2 , data = barro)
5fit1 <- rq(y.net ~  lgdp2 + fse2 + gedy2 + Iy2 + gcony2, data = barro)
6
7a01 <- anova(fit1,fit0)
8a01
9
10fit2 <- rq(y.net ~  lgdp2 + fse2 + gedy2 + Iy2 + gcony2, data = barro,
11	   tau = 0.75)
12fit3 <- rq(y.net ~  lgdp2 + fse2 + gedy2 + Iy2 + gcony2, data = barro,
13	   tau = 0.25)
14
15a123  <- anova(fit1,fit2, fit3)
16a.123 <- anova(fit1,fit2, fit3, joint=FALSE)
17a.123
18
19AE <- function(x,y) all.equal(x, y, tol = 1e-5)
20##                                  ---------- {giving a bit more digits below}
21stopifnot(
22	  AE(100 * unname(coef(fit0)),
23	     c(-0.74679759, 0.46539963, 0.15902838, -36.619915))
24	  ,
25	  AE(unlist(a01$table),
26	     c(ndf=2, ddf=155, Tn = 18.878717, pvalue= 4.6e-08))
27	  ,
28	  AE(100* unname(coef( fit2 )),
29	     c(13.103018, -1.4885239, -0.026452369,
30	       0.3999839, 14.526663, -13.504643))
31	  ,
32	  AE(100* unname(coef( fit3 )),
33	     c(6.0860719, -0.88350554, 0.24596781,
34	       -14.962498, 15.592489, -15.861804))
35	  ,
36	  AE(unlist(a123$table),
37	     c(ndf = 10, ddf = 473, Tn = 1.80385526, pvalue=0.0575117558))
38	  ,
39	  AE(a.123$table[,"Tn"],
40	     c(1.0655561, 2.6398508, 0.78623238, 0.04467014, 0.065344348))
41	  )
42