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