tests/newquantile.R

In survey: Analysis of Complex Survey Samples

## quantiles with equal weights

library(survey)
data(api)
dclus1<-svydesign(id=~dnum, weights=~pw, data=apiclus1, fpc=~fpc)


for(i in 1:9){
    print(i)
    all.equal(
        as.vector(coef(svyquantile(~ell, dclus1, c(0.2,0.5,0.9), qrule=paste0("hf",i)))),
        as.vector(quantile(apiclus1$ell,  c(0.2,0.5,0.9), type=i))
    )
    
}