pctiles <- seq( 0 , 1 , .01 )

library(DBI)            # load the DBI package (implements the R-database coding)
library(MonetDB.R)        # load the MonetDB.R package (connects r to a monet database)
library(MonetDBLite)    # load MonetDBLite package (creates database files in R)

db <- dbConnect( MonetDBLite() )

# requires a unique identifier because i could not get row_number() working
# https://www.monetdb.org/bugzilla/show_bug.cgi?id=4000
mtcars$id <- seq( nrow( mtcars ) )

# store mtcars
dbWriteTable( db , 'mtcars' , mtcars )


sql <-
    "SELECT
        t1.id,
        t1.hp ,
        t1.wt ,
        SUM( t2.wt ) AS cumwt
    FROM mtcars t1 INNER JOIN mtcars t2 ON t1.hp >= t2.hp
    GROUP BY t1.id , t1.hp , t1.wt"

dbGetQuery( db , sql )

sql2 <-
    paste("( SELECT
            row_number() over () AS id ,
            wt ,
            cumwt,
            totwt,
            hp
    FROM
            mtcars AS td
            NATURAL JOIN
            (
                    " , sql , "
            ) AS cj
            ,
            (
                    SELECT
                            SUM(wt) AS totwt
                    FROM
                            mtcars
            ) AS tj ) ")


dbGetQuery( db , sql2 )

sql3 <-
    paste(
        "select pcs.percentile ,
            min(case when cumwt >= totwt * pcs.percentile then hp end)
            from (" , sql2 , ") AS v ,
     (" ,
     
     paste( "select" , pctiles , " as percentile" , collapse = " union all " ) ,
     " ) pcs group by pcs.percentile" )


library(Hmisc)
library(survey)
ms <- svydesign(~1,data=mtcars,weight=~wt)

cbind( dbGetQuery( db , sql3 ),
wtd.quantile( mtcars$hp , mtcars$wt , pctiles ) ,
as.numeric(svyquantile(~ hp,ms,pctiles)))